The UCSC Genome Browser database: 2019 update

Multifaceted investigations of PSMB8 provides insights into prognostic prediction and immunological target in thyroid carcinoma

The Proteasome 20S subunit beta 8 (PSMB8) is an integral element of the immunoproteasome complex, playing a pivotal role in antigen processing. Despite its significance, the contributory role of PSMB8 in oncogenesis, particularly in thyroid carcinoma (THCA), has not been well-characterized. To address this gap in knowledge, our study endeavored to delineate the potential associations between PSMB8 and THCA. Transcriptomic profiles and clinical data of patients with THCA were retrieved from The Cancer Genome Atlas (TCGA) database to facilitate comprehensive analysis. Complementary resources from additional online databases were utilized to augment the study. Logistic regression analysis was employed to elucidate the relationship between PSMB8 and various clinicopathological parameters. Uni/multivariate Cox regression analyses were conducted to ascertain the independent prognostic factors for THCA patient outcomes. Quantitative polymerase chain reaction (qPCR) and western blot assays were employed to verify the expression level of PSMB8 in vitro. Our study demonstrated that PSMB8 was significantly upregulated in THCA, with its overexpression correlating with lymph node metastasis, extrathyroidal extension, and favorable prognosis. Immunohistochemistry substantiated a higher PSMB8 protein presence in THCA tissue compared to the normal, supporting its potential role as a moderately accurate diagnostic biomarker. Logistic regression analysis identified PSMB8 as a significant indicator of the N1 stage, classical histological subtype, and extrathyroidal extension. Age, T stage, and PSMB8 were further determined as independent prognostic factors for THCA. Functional investigations linked PSMB8 to immune processes, evidenced by its association with increased immune cell infiltration and higher stromal/immune scores, as well as a positive co-expression with several immune checkpoints. A constructed predicted competing endogenous RNA (ceRNA) network implicated PSMB8 in complex post-transcriptional regulation. Finally, in vitro assays confirmed the upregulation of PSMB8, underscoring its relevance in THCA and as a target for future research. Our work has preliminarily appraised PSMB8 as a biomarker with certain prognostic and diagnostic significance, and as a potential target for immunotherapy in THCA.

TRAPT: a multi-stage fused deep learning framework for predicting transcriptional regulators based on large-scale epigenomic data

It is challenging to identify regulatory transcriptional regulators (TRs), which control gene expression via regulatory elements and epigenomic signals, in context-specific studies on the onset and progression of diseases. The use of large-scale multi-omics epigenomic data enables the representation of the complex epigenomic patterns of control of the regulatory elements and the regulators. Herein, we propose Transcription Regulator Activity Prediction Tool (TRAPT), a multi-modality deep learning framework, which infers regulator activity by learning and integrating the regulatory potentials of target gene cis-regulatory elements and genome-wide binding sites. The results of experiments on 570 TR-related datasets show that TRAPT outperformed state-of-the-art methods in predicting the TRs, especially in terms of forecasting transcription co-factors and chromatin regulators. Moreover, we successfully identify key TRs associated with diseases, genetic variations, cell-fate decisions, and tissues. Our method provides an innovative perspective on identifying TRs by using epigenomic data.

Identification and validation of pyroptosis patterns in AML via comprehensive bioinformatics analysis

Pyroptosis, a lytic inflammatory cell death mechanism, plays dual roles in tumorigenesis, but its clinical relevance in acute myeloid leukemia (AML) remains poorly understood. Through an integrative analysis of 40 pyroptosis-related genes in newly diagnosed AML patients (TCGA, n = 151) and healthy controls (GTEx, n = 386), we identified 32 genes with aberrant expression. Among these, 9 genes were found to be significant prognostic markers, including ELANE (protective), and CASP1, CHMP4B, BAK1, and CHMP2A (risk), which retained their prognostic significance after adjusting for age and gender. Using unsupervised nonnegative matrix factorization (NMF) on TCGA data, we classified AML into two pyroptosis patterns: the ELANEhigh subtype, associated with favorable survival, and the ELANElow subtype, which was enriched in poor karyotypes and adverse outcomes. This classification was validated in an independent cohort (GSE10358, n = 91). Single-cell RNA sequencing data (GSE116256, n = 15) revealed that the ELANElow subtype is characterized by an immunologically active microenvironment, marked by an expansion of cytotoxic T cells and naive CD4 + /CD8 + T cells. Factor analysis revealed associations between pyroptosis patterns and other forms of cell death, including ferroptosis, autophagy, and apoptosis, as well as with karyotype, leukemia stemness, and TP53/FLT3-ITD mutations. Prognostic immune gene sets enriched in the ELANElow subtype were associated with interferon signaling and ubiquitin-mediated degradation pathways. Furthermore, protein-protein interaction (PPI) network analysis identified three sub-networks and nine key hub genes. This study integrates gene expression data from newly diagnosed AML patients, revealing the heterogeneity of pyroptosis patterns within the population. It highlights the potential links between distinct pyroptosis patterns, the immune microenvironment, various cell death pathways, leukemia stemness, and genomic alterations, offering novel biomarkers and therapeutic targets for risk stratification and immunomodulatory interventions in AML.

Genomic analysis of 11,555 probands identifies 60 dominant congenital heart disease genes

Congenital heart disease (CHD) is a leading cause of infant mortality. We analyzed de novo mutations (DNMs) and very rare transmitted/unphased damaging variants in 248 prespecified genes in 11,555 CHD probands. The results identified 60 genes with a significant burden of heterozygous damaging variants. Variants in these genes accounted for CHD in 10.1% of probands with similar contributions from de novo and transmitted variants in parent-offspring trios that showed incomplete penetrance. DNMs in these genes accounted for 58% of the signal from DNMs. Thirty-three genes were linked to a single CHD subtype while 12 genes were associated with 2 to 4 subtypes. Seven genes were only associated with isolated CHD, while 37 were associated with 1 or more extracardiac abnormalities. Genes selectively expressed in the cardiomyocyte lineage were associated with isolated CHD, while those widely expressed in the brain were also associated with neurodevelopmental delay (NDD). Missense variants introducing or removing cysteines in epidermal growth factor (EGF)-like domains of NOTCH1 were enriched in tetralogy of Fallot and conotruncal defects, unlike the broader CHD spectrum seen with loss of function variants. Transmitted damaging missense variants in MYH6 were enriched in multiple CHD phenotypes and account for ~1% of all probands. Probands with characteristic mutations causing syndromic CHD were frequently not diagnosed clinically, often due to missing cardinal phenotypes. CHD genes that were positively or negatively associated with development of NDD suggest clinical value of genetic testing. These findings expand the understanding of CHD genetics and support the use of molecular diagnostics in CHD.

Widespread impact of transposable elements on the evolution of post-transcriptional regulation in the cotton genus Gossypium

BACKGROUND

Transposable element (TE) expansion has long been known to mediate genome evolution and phenotypic diversity in organisms, but its impact on the evolution of post-transcriptional regulation following species divergence remains unclear.

RESULTS

To address this issue, we perform long-read direct RNA sequencing, polysome profiling sequencing, and small RNA sequencing in the cotton genus Gossypium, the species of which range more than three folds in genome size. We find that TE expansion contributes to the turnover of transcription splicing sites and regulatory sequences, leading to changes in alternative splicing patterns and the expression levels of orthologous genes. We also find that TE-derived upstream open reading frames and microRNAs serve as regulatory elements mediating differences in the translation levels of orthologous genes. We further identify genes that exhibit lineage-specific divergence at the transcriptional, splicing, and translational levels, and showcase the high flexibility of gene expression regulation in the evolutionary process.

CONCLUSIONS

Our work highlights the significant role of TE in driving post-transcriptional regulation divergence in the cotton genus. It offers insights for deciphering the evolutionary mechanisms of cotton species and the formation of biological diversity.

Glutamine missense suppressor transfer RNAs inhibit polyglutamine aggregation

Huntington's disease (HD) is caused by polyglutamine (polyQ) repeat expansions in the huntingtin gene. HD-causative polyQ alleles lead to protein aggregation, which is a prerequisite for disease. Translation fidelity modifies protein aggregation, and several studies suggest that mutating one or two glutamine (Gln) residues in polyQ reduces aggregation. Thus, we hypothesized that missense suppression of Gln codons with other amino acids will reduce polyQ aggregate formation in cells. In neuroblastoma cells, we assessed tRNA variants that misread Gln codons with serine (tRNASer C/UUG) or alanine (tRNAAla C/UUG). The tRNAs with the CUG anticodon were more effective at suppressing the CAG repeats in polyQ, and serine and alanine mis-incorporation had differential impacts on polyQ. The expression of tRNASer CUG reduced polyQ protein production as well as both soluble and insoluble aggregate formation. In contrast, cells expressing tRNAAla CUG selectively decreased insoluble polyQ aggregate formation by 2-fold. Mass spectrometry confirmed Ala mis-incorporation at an average level of ∼20% per Gln codon. Cells expressing the missense suppressor tRNAs showed no cytotoxic effects and no defects in growth or global protein synthesis levels. Our findings demonstrate that tRNA-dependent missense suppression of Gln codons is well tolerated in mammalian cells and significantly reduces polyQ levels and aggregates that cause HD.

Optical genome and epigenome mapping of clear cell renal cell carcinoma

Cancer cells display complex genomic aberrations that include large-scale genetic rearrangements and epigenetic modulation that are not easily captured by short-read sequencing. This study presents a novel approach for simultaneous profiling of long-range genetic and epigenetic changes in matched cancer samples, focusing on clear cell renal cell carcinoma (ccRCC). ccRCC is a common kidney cancer subtype frequently characterized by a 3p deletion and the inactivation of the von Hippel-Lindau (VHL) gene. We performed integrated genetic, cytogenetic, and epigenetic analyses on paired tumor and adjacent nontumorous tissue samples. Optical genome mapping identified genomic aberrations as structural and copy number variations, complementing exome-sequencing findings. Single-molecule methylome and hydroxymethylome mapping revealed a significant global reduction in 5hmC level in both sample pairs, and a correlation between both epigenetic signals and gene expression was observed. The single-molecule epigenetic analysis identified numerous differentially modified regions, some implicated in ccRCC pathogenesis, including the genes VHL, PRCC, and PBRM1. Notably, pathways related to metabolism and cancer development were significantly enriched among these differential regions. This study demonstrates the feasibility of integrating optical genome and epigenome mapping for comprehensive characterization of matched tumor and adjacent tissue, uncovering both established and novel somatic aberrations.

Investigation of a pathogenic inversion in UNC13D and comprehensive analysis of chromosomal inversions across diverse datasets

Inversions are known contributors to the pathogenesis of genetic diseases. Identifying inversions poses significant challenges, making it one of the most demanding structural variants (SVs) to detect and interpret. Recent advancements in sequencing technologies and the development of publicly available SV datasets have substantially enhanced our capability to explore inversions. However, a cross-comparison in those datasets remains unexplored. In this study, we reported a proband with familial hemophagocytic lymphohistiocytosis type-3 carrying a splicing variant (c.1389+1G>A) in trans with an inversion present in 0.006345% of individuals in gnomAD (v4.0) that disrupts UNC13D. Based on this result, we investigate the features of potentially pathogenic inversions in gnomAD which revealed 98.9% of them are rare and disrupt 5% of protein-coding genes associated with a phenotype in OMIM. We then conducted a comparative analysis of additional public datasets, including DGV, 1KGP, and two recent studies from the Human Genome Structural Variation Consortium which revealed common and dataset-specific inversion characteristics suggesting methodology detection biases. Next, we investigated the genetic features of inversions disrupting the protein-coding genes. Notably, we found that the majority of protein-coding genes in OMIM disrupted by inversions are associated with autosomal recessive phenotypes supporting the hypothesis that inversions in trans with other variants are potential hidden causes of monogenic diseases. This effort aims to fill the gap in our understanding of the molecular characteristics of inversions with low frequency in the population and highlight the importance of identifying them in rare disease studies.

Deep5mC: Predicting 5-methylcytosine (5mC) methylation status using a deep learning transformer approach

DNA methylations, such as 5-methylcytosine (5mC), are crucial in biological processes, and aberrant methylations are strongly linked to various human diseases. Genomic 5mC is not randomly distributed but exhibits a strong association with genomic sequences. Thus, various computational methods were developed to predict 5mC status based on DNA sequences. These methods generated promising achievements and overcome the limitations of experimental approaches. However, few studies have comprehensively investigated the dependency of 5mC on genomic sequences, and most existing methods focus on specific genomic regions. In this work, we introduce Deep5mC, a deep learning transformer-based method designed to predict 5mC methylations. Deep5mC leverages long-range dependencies within genomic sequences to estimate the probability of cytosine methylations. Through cross-chromosome evaluation, Deep5mC achieves Matthew's correlation coefficient over 0.86 and F1-score over 0.93, substantially outperforming state-of-the-art methods. Deep5mC not only confirms the influence of long-range sequence context on 5mC prediction but also paves the way for further studying 5mC-sequence dependency across species and in human diseases.

LncPepAtlas: a comprehensive resource for exploring the translational landscape of long non-coding RNAs

Long non-coding RNAs were commonly viewed as non-coding elements. However, they are increasingly recognized for their ability to be translated into proteins, thereby playing a significant role in various cellular processes and diseases. With developments in biotechnology and computational algorithms, a range of novel approaches are being applied to investigate the translation of long non-coding RNA (lncRNAs). Herein, we developed the LncPepAtlas database (http://www.cnitbiotool.net/LncPepAtlas/), which aims to compile multiple evidences for the translation of lncRNAs and annotations for the upstream regulation of lncRNAs across various species. LncPepAtlas integrated compelling evidence from nine distinct sources for the translation of lncRNAs. These include a dataset comprising 2631 publicly available Ribo-seq samples from nine species, which has been collected and analysed. LncPepAtlas offers extensive annotation for lncRNA upstream regulation and expression profiles across various cancers, tissues or cell lines at transcriptional and translational levels. Importantly, it enables novel antigen predictions for lncRNA-encoded peptides. By identifying numerous peptide candidates that could potentially bind to major histocompatibility complex class I and II molecules, this work may provide new insights into cancer immunotherapy. The function of peptides were inferred by aligning them with experimentally detected proteins. LncPepAtlas aims to become a convenient resource for exploring translatable lncRNAs.

Aortic Root Dilation and Genotype Associations in Phelan-McDermid Syndrome

Phelan-McDermid syndrome (PMS) is a rare genetic neurodevelopmental disorder that results from the loss of one functional copy of the SHANK3 gene. While many clinical features of PMS are well-understood, there is currently limited literature on cardiovascular abnormalities in PMS. This report aims to evaluate the prevalence of aortic root dilation (ARD) among individuals with PMS and to understand if underlying genetic variation relates to risk for ARD. We present findings from 59 participants collected from a multisite observational study evaluating the phenotype and natural history of PMS. Individual echocardiographic and genetic reports were analyzed for aortic root measurements and genetic variant data, respectively. Our a priori hypothesis was that participants with chromosome 22 deletions with hg19 start coordinates on or before 49,900,000 (larger deletions) would have more instances of ARD than participants with deletion start coordinates after 49,900,000 (smaller deletions). Eight participants (14%) had ARD, and its presence was statistically significantly associated with large deletions (p = 0.047). Relatedly, participants with ARD had significantly more genes deleted on chromosome 22 than participants without ARD (p = 0.013). These results could aid in the identification of individuals with PMS who are at higher risk for ARD.

DSN1 may predict poor prognosis of lower-grade glioma patients and be a potential target for immunotherapy

DSN1 has been previously found to be positively correlated with various cancers. However, the effect of DSN1 or its methylation on the prognosis, molecular characteristics, and immune cell infiltration of low-grade glioma (LGG) has not yet been studied. We obtained 1046 LGG samples from the The Cancer Genome Atlas, The Chinese Glioma Genome Atlas (CGGA) microarray, and CGGA RNA-Seq databases. Bioinformatic methods (gene set enrichment analysis (GSEA), chi-square test, multivariate), and laboratory validation were used to investigate DSN1 in LGG. The expression levels of DSN1 mRNA and protein in LGG were substantially higher than those in normal brain tissue, and their expression was negatively regulated by methylation. The survival time of patients with low expression of DSN1 and cg12601032 hypermethylation was considerably prolonged. DSN1 was a risk factor, and of good diagnostic and prognostic value for LGG. Importantly, the expression of DSN1 is related to many types of tumor-infiltrating immune cells and has a positive correlation with PDL1. DSN1 promoted the activation of multiple cancer-related pathways, such as the cell cycle. Additionally, knockdown of DSN1 substantially inhibited the proliferation and invasion of LGG cells. To the best of our knowledge, this study is the first comprehensive analysis of the mechanism of DSN1 leading to poor prognosis of LGG, which provides a new perspective for revealing the pathogenesis of LGG. DSN1 or its methylation has diagnostic value for the prognosis of glioma, and may become a new biological target of anti-tumor immunotherapy.

DCLRE1B as a novel prognostic biomarker associated with immune infiltration: a pancancer analysis

The DNA cross-link repair 1B (DCLRE1B) gene is involved in repairing cross-links between DNA strands, including those associated with Hoyeraal-Hreidarsson syndrome and congenital dyskeratosis. However, its role in tumours is not well understood. DCLRE1B expression profiles were examined in tumour tissues and normal tissues using TCGA, GTEx, and TARGET datasets. Additionally, we performed experiments with clinical melanoma samples to verify DCLRE1B expression patterns. We also performed pancancer analyses to investigate the diverse roles of DCLRE1B in the biological functions of various cancers. DCLRE1B exhibited distinct expression patterns and played crucial prognostic roles in most tumours. In particular, high expression of DCLRE1B in melanoma was significantly correlated with a poor prognosis and increased malignancy. DCLRE1B was also found to be associated with the immune landscape and various immune biomarkers and regulators. Furthermore, our analysis identified potential small molecules that could target DCLRE1B in different cancer types. The DCLRE1B gene may be involved in the development and occurrence of a variety of cancers. Additionally, DCLRE1B affects various tumour types not only by mediating DNA repair but also by shaping the differential immune microenvironment. In conclusion, our research offers fresh perspectives on the diagnosis and treatment of different types of cancers.

Regulation of sarcomere formation and function in the healthy heart requires a titin intronic enhancer

Heterozygous truncating variants in the sarcomere protein titin (TTN) are the most common genetic cause of heart failure. To understand mechanisms that regulate abundant cardiomyocyte (CM) TTN expression, we characterized highly conserved intron 1 sequences that exhibited dynamic changes in chromatin accessibility during differentiation of human CMs from induced pluripotent stem cells (hiPSC-CMs). Homozygous deletion of these sequences in mice caused embryonic lethality, whereas heterozygous mice showed an allele-specific reduction in Ttn expression. A 296 bp fragment of this element, denoted E1, was sufficient to drive expression of a reporter gene in hiPSC-CMs. Deletion of E1 downregulated TTN expression, impaired sarcomerogenesis, and decreased contractility in hiPSC-CMs. Site-directed mutagenesis of predicted binding sites of NK2 homeobox 5 (NKX2-5) and myocyte enhancer factor 2 (MEF2) within E1 abolished its transcriptional activity. In embryonic mice expressing E1 reporter gene constructs, we validated in vivo cardiac-specific activity of E1 and the requirement for NKX2-5- and MEF2-binding sequences. Moreover, isogenic hiPSC-CMs containing a rare E1 variant in the predicted MEF2-binding motif that was identified in a patient with unexplained dilated cardiomyopathy (DCM) showed reduced TTN expression. Together, these discoveries define an essential, functional enhancer that regulates TTN expression. Manipulation of this element may advance therapeutic strategies to treat DCM caused by TTN haploinsufficiency.

TFTG: A comprehensive database for human transcription factors and their targets

Transcription factors (TFs) are major contributors to gene transcription, especially in controlling cell-specific gene expression and disease occurrence and development. Uncovering the relationship between TFs and their target genes is critical to understanding the mechanism of action of TFs. With the development of high-throughput sequencing techniques, a large amount of TF-related data has accumulated, which can be used to identify their target genes. In this study, we developed TFTG (Transcription Factor and Target Genes) database (http://tf.liclab.net/TFTG), which aimed to provide a large number of available human TF-target gene resources by multiple strategies, besides performing a comprehensive functional and epigenetic annotations and regulatory analyses of TFs. We identified extensive available TF-target genes by collecting and processing TF-associated ChIP-seq datasets, perturbation RNA-seq datasets and motifs. We also obtained experimentally confirmed relationships between TF and target genes from available resources. Overall, the target genes of TFs were obtained through integrating the relevant data of various TFs as well as fourteen identification strategies. Meanwhile, TFTG was embedded with user-friendly search, analysis, browsing, downloading and visualization functions. TFTG is designed to be a convenient resource for exploring human TF-target gene regulations, which will be useful for most users in the TF and gene expression regulation research.

BICEP: Bayesian inference for rare genomic variant causality evaluation in pedigrees

Next-generation sequencing is widely applied to the investigation of pedigree data for gene discovery. However, identifying plausible disease-causing variants within a robust statistical framework is challenging. Here, we introduce BICEP: a Bayesian inference tool for rare variant causality evaluation in pedigree-based cohorts. BICEP calculates the posterior odds that a genomic variant is causal for a phenotype based on the variant cosegregation as well as a priori evidence such as deleteriousness and functional consequence. BICEP can correctly identify causal variants for phenotypes with both Mendelian and complex genetic architectures, outperforming existing methodologies. Additionally, BICEP can correctly down-weight common variants that are unlikely to be involved in phenotypic liability in the context of a pedigree, even if they have reasonable cosegregation patterns. The output metrics from BICEP allow for the quantitative comparison of variant causality within and across pedigrees, which is not possible with existing approaches.

Integrated analysis of circRNA regulation with ADARB2 enrichment in inhibitory neurons

This study investigates the regulation of circular RNAs (circRNAs) with Adenosine Deaminase RNA Specific B2 (ADARB2) enrichment specifically in inhibitory neurons. Using an integrated analysis combining high-throughput sequencing and bioinformatics approaches, we identified a group of circRNAs that are potentially enhanced by ADARB2. Our findings highlight the pivotal role of ADARB2 in circRNA synthesis within inhibitory neurons, likely through its specific binding to precursor RNAs, which facilitates circRNA biogenesis.

Expression of the non-coding RNA nc886 facilitates the development of tyrosine kinase inhibitor resistance in EGFR-mutated non-small-cell lung cancer cells

Treatment of non-small-cell lung cancer (NSCLC) patients possessing EGFR-activating mutations with tyrosine kinase inhibitors (TKIs) can confer an initial promising response. However, TKI resistance inevitably arises. Numerous TKI resistance mechanisms are identified including EGFR secondary mutations, bypass receptor tyrosine kinase (RTK) signaling, and cellular transition e.g. epithelial-mesenchymal transition (EMT). To increase the knowledge of TKI resistance we performed an epigenetic screen to identify small non-coding (nc) genes with DNA methylation alterations in HCC827 NSCLC EGFR-mutated cells with acquired TKI resistance. We analyzed Infinium Methylation EPIC 850K Array data for DNA methylation changes present in both TKI-resistant HCC827 cells with EMT and MET-amplification. Hereby, we identified that the polymorphic maternal imprinted gene nc886 (vtRNA2-1) has a decrease in promoter DNA methylation in TKI-resistant cells. This epigenetic change was associated with an increase in the expression of nc886. The induction of EMT did not affect nc886 expression. CRISPR/Cas9-mediated distortion of the nc886 sequence increased the sensitivity of HCC827 cells towards TKI. Finally, nc886 sequence distortion hindered MET RTK activation and instead was EMT the endpoint TKI resistance mechanism. In conclusion, the expression of nc886 contributes to TKI resistance in the HCC827 NSCLC cell line by supporting cell survival and selection of the endpoint TKI resistance mechanism. We propose DNA methylation and expression changes for nc886 to constitute a novel TKI resistance contributing mechanism in NSCLC.

Investigation of a Pathogenic Inversion in UNC13D and Comprehensive Analysis of Chromosomal Inversions Across Diverse Datasets

Inversions are known contributors to the pathogenesis of genetic diseases. Identifying inversions poses significant challenges, making it one of the most demanding structural variants (SVs) to detect and interpret. Recent advancements in sequencing technologies and the development of publicly available SV datasets have substantially enhanced our capability to explore inversions. However, a cross-comparison in those datasets remains unexplored. In this study, we reported a proband with familial hemophagocytic lymphohistiocytosis type-3 carrying c.1389+1G>A in trans with NC_000017.11:75576992_75829587inv disrupting UNC13D , an inversion present in 0.006345% of individuals in gnomAD(v4.0). Based on this result, we investigate the features of potentially pathogenic inversions in public datasets. 98.9% of inversions are rare in gnomAD, and they disrupt 5% of protein-coding genes associated with a phenotype in OMIM. We then conducted a comparative analysis of the datasets, including gnomAD, DGV, and 1KGP, and two recent studies from the Human Genome Structural Variation Consortium revealed common and dataset-specific inversion characteristics suggesting methodology detection biases. Next, we investigated the genetic features of inversions disrupting the protein-coding genes by classifying the intersections between them into three categories. We found that most of the protein-coding genes in OMIM disrupted by inversions are associated with autosomal recessive phenotypes regardless of categories supporting the hypothesis that inversions in trans with other variants are hidden causes of monogenic diseases. This effort aims to fill the gap in our understanding of the molecular characteristics of inversions with low frequency in the population and highlight the importance of identifying them in rare disease studies.

Copy-number dosage regulates telomere maintenance and disease-associated pathways in neuroblastoma

Telomere maintenance in neuroblastoma is linked to poor outcome and caused by either telomerase reverse transcriptase (TERT) activation or through alternative lengthening of telomeres (ALT). In contrast to TERT activation, commonly caused by genomic rearrangements or MYCN amplification, ALT is less well understood. Alterations at the ATRX locus are key drivers of ALT but only present in ∼50% of ALT tumors. To identify potential new pathways to telomere maintenance, we investigate allele-specific gene dosage effects from whole genomes and transcriptomes in 115 primary neuroblastomas. We show that copy-number dosage deregulates telomere maintenance, genomic stability, and neuronal pathways and identify upregulation of variants of histone H3 and H2A as a potential alternative pathway to ALT. We investigate the interplay between TERT activation, overexpression and copy-number dosage and reveal loss of imprinting at the RTL1 gene associated with poor clinical outcome. These results highlight the importance of gene dosage in key oncogenic mechanisms in neuroblastoma.

[Single nucleotide polymorphism heritability of non-syndromic cleft lip with or without cleft palate in Chinese population]

OBJECTIVE

To delve into the intricate relationship between common genetic variations across the entire genome and the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P).

METHODS

Utilizing summary statistics data from genome-wide association studies (GWAS), a thorough investigation to evaluate the impact of common variations on the genome were undertook. This involved assessing single nucleotide polymorphism (SNP) heritability across the entire genome, as well as within specific genomic regions. To ensure the robustness of our analysis, stringent quality control measures were applied to the GWAS summary statistics data. Criteria for inclusion encompassed the absence of missing values, a minor allele frequency ≥1%, P-values falling within the range of 0 to 1, and clear SNP strand orientation. SNP meeting these stringent criteria were then meticulously included in our analysis. The SNP heritability of NSCL/P was calculated using linkage disequilibrium score regression. Additionally, hierarchical linkage disequilibrium score regression to partition SNP heritability within coding regions, promoters, introns, enhancers, and super enhancers were employed, and the enrichment levels within different genomic regions using LDSC (v1.0.1) software were further elucidated.

RESULTS

Our study drew upon GWAS summary statistics data obtained from 806 NSCL/P trios, comprising a total of 2 418 individuals from the Chinese population. Following rigorous quality control procedures, 490 593 out of 492 993 SNP were deemed suitable for inclusion in SNP heritability calculations. The observed SNP heritability of NSCL/P was 0.55 (95%CI: 0.28-0.82). Adjusting for the elevated disease pre-valence within our sample, the SNP heritability scaled down to 0.37 (95%CI: 0.19-0.55) based on the prevalence observed in the general Chinese population. Notably, our enrichment analysis unveiled significant enrichment of SNP heritability within enhancer regions (15.70, P=0.04) and super enhancer regions (3.18, P=0.03).

CONCLUSION

Our study sheds light on the intricate interplay between common genetic variations and the risk of NSCL/P in the Chinese population. By elucidating the SNP heritability landscape across different genomic regions, we contribute valuable insights into the genetic basis of NSCL/P. The significant enrichment of SNP heritability within enhancer and super enhancer regions underscores the potential role of these regulatory elements in shaping the genetic susceptibility to NSCL/P. This paves the way for further research aimed at uncovering novel genetic pathogenic factors underlying NSCL/P pathogenesis.

Association of Unc-51-like Kinase 4 (ULK4) with the reactivity of the extended reward system in response to conditioned stimuli

OBJECTIVES

ULK4 is an established candidate gene for mental disorders and antipsychotic treatment response. We investigated the association of functional genetic variation at the ULK4 locus with the human extended dopaminergic reward system using fMRI during the performance of a well-established reward paradigm.

METHODS

Two hundred and thirty-four patients were included in this study. Association of genetic variation in the ULK4 gene with reward system functioning were determined using the Desire-Reason-Dilemma (DRD) paradigm which allows to assess brain activation in response to conditioned reward stimuli.

RESULTS

Variant prioritisation revealed the strongest functional signatures for the ULK4 variant rs17215589, coding for amino acid exchange Ala715Thr. For rs17215589 minor allele carriers, we detected increased activation responses to conditioned reward stimuli in the ventral tegmental area, nucleus accumbens and several cortical brain regions of the extended reward system.

CONCLUSIONS

Our findings provide further evidence in humans that genetic variation in ULK4 may increase the vulnerability to mental disorders, by modulating the extended reward system function. Future studies are needed to confirm the modulation of the extended reward system by ULK4 and to specify the role of this mechanism in the pathogenesis of psychiatric disorders.

Unraveling the genomic diversity and admixture history of captive tigers in the United States

Genomic studies of endangered species have primarily focused on describing diversity patterns and resolving phylogenetic relationships, with the overarching goal of informing conservation efforts. However, few studies have investigated genomic diversity housed in captive populations. For tigers (Panthera tigris), captive individuals vastly outnumber those in the wild, but their diversity remains largely unexplored. Privately owned captive tiger populations have remained an enigma in the conservation community, with some believing that these individuals are severely inbred, while others believe they may be a source of now-extinct diversity. Here, we present a large-scale genetic study of the private (non-zoo) captive tiger population in the United States, also known as "Generic" tigers. We find that the Generic tiger population has an admixture fingerprint comprising all six extant wild tiger subspecies. Of the 138 Generic individuals sequenced for the purpose of this study, no individual had ancestry from only one subspecies. We show that the Generic tiger population has a comparable amount of genetic diversity relative to most wild subspecies, few private variants, and fewer deleterious mutations. We observe inbreeding coefficients similar to wild populations, although there are some individuals within both the Generic and wild populations that are substantially inbred. Additionally, we develop a reference panel for tigers that can be used with imputation to accurately distinguish individuals and assign ancestry with ultralow coverage (0.25×) data. By providing a cost-effective alternative to whole-genome sequencing (WGS), the reference panel provides a resource to assist in tiger conservation efforts for both ex- and in situ populations.

Biophysical characterization of high-confidence, small human proteins

Significant efforts have been made to characterize the biophysical properties of proteins. Small proteins have received less attention because their annotation has historically been less reliable. However, recent improvements in sequencing, proteomics, and bioinformatics techniques have led to the high-confidence annotation of small open reading frames (smORFs) that encode for functional proteins, producing smORF-encoded proteins (SEPs). SEPs have been found to perform critical functions in several species, including humans. While significant efforts have been made to annotate SEPs, less attention has been given to the biophysical properties of these proteins. We characterized the distributions of predicted and curated biophysical properties, including sequence composition, structure, localization, function, and disease association of a conservative list of previously identified human SEPs. We found significant differences between SEPs and both larger proteins and control sets. In addition, we provide an example of how our characterization of biophysical properties can contribute to distinguishing protein-coding smORFs from noncoding ones in otherwise ambiguous cases.

The genomic landscape of 2,023 colorectal cancers

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

MTR3D-AF2: Expanding the coverage of spatially derived missense tolerance scores across the human proteome using AlphaFold2

The missense tolerance ratio (MTR) was developed as a novel approach to assess the deleteriousness of variants. Its three-dimensional successor, MTR3D, was demonstrated powerful at discriminating pathogenic from benign variants. However, its reliance on experimental structures and homologs limited its coverage of the proteome. We have now utilized AlphaFold2 models to develop MTR3D-AF2, which covers 89.31% of proteins and 85.39% of residues across the human proteome. This work has improved MTR3D's ability to distinguish clinically established pathogenic from benign variants. MTR3D-AF2 is freely available as an interactive web server at https://biosig.lab.uq.edu.au/mtr3daf2/.

Investigating copy number variants in schizophrenia pedigrees using a new consensus pipeline called PECAN

Copy number variants (CNVs) have been implicated in many human diseases, including psychiatric disorders. Whole genome sequencing offers advantages in CNV calling compared to previous array-based methods. Here we present a robust and transparent CNV calling pipeline, PECAN (PEdigree Copy number vAriaNt calling), for short-read, whole genome sequencing data, comprised of a novel combination of four calling methods and structural variant genotyping. This method is scalable and can incorporate pedigree information to retain lower-confidence CNVs that would otherwise be discarded. We have robustly benchmarked PECAN using gold-standard CNV calls for two well-established evaluation samples, NA12878 and HG002, showing that PECAN performs with high precision and recall on both datasets, outperforming another pedigree-based CNV calling pipeline. As part of this work, we provide a list of high-confidence gold standard CNVs for the NA12878 reference sample, curated from multiple studies. We applied PECAN to a collection of pedigrees multiply affected with schizophrenia and identified a rare deletion that perfectly co-segregates with schizophrenia in one of the pedigrees. The CNV overlaps the gene PITRM1, which has been implicated in a complex phenotype including ataxia, developmental delay, and schizophrenia-like episodes in affected adults.

Specificity of Key Sex Determination Genes in a Mammal with Ovotestes: The European Mole Talpa europaea

Here, for the first time, the structure of genes involved in sex determination in mammals (full Sry and partial Rspo1, Eif2s3x, and Eif2s3y) was analyzed for the European mole Talpa europaea with ovotestes in females. We confirmed male-specificity for Eif2s3y and Sry. Five exons were revealed for Rspo1 and the deep similarity with the structure of this gene in T. occidentalis was proved. The most intriguing result was obtained for the Sry gene, which, in placental mammals, initiates male development. We described two exons for this canonically single-exon gene: the first (initial) exon is only 15 bp while the second exon includes 450 bp. The exons are divided by an extended intron of about 1894 bp, including the fragment of the LINE retroposon. Moreover, in chromatogram fragments, which correspond to intron and DNA areas, flanking both exons, we revealed double peaks, similar to heterozygous nucleotide sites of autosomal genes. This may indicate the existence of two or more copies of the Sry gene. Proof of copies requires an additional in-depth study. We hypothesize that unusual structure and possible supernumerary copies of Sry may be involved in ovotestes formation.

RNA editing regulates glutamatergic synapses in the frontal cortex of a molecular subtype of Amyotrophic Lateral Sclerosis

BACKGROUND

Amyotrophic Lateral Sclerosis (ALS) is a highly heterogenous neurodegenerative disorder that primarily affects upper and lower motor neurons, affecting additional cell types and brain regions. Underlying molecular mechanisms are still elusive, in part due to disease heterogeneity. Molecular disease subtyping through integrative analyses including RNA editing profiling is a novel approach for identification of molecular networks involved in pathogenesis.

METHODS

We aimed to highlight the role of RNA editing in ALS, focusing on the frontal cortex and the prevalent molecular disease subtype (ALS-Ox), previously determined by transcriptomic profile stratification. We established global RNA editing (editome) and gene expression (transcriptome) profiles in control and ALS-Ox cases, utilizing publicly available RNA-seq data (GSE153960) and an in-house analysis pipeline. Functional annotation and pathway analyses identified molecular processes affected by RNA editing alterations. Pearson correlation analyses assessed RNA editing effects on expression. Similar analyses on additional ALS-Ox and control samples (GSE124439) were performed for verification. Targeted re-sequencing and qRT-PCR analysis targeting CACNA1C, were performed using frontal cortex tissue from ALS and control samples (n = 3 samples/group).

RESULTS

We identified reduced global RNA editing in the frontal cortex of ALS-Ox cases. Differentially edited transcripts are enriched in synapses, particularly in the glutamatergic synapse pathway. Bioinformatic analyses on additional ALS-Ox and control RNA-seq data verified these findings. We identified increased recoding at the Q621R site in the GRIK2 transcript and determined positive correlations between RNA editing and gene expression alterations in ionotropic receptor subunits GRIA2, GRIA3 and the CACNA1C transcript, which encodes the pore forming subunit of a post-synaptic L-type calcium channel. Experimental data verified RNA editing alterations and editing-expression correlation in CACNA1C, highlighting CACNA1C as a target for further study.

CONCLUSIONS

We provide evidence on the involvement of RNA editing in the frontal cortex of an ALS molecular subtype, highlighting a modulatory role mediated though recoding and gene expression regulation on glutamatergic synapse related transcripts. We report RNA editing effects in disease-related transcripts and validated editing alterations in CACNA1C. Our study provides targets for further functional studies that could shed light in underlying disease mechanisms enabling novel therapeutic approaches.

A single-cell atlas of chromatin accessibility in mouse organogenesis

Organogenesis is a highly complex and precisely regulated process. Here we profiled the chromatin accessibility in >350,000 cells derived from 13 mouse embryos at four developmental stages from embryonic day (E) 10.5 to E13.5 by SPATAC-seq in a single experiment. The resulting atlas revealed the status of 830,873 candidate cis-regulatory elements in 43 major cell types. By integrating the chromatin accessibility atlas with the previous transcriptomic dataset, we characterized cis-regulatory sequences and transcription factors associated with cell fate commitment, such as Nr5a2 in the development of gastrointestinal tract, which was preliminarily supported by the in vivo experiment in zebrafish. Finally, we integrated this atlas with the previous single-cell chromatin accessibility dataset from 13 adult mouse tissues to delineate the developmental stage-specific gene regulatory programmes within and across different cell types and identify potential molecular switches throughout lineage development. This comprehensive dataset provides a foundation for exploring transcriptional regulation in organogenesis.

Gestational exposure to environmental chemicals and epigenetic alterations in the placenta and cord blood mononuclear cells

Background

Exposure to environmental chemicals such as phthalates, phenols, and polycyclic aromatic hydrocarbons (PAHs) during pregnancy can increase the risk of adverse newborn outcomes. We explored the associations between maternal exposure to select environmental chemicals and DNA methylation in cord blood mononuclear cells (CBMC) and placental tissue (maternal and fetal sides) to identify potential mechanisms underlying these associations.

Method

This study included 75 pregnant individuals who planned to give birth at the University of Cincinnati Hospital between 2014 and 2017. Maternal urine samples during the delivery visit were collected and analyzed for 37 biomarkers of phenols (12), phthalates (13), phthalate replacements (4), and PAHs (8). Cord blood and placenta tissue (maternal and fetal sides) were also collected to measure the DNA methylation intensities using the Infinium HumanMethylation450K BeadChip. We used linear regression, adjusting for potential confounders, to assess CpG-specific methylation changes in CBMC (n = 54) and placenta [fetal (n = 67) and maternal (n = 68) sides] associated with gestational chemical exposures (29 of 37 biomarkers measured in this study). To account for multiple testing, we used a false discovery rate q-values < 0.05 and presented results by limiting results with a genomic inflation factor of 1±0.5. Additionally, gene set enrichment analysis was conducted using the Kyoto Encyclopedia of Genes and Genomics pathways.

Results

Among the 29 chemical biomarkers assessed for differential methylation, maternal concentrations of PAH metabolites (1-hydroxynaphthalene, 2-hydroxyfluorene, 4-hydroxyphenanthrene, 1-hydroxypyrene), monocarboxyisononyl phthalate, mono-3-carboxypropyl phthalate, and bisphenol A were associated with altered methylation in placenta (maternal or fetal side). Among exposure biomarkers associated with epigenetic changes, 1-hydroxynaphthalene, and mono-3-carboxypropyl phthalate were consistently associated with differential CpG methylation in the placenta. Gene enrichment analysis indicated that maternal 1-hydroxynaphthalene was associated with lipid metabolism and cellular processes of the placenta. Additionally, mono-3-carboxypropyl phthalate was associated with organismal systems and genetic information processing of the placenta.

Conclusion

Among the 29 chemical biomarkers assessed during delivery, 1-hydroxynaphthalene and mono-3-carboxypropyl phthalate were associated with DNA methylation in the placenta.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43682-024-00027-7.

A combination of long- and short-read genomics reveals frequent p-arm breakpoints within chromosome 21 complex genomic rearrangements

Purpose

Although chromosome 21 is the smallest human chromosome, it is highly relevant in the pathogenicity of both cancer and congenital diseases, including Alzheimer disease and trisomy 21 (Down syndrome). In addition, cases with rare structural variants (SVs) of chromosome 21 have been reported. These events vary in size and include large chromosomal events, such as ring chromosomes and small partial aneuploidies. The p-arm of the acrocentric chromosome 21 was devoid of reference genomic sequence in GRCh37 and GRCh38, which hampered our ability to solve genomic rearrangements and find the mechanism of formation of disease-causing SVs. We hypothesize that conserved satellite structures and segmental duplications located on the p-arm play an important role in the formation of complex SVs involving chromosome 21.

Methods

Three cases with complex chromosome 21 rearrangements were studied with a combination of short-read and long-read genome sequencing, as well as optical genome mapping. The data were aligned to the T2T-CHM13 assembly.

Results

We were able to resolve all 3 complex chromosome 21 rearrangements in which 15, 8, and 26 breakpoints were identified, respectively. By comparing the identified SV breakpoints, we were able to pinpoint a region between 21p13 and 21p12 that appears to be frequently involved in chromosome 21 rearrangements. Importantly, we observed acrocentric satellite DNA at several breakpoint junctions suggesting an important role for those elements in the formation of complex SVs.

Conclusion

Taken together, our results provide further insights into the architecture and underlying mechanisms of complex rearrangements on acrocentric chromosomes.

Systematic single-cell analysis reveals dynamic control of transposable element activity orchestrating the endothelial-to-hematopoietic transition

BACKGROUND

The endothelial-to-hematopoietic transition (EHT) process during definitive hematopoiesis is highly conserved in vertebrates. Stage-specific expression of transposable elements (TEs) has been detected during zebrafish EHT and may promote hematopoietic stem cell (HSC) formation by activating inflammatory signaling. However, little is known about how TEs contribute to the EHT process in human and mouse.

RESULTS

We reconstructed the single-cell EHT trajectories of human and mouse and resolved the dynamic expression patterns of TEs during EHT. Most TEs presented a transient co-upregulation pattern along the conserved EHT trajectories, coinciding with the temporal relaxation of epigenetic silencing systems. TE products can be sensed by multiple pattern recognition receptors, triggering inflammatory signaling to facilitate HSC emergence. Interestingly, we observed that hypoxia-related signals were enriched in cells with higher TE expression. Furthermore, we constructed the hematopoietic cis-regulatory network of accessible TEs and identified potential TE-derived enhancers that may boost the expression of specific EHT marker genes.

CONCLUSIONS

Our study provides a systematic vision of how TEs are dynamically controlled to promote the hematopoietic fate decisions through transcriptional and cis-regulatory networks, and pre-train the immunity of nascent HSCs.

Repurposing homoharringtonine for thyroid cancer treatment through TIMP1/FAK/PI3K/AKT signaling pathway

Homoharringtonine (HHT), an alkaloid isolated from Cephalotaxus, is an effective anti-leukemia agent and exhibits inhibitory effects in various solid tumors. However, the impacts of HHT treatment on thyroid cancer (TC) remain unclear. Our findings demonstrated that HHT exhibited remarkable anti-TC activity that involved inhibiting cell proliferation, invasion, and migration, as well as inducing apoptosis. Proteomics analysis revealed that the expression of the tissue inhibitor of metalloproteinase 1 (TIMP1) was downregulated in TC cells after HHT treatment. TIMP1 overexpression promoted TC progression and partially reversed the anti-TC effects of HHT, while TIMP1 downregulation inhibited TC progression and enhanced the anti-TC effects of HHT. Furthermore, TIMP1 re-expression attenuated the enhancement of anti-TC effects of HHT induced by TIMP1 knockdown. Mechanistically, HHT exerted anti-TC effects by downregulating TIMP1 expression and then inactivating the FAK/PI3K/AKT signaling pathway. Taken together, our study demonstrated that HHT could inhibit TC progression by inhibiting the TIMP1/FAK/PI3K/AKT signaling pathway.

Osteocytes support bone metastasis of melanoma cells by CXCL5

Bone metastasis is a common complication of certain cancers such as melanoma. The spreading of cancer cells into the bone is supported by changes in the bone marrow environment. The specific role of osteocytes in this process is yet to be defined. By RNA-seq and chemokines screening we show that osteocytes release the chemokine CXCL5 when they are exposed to melanoma cells. Osteocytes-mediated CXCL5 secretion enhanced the migratory and invasive behaviour of melanoma cells. When the expression of the CXCL5 receptor, CXCR2, was down-regulated in melanoma cells in vitro, we observed a significant decrease in melanoma cell migration in response to osteocytes. Furthermore, melanoma cells with down-regulated CXCR2 expression showed less bone metastasis and less bone loss in the bone metastasis model in vivo. Furthermore, when simultaneously down-regulating CXCL5 in osteocytes and CXCR2 in melanoma cells, melanoma progression was abrogated in vivo. In summary, these data suggest a significant role of osteocytes in bone metastasis of melanoma, which is mediated through the CXCL5-CXCR2 pathway.

Estrogen receptor 1 chromatin profiling in human breast tumors reveals high inter-patient heterogeneity with enrichment of risk SNPs and enhancer activity at most-conserved regions

Estrogen Receptor 1 (ESR1; also known as ERα, encoded by ESR1 gene) is the main driver and prime drug target in luminal breast cancer. ESR1 chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ESR1 chromatin action, along with its biological implications. Here, we use a large set of ESR1 ChIP-seq data from 70 ESR1+ breast cancers to explore inter-patient heterogeneity in ESR1 DNA binding to reveal a striking inter-tumor heterogeneity of ESR1 action. Of note, commonly shared ESR1 sites show the highest estrogen-driven enhancer activity and are most engaged in long-range chromatin interactions. In addition, the most commonly shared ESR1-occupied enhancers are enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ESR1 and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we can confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ESR1-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ESR1 landscape, with the most common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

Methylation-based markers for the estimation of age in African cheetah, Acinonyx jubatus

Age is a key demographic in conservation where age classes show differences in important population metrics such as morbidity and mortality. Several traits, including reproductive potential, also show senescence with ageing. Thus, the ability to estimate age of individuals in a population is critical in understanding the current structure as well as their future fitness. Many methods exist to determine age in wildlife, with most using morphological features that show inherent variability with age. These methods require significant expertise and become less accurate in adult age classes, often the most critical groups to model. Molecular methods have been applied to measuring key population attributes, and more recently epigenetic attributes such as methylation have been explored as biomarkers for age. There are, however, several factors such as permits, sample sovereignty, and costs that may preclude the use of extant methods in a conservation context. This study explored the utility of measuring age-related changes in methylation in candidate genes using mass array technology. Novel methods are described for using gene orthologues to identify and assay regions for differential methylation. To illustrate the potential application, African cheetah was used as a case study. Correlation analyses identified six methylation sites with an age relationship, used to develop a model with sufficient predictive power for most conservation contexts. This model was more accurate than previous attempts using PCR and performed similarly to candidate gene studies in other mammal species. Mass array presents an accurate and cost-effective method for age estimation in wildlife of conservation concern.

Solute carrier family 16 member 1 as a potential prognostic factor for pancreatic ductal adenocarcinoma and its influence on tumor immunity

Background

Solute carrier family 16 member 1 (SLC16A1) serves as a biomarker in numerous types of cancer. Tumor immune infiltration has drawn increasing attention in cancer progression and treatment. The objective of our study was to explore the association between SLC16A1 and the tumor immune microenvironment in pancreatic ductal adenocarcinoma (PDAC).

Methods

Data were obtained from The Cancer Genome Atlas. The xCell web tool was used to calculate the proportion of immune cells according to SLC16A1 expression. To further explore the mechanism of SLC16A1, immunity-related genes were screened from differentially expressed genes through weighted gene coexpression network analysis, examined via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, and filtrated using univariate Cox regression and least absolute shrinkage and selection operator regression model combined correlation analysis (P<0.05). Next, CIBERSORT was used to analyze the correlation between immune cells and five important genes. SLC16A1 expression and its clinical role in pancreatic cancer was clarified via immunohistochemical staining experiments. Finally, the effects of SLC16A1 on the results of cancer immunity were evaluated by in vitro experiments.

Results

SLC16A1 was overexpressed in PDAC tissues and could be an independent prognostic factor. SLC16A1 was significantly negatively correlated with overall survival and suppressed the tumor immunity of PDAC. In clinic, SLC16A1 expression was significantly positively correlated with tumor progression and poor prognosis. We also found that SLC16A1 could suppress the antitumor ability of CD8+ T cells.

Conclusions

SLC16A1 is a biomarker for the prognosis of PDAC and can influence the immune environment of PDAC. These findings provide new insights into the treatment of PDAC.

Genome- and epigenome-wide association studies identify susceptibility of CpG sites and regions for metabolic syndrome in a Korean population

BACKGROUND

While multiple studies have investigated the relationship between metabolic syndrome (MetS) and its related traits (fasting glucose, triglyceride, HDL cholesterol, blood pressure, waist circumference) and DNA methylation, our understanding of the epigenetic mechanisms in MetS remains limited. Therefore, we performed an epigenome-wide meta-analysis of blood DNA methylation to identify differentially methylated probes (DMPs) and differentially methylated regions (DMRs) associated with MetS and its components using two independent cohorts comprising a total of 2,334 participants. We also investigated the specific genetic effects on DNA methylation, identified methylation quantitative trait loci (meQTLs) through genome-wide association studies and further utilized Mendelian randomization (MR) to assess how these meQTLs subsequently influence MetS status.

RESULTS

We identified 40 DMPs and 27 DMRs that are significantly associated with MetS. In addition, we identified many novel DMPs and DMRs underlying inflammatory and steroid hormonal processes. The most significant associations were observed in 3 DMPs (cg19693031, cg26974062, cg02988288) and a DMR (chr1:145440444-145441553) at the TXNIP, which are involved in lipid metabolism. These CpG sites were identified as coregulators of DNA methylation in MetS, TG and FAG levels. We identified a total of 144 cis-meQTLs, out of which only 13 were found to be associated with DMPs for MetS. Among these, we confirmed the identified causal mediators of genetic effects at CpG sites cg01881899 at ABCG1 and cg00021659 at the TANK genes for MetS.

CONCLUSIONS

This study observed whether specific CpGs and methylated regions act independently or are influenced by genetic effects for MetS and its components in the Korean population. These associations between the identified DNA methylation and MetS, along with its individual components, may serve as promising targets for the development of preventive interventions for MetS.

Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease

BACKGROUND

X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique.

RESULTS

In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis.

CONCLUSIONS

This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Biophysical characterization of high-confidence, small human proteins

Significant efforts have been made to characterize the biophysical properties of proteins. Small proteins have received less attention because their annotation has historically been less reliable. However, recent improvements in sequencing, proteomics, and bioinformatics techniques have led to the high-confidence annotation of small open reading frames (smORFs) that encode for functional proteins, producing smORF-encoded proteins (SEPs). SEPs have been found to perform critical functions in several species, including humans. While significant efforts have been made to annotate SEPs, less attention has been given to the biophysical properties of these proteins. We characterized the distributions of predicted and curated biophysical properties, including sequence composition, structure, localization, function, and disease association of a conservative list of previously identified human SEPs. We found significant differences between SEPs and both larger proteins and control sets. Additionally, we provide an example of how our characterization of biophysical properties can contribute to distinguishing protein-coding smORFs from non-coding ones in otherwise ambiguous cases.

A paradigm shift in cancer research based on integrative multi-omics approaches: glutaminase serves as a pioneering cuproptosis-related gene in pan-cancer

BACKGROUND

Cuproptosis is a newly identified form of unprogrammed cell death. As a pivotal metabolic regulator, glutaminase (GLS) has recently been discovered to be linked to cuproptosis. Despite this discovery, the oncogenic functions and mechanisms of GLS in various cancers are still not fully understood.

METHODS

In this study, a comprehensive omics analysis was performed to investigate the differential expression levels, diagnostic and prognostic potential, correlation with tumor immune infiltration, genetic alterations, and drug sensitivity of GLS across multiple malignancies.

RESULTS

Our findings revealed unique expression patterns of GLS across various cancer types and molecular subtypes of carcinomas, underscoring its pivotal role primarily in energy and nutrition metabolism. Additionally, GLS showed remarkable diagnostic and prognostic performance in specific cancers, suggesting its potential as a promising biomarker for cancer detection and prognosis. Furthermore, we focused on uterine corpus endometrial carcinoma (UCEC) and developed a novel prognostic model associated with GLS, indicating a close correlation between GLS and UCEC. Moreover, our exploration into immune infiltration, genetic heterogeneity, tumor stemness, and drug sensitivity provided novel insights and directions for future research and laid the foundation for high-quality verification.

CONCLUSION

Collectively, our study is the first comprehensive investigation of the biological and clinical significance of GLS in pan-cancer. In our study, GLS was identified as a promising biomarker for UCEC, providing valuable evidence and a potential target for anti-tumor therapy. Overall, our findings shed light on the multifaceted functions of GLS in cancer and offer new avenues for further research.

The prognostic biomarker TPGS2 is correlated with immune infiltrates in pan-cancer: a bioinformatics analysis

Background

Tubulin polyglutamylase complex subunit 2 (TPGS2) is an element of the neuronal polyglutamylase complex that plays a role in the post-translational addition of glutamate residues to C-terminal tubulin tails. Recent research has shown that TPGS2 is associated with some tumors, but the roles of TPGS2 in tumor immunity remain unclear.

Methods

The research data were mainly sourced from The Cancer Genome Atlas. The data were analyzed to identify potential correlations between TPGS2 expression and survival, gene alterations, the tumor mutational burden (TMB), microsatellite instability (MSI), immune infiltration, and various immune-related genes across various cancers. The Wilcoxon rank-sum test was used to identify the significance. A log-rank test and univariate Cox regression analysis were performed to assess the survival state of the patients. Spearman's correlation coefficients were used to show the correlations.

Results

TPGS2 exhibited abnormal expression patterns in most types of cancers, and has promising prognostic potential in adrenocortical carcinoma and liver hepatocellular carcinoma. Further, TPGS2 expression was significantly correlated with molecular and immune subtypes. Moreover, the single-cell analyses showed that the expression of TPGS2 was associated with the cell cycle, metastasis, invasion, inflammation, and DNA damage. In addition, the immune cell infiltration analysis and gene-set enrichment analysis demonstrated that a variety of immune cells and immune processes were associated with TPGS2 expression in various cancers. Further, immune regulators, including immunoinhibitors, immunostimulators, the major histocompatibility complex, chemokines, and chemokine receptors, were correlated with TPGS2 expression in different cancer types. Finally, the TMB and MSI, which have been identified as powerful predictors of immunotherapy, were shown to be correlated with the expression of TPGS2 across human cancers.

Conclusions

TPGS2 is aberrantly expressed in most cancer tissues and might be associated with immune cell infiltration in the tumor microenvironment. TPGS2 could serve not only as a biomarker for predicting clinical outcomes, but also as a promising biomarker for evaluating and developing new approaches to immunotherapy in many types of cancers, especially colon adenocarcinoma and stomach adenocarcinoma.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

Self-supervised learning on millions of primary RNA sequences from 72 vertebrates improves sequence-based RNA splicing prediction

Language models pretrained by self-supervised learning (SSL) have been widely utilized to study protein sequences, while few models were developed for genomic sequences and were limited to single species. Due to the lack of genomes from different species, these models cannot effectively leverage evolutionary information. In this study, we have developed SpliceBERT, a language model pretrained on primary ribonucleic acids (RNA) sequences from 72 vertebrates by masked language modeling, and applied it to sequence-based modeling of RNA splicing. Pretraining SpliceBERT on diverse species enables effective identification of evolutionarily conserved elements. Meanwhile, the learned hidden states and attention weights can characterize the biological properties of splice sites. As a result, SpliceBERT was shown effective on several downstream tasks: zero-shot prediction of variant effects on splicing, prediction of branchpoints in humans, and cross-species prediction of splice sites. Our study highlighted the importance of pretraining genomic language models on a diverse range of species and suggested that SSL is a promising approach to enhance our understanding of the regulatory logic underlying genomic sequences.

Comparative study on chromatin loop callers using Hi-C data reveals their effectiveness

BACKGROUND

Chromosome is one of the most fundamental part of cell biology where DNA holds the hierarchical information. DNA compacts its size by forming loops, and these regions house various protein particles, including CTCF, SMC3, H3 histone. Numerous sequencing methods, such as Hi-C, ChIP-seq, and Micro-C, have been developed to investigate these properties. Utilizing these data, scientists have developed a variety of loop prediction techniques that have greatly improved their methods for characterizing loop prediction and related aspects.

RESULTS

In this study, we categorized 22 loop calling methods and conducted a comprehensive study of 11 of them. Additionally, we have provided detailed insights into the methodologies underlying these algorithms for loop detection, categorizing them into five distinct groups based on their fundamental approaches. Furthermore, we have included critical information such as resolution, input and output formats, and parameters. For this analysis, we utilized the GM12878 Hi-C datasets at 5 KB, 10 KB, 100 KB and 250 KB resolutions. Our evaluation criteria encompassed various factors, including memory usages, running time, sequencing depth, and recovery of protein-specific sites such as CTCF, H3K27ac, and RNAPII.

CONCLUSION

This analysis offers insights into the loop detection processes of each method, along with the strengths and weaknesses of each, enabling readers to effectively choose suitable methods for their datasets. We evaluate the capabilities of these tools and introduce a novel Biological, Consistency, and Computational robustness score ( B C C score ) to measure their overall robustness ensuring a comprehensive evaluation of their performance.

4-Hydroxyphenylpyruvate Dioxygenase-Like predicts the prognosis and the immunotherapy response of cancers: a pan-cancer analysis

The 4-Hydroxyphenylpyruvate Dioxygenase-Like (HPDL) protein plays a crucial role in safeguarding cells from oxidative stress by orchestrating metabolic reprogramming. New research suggests that HPDL is considerably increased in pancreatic ductal adenocarcinoma, although its impact on cancer immunotherapy is still unclear. Pancancer transcriptional data were obtained from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression datasets. The cBioPortal webtool was utilized to examine genomic changes in different cancer types. The prognostic significance of HPDL in pancancer was evaluated using univariate Cox regression analysis. Extensive utilization of the CTRP and PRISM databases was performed to forecast potential medications that specifically target HPDL in LUAD. In summary, studies were conducted to evaluate the impact of HPDL on the proliferation and movement of LUAD cells using loss-of-function experiments. HPDL is expressed excessively in a wide variety of cancer types, indicating its prognostic and predictive value. Moreover, we emphasized the strong correlation between HPDL and indicators of immune stimulation, infiltration of immune cells, and expression of immunoregulators. The remarkable finding of the HPDL was its capacity to precisely anticipate responses to cancer therapies using anti-PDL1 and anti-PD1 antibodies among individuals. Moreover, HPDL can function as a predictive marker for specific inhibitors in instances of cancer. Suppression of HPDL resulted in reduced growth and movement of LUAD cells. To summarize, our results suggest that HPDL acts as a prospective predictor of outcomes and a positive indication of response to immunotherapy in patients undergoing treatment with immune checkpoint inhibitors (ICIs).

Genetics of kidney disorders in Phelan-McDermid syndrome: evidence from 357 registry participants

BACKGROUND

Phelan-McDermid syndrome (PMS) is a rare genetic disorder caused by SHANK3 pathogenic variants or chromosomal rearrangements affecting the chromosome 22q13 region. Previous research found that kidney disorders, primarily congenital anomalies of the kidney and urinary tract, are common in people with PMS, yet research into candidate genes has been hampered by small study sizes and lack of attention to these problems.

METHODS

We used a cohort of 357 people from the Phelan-McDermid Syndrome Foundation International Registry to investigate the prevalence of kidney disorders in PMS using a cross-sectional design and to identify 22q13 genes contributing to these disorders.

RESULTS

Kidney disorders reported included vesicoureteral reflux (n = 37), hydronephrosis (n = 36), dysplastic kidneys (n = 19), increased kidney size (n = 19), polycystic kidneys (15 cases), and kidney stones (n = 4). Out of 315 subjects with a 22q13 deletion, 101 (32%) had at least one kidney disorder, while only one out of 42 (2%) individuals with a SHANK3 pathogenic variant had a kidney disorder (increased kidney size). We identified two genomic regions that were significantly associated with having a kidney disorder with the peak associations observed near positions approximately 5 Mb and 400 Kb from the telomere.

CONCLUSIONS

The candidate genes for kidney disorders include FBLN1, WNT7B, UPK3A, CELSR1, and PLXNB2. This study demonstrates the utility of patient registries for uncovering genetic contributions to rare diseases. Future work should focus on functional studies for these genes to assess their potential pathogenic contribution to the different subsets of kidney disorders.

Convergence of the dysregulated regulome in schizophrenia with polygenic risk and evolutionarily constrained enhancers

Enhancers play an essential role in the etiology of schizophrenia; however, the dysregulation of enhancer activity and its impact on the regulome in schizophrenia remains understudied. To address this gap in our knowledge, we assessed enhancer and gene expression in 1,382 brain samples comprising cases with schizophrenia and unaffected controls. Dysregulation of enhancer expression was concordant with changes in gene expression, and was more closely associated with schizophrenia polygenic risk, suggesting that enhancer dysregulation is proximal to the genetic etiology of the disease. Modeling the shared variance of cis-coordinated genes and enhancers revealed a gene regulatory program that was highly associated with genetic vulnerability to schizophrenia. By integrating coordinated factors with evolutionary constraints, we found that enhancers acquired during human evolution are more likely to regulate genes that are implicated in neuropsychiatric disorders and, thus, hold potential as therapeutic targets. Our analysis provides a systematic view of regulome dysregulation in schizophrenia and highlights its convergence with schizophrenia polygenic risk and human-gained enhancers.