The UCSC Genome Browser database: 2019 update

DPM2 serve as novel oncogene and prognostic marker transactivated by ESR1 in breast cancer

Breast cancer (BRCA) is the most common malignancies worldwide with increasing rate. Dolichol phosphate mannose synthase (DPMS) is a critical mannosyltransferase involved in the posttranslational modification of proteins. At present, there is limited knowledge regarding the function of DPMS in breast cancer. In this study, silica analysis in multiple datasets found that dolichyl-phosphate mannosyltransferase subunit 2 (DPM2) is an unfavorable prognostic marker, suggesting its oncogenic role. Cell counting kit-8 and apoptosis assays show that DPM2-silenced cancer cells exhibit decreased growth potential and enhanced cell death rate. Further, transwell and wound healing assays show reduced invasion and migration capabilities in DPM2 knockdown groups, xenograft nude mice model demonstrated smaller tumor volume in DPM2 silenced BC cells. Then, the underlying downstream mechanism of DPM2 in BC was predicted and analyzed, highlighting classical tumorigenic pathways like JAK/STAT signaling pathway and oxidative phosphorylation activated in the cancer group. Finally, ChIP-seq analysis, expression correlation analysis, inhibitor treatment, and dual luciferase assays show that DPM2 is transcriptionally activated by estrogen receptor1 (ESR1). The results show that high expression of DPM2 mRNA is significantly correlated with shorter overall survival (OS) and disease-free survival (DFS) in breast cancer patients, and in vitro knockdown of DPM2 can significantly inhibit the malignant phenotypes of cells, including proliferation, invasion, migration, and apoptosis. These results suggest that DPM2 may play an important role in breast cancer. Altogether, we first uncovered the tumorigenic and prognostic role of DPM2 in breast cancer, cellular assays, and bioinformatics analysis highlighted DPM2 as oncogene via inhibited cancer-related signaling pathways in breast cancer. Besides, DPM2 is transcriptionally activated by ESR1, the signaling axis of ESR1/DPM2 provides a new strategy for BC-targeted therapy.

Hypoxia-induced immortalization of primary cells depends on Tfcp2L1 expression

Cellular senescence is a stress response mechanism that induces proliferative arrest. Hypoxia can bypass senescence and extend the lifespan of primary cells, mainly by decreasing oxidative damage. However, how hypoxia promotes these effects prior to malignant transformation is unknown. Here we observed that the lifespan of mouse embryonic fibroblasts (MEFs) is increased when they are cultured in hypoxia by reducing the expression of p16INK4a, p15INK4b and p21Cip1. We found that proliferating MEFs in hypoxia overexpress Tfcp2l1, which is a main regulator of pluripotency and self-renewal in embryonic stem cells, as well as stemness genes including Oct3/4, Sox2 and Nanog. Tfcp2l1 expression is lost during culture in normoxia, and its expression in hypoxia is regulated by Hif1α. Consistently, its overexpression in hypoxic levels increases the lifespan of MEFs and promotes the overexpression of stemness genes. ATAC-seq and Chip-seq experiments showed that Tfcp2l1 regulates genes that control proliferation and stemness such as Sox2, Sox9, Jarid2 and Ezh2. Additionally, Tfcp2l1 can replicate the hypoxic effect of increasing cellular reprogramming. Altogether, our data suggest that the activation of Tfcp2l1 by hypoxia contributes to immortalization prior to malignant transformation, facilitating tumorigenesis and dedifferentiation by regulating Sox2, Sox9, and Jarid2.

Essential role of PLD2 in hypoxia-induced stemness and therapy resistance in ovarian tumors

BACKGROUND

Hypoxia in solid tumors is an important source of chemoresistance that can determine poor patient prognosis. Such chemoresistance relies on the presence of cancer stem cells (CSCs), and hypoxia promotes their generation through transcriptional activation by HIF transcription factors.

METHODS

We used ovarian cancer (OC) cell lines, xenograft models, OC patient samples, transcriptional databases, induced pluripotent stem cells (iPSCs) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq).

RESULTS

Here, we show that hypoxia induces CSC formation and chemoresistance in ovarian cancer through transcriptional activation of the PLD2 gene. Mechanistically, HIF-1α activates PLD2 transcription through hypoxia response elements, and both hypoxia and PLD2 overexpression lead to increased accessibility around stemness genes, detected by ATAC-seq, at sites bound by AP-1 transcription factors. This in turn provokes a rewiring of stemness genes, including the overexpression of SOX2, SOX9 or NOTCH1. PLD2 overexpression also leads to decreased patient survival, enhanced tumor growth and CSC formation, and increased iPSCs reprograming, confirming its role in dedifferentiation to a stem-like phenotype. Importantly, hypoxia-induced stemness is dependent on PLD2 expression, demonstrating that PLD2 is a major determinant of de-differentiation of ovarian cancer cells to stem-like cells in hypoxic conditions. Finally, we demonstrate that high PLD2 expression increases chemoresistance to cisplatin and carboplatin treatments, both in vitro and in vivo, while its pharmacological inhibition restores sensitivity.

CONCLUSIONS

Altogether, our work highlights the importance of the HIF-1α-PLD2 axis for CSC generation and chemoresistance in OC and proposes an alternative treatment for patients with high PLD2 expression.

Dissecting the sequence and structural determinants guiding m6A deposition and evolution via inter- and intra-species hybrids

BACKGROUND

N6-methyladenosine (m6A) is the most abundant mRNA modification, and controls mRNA stability. m6A distribution varies considerably between and within species. Yet, it is unclear to what extent this variability is driven by changes in genetic sequences ('cis') or cellular environments ('trans') and via which mechanisms.

RESULTS

Here we dissect the determinants governing RNA methylation via interspecies and intraspecies hybrids in yeast and mammalian systems, coupled with massively parallel reporter assays and m6A-QTL reanalysis. We find that m6A evolution and variability is driven primarily in 'cis', via two mechanisms: (1) variations altering m6A consensus motifs, and (2) variation impacting mRNA secondary structure. We establish that mutations impacting RNA structure - even when distant from an m6A consensus motif - causally dictate methylation propensity. Finally, we demonstrate that allele-specific differences in m6A levels lead to allele-specific changes in gene expression.

CONCLUSIONS

Our findings define the determinants governing m6A evolution and diversity and characterize the consequences thereof on gene expression regulation.

Zika virus co-opts microRNA networks to persist in placental niches detected by spatial transcriptomics

BACKGROUND

Zika virus congenital infection evades double-stranded RNA detection and may persist in the placenta for the duration of pregnancy without accompanying overt histopathologic inflammation. Understanding how viruses can persist and replicate in the placenta without causing overt cellular or tissue damage is fundamental to deciphering mechanisms of maternal-fetal vertical transmission.

OBJECTIVE

Placenta-specific microRNAs are believed to be a tenet of viral resistance at the maternal-fetal interface. We aimed to test the hypothesis that the Zika virus functionally disrupts placental microRNAs, enabling viral persistence and fetal pathogenesis.

STUDY DESIGN

To test this hypothesis, we used orthogonal approaches in human and murine experimental models. In primary human trophoblast cultures (n=5 donor placentae), we performed Argonaute high-throughput sequencing ultraviolet-crosslinking and immunoprecipitation to identify any significant alterations in the functional loading of microRNAs and their targets onto the RNA-induced silencing complex. Trophoblasts from same-donors were split and infected with a contemporary first-passage Zika virus strain HN16 (multiplicity of infection=1 plaque forming unit per cell) or mock infected. To functionally cross-validate microRNA-messenger RNA interactions, we compared our Argonaute high-throughput sequencing ultraviolet-crosslinking and immunoprecipitation results with an independent analysis of published bulk RNA-sequencing data from human placental disk specimens (n=3 subjects; Zika virus positive in first, second, or third trimester, CD45- cells sorted by flow cytometry) and compared it with uninfected controls (n=2 subjects). To investigate the importance of these microRNA and RNA interference networks in Zika virus pathogenesis, we used a gnotobiotic mouse model uniquely susceptible to the Zika virus. We evaluated if small-molecule enhancement of microRNA and RNA interference pathways with enoxacin influenced Zika virus pathogenesis (n=20 dams total yielding 187 fetal specimens). Lastly, placentae (n=14 total) from this mouse model were analyzed with Visium spatial transcriptomics (9743 spatial transcriptomes) to identify potential Zika virus-associated alterations in immune microenvironments.

RESULTS

We found that Zika virus infection of primary human trophoblast cells led to an unexpected disruption of placental microRNA regulation networks. When compared with uninfected controls, Zika virus-infected placentae had significantly altered SLC12A8, SDK1, and VLDLR RNA-induced silencing complex loading and transcript levels (-22; adjusted P value <.05; Wald-test with false discovery rate correction q<0.05). In silico microRNA target analyses revealed that 26 of 119 transcripts (22%) in the transforming growth factor-β signaling pathway were targeted by microRNAs that were found to be dysregulated following Zika virus infection in trophoblasts. In gnotobiotic mice, relative to mock controls, Zika virus-associated fetal pathogenesis included fetal growth restriction (P=.036) and viral persistence in placental tissue (P=.011). Moreover, spatial transcriptomics of murine placentae revealed that Zika virus-specific placental niches were defined by significant up-regulation of complement cascade components and coordinated changes in transforming growth factor-β gene expression. Finally, treatment of Zika virus-infected mice with enoxacin abolished placental Zika virus persistence, rescued the associated fetal growth restriction, and the Zika virus-associated transcriptional changes in placental immune microenvironments were no longer observed.

CONCLUSION

These results collectively suggest that (1) Zika virus infection and persistence is associated with functionally perturbed microRNA and RNA interference pathways specifically related to immune regulation in placental microenvironments and (2) enhancement of placental microRNA and RNA interference pathways in mice rescued Zika virus-associated pathogenesis, specifically persistence of viral transcripts in placental microenvironments and fetal growth restriction.

Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

Disparate macrophage responses are linked to infection outcome of Hantan virus in humans or rodents

Hantaan virus (HTNV) is asymptomatically carried by rodents, yet causes lethal hemorrhagic fever with renal syndrome in humans, the underlying mechanisms of which remain to be elucidated. Here, we show that differential macrophage responses may determine disparate infection outcomes. In mice, late-phase inactivation of inflammatory macrophage prevents cytokine storm syndrome that usually occurs in HTNV-infected patients. This is attained by elaborate crosstalk between Notch and NF-κB pathways. Mechanistically, Notch receptors activated by HTNV enhance NF-κB signaling by recruiting IKKβ and p65, promoting inflammatory macrophage polarization in both species. However, in mice rather than humans, Notch-mediated inflammation is timely restrained by a series of murine-specific long noncoding RNAs transcribed by the Notch pathway in a negative feedback manner. Among them, the lnc-ip65 detaches p65 from the Notch receptor and inhibits p65 phosphorylation, rewiring macrophages from the pro-inflammation to the pro-resolution phenotype. Genetic ablation of lnc-ip65 leads to destructive HTNV infection in mice. Thus, our findings reveal an immune-braking function of murine noncoding RNAs, offering a special therapeutic strategy for HTNV infection.

Candidate biomarkers of EV-microRNA in detecting REM sleep behavior disorder and Parkinson's disease

Parkinson's disease (PD) lacks reliable, non-invasive biomarker tests for early intervention and management. Thus, a minimally invasive test for the early detection and monitoring of PD and REM sleep behavior disorder (iRBD) is a highly unmet need for developing drugs and planning patient care. Extracellular vehicles (EVs) are found in a wide variety of biofluids, including plasma. EV-mediated functional transfer of microRNAs (miRNAs) may be viable candidates as biomarkers for PD and iRBD. Next-generation sequencing (NGS) of EV-derived small RNAs was performed in 60 normal controls, 56 iRBD patients and 53 PD patients to profile small non-coding RNAs (sncRNAs). Moreover, prospective follow-up was performed for these 56 iRBD patients for an average of 3.3 years. Full-scale miRNA profiles of plasma EVs were evaluated by machine-learning methods. After optimizing the library construction method for low RNA inputs (named EVsmall-seq), we built a machine learning algorithm that identified diagnostic miRNA signatures for distinguishing iRBD patients (AUC 0.969) and PD patients (AUC 0.916) from healthy individuals; and PD patients (AUC 0.929) from iRBD patients. We illustrated all the possible expression patterns across healthy-iRBD-PD hierarchy. We also showed 20 examples of miRNAs with consistently increasing or decreasing expression levels from controls to iRBD to PD. In addition, four miRNAs were found to be correlated with iRBD conversion. Distinct characteristics of the miRNA profiles among normal, iRBD and PD samples were discovered, which provides a panel of promising biomarkers for the identification of PD patients and those in the prodromal stage iRBD.

Integration of RRBS and RNA-seq unravels the regulatory role of DNMT3A in porcine Sertoli cell proliferation

DNMT3A participates in de novo methylation, yet its impact on the proliferation of testicular Sertoli cells remains unclear. Development-specific methylation has been proven to be associated with cellular development. Therefore, in this study, we simulated DNMT3A expression pattern during testicular development by DNMT3A interference. Then, RRBS and RNA-seq were used to decipher DNMT3A regulatory mechanisms on Sertoli cell proliferation. Immunofluorescence staining revealed the expression of DNMT3A in the Sertoli cells of the prepubertal testis. DNMT3A was demonstrated to inhibit the cell cycle and proliferation of Sertoli cells, while promoting cell apoptosis. After transfected with DNMT3A interference, a total of 560 DEGs and 2,091 DMGs produced by DNMT3A interference were identified between two treated groups, respectively. Integrating the results from RRBS and RNA-seq, the overlapping genes between DMGs and DEGs were found to be enriched in the Gene Ontology (GO) terms related to cellular development and the Apelin signaling pathway. The present study demonstrated the impact of DNMT3A on the proliferation of porcine testicular Sertoli cells, suggesting that DNMT3A primarily acts through the Apelin signaling pathway. These findings provide valuable insights into how DNMT3A influences testicular development and health, offering new perspectives.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

RIP-PEN-seq identifies a class of kink-turn RNAs as splicing regulators

A kink-turn (K-turn) is a three-dimensional RNA structure that exists in all three primary phylogenetic domains. In this study, we developed the RIP-PEN-seq method to identify the full-length sequences of RNAs bound by the K-turn binding protein 15.5K and discovered a previously uncharacterized class of RNAs with backward K-turn motifs (bktRNAs) in humans and mice. All bktRNAs share two consensus sequence motifs at their fixed terminal position and have complex folding properties, expression and evolution patterns. We found that a highly conserved bktRNA1 guides the methyltransferase fibrillarin to install RNA methylation of U12 small nuclear RNA in humans. Depletion of bktRNA1 causes global splicing dysregulation of U12-type introns by impairing the recruitment of ZCRB1 to the minor spliceosome. Most bktRNAs regulate the splicing of local introns by interacting with the 15.5K protein. Taken together, our findings characterize a class of small RNAs and uncover another layer of gene expression regulation that involves crosstalk among bktRNAs, RNA splicing and RNA methylation.

Reprogramming of VEGF-mediated extracellular matrix changes through autocrine signaling

Vascular endothelial growth factor (VEGF) plays key roles in angiogenesis, vasculogenesis, and wound healing. In cancers, including triple negative breast cancer (TNBC), VEGF has been associated with increased invasion and metastasis, processes that require cancer cells to traverse through the extracellular matrix (ECM) and establish angiogenesis at distant sites. To further understand the role of VEGF in modifying the ECM, we characterized VEGF-mediated changes in the ECM of tumors derived from TNBC MDA-MB-231 cells engineered to overexpress VEGF. We established that increased VEGF expression by these cells resulted in tumors with reduced collagen 1 (Col1) fibers, fibronectin, and hyaluronan. Molecular characterization of tumors identified an increase of MMP1, uPAR, and LOX, and a decrease of MMP2, and ADAMTS1. α-SMA, a marker of cancer associated fibroblasts (CAFs), increased, and FAP-α, a marker of a subset of CAFs associated with immune suppression, decreased with VEGF overexpression. Analysis of human data from The Cancer Genome Atlas Program confirmed mRNA differences for several molecules when comparing TNBC with high and low VEGF expression. We additionally characterized enzymatic changes induced by VEGF overexpression in three different cancer cell lines that clearly identified autocrine-mediated changes, specifically uPAR, in these enzymes. Unlike the increase of Col1 fibers and fibronectin mediated by VEGF during wound healing, in the TNBC model, VEGF significantly reduced key protein components of the ECM. These results further expand our understanding of the role of VEGF in cancer progression and identify potential ECM-related targets to disrupt this progression.

New insights into COL26A1 in thyroid carcinoma: prognostic prediction, functional characterization, immunological drug target and ceRNA network

Background

Thyroid carcinoma (THCA) is one of the most commonly diagnosed malignancies. Collagen is the main component in extracellular matrix. Rising studies have determined the oncogenic effect of collagen in cancer progression, which is intriguing to be further explored. Collagen type XXVI alpha 1 chain (COL26A1) is a newly discovered collagen subtype, functions of which still remain poorly demonstrated in THCA.

Methods

Based on the transcriptome data from The Cancer Genome Atlas (TCGA) and other public databases, we conducted investigations of COL26A1 in THCA with respects to diagnostic/prognostic prediction, functional characterization, immune infiltration, chemical drug target and non-coding RNA regulatory network. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) and western blot were used to verify the expression of COL26A1 in THCA.

Results

COL26A1 was significantly upregulated in THCA, and the high COL26A1 expression inferred poor prognosis [hazard ratio (HR) =4.76; 95% confidence interval (CI): 1.36-16.73; P=0.015]. The diagnostic area under the curve (AUC) of COL26A1 achieved 0.736 (95% CI: 0.669-0.802). COL26A1 was also identified as an independent prognostic predictor for THCA (HR =3.928; 95% CI: 3.716-4.151; P<0.001). Besides, logistic regression analysis indicated that age >45 years [odds ratio (OR) =1.532; 95% CI: 1.081-2.176; P=0.017], pathological stage III (OR =2.055; 95% CI: 1.314-3.184; P=0.001), tall cell subtype (OR =5.533; 95% CI: 2.420-14.957; P<0.001), residual tumor R1 (OR =1.844; 95% CI: 1.035-3.365; P=0.041) and extrathyroidal extension (OR =1.800; 95% CI: 1.225-1.660; P=0.003) were risk factors associated with high COL26A1 expression in THCA. Functional characterizations implied that COL26A1 was associated with immunological processes and oncogenic signaling pathways. High COL26A1 expression was accompanied by more abundant infiltration of immune cells and higher stromal/immune score. In addition, most immune checkpoints were significantly positively co-expressed with COL26A1, including PD-1, PD-L1 and CTLA4. Drugs including trichloroethylene, acetamide and thioacetamide etc. that can decrease the expression of COL26A1 were also identified. The predicted long noncoding RNA (lncRNA)-microRNA (miRNA)-COL26A1 regulatory axes were successfully deciphered. qRT-PCR and western blot verified the upregulation of COL26A1 in THCA.

Conclusions

Our work has primarily appraised COL26A1 as a promising biomarker for diagnosis/prognosis and immuno/therapeutic target in THCA.

A naturally occurring variant of MBD4 causes maternal germline hypermutation in primates

As part of an ongoing genome sequencing project at the Oregon National Primate Research Center, we identified a rhesus macaque with a rare homozygous frameshift mutation in the gene methyl-CpG binding domain 4, DNA glycosylase (MBD4). MBD4 is responsible for the repair of C > T deamination mutations at CpG dinucleotides and has been linked to somatic hypermutation and cancer predisposition in humans. We show here that MBD4-associated hypermutation also affects the germline: The six offspring of the MBD4-null dam have a fourfold to sixfold increase in de novo mutation burden. This excess burden was predominantly C > T mutations at CpG dinucleotides consistent with MBD4 loss of function in the dam. There was also a significant excess of C > T at CpA sites, indicating an important, unappreciated role for MBD4 to repair deamination in CpA contexts. The MBD4-null dam developed sustained eosinophilia later in life, but we saw no other signs of neoplastic processes associated with MBD4 loss of function in humans nor any obvious disease in the hypermutated offspring. This work provides the first evidence for a genetic factor causing hypermutation in the maternal germline of a mammal and adds to the very small list of naturally occurring variants known to modulate germline mutation rates in mammals.

K6-linked ubiquitylation marks formaldehyde-induced RNA-protein crosslinks for resolution

Reactive aldehydes are produced by normal cellular metabolism or after alcohol consumption, and they accumulate in human tissues if aldehyde clearance mechanisms are impaired. Their toxicity has been attributed to the damage they cause to genomic DNA and the subsequent inhibition of transcription and replication. However, whether interference with other cellular processes contributes to aldehyde toxicity has not been investigated. We demonstrate that formaldehyde induces RNA-protein crosslinks (RPCs) that stall the ribosome and inhibit translation in human cells. RPCs in the messenger RNA (mRNA) are recognized by the translating ribosomes, marked by atypical K6-linked ubiquitylation catalyzed by the RING-in-between-RING (RBR) E3 ligase RNF14, and subsequently resolved by the ubiquitin- and ATP-dependent unfoldase VCP. Our findings uncover an evolutionary conserved formaldehyde-induced stress response pathway that protects cells against RPC accumulation in the cytoplasm, and they suggest that RPCs contribute to the cellular and tissue toxicity of reactive aldehydes.

Chemoproteomic mapping of the glycolytic targetome in cancer cells

Hyperactivated glycolysis is a metabolic hallmark of most cancer cells. Although sporadic information has revealed that glycolytic metabolites possess nonmetabolic functions as signaling molecules, how these metabolites interact with and functionally regulate their binding targets remains largely elusive. Here, we introduce a target-responsive accessibility profiling (TRAP) approach that measures changes in ligand binding-induced accessibility for target identification by globally labeling reactive proteinaceous lysines. With TRAP, we mapped 913 responsive target candidates and 2,487 interactions for 10 major glycolytic metabolites in a model cancer cell line. The wide targetome depicted by TRAP unveils diverse regulatory modalities of glycolytic metabolites, and these modalities involve direct perturbation of enzymes in carbohydrate metabolism, intervention of an orphan transcriptional protein's activity and modulation of targetome-level acetylation. These results further our knowledge of how glycolysis orchestrates signaling pathways in cancer cells to support their survival, and inspire exploitation of the glycolytic targetome for cancer therapy.

Comparing methods for constructing and representing human pangenome graphs

BACKGROUND

As a single reference genome cannot possibly represent all the variation present across human individuals, pangenome graphs have been introduced to incorporate population diversity within a wide range of genomic analyses. Several data structures have been proposed for representing collections of genomes as pangenomes, in particular graphs.

RESULTS

In this work, we collect all publicly available high-quality human haplotypes and construct the largest human pangenome graphs to date, incorporating 52 individuals in addition to two synthetic references (CHM13 and GRCh38). We build variation graphs and de Bruijn graphs of this collection using five of the state-of-the-art tools: Bifrost, mdbg, Minigraph, Minigraph-Cactus and pggb. We examine differences in the way each of these tools represents variations between input sequences, both in terms of overall graph structure and representation of specific genetic loci.

CONCLUSION

This work sheds light on key differences between pangenome graph representations, informing end-users on how to select the most appropriate graph type for their application.

The ferroptosis landscape in acute myeloid leukemia

Ferroptosis induction through the suppression of glutathione peroxidase 4 (GPX4) and apoptosis-inducing factor mitochondria-associated 2 (AIFM2) has proven to be an effective approach in eliminating chemotherapy-resistant cells of various types. However, a comprehensive understanding of the roles of GPX4 and AIFM2 in acute myeloid leukemia (AML) has not yet been achieved. Using cBioPortal, DepMap, GEPIA, Metascape, and ONCOMINE, we compared the transcriptional expression, survival data, gene mutation, methylation, and network analyses of GPX4- and AIFM2-associated signaling pathways in AML. The results revealed that high expression levels of GPX4 and AIFM2 are associated with an adverse prognosis for AML patients. Overexpression of AIFM2 correlated with elevated mutation frequencies in NPM1 and DNMT3A. GPX4 upregulation modulated the following pathways: GO:0045333, cellular respiration; R-HSA-5389840, mitochondrial translation elongation; GO:0009060, aerobic respiration; R-HSA-9609507, protein localization; and R-HSA-8953854, metabolism of RNA. On the other hand, the overexpression of AIFM2 influenced the following processes: GO:0048704, embryonic skeletal system morphogenesis; GO:0021546, rhombomere development; GO:0009954, proximal/distal pattern formation; and GO:0048732, gland development. This study identifies the high expression of GPX4 and AIFM2 as novel biomarkers predicting a poor prognosis for AML patients. Furthermore, ferroptosis induction may improve the stratified treatment of AML.

Distinct disease mutations in DNMT3A result in a spectrum of behavioral, epigenetic, and transcriptional deficits

Phenotypic heterogeneity in monogenic neurodevelopmental disorders can arise from differential severity of variants underlying disease, but how distinct alleles drive variable disease presentation is not well understood. Here, we investigate missense mutations in DNA methyltransferase 3A (DNMT3A), a DNA methyltransferase associated with overgrowth, intellectual disability, and autism, to uncover molecular correlates of phenotypic heterogeneity. We generate a Dnmt3aP900L/+ mouse mimicking a mutation with mild to moderate severity and compare phenotypic and epigenomic effects with a severe R878H mutation. P900L mutants exhibit core growth and behavioral phenotypes shared across models but show subtle epigenomic changes, while R878H mutants display extensive disruptions. We identify mutation-specific dysregulated genes that may contribute to variable disease severity. Shared transcriptomic disruption identified across mutations overlaps dysregulation observed in other developmental disorder models and likely drives common phenotypes. Together, our findings define central drivers of DNMT3A disorders and illustrate how variable epigenomic disruption contributes to phenotypic heterogeneity in neurodevelopmental disease.

Breast cancer risk SNPs converge on estrogen receptor binding sites commonly shared between breast tumors to locally alter estrogen signalling output

Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα 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 ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

hipFG: high-throughput harmonization and integration pipeline for functional genomics data

SUMMARY

Preparing functional genomic (FG) data with diverse assay types and file formats for integration into analysis workflows that interpret genome-wide association and other studies is a significant and time-consuming challenge. Here we introduce hipFG (Harmonization and Integration Pipeline for Functional Genomics), an automatically customized pipeline for efficient and scalable normalization of heterogenous FG data collections into standardized, indexed, rapidly searchable analysis-ready datasets while accounting for FG datatypes (e.g. chromatin interactions, genomic intervals, quantitative trait loci).

AVAILABILITY AND IMPLEMENTATION

hipFG is freely available at https://bitbucket.org/wanglab-upenn/hipFG. A Docker container is available at https://hub.docker.com/r/wanglab/hipfg.

RGMb impacts partial epithelial-mesenchymal transition and BMP2-Induced ID mRNA expression independent of PD-L2 in nonsmall cell lung cancer cells

PD-1/PD-ligand-axis immunotherapy-mediated activation of T-cells for cancer cell elimination is a promising treatment of nonsmall cell lung cancer (NSCLC). However, the effect of immunotherapy on intracellular signaling pathways in cancer cells still needs further delineation. Repulsive Guidance Molecule b (RGMb), a regulator of Bone Morphogenetic Proteins (BMPs) signaling, interacts with the PD-ligand, PD-L2, at cancer cell membranes. Accordingly, a clarification of the functions of RGMb and its relation to PD-L2 might provide insight into NSCLC cell signaling responses to PD-1/PD-ligand-axis immunotherapy. In this study, the functions of RGMb and PD-L2 were examined using the two NSCLC cell lines HCC827 and A549. CRISPR/Cas9 was used to decrease the expression of RGMb and PD-L2, while lentiviral vectors were used to increase their expression. Downstream effects were examined by RT-qPCR and immunoassays. Ectopic expression of RGMb impacted BMP2-induced expression of ID1 and ID2 messenger RNA (mRNA) independently of PD-L2, while RGMb depletion by CRISPR/Cas9 did not affect the BMP2-mediated induction of ID1, ID2, and ID3 mRNA. However, depletion of RGMb resulted in a partial epithelial-mesenchymal transition (EMT) gene expression profile in HCC827 cells, which was not mimicked by PD-L2 depletion. The results show that RGMb is a coregulator of BMP signaling and hence, ID mRNA expression and that RGMb can control the EMT balance in NSCLC cells. However, RGMb appears to exert these functions independently of PD-L2, and accordingly, the PD-1/PD-ligand axis for immune surveillance in NSCLC cells.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Splicing across adipocyte differentiation is highly dynamic and impacted by metabolic phenotype

Adipose tissue dysfunction underlies many of the metabolic complications associated with obesity. A better understanding of the gene regulation differences present in metabolically unhealthy adipose tissue can provide insights into the mechanisms underlying adipose tissue dysfunction. Here, we used RNA-seq data collected from a differentiation time course of lean, obese, and obese with type 2 diabetes (T2D) individuals to characterize the role of alterative splicing in adipocyte differentiation and function. We found that splicing was highly dynamic across adipocyte differentiation in all three cohorts, and that the dynamics of splicing were significantly impacted by metabolic phenotype. We also found that there was very little overlap between genes that were differentially spliced in adipocyte differentiation and those that were differentially expressed, positioning alternative splicing as a largely independent gene regulatory mechanism whose impact would be missed when looking at gene expression changes alone. To assess the impact of alternative splicing across adipocyte differentiation on genetic risk for metabolic diseases, we integrated the differential splicing results generated here with genome-wide association study results for body mass index and T2D, and found that variants associated with T2D were enriched in regions that were differentially spliced in early differentiation. These findings provide insight into the role of alternative splicing in adipocyte differentiation and can serve as a resource to guide future variant-to-function studies.

Super enhancer-driven core transcriptional regulatory circuitry crosstalk with cancer plasticity and patient mortality in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is a clinically aggressive subtype of breast cancer. Core transcriptional regulatory circuitry (CRC) consists of autoregulated transcription factors (TFs) and their enhancers, which dominate gene expression programs and control cell fate. However, there is limited knowledge of CRC in TNBC. Herein, we systemically characterized the activated super-enhancers (SEs) and interrogated 14 CRCs in breast cancer. We found that CRCs could be broadly involved in DNA conformation change, metabolism process, and signaling response affecting the gene expression reprogramming. Furthermore, these CRC TFs are capable of coordinating with partner TFs bridging the enhancer-promoter loops. Notably, the CRC TF and partner pairs show remarkable specificity for molecular subtypes of breast cancer, especially in TNBC. USF1, SOX4, and MYBL2 were identified as the TNBC-specific CRC TFs. We further demonstrated that USF1 was a TNBC immunophenotype-related TF. Our findings that the rewiring of enhancer-driven CRCs was related to cancer immune and mortality, will facilitate the development of epigenetic anti-cancer treatment strategies.

Lifetime Postnatal Exposure to Perfluoroalkyl Substance Mixture and DNA Methylation at Twelve Years of Age

Per- and polyfluoroalkyl substance (PFAS) exposure has been linked to DNA methylation changes in neonates and adults. We previously reported that prenatal PFAS exposure may have a durable impact on DNA methylation from birth to adolescence. However, few studies have examined the association of postnatal PFAS exposure with alterations in DNA methylation. We examined the associations of lifetime postnatal PFAS mixture exposure with leukocyte DNA methylation in 154 adolescents from the HOME Study (2003-2006; Cincinnati, Ohio). Lifetime postnatal PFAS mixture exposure was estimated using latent profile analysis of four PFAS concentrations measured at birth, and ages 3, 8, and 12 years. We measured DNA methylation in peripheral leukocytes at 12 years using the Illumina HumanMethylation EPIC BeadChip. We estimated covariate-adjusted associations between postnatal PFAS mixture concentrations and DNA methylation measures using linear regression, and used KEGG enrichment analysis to identify molecular pathways. Four significant differentially methylated positions were observed in the higher vs. lower PFAS profile (FDR p-value <0.05). These PFAS-associated CpG sites annotated to gene regions related to various cancers, cognition, and cardiometabolic health. We identified 17 pathways (FDR p-value <0.05), which indicates possible mechanism linking PFAS exposure to several health effects.

Ultrarare Missense Variants Implicated in Utah Pedigrees Multiply Affected With Schizophrenia

Background

Recent work from the Schizophrenia Exome Sequencing Meta-analysis (SCHEMA) consortium showed significant enrichment of ultrarare variants in schizophrenia cases. Family-based studies offer a unique opportunity to evaluate rare variants because risk in multiplex pedigrees is more likely to be influenced by the same collection of variants than an unrelated cohort.

Methods

Here, we examine whole genome sequencing data from 35 individuals across 6 pedigrees multiply affected by schizophrenia. We applied a rigorous filtering pipeline to search for classes of protein-coding variants that cosegregated with disease status, and we examined these for evidence of enrichment in the SCHEMA dataset. Additionally, we applied a family-based consensus approach to call copy number variants and screen against a list of schizophrenia-associated risk variants.

Results

We identified deleterious missense variants in 3 genes (ATP2B2, SLC25A28, and GSK3A) that cosegregated with disease in 3 of the pedigrees. In the SCHEMA, the gene ATP2B2 shows highly suggestive evidence for deleterious missense variants in schizophrenia cases (p = .000072). ATP2B2 is involved in intracellular calcium homeostasis, expressed in multiple brain tissue types, and predicted to be intolerant to loss-of-function and missense variants.

Conclusions

We have identified genes that are likely to increase schizophrenia risk in 3 of the 6 pedigrees examined, the strongest evidence being for a gene involved in calcium homeostasis. Further work is required to examine other classes of variants that may be contributing to disease burden.

Identification of epigenetically active L1 promoters in the human brain and their relationship with psychiatric disorders

Long interspersed nuclear element-1 (LINE-1, L1) affects the transcriptome landscape in multiple ways. Promoter activity within its 5'UTR plays a critical role in regulating diverse L1 activities. However, the epigenetic status of L1 promoters in adult brain cells and their relationship with psychiatric disorders remain poorly understood. Here, we examined DNA methylation and hydroxymethylation of the full-length L1s in neurons and nonneurons and identified "epigenetically active" L1s. Notably, some of epigenetically active L1s were retrotransposition competent, which even had chimeric transcripts from the antisense promoters at their 5'UTRs. We also identified differentially methylated L1s in the prefrontal cortices of patients with psychiatric disorders. In nonneurons of bipolar disorder patients, one L1 was significantly hypomethylated and showed an inverse correlation with the expression level of the overlapping gene NREP. Finally, we observed that altered DNA methylation levels of L1 in patients with psychiatric disorders were not affected by the surrounding genomic regions but originated from the L1 sequences. These results suggested that altered epigenetic regulation of the L1 5'UTR in the brain was involved in the pathophysiology of psychiatric disorders.

Methamphetamine-induced region-specific transcriptomic and epigenetic changes in the brain of male rats

Psychostimulant methamphetamine (METH) is neurotoxic to the brain and, therefore, its misuse leads to neurological and psychiatric disorders. The gene regulatory network (GRN) response to neurotoxic METH binge remains unclear in most brain regions. Here we examined the effects of binge METH on the GRN in the nucleus accumbens, dentate gyrus, Ammon's horn, and subventricular zone in male rats. At 24 h after METH, ~16% of genes displayed altered expression and over a quarter of previously open chromatin regions - parts of the genome where genes are typically active - showed shifts in their accessibility. Intriguingly, most changes were unique to each area studied, and independent regulation between transcriptome and chromatin accessibility was observed. Unexpectedly, METH differentially impacted gene activity and chromatin accessibility within the dentate gyrus and Ammon's horn. Around 70% of the affected chromatin-accessible regions in the rat brain have conserved DNA sequences in the human genome. These regions frequently act as enhancers, ramping up the activity of nearby genes, and contain mutations linked to various neurological conditions. By sketching out the gene regulatory networks associated with binge METH in specific brain regions, our study offers fresh insights into how METH can trigger profound, region-specific molecular shifts.

A pan-cancer analysis of STAT3 expression and genetic alterations in human tumors

Combined cancer immunotherapy and targeted therapy have proven to be effective against various cancers and therefore have recently become the focus of cancer research. Signal transducer and activator of transcription 3 (STAT3) is a member of the STAT protein family of transcription factors. Several studies have shown that STAT3 can affect the prognosis of cancer patients by regulating immune microenvironment (IME). Therefore, STAT3 may have high research value for the development of combined immunotherapy/targeted therapy approaches for the treatment of cancer patients. We found differences in STAT3 expression between tumor and normal tissues. Kaplan-Meier survival and Cox regression analyses showed that high expression of STAT3 is associated with poor prognosis in low-grade glioma (LGG) patients. The results of the analysis of the area under the curve of the receiver operating characteristic curve further suggested that the expression of STAT3 is an effective way to evaluate the prognosis of patients with glioma. The results of the IME analysis revealed that the immune and matrix scores of LGGs were positively correlated with the expression of STAT3 (P < 0.05). The results of immune cell infiltration analysis showed that STAT3 was positively correlated with resting dendritic cells, eosinophils, neutrophils, M0 macrophages, M1 macrophages, CD4 memory resting T cells, and CD8 T cells in LGG patients, but negatively correlated with activated mast cells and M2 macrophages (P < 0.05). Our gene set enrichment analysis identified 384 enriched pathways. According to the enrichment scores, the top ten most significantly upregulated pathways were related to immune response. The top ten most significantly downregulated pathways were related to cell signal transduction and the regulation of cell survival, proliferation, and metabolism. Genetic alteration analysis showed that missense mutations in STAT3 account for the majority of mutations, and STAT3 mutations mostly occur in the Src homology domain. In conclusion overexpression of STAT3 can promote the development and growth of tumors by regulating IME, which is significantly related to the poor prognosis of cancer patients. Therefore, targeted inhibition of STAT3 expression may have high research value for the development of combined immunotherapy/targeted therapy approaches for the treatment of cancer patients.

Vascular inflammation and endothelial injury in SARS-CoV-2 infection: the overlooked regulatory cascades implicated by the ACE2 gene cluster

Coronavirus disease 2019 (COVID-19) has presented physicians with an unprecedented number of challenges and mortality. The basic question is why, in contrast to other 'respiratory' viruses, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection can result in such multi-systemic, life-threatening complications and a severe pulmonary vasculopathy. It is widely known that SARS-CoV-2 uses membrane-bound angiotensin-converting enzyme 2 (ACE2) as a receptor, resulting in internalization of the complex by the host cell. We discuss the evidence that failure to suppress coronaviral replication within 5 days results in sustained downregulation of ACE2 protein expression and that ACE2 is under negative-feedback regulation. We then expose openly available experimental repository data that demonstrate the gene for ACE2 lies in a novel cluster of inter-regulated genes on the X chromosome including PIR encoding pirin (quercetin 2,3-dioxygenase), and VEGFD encoding the predominantly lung-expressed vascular endothelial growth factor D. The five double-elite enhancer/promoters pairs that are known to be operational, and shared read-through lncRNA transcripts, imply that ongoing SARS-CoV-2 infection will reduce host defences to reactive oxygen species, directly generate superoxide O2·- and H2O2 (a ' ROS storm'), and impair pulmonary endothelial homeostasis. Published cellular responses to oxidative stress complete the loop to pathophysiology observed in severe COVID-19. Thus, for patients who fail to rapidly suppress viral replication, the newly appreciated ACE2 co-regulated gene cluster predicts delayed responses that would account for catastrophic deteriorations. We conclude that ACE2 homeostatic drives provide a unified understanding that should help optimize therapeutic approaches during the wait until safe, effective vaccines and antiviral therapies for SARS-CoV-2 are delivered.

Pan-cancer analysis of IFN-γ with possible immunotherapeutic significance: a verification of single-cell sequencing and bulk omics research

Background

Interferon-gamma (IFN-γ), commonly referred to as type II interferon, is a crucial cytokine that coordinates the tumor immune process and has received considerable attention in tumor immunotherapy research. Previous studies have discussed the role and mechanisms associated with IFN-γ in specific tumors or diseases, but the relevant role of IFN-γ in pan-cancer remains uncertain.

Methods

TCGA and GTEx RNA expression data and clinical data were downloaded. Additionally, we analyzed the role of IFN-γ on tumors by using a bioinformatic approach, which included the analysis of the correlation between IFN-γ in different tumors and expression, prognosis, functional status, TMB, MSI, immune cell infiltration, and TIDE. We also developed a PPI network for topological analysis of the network, identifying hub genes as those having a degree greater than IFN-γ levels.

Result

IFN-γ was differentially expressed and predicted different survival statuses in a majority of tumor types in TCGA. Additionally, IFN-γ expression was strongly linked to factors like infiltration of T cells, immune checkpoints, immune-activating genes, immunosuppressive genes, chemokines, and chemokine receptors, as well as tumor purity, functional statuses, and prognostic value. Also, prognosis, CNV, and treatment response were all substantially correlated with IFN-γ-related gene expression. Particularly, the IFN-γ-related gene STAT1 exhibited the greatest percentage of SNVs and the largest percentage of SNPs in UCEC. Elevated expression levels of IFN-γ-related genes were found in a wide variety of tumor types, and this was shown to be positively linked to drug sensitivity for 20 different types of drugs.

Conclusion

IFN-γ is a good indicator of response to tumor immunotherapy and is likely to limit tumor progression, offering a novel approach for immunotherapy's future development.

Genetic features and genomic targets of human KRAB-zinc finger proteins

Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KZFPs) are one of the largest groups of transcription factors encoded by tetrapods, with 378 members in human alone. KZFP genes are often grouped in clusters reflecting amplification by gene and segment duplication since the gene family first emerged more than 400 million years ago. Previous work has revealed that many KZFPs recognize transposable element (TE)-embedded sequences as genomic targets, and that KZFPs facilitate the co-option of the regulatory potential of TEs for the benefit of the host. Here, we present a comprehensive survey of the genetic features and genomic targets of human KZFPs, notably completing past analyses by adding data on close to a hundred family members. General principles emerge from our study of the TE-KZFP regulatory system, which point to multipronged evolutionary mechanisms underlaid by highly complex and combinatorial modes of action with strong influences on human speciation.

Lack of Stim2 Affects Vision-Dependent Behavior and Sensitivity to Hypoxia

Stromal interaction molecules (STIMs) are endoplasmic reticulum-resident proteins that regulate Ca2+ homeostasis and signaling by store-operated calcium entry (SOCE). The different properties and functions of STIM1 and STIM2 have been described mostly based on work in vitro. STIM2 knockout mice do not survive until adulthood. Therefore, we generated and characterized stim2a and stim2b double-knockout zebrafish. The (stim2a;stim2b)-/- zebrafish did not have any apparent morphological phenotype. However, RNA sequencing revealed 1424 differentially expressed genes. One of the most upregulated genes was annexin A3a, which is a marker of activated microglia. This corresponded well to an increase in Neutral Red staining in the in vivo imaging of the (stim2a;stim2b)-/- zebrafish brain. The lack of Stim2 decreased zebrafish survival under low oxygen conditions. Behavioral tests, such as the visual-motor response test and dark-light preference test, indicated that (stim2a;stim2b)-/- larvae might have problems with vision. This was consistent with the downregulation of many genes that are related to light perception. The periodic acid-Schiff staining of retina sections from adult zebrafish revealed alterations of the stratum pigmentosum, suggesting the involvement of a Stim2-dependent process in visual perception. Altogether, these data reveal new functions for Stim2 in the nervous system.

RNA covariation at helix-level resolution for the identification of evolutionarily conserved RNA structure

Many biologically important RNAs fold into specific 3D structures conserved through evolution. Knowing when an RNA sequence includes a conserved RNA structure that could lead to new biology is not trivial and depends on clues left behind by conservation in the form of covariation and variation. For that purpose, the R-scape statistical test was created to identify from alignments of RNA sequences, the base pairs that significantly covary above phylogenetic expectation. R-scape treats base pairs as independent units. However, RNA base pairs do not occur in isolation. The Watson-Crick (WC) base pairs stack together forming helices that constitute the scaffold that facilitates the formation of the non-WC base pairs, and ultimately the complete 3D structure. The helix-forming WC base pairs carry most of the covariation signal in an RNA structure. Here, I introduce a new measure of statistically significant covariation at helix-level by aggregation of the covariation significance and covariation power calculated at base-pair-level resolution. Performance benchmarks show that helix-level aggregated covariation increases sensitivity in the detection of evolutionarily conserved RNA structure without sacrificing specificity. This additional helix-level sensitivity reveals an artifact that results from using covariation to build an alignment for a hypothetical structure and then testing the alignment for whether its covariation significantly supports the structure. Helix-level reanalysis of the evolutionary evidence for a selection of long non-coding RNAs (lncRNAs) reinforces the evidence against these lncRNAs having a conserved secondary structure.

PET/MRI and Bioluminescent Imaging Identify Hypoxia as a Cause of Programmed Cell Death Ligand 1 Image Heterogeneity

Purpose To examine the association between hypoxia and programmed cell death ligand 1 (PD-L1) expression using bioluminescence imaging (BLI) and PET/MRI in a syngeneic mouse model of triple-negative breast cancer (TNBC). Materials and Methods PET/MRI and optical imaging were used to determine the role of hypoxia in altering PD-L1 expression using a syngeneic TNBC model engineered to express luciferase under hypoxia. Results Imaging showed a close spatial association between areas of hypoxia and increased PD-L1 expression in the syngeneic murine (4T1) tumor model. Mouse and human TNBC cells exposed to hypoxia exhibited a significant increase in PD-L1 expression, consistent with the in vivo imaging data. The role of hypoxia in increasing PD-L1 expression was further confirmed by using The Cancer Genome Atlas analyses of different human TNBCs. Conclusion These results have identified the potential role of hypoxia in contributing to PD-L1 heterogeneity in tumors by increasing cancer cell PD-L1 expression. Keywords: Hypoxia, PD-L1, Triple-Negative Breast Cancer, PET/MRI, Bioluminescence Imaging Supplemental material is available for this article. © RSNA, 2023.

Tumour mutations in long noncoding RNAs enhance cell fitness

Long noncoding RNAs (lncRNAs) are linked to cancer via pathogenic changes in their expression levels. Yet, it remains unclear whether lncRNAs can also impact tumour cell fitness via function-altering somatic "driver" mutations. To search for such driver-lncRNAs, we here perform a genome-wide analysis of fitness-altering single nucleotide variants (SNVs) across a cohort of 2583 primary and 3527 metastatic tumours. The resulting 54 mutated and positively-selected lncRNAs are significantly enriched for previously-reported cancer genes and a range of clinical and genomic features. A number of these lncRNAs promote tumour cell proliferation when overexpressed in in vitro models. Our results also highlight a dense SNV hotspot in the widely-studied NEAT1 oncogene. To directly evaluate the functional significance of NEAT1 SNVs, we use in cellulo mutagenesis to introduce tumour-like mutations in the gene and observe a significant and reproducible increase in cell fitness, both in vitro and in a mouse model. Mechanistic studies reveal that SNVs remodel the NEAT1 ribonucleoprotein and boost subnuclear paraspeckles. In summary, this work demonstrates the utility of driver analysis for mapping cancer-promoting lncRNAs, and provides experimental evidence that somatic mutations can act through lncRNAs to enhance pathological cancer cell fitness.

Comprehensive analysis of the coding and non-coding RNA transcriptome expression profiles of hippocampus tissue in tx-J animal model of Wilson's disease

Wilson's disease (WD) is an autosomal recessive disorder with a genetic basis. The predominant non-motor symptom of WD is cognitive dysfunction, although the specific genetic regulatory mechanism remains unclear. Tx-J mice, with an 82% sequence homology of the ATP7B gene to the human gene, are considered the most suitable model for WD. This study employs deep sequencing to investigate the differences in RNA transcript profiles, both coding and non-coding, as well as the functional characteristics of the regulatory network involved in WD cognitive impairment. The cognitive function of tx-J mice was evaluated using the Water Maze Test (WMT). Long non-coding RNA (lncRNA), circular RNA (circRNA), and messenger RNA (mRNA) profiles were analyzed in the hippocampal tissue of tx-J mice to identify differentially expressed RNAs (DE-RNAs). Subsequently, the DE-RNAs were used to construct protein-protein interaction (PPI) networks, as well as DE-circRNAs and lncRNAs-associated competing endogenous RNA (ceRNA) expression networks, and coding-noncoding co-expression (CNC) networks. To elucidate their biological functions and pathways, the PPI and ceRNA networks were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. A total of 361 differentially expressed mRNAs (DE-mRNAs), comprising 193 up-regulated and 168 down-regulated mRNAs, 2627 differentially expressed long non-coding RNAs (DE-lncRNAs), consisting of 1270 up-regulated and 1357 down-regulated lncRNAs, and 99 differentially expressed circular RNAs (DE-circRNAs), consisting of 68 up-regulated and 31 down-regulated circRNAs, were observed in the tx-J mice group when compared to the control mice group. Gene Ontology (GO) and pathway analyses revealed that DE-mRNAs were enriched in cellular processes, calcium signaling pathways, and mRNA surveillance pathways. In contrast, the DE-circRNAs-associated competing endogenous RNA (ceRNA) network was enriched for covalent chromatin modification, histone modification, and axon guidance, whereas the DE-lncRNAs-associated ceRNA network was enriched for dendritic spine, regulation of cell morphogenesis involved in differentiation, and mRNA surveillance pathway. The study presented the expression profiles of lncRNA, circRNA, and mRNA in the hippocampal tissue of tx-J mice. Furthermore, the study constructed PPI, ceRNA, and CNC expression networks. The findings are significant in comprehending the function of regulatory genes in WD associated with cognitive impairment. These results also offer valuable information for the diagnosis and treatment of WD.

Thirteen dubious ways to detect conserved structural RNAs

Covariation induced by compensatory base substitutions in RNA alignments is a great way to deduce conserved RNA structure, in principle. In practice, success depends on many factors, importantly the quality and depth of the alignment and the choice of covariation statistic. Measuring covariation between pairs of aligned positions is easy. However, using covariation to infer evolutionarily conserved RNA structure is complicated by other extraneous sources of covariation such as that resulting from homologous sequences having evolved from a common ancestor. In order to provide evidence of evolutionarily conserved RNA structure, a method to distinguish covariation due to sources other than RNA structure is necessary. Moreover, there are several sorts of artifactually generated covariation signals that can further confound the analysis. Additionally, some covariation signal is difficult to detect due to incomplete comparative data. Here, we investigate and critically discuss the practice of inferring conserved RNA structure by comparative sequence analysis. We provide new methods on how to approach and decide which of the numerous long non-coding RNAs (lncRNAs) have biologically relevant structures.

Pan-cancer analysis identifies the correlations of Thymosin Beta 10 with predicting prognosis and immunotherapy response

INTRODUCTION

The biological function and prognosis roles of thymosin β(TMSB) 10 are still unclear in pan-cancer.

METHODS

We retrieved The Cancer Genome Atlas and Genotype-tissue expression datasets to obtain the difference of TMSB10 expression between pan-cancer and normal tissues, and analyzed the biological function and prognosis role of TMSB10 in pan-cancer by using cBioPortal Webtool.

RESULTS

The expression of TMSB10 in tumor tissues was significantly higher than normal tissues, and showed the potential ability to predict the prognosis of patients in Pan-cancer. It was found that TMSB10 was significantly correlated with tumor microenvironment, immune cell infiltration and immune regulatory factor expression. TMSB10 is involved in the regulation of cellular signal transduction pathways in a variety of tumors, thereby mediating the occurrence of tumor cell invasion and metastasis. Finally, TMSB10 can not only effectively predict the anti-PD-L1 treatment response of cancer patients, but also be used as an important indicator to evaluate the sensitivity of chemotherapy. In vitro, low expression of TMSB10 inhibited clonogenic formation ability, invasion, and migration in glioma cells. Furthermore, TMSB10 may involve glioma immune regulation progression by promoting PD-L1 expression levels via activating STAT3 signaling pathway.

CONCLUSIONS

Our results show that TMSB10 is abnormally expressed in tumor tissues, which may be related to the infiltration of immune cells in the tumor microenvironment. Clinically, TMSB10 is not only an effective prognostic factor for predicting the clinical treatment outcome of cancer patients, but also a promising biomarker for predicting the effect of tumor immune checkpoint inhibitors (ICIs) and chemotherapy in some cancers.

hipFG: High-throughput harmonization and integration pipeline for functional genomics data

Preparing functional genomic (FG) data with diverse assay types and file formats for integration into analysis workflows that interpret genome-wide association and other studies is a significant and time-consuming challenge. Here we introduce hipFG, an automatically customized pipeline for efficient and scalable normalization of heterogenous FG data collections into standardized, indexed, rapidly searchable analysis-ready datasets while accounting for FG datatypes (e.g., chromatin interactions, genomic intervals, quantitative trait loci).

RNA covariation at helix-level resolution for the identification of evolutionarily conserved RNA structure

Many biologically important RNAs fold into specific 3D structures conserved through evolution. Knowing when an RNA sequence includes a conserved RNA structure that could lead to new biology is not trivial and depends on clues left behind by conservation in the form of covariation and variation. For that purpose, the R-scape statistical test was created to identify from alignments of RNA sequences, the base pairs that significantly covary above phylogenetic expectation. R-scape treats base pairs as independent units. However, RNA base pairs do not occur in isolation. The Watson-Crick (WC) base pairs stack together forming helices that constitute the scaffold that facilitates the formation of the non-WC base pairs, and ultimately the complete 3D structure. The helix-forming WC base pairs carry most of the covariation signal in an RNA structure. Here, I introduce a new measure of statistically significant covariation at helix-level by aggregation of the covariation significance and covariation power calculated at base-pair-level resolution. Performance benchmarks show that helix-level aggregated covariation increases sensitivity in the detection of evolutionarily conserved RNA structure without sacrificing specificity. This additional helix-level sensitivity reveals an artifact that results from using covariation to build an alignment for a hypothetical structure and then testing the alignment for whether its covariation significantly supports the structure. Helix-level reanalysis of the evolutionary evidence for a selection of long non-coding RNAs (lncRNAs) reinforces the evidence against these lncRNAs having a conserved secondary structure.

Availability

Helix aggregated E-values are integrated in the R-scape software package (version 2.0.0.p and higher). The R-scape web server eddylab.org/R-scape includes a link to download the source code.

Contact

elenarivas@fas.harvard.edu.

Supplementary information

Supplementary data and code are provided with this manuscript at rivaslab.org .

The origins and functional effects of postzygotic mutations throughout the human life span

Postzygotic mutations (PZMs) begin to accrue in the human genome immediately after fertilization, but how and when PZMs affect development and lifetime health remain unclear. To study the origins and functional consequences of PZMs, we generated a multitissue atlas of PZMs spanning 54 tissue and cell types from 948 donors. Nearly half the variation in mutation burden among tissue samples can be explained by measured technical and biological effects, and 9% can be attributed to donor-specific effects. Through phylogenetic reconstruction of PZMs, we found that their type and predicted functional impact vary during prenatal development, across tissues, and through the germ cell life cycle. Thus, methods for interpreting effects across the body and the life span are needed to fully understand the consequences of genetic variants.

Application of single-cell RNA sequencing on human testicular samples: a comprehensive review

So far there has been no comprehensive review using systematic literature search strategies to show the application of single-cell RNA sequencing (scRNA-seq) in the human testis of the whole life cycle (from embryos to aging males). Here, we summarized the application of scRNA-seq analyses on various human testicular biological samples. A systematic search was conducted in PubMed and Gene Expression Omnibus (GEO), focusing on English researches published after 2009. Articles related to GEO data-series were also retrieved in PubMed or BioRxiv. 81 full-length studies were finally included in the review. ScRNA-seq has been widely used on different human testicular samples with various library strategies, and new cell subtypes such as State 0 spermatogonial stem cells (SSC) and stage_a/b/c Sertoli cells (SC) were identified. For the development of normal testes, scRNA-seq-based evidence showed dynamic transcriptional changes of both germ cells and somatic cells from embryos to adults. And dysregulated metabolic signaling or hedgehog signaling were revealed by scRNA-seq in aged SC or Leydig cells (LC), respectively. For infertile males, scRNA-seq studies revealed profound changes of testes, such as the increased proportion of immature SC/LC of Klinefelter syndrome, the somatic immaturity and altered germline autophagy of patients with non-obstructive azoospermia, and the repressed differentiation of SSC in trans-females receiving testosterone inhibition therapy. Besides, the re-analyzing of public scRNA-seq data made further discoveries such as the potential vulnerability of testicular SARS-CoV-2 infection, and both evolutionary conservatism and divergence among species. ScRNA-seq analyses would unveil mechanisms of testes' development and changes so as to help developing novel treatments for male infertility.

Multiallelic Copy Number Variation in ORM1 is Associated with Plasma Cell-Free DNA Levels as an Intermediate Phenotype for Venous Thromboembolism

Venous thromboembolism (VTE) is a common disease with high heritability. However, only a small portion of the genetic variance of VTE can be explained by known genetic risk factors. Neutrophil extracellular traps (NETs) have been associated with prothrombotic activity. Therefore, the genetic basis of NETs could reveal novel risk factors for VTE. A recent genome-wide association study of plasma cell-free DNA (cfDNA) levels in the Genetic Analysis of Idiopathic Thrombophilia 2 (GAIT-2) Project showed a significant associated locus near ORM1. We aimed to further explore this candidate region by next-generation sequencing, copy number variation (CNV) quantification, and expression analysis using an extreme phenotype sampling design involving 80 individuals from the GAIT-2 Project. The RETROVE study with 400 VTE cases and 400 controls was used to replicate the results. A total of 105 genetic variants and a multiallelic CNV (mCNV) spanning ORM1 were identified in GAIT-2. Of these, 17 independent common variants, a region of 22 rare variants, and the mCNV were significantly associated with cfDNA levels. In addition, eight of these common variants and the mCNV influenced ORM1 expression. The association of the mCNV and cfDNA levels was replicated in RETROVE (p-value = 1.19 × 10-6). Additional associations between the mCNV and thrombin generation parameters were identified. Our results reveal that increased mCNV dosages in ORM1 decreased gene expression and upregulated cfDNA levels. Therefore, the mCNV in ORM1 appears to be a novel marker for cfDNA levels, which could contribute to VTE risk.

FZD7, Regulated by Non-CpG Methylation, Plays an Important Role in Immature Porcine Sertoli Cell Proliferation

The regulatory role of non-CpG methylation in mammals has been important in whole-genome bisulfite sequencing. It has also been suggested that non-CpG methylation regulates gene expression to affect the development and health of mammals. However, the dynamic regulatory mechanisms of genome-wide, non-CpG methylation during testicular development still require intensive study. In this study, we analyzed the dataset from the whole-genome bisulfite sequencing (WGBS) and the RNA-seq of precocious porcine testicular tissues across two developmental stages (1 and 75 days old) in order to explore the regulatory roles of non-CpG methylation. Our results showed that genes regulated by non-CpG methylation affect the development of testes in multiple pathways. Furthermore, several hub genes that are regulated by non-CpG methylation during testicular development-such as VEGFA, PECAM1, and FZD7-were also identified. We also found that the relative expression of FZD7 was downregulated by the zebularine-induced demethylation of the first exon of FZD7. This regulatory relationship was consistent with the results of the WGBS and RNA-seq analysis. The immature porcine Sertoli cells were transfected with RNAi to mimic the expression patterns of FZD7 during testicular development. The results of the simulation test showed that cell proliferation was significantly impeded and that cell cycle arrest at the G2 phase was caused by the siRNA-induced FZD7 inhibition. We also found that the percentage of early apoptotic Sertoli cells was decreased by transfecting them with the RNAi for FZD7. This indicates that FZD7 is an important factor in linking the proliferation and apoptosis of Sertoli cells. We further demonstrated that Sertoli cells that were treated with the medium collected from apoptotic cells could stimulate proliferation. These findings will contribute to the exploration of the regulatory mechanisms of non-CpG methylation in testicular development and of the relationship between the proliferation and apoptosis of normal somatic cells.

Identification of ZBTB4 as an immunological biomarker that can inhibit the proliferation and invasion of pancreatic cancer

BACKGROUND

Zinc finger and BTB domain-containing protein 4 (ZBTB4) belongs to the zinc finger protein family, which has a role in regulating epigenetic inheritance and is associated with cell differentiation and proliferation. Previous studies have identified aberrant ZBTB4 expression in cancer and its ability to modulate disease progression, but studies on the immune microenvironment, immunotherapy and its role in cancer are still lacking.

METHODS

Human pan-cancer and normal tissue transcriptome data were obtained from The Cancer Genome Atlas. The pan-cancer genomic alteration landscape of ZBTB4 was investigated with the online tool. The Kaplan-Meier method was used to evaluate the prognostic significance of ZBTB4 in pancreatic cancer. In parallel, ZBTB4 interacting molecules and potential functions were analyzed by co-expression and the correlation between ZBTB4 and immune cell infiltration, immune modulatory cells and efficacy of immune checkpoint therapy was explored. Next, we retrieved the Gene Expression Omnibus database expression datasets of ZBTB4 and investigated ZBTB4 expression and clinical significance in pancreatic cancer by immunohistochemical staining experiments. Finally, cell experiments were performed to investigate changes in pancreatic cancer cell proliferation, migration and invasion following overexpression and knockdown of ZBTB4.

FINDINGS

ZBTB4 showed loss of expression in the majority of tumors and possessed the ability to predict cancer prognosis. ZBTB4 was closely related to the tumor immune microenvironment, immune cell infiltration and immunotherapy efficacy. ZBTB4 had good diagnostic performance for pancreatic cancer in the clinic, and ZBTB4 protein expression was lost in pancreatic cancer tumor tissues. Cell experiments revealed that overexpression of ZBTB4 inhibited the proliferation, migration and invasion of pancreatic cancer cells, while silencing ZBTB4 showed the opposite effect.

CONCLUSIONS

According to our results, ZBTB4 is present in pancreatic cancer with aberrant expression and is associated with an altered immune microenvironment. We show that ZBTB4 is a promising marker for cancer immunotherapy and cancer prognosis and has the potential to influence pancreatic cancer progression.

dbCNV: deleteriousness-based model to predict pathogenicity of copy number variations

BACKGROUND

Copy number variation (CNV) is a type of structural variation, which is a gain or loss event with abnormal changes in copy number. Methods to predict the pathogenicity of CNVs are required to realize the relationship between these variants and clinical phenotypes. ClassifyCNV, X-CNV, StrVCTVRE, etc. have been trained to predict the pathogenicity of CNVs, but few studies have been reported based on the deleterious significance of features.

RESULTS

From single nucleotide polymorphism (SNP), gene and region dimensions, we collected 79 informative features that quantitatively describe the characteristics of CNV, such as CNV length, the number of protein genes, the number of three prime untranslated region. Then, according to the deleterious significance, we formulated quantitative methods for features, which fall into two categories: the first is variable type, including maximum, minimum and mean; the second is attribute type, which is measured by numerical sum. We used Gradient Boosted Trees (GBT) algorithm to construct dbCNV, which can be used to predict pathogenicity for five-tier classification and binary classification of CNVs. We demonstrated that the distribution of most feature values was consistent with the deleterious significance. The five-tier classification model accuracy for 0.85 and 0.79 in loss and gain CNVs, which proved that it has high discrimination power in predicting the pathogenicity of five-tier classification CNVs. The binary model achieved area under curve (AUC) values of 0.96 and 0.81 in the validation set, respectively, in gain and loss CNVs.

CONCLUSION

The performance of the dbCNV suggest that functional deleteriousness-based model of CNV is a promising approach to support the classification prediction and to further understand the pathogenic mechanism.

Investigation of the clinical utility of adhesion molecules in the management of thyroid nodules

To better understand the relationship among cell adhesion molecules (CAM) and investigate the clinical diagnostic and prognostic application of ICAM-1 (ICAM1), LFA-1 (ITGAL), and L-selectin (SELL) proteins and mRNA corresponding expression in thyroid cancer. Gene expression was evaluated by RT-qPCR, and protein expression was evaluated by immunohistochemistry. We evaluated 275 patients (218 women, 57 men, 48.4 ± 14.5 years old), including 102 benign and 173 malignant nodules. The 143 papillary thyroid carcinoma (PTC) and 30 follicular thyroid carcinoma (FTC) patients were managed according to current guidelines and followed-up for 78.7 ± 54.2 months. Malignant and benign nodules differed concerning mRNA (p = 0.0027) and protein (p = 0.0020 for nuclear) expression of L-selectin and ICAM-1 (mRNA: p = 0.0001 and protein: p = 0.0014) and protein expression of LFA-1 (p = 0.0168), but not mRNA expression of LFA-1 (p = 0.2131). SELL expression was more intense in malignant tumors (p = 0.0027). ICAM1 (p = 0.0064) and ITGAL (p = 0.0244) mRNA expression was higher in tumors with lymphocyte infiltrate. ICAM-1 expression correlated with younger age at diagnosis (p = 0.0312) and smaller tumor size (p = 0.0443). Also, LFA-1 expression correlated with higher age at diagnosis (p = 0.0376) and was more intense at stage III and IV (p = 0.0077). In general, the protein expression of the 3 CAM decreased as the process of cellular dedifferentiation occurred. We suggest that the SELL and ICAM1 genes and L-selectin and LFA-1 protein expression may help confirm malignancy and assist in the histological characterization of follicular patterned lesions, but we were unable to correlate these CAMs with patient outcomes.

Multiomics reveals glutathione metabolism as a driver of bimodality during stem cell aging

With age, skeletal muscle stem cells (MuSCs) activate out of quiescence more slowly and with increased death, leading to defective muscle repair. To explore the molecular underpinnings of these defects, we combined multiomics, single-cell measurements, and functional testing of MuSCs from young and old mice. The multiomics approach allowed us to assess which changes are causal, which are compensatory, and which are simply correlative. We identified glutathione (GSH) metabolism as perturbed in old MuSCs, with both causal and compensatory components. Contrary to young MuSCs, old MuSCs exhibit a population dichotomy composed of GSHhigh cells (comparable with young MuSCs) and GSHlow cells with impaired functionality. Mechanistically, we show that antagonism between NRF2 and NF-κB maintains this bimodality. Experimental manipulation of GSH levels altered the functional dichotomy of aged MuSCs. These findings identify a novel mechanism of stem cell aging and highlight glutathione metabolism as an accessible target for reversing MuSC aging.