Catalog of 5' Fusion Partners in ALK-positive NSCLC Circa 2020

Since the discovery of anaplastic lymphoma kinase fusion-positive (ALK+) NSCLC in 2007, the methods to detect ALK+ NSCLC have evolved and expanded from fluorescence in situ hybridization and immunohistochemistry to next-generation DNA sequencing, targeted RNA sequencing, and whole transcriptome sequencing. As such, the deep sequencing methods have resulted in the expansion of distinct fusion partners identified in ALK+ NSCLC to 90 (one variant PLEKHM2-ALK is found in small cell lung cancer but included in this catalog) by the end of January 2020; about 65 of them (since 2018) and most of the recent novel fusion partners were reported from China. Thirty-four of the distinct fusion partners are located on the short arm of chromosome 2; 28 of these 34 fusion partners are located on 2p21-25, in which ALK is located on 2p23.2-p23.1. Many of these new ALK+ NSCLC fusion variants have responded to ALK tyrosine kinase inhibitors (TKIs). Several of these novel ALK fusion variants were identified as being resistant to EGFR TKIs or as dual 3'ALK fusions. In addition, at least 28 intergenic ALK rearrangements have also been reported, with three of them reported as responding to crizotinib. This review aims to serve as a central source of reference of fusion partners in ALK+ NSCLC for clinicians and scientists. We aim to update and improve the list going forward.

The comparison of clinical and biological characteristics between IDH1 and IDH2 mutations in gliomas

Background

Mutations in isocitrate dehydrogenase 1 (IDH1) and isocitrate dehydrogenase 2 (IDH2) are frequent in low-grade gliomas and secondary glioblastomas (sGBM). Because they yield the same oncometabolite, D-2-hydroxyglutarate, they are often treated as equivalent and pooled. The objective of this study was to provide insight into the differences between IDH1 and IDH2 mutant gliomas.

Methods

To investigate the different clinical and molecular characterization between IDH1 mutant and IDH2 mutant gliomas, we studied 811 patients with IDH1 mutations, IDH2 mutations and IDH1/2 wild-type. In addition, whole-transcriptome sequencing and DNA methylation data were used to assess the distribution of genetic changes in IDH1 and IDH2 mutant gliomas in a Chinese population-based cohort.

Results

Among 811 gliomas in our cohort, 448 cases (55.2 %) harbored an IDH1 mutation, 18 cases (2.2 %) harbored an IDH2 mutation and 345 cases (42.6 %) harbored an IDH1/2 wild-type. We found that IDH1 and IDH2 are mutually exclusive in gliomas, and IDH2 mutations are mutually exclusive with PTEN, P53 and ATRX mutations. Patients with IDH2 mutations had a higher frequency of 1p/19q co-deletion (p < 0.05) than IDH1 mutant patients. In addition, a Gene Set Enrichment Analysis (GSEA) showed that IDH2 mutant gliomas were associated with the oxidative phosphorylation gene set, and the four most representative biological processes for genes commonly altered by hypermethylation in IDH2 mutant gliomas were the regulation of cell proliferation, cell motion, cell migration and response to hypoxia. Patients with IDH2 mutant gliomas exhibited longer Overall survival (OS) (p < 0.05) and longer Progression-free survival (PFS) (p < 0.05) than patients with IDH1/2 wild-type gliomas. However, their OS and PFS did not differ from that of IDH1 mutant patients.

Conclusions

Our study revealed an intrinsic distinction between IDH1 and IDH2 mutant gliomas, and these mutations should be considered separately because their differences could have implications for the diagnosis and treatment of IDH1/2 mutant gliomas.

Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots

BACKGROUND

Environmental stimuli can activate a series of physiological and biochemical responses in plants accompanied by extensive transcriptional reprogramming. Long non-coding RNAs (lncRNAs), as versatile regulators, control gene expression in multiple ways and participate in the adaptation to biotic and abiotic stresses.

RESULTS

In this study, soybean seedlings were continuously cultured for 15 days with high salinity solutions started from seed germination. Strand-specific whole transcriptome sequencing and stringent bioinformatic analysis led to the identification of 3030 long intergenic non-coding RNAs (lincRNAs) and 275 natural antisense transcripts (lncNATs) in soybean roots. In contrast to mRNAs, newly identified lncRNAs exhibited less exons, similar AU content to UTRs, even distribution across the genome and low evolutionary conservation. Remarkably, more than 75% of discovered lncRNAs that were activated or up-regulated by continuous salt stress mainly targeted proteins with binding and catalytic activities. Furthermore, two DNA methylation maps with single-base resolution were generated by using reduced representation bisulfite sequencing, offering a genome-wide perspective and important clues for epigenetic regulation of stress-associated lncRNAs and protein-coding genes.

CONCLUSIONS

Taken together, our findings systematically demonstrated the characteristics of continuous salt stress-induced lncRNAs and extended the knowledge of corresponding methylation profiling, providing valuable evidence for a better understanding of how plants cope with long-term salt stress circumstances.

Reproducible bioinformatics project: a community for reproducible bioinformatics analysis pipelines

Background

Reproducibility of a research is a key element in the modern science and it is mandatory for any industrial application. It represents the ability of replicating an experiment independently by the location and the operator. Therefore, a study can be considered reproducible only if all used data are available and the exploited computational analysis workflow is clearly described. However, today for reproducing a complex bioinformatics analysis, the raw data and the list of tools used in the workflow could be not enough to guarantee the reproducibility of the results obtained. Indeed, different releases of the same tools and/or of the system libraries (exploited by such tools) might lead to sneaky reproducibility issues.

Results

To address this challenge, we established the Reproducible Bioinformatics Project (RBP), which is a non-profit and open-source project, whose aim is to provide a schema and an infrastructure, based on docker images and R package, to provide reproducible results in Bioinformatics. One or more Docker images are then defined for a workflow (typically one for each task), while the workflow implementation is handled via R-functions embedded in a package available at github repository. Thus, a bioinformatician participating to the project has firstly to integrate her/his workflow modules into Docker image(s) exploiting an Ubuntu docker image developed ad hoc by RPB to make easier this task. Secondly, the workflow implementation must be realized in R according to an R-skeleton function made available by RPB to guarantee homogeneity and reusability among different RPB functions. Moreover she/he has to provide the R vignette explaining the package functionality together with an example dataset which can be used to improve the user confidence in the workflow utilization.

Conclusions

Reproducible Bioinformatics Project provides a general schema and an infrastructure to distribute robust and reproducible workflows. Thus, it guarantees to final users the ability to repeat consistently any analysis independently by the used UNIX-like architecture.

Dibutyl phthalate disrupts conserved circadian rhythm in Drosophila and human cells

People are constantly exposed to phthalates, due to their common use in the production of plastics, pharmaceuticals, cosmetics and skin care products. The ability of phthalates to disrupt endocrine signaling, leading to developmental, reproductive and metabolic defects, has been studied, yet how phthalates interfere with these biological functions is still unclear. To uncover DBP interacting molecular pathways, we raised Drosophila melanogaster on food containing dibutyl phthalate (DBP) at various concentrations. Whole transcriptome analysis of adult Drosophila reveals that DBP exposure throughout development disrupts the expression of genes central to circadian rhythm regulation, including increased expression of vrille (vri, human NFIL3), timeless (tim, human TIMELESS) and period (per, human PER3), with decreased expression of Pigment-dispersing factor (Pdf). DBP exposure also alters the expression of the evolutionarily conserved nuclear receptor Hormone receptor-like in 38 (Hr38, human NR4A2), which is known to regulate Pdf expression. Furthermore, behavioral assays determined that exposing Drosophila to DBP throughout development modifies the circadian rhythm of adults. Although DBP inhibits the expression of signaling systems regulating vision, including Rh5 and Rh6, two light-sensing G-protein coupled receptors involved in the daily resetting of circadian rhythm, it does not influence eye development. Circadian rhythm genes are well conserved from flies to humans; therefore, we tested the effect of DBP exposure on human breast cells (MCF10A) and demonstrate that, similar to the fruit fly model, this exposure disrupts circadian rhythm (BMAL1 expression) at doses that promote the proliferation and migration ability of MCF10A cells. Our results are the first to provide comprehensive evidence that DBP interferes with circadian rhythm in both adult Drosophila and human cells, which may help to explain the broad physiological action of phthalates.

Clinical application of whole transcriptome sequencing for the classification of patients with acute lymphoblastic leukemia

BACKGROUND

Considering the clinical and genetic characteristics, acute lymphoblastic leukemia (ALL) is a rather heterogeneous hematological neoplasm for which current standard diagnostics require various analyses encompassing morphology, immunophenotyping, cytogenetics, and molecular analysis of gene fusions and mutations. Hence, it would be desirable to rely on a technique and an analytical workflow that allows the simultaneous analysis and identification of all the genetic alterations in a single approach. Moreover, based on the results with standard methods, a significant amount of patients have no established abnormalities and hence, cannot further be stratified.

METHODS

We performed WTS and WGS in 279 acute lymphoblastic leukemia (ALL) patients (B-cell: n = 211; T-cell: n = 68) to assess the accuracy of WTS, to detect relevant genetic markers, and to classify ALL patients.

RESULTS

DNA and RNA-based genotyping was used to ensure correct WTS-WGS pairing. Gene expression analysis reliably assigned samples to the B Cell Precursor (BCP)-ALL or the T-ALL group. Subclassification of BCP-ALL samples was done progressively, assessing first the presence of chromosomal rearrangements by the means of fusion detection. Compared to the standard methods, 97% of the recurrent risk-stratifying fusions could be identified by WTS, assigning 76 samples to their respective entities. Additionally, read-through fusions (indicative of CDKN2A and RB1 gene deletions) were recurrently detected in the cohort along with 57 putative novel fusions, with yet untouched diagnostic potentials. Next, copy number variations were inferred from WTS data to identify relevant ploidy groups, classifying an additional of 31 samples. Lastly, gene expression profiling detected a BCR-ABL1-like signature in 27% of the remaining samples.

CONCLUSION

As a single assay, WTS allowed a precise genetic classification for the majority of BCP-ALL patients, and is superior to conventional methods in the cases which lack entity defining genetic abnormalities.

Transcriptomics paving the way for improved diagnostics and precision medicine of acute leukemia

Transcriptional profiling of acute leukemia, specifically by RNA-sequencing or whole transcriptome sequencing (WTS), has provided fundamental insights into its underlying disease biology and allows unbiased detection of oncogenic gene fusions, as well as of gene expression signatures that can be used for improved disease classification. While used as a research tool for many years, RNA-sequencing is becoming increasingly used in clinical diagnostics. Here, we highlight key transcriptomic studies of acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML) that have improved our biological understanding of these heterogeneous malignant disorders and have paved the way for translation into clinical diagnostics. Recent single-cell transcriptomic studies of ALL and AML, which provide new insights into the cellular ecosystem of acute leukemia and point to future clinical utility, are also reviewed. Finally, we discuss current challenges that need to be overcome for a more wide-spread adoption of RNA-sequencing in clinical diagnostics and how this technology significantly can aid the identification of genetic alterations in current guidelines and of newly emerging disease entities, some of which are critical to identify because of the availability of targeted therapies, thereby paving the way for improved precision medicine of acute leukemia.

Whole transcriptome RNA-Seq analysis reveals extensive cell type-specific compartmentalization in Volvox carteri

Background

One of evolution’s most important achievements is the development and radiation of multicellular organisms with different types of cells. Complex multicellularity has evolved several times in eukaryotes; yet, in most lineages, an investigation of its molecular background is considerably challenging since the transition occurred too far in the past and, in addition, these lineages evolved a large number of cell types. However, for volvocine green algae, such as Volvox carteri, multicellularity is a relatively recent innovation. Furthermore, V. carteri shows a complete division of labor between only two cell types – small, flagellated somatic cells and large, immotile reproductive cells. Thus, V. carteri provides a unique opportunity to study multicellularity and cellular differentiation at the molecular level.

Results

This study provides a whole transcriptome RNA-Seq analysis of separated cell types of the multicellular green alga V. carteri f. nagariensis to reveal cell type-specific components and functions. To this end, 246 million quality filtered reads were mapped to the genome and valid expression data were obtained for 93% of the 14,247 gene loci. In the subsequent search for protein domains with assigned molecular function, we identified 9435 previously classified domains in 44% of all gene loci. Furthermore, in 43% of all gene loci we identified 15,254 domains that are involved in biological processes. All identified domains were investigated regarding cell type-specific expression. Moreover, we provide further insight into the expression pattern of previously described gene families (e.g., pherophorin, extracellular matrix metalloprotease, and VARL families). Our results demonstrate an extensive compartmentalization of the transcriptome between cell types: More than half of all genes show a clear difference in expression between somatic and reproductive cells.

Conclusions

This study constitutes the first transcriptome-wide RNA-Seq analysis of separated cell types of V. carteri focusing on gene expression. The high degree of differential expression indicates a strong differentiation of cell types despite the fact that V. carteri diverged relatively recently from its unicellular relatives. Our expression dataset and the bioinformatic analyses provide the opportunity to further investigate and understand the mechanisms of cell type-specific expression and its transcriptional regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0450-y) contains supplementary material, which is available to authorized users.

Exposure to the environmental pollutant bisphenol A diglycidyl ether (BADGE) causes cell over-proliferation in Drosophila

Bisphenol A diglycidyl ether (BADGE), a derivative of bisphenol A (BPA), is widely used in the manufacture of epoxy resins as well as a coating on food containers. Recent studies have demonstrated the adverse effects of BADGE on reproduction and development in rodents and amphibians, but how BADGE affects biological activity is not understood. To gain a better understanding of the biological effects of BADGE exposure during development, we used the model organism Drosophila melanogaster and performed whole transcriptome sequencing. Interestingly, when Drosophila are raised on food containing BADGE, genes having significantly increased transcript numbers are enriched for those involved in regulating cell proliferation, including DNA replication and cell cycle control. Furthermore, raising larvae on BADGE-containing food induces hemocyte (blood cell) over-proliferation. This effect can be stimulated with even lower concentrations of BADGE if the hemocytes are already primed for cell proliferation by the expression of dominant active Ras GTPase. We conclude that chronic exposure to the xenobiotic BADGE throughout development can induce cell proliferation.

Overexpression of protocadherin 7 inhibits neuronal survival by downregulating BIRC5 in vitro

Protocadherins (Pcdhs) are widely-expressed transmembrane proteins in the nervous system. Recent studies suggest that Pcdhs play multiple critical roles during neuronal development. However, the cellular mechanisms of Pcdh7 in neurons are still largely unknown. In the current study, we demonstrated that the expression of Pcdh7 during mouse brain development was regulated spatiotemporally. We observed that the elevated expression of Pcdh7 led to activation of the intrinsic apoptotic pathway in primary cortical neurons. Whole transcriptome sequencing revealed that 12 genes were involved in the apoptotic pathway including baculoviral inhibitor of apoptosis (IAP) repeat containing 5 (BIRC5). The neuronal apoptosis caused by Pcdh7 overexpression could be significantly inhibited by either a missense mutation in the conserved motif CM2 domain of Pcdh7 or BIRC5 overexpression. These results suggest the existence of Pcdh7-BIRC5 signaling cascade in the cortical neurons and represent a potential therapeutic area for further investigation.

Whole transcriptome analyses and comparison reveal the metabolic differences between oxidative and glycolytic skeletal muscles of yak

Mammalian skeletal muscle is composed of various muscle fibers that exhibit different physiological and metabolic features. Muscle fiber type composition has significant influences on the meat quality of livestock. In this study, we comprehensively analyzed the whole transcriptome profiles of the oxidative muscle biceps femoris (BF) and the glycolytic muscle obliquus externus abdominis (OEA) of yak. A total of 1436 mRNAs, 1172 lncRNAs, and 218 circRNAs were differentially expressed in the oxidative muscles compared with the glycolytic muscles. KEGG annotation showed that differentially expressed mRNAs regulated by lncRNA and circRNA were mainly involved in PPAR signaling pathway, citrate cycle (TCA cycle), and PI3K-Akt signaling pathway, which reflect the different metabolic properties between oxidative and glycolytic muscles. In addition, regulatory networks associated with muscle fiber type conversion and mitochondria energy metabolism in muscles were constructed. Our study provides new evidence for a better understanding of the molecular mechanisms underlying skeletal muscle fiber determination and meat quality traits of yak.

Successful multidisciplinary clinical approach and molecular characterization by whole transcriptome sequencing of a cardiac myxofibrosarcoma: A case report

BACKGROUND

Cardiac tumors are rare and complex entities. Surgery represents the cornerstone of therapy, while the role of adjuvant treatment remains unclear and, in case of relapse or metastatic disease, the prognosis is very poor. Lack of prospective, randomized clinical trials hinders the generation of high level evidence for the optimal diagnostic workup and multimodal treatment of cardiac sarcomas. Herein, we describe the multidisciplinary clinical management and molecular characterization of a rare case of cardiac myxofibrosarcoma in an elderly woman.

CASE SUMMARY

A 73-year-old woman presented signs and symptoms of acute left-sided heart failure. Imaging examination revealed a large, left atrial mass. With suspicion of a myxoma, she underwent surgery, and symptoms were promptly relieved. Histology showed a cardiac myxofibrosarcoma, a rare histotype of cardiac sarcoma. Eight months later, disease unfortunately relapsed, and after a multidisciplinary discussion, a chemotherapy with doxorubicin and then gemcitabine was started, achieving partial radiologic and complete metabolic response, which was maintained up to 2 years and is still present. This report is focused on the entire clinical path of our patient from diagnosis to follow-up, through surgery and strategies adopted at relapse. Moreover, due to their rarity, very little is known about the molecular landscape of myxofibrosarcomas. Thus, we also performed and described preliminary genome analysis of the tumor tissue to get further insight on mechanisms involved in tumor growth, and to possibly unveil new clinically actionable targets.

CONCLUSION

We report a case of cardiac myxofibrosarcoma that achieved a very good prognosis due to an integrated surgical, cardiac and oncologic treatment strategy.

Novel ATXN1/ATXN1L::NUTM2A fusions identified in aggressive infant sarcomas with gene expression and methylation patterns similar to CIC-rearranged sarcoma

CIC-rearranged sarcomas are newly defined undifferentiated soft tissue tumors with CIC-associated fusions, and dismal prognosis. CIC fusions activate PEA3 family genes, ETV1/4/5, leading to tumorigenesis and progression. We report two high-grade CNS sarcomas of unclear histological diagnosis and one disseminated tumor of unknown origin with novel fusions and similar gene-expression/methylation patterns without CIC rearrangement. All three patients were infants with aggressive diseases, and two experienced rapid disease deterioration and death. Whole-transcriptome sequencing identified an ATXN1-NUTM2A fusion in the two CNS tumors and an ATXN1L-NUTM2A fusion in case 3. ETV1/4/5 and WT1 overexpression were observed in all three cases. Methylation analyses predicted CIC-rearranged sarcoma for all cases. Retrospective IHC staining on case 2 demonstrated ETV4 and WT1 overexpression. ATXN1 and ATXN1L interact with CIC forming a transcription repressor complex. We propose that ATXN1/ATXN1L-associated fusions disrupt their interaction with CIC and decrease the transcription repressor complex, leading to downstream PEA3 family gene overexpression. These three cases with novel ATXN1/ATXN1L-associated fusions and features of CIC-rearranged sarcomas may further expand the scope of "CIC-rearranged" sarcomas to include non-CIC rearrangements. Additional cases are needed to demonstrate if ATXN1/ATXN1L-NUTM2A fusions are associated with younger age and more aggressive diseases.

Identification and comparison of novel circular RNAs with associated co-expression and competing endogenous RNA networks in postmenopausal osteoporosis

BACKGROUND

Circular RNAs (circRNAs) are emerging as crucial regulators in various human diseases. So far, the expression profile and regulatory mechanism of circRNAs in postmenopausal osteoporosis (PMOP) are less studied and should be deciphered urgently. Herein, we aimed to reveal key circRNAs affecting PMOP and clarify their compounding regulatory actions.

METHODS

To reveal key circRNAs affecting PMOP and clarify their compounding regulatory actions, whole transcriptome sequencing and bioinformatics analysis were performed to identify differentially expressed circRNAs (DECs). The expression pattern and regulatory networks of DECs in peripheral blood mononuclear cells (PBMCs) were unearthed.

RESULTS

A total of 373 DECs comprising 123 intronic, 100 antisense, 70 exonic, 55 intergenic, and 25 sense-overlapping circRNAs were identified. Among these, 73 circRNAs were upregulated and 300 were downregulated. These DECs exerted pivotal functions in the pathogenesis of PMOP as demonstrated by Gene Ontology (GO) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The circRNA-miRNA-mRNA co-expression network comprising 28 DECs, 145 miRNAs, and 175 differentially expressed mRNAs predicted the possible mechanism of the pathogenesis and progression of PMOP.

CONCLUSION

The results of the present study provided a further comprehension of circRNA-associated competing endogenous RNA regulatory mechanism in PMOP. The steadily expressed and disease-specific DECs may serve as promising diagnostic and prognostic biomarkers for PMOP.

Gene Fusion Discovery with INTEGRATE

Next-generation sequencing (NGS) has become the primary technology for discovering gene fusions. Decreasing NGS costs have resulted in a growing quantity of patients with whole transcriptome sequencing (RNA-seq) and whole genome sequencing (WGS) data. We developed a gene fusion discovery tool, INTEGRATE, that leverages both RNA-seq and WGS data to reconstruct gene fusion junctions and genomic breakpoints by split-read alignment. INTEGRATE has become widely adopted by the larger cancer research community to discover biologically and clinically relevant gene fusions. Here we explain the rationale driving the development of the INTEGRATE tool and describe the detailed practical procedures for applying INTEGRATE to discover gene fusions using NGS data. INTEGRATE can be applied to both combined data and RNA-seq only data.

Identification of Novel CDH1-NRG2α and F11R-NRG2α Fusions in NSCLC Plus Additional Novel NRG2α Fusions in Other Solid Tumors by Whole Transcriptome Sequencing

Introduction

A novel CD74-NRG2α fusion has recently been identified in NSCLC. We surveyed a large tumor database comprehensively profiled by whole transcriptome sequencing to investigate the incidence and distribution of NRG2 fusions among various solid tumors.

Methods

Tumor samples submitted for clinical molecular profiling at Caris Life Sciences (Phoenix, AZ) that underwent whole transcriptome sequencing (NovaSeq [Illumina, San Diego, CA]) were retrospectively analyzed for NRG2 fusion events. All NRG2 fusions with sufficient reads (> three junctional reads spanning ≥ seven nucleotides) were identified for manual review, characterization of fusion class, intact functional domains, EGF-like domain isoforms, breakpoints, frame retention, and co-occurring alterations by next-generation sequencing (NextSeq [Illumina, San Diego, CA], 592 genes).

Results

Seven inframe functional (containing the intact EGF-like domain) NRG2α fusions were identified, namely, the following: (1) NSCLC (two of 9600, 0.02%: CDH1-NRG2α [C11, N2], F11R-NRG2α [F1, N4]); (2) endometrial (two of 3060, 0.065%: CPM-NRG2α [C2, N2], OPA3-NRG2α [O1, N2]); (3) ovarian (one of 5030, 0.02%: SPON1-NRG2α [S6, N2]); (4) prostate (one of 1600, 0.063%: PLPP1-NRG2α [P1, N2]); and (5) carcinoma of unknown origin (one of 1400, 0.07%: CYSTM1-NRG2α [C2, N2]). No NRG2β fusions were identified. Both NSCLC samples contained the reciprocal NRG2 fusions (NRG2-CDH1, NRG2-F11R). Almost all inframe NRG2α fusions have no (N = 6, 85.7%) or low (N = 1, 14.3%) programmed death-ligand 1 expression. No additional known driver mutations were identified in these seven NRG2α fusion-positive tumor samples.

Conclusions

Similar to NRG1 fusions, NRG2α fusions are recurrent and rare ligand-fusions in NSCLC and other multiple tumor types, especially gynecologic malignancies.

Whole transcriptome sequencing identifies a competitive endogenous RNA network that regulates the immunity of bladder cancer

Several types of non-coding RNAs such as circRNAs, lncRNAs, and miRNAs have been identified to regulate mRNAs through the mechanism known as the competitive endogenous RNA (ceRNA) network. To explore the role of the ceRNA regulatory network in the immune microenvironment of bladder cancer, whole-transcriptome sequencing of bladder tumor and its peritumoral tissues from 38 bladder cancer patients, with a total of 63 samples, was performed to screen differentially expressed circ-, lnc-, mi-, and mRNAs to construct a circ/lnc-mi-mRNA regulatory network with pruning algorithms. We excavated a key immune-related gene BDNF to build the final ceRNA network as hsa-miR-107 sponged by hsa-circ-000211, AC108488.1, and LINC00163. Finally, a meta-analysis of 7 public datasets demonstrated that low expression of BDNF and high expression of hsa-miR-107 were associated with longer survival. Our study identified a ceRNA regulatory network as a potentially new prognostic marker and molecular therapeutic target of bladder cancer.

Key circular RNAs identified in male osteoporosis patients by whole transcriptome sequencing

Background

Osteoporosis (OP) is a systemic disease with bone loss and microstructural deterioration. Numerous noncoding RNAs (ncRNAs) have been proved to participate in various diseases, especially circular RNAs (circRNAs). However, the expression profile and mechanisms underlying circRNAs in male osteoporosis have not yet been explored.

Methods

The whole transcriptome expression profile and differences in mRNAs, circRNAs, and microRNAs (miRNAs) were investigated in peripheral blood samples of patients with osteoporosis and healthy controls consisting of males ≥ 60-years-old.

Results

A total of 398 circRNAs, 51 miRNAs, and 642 mRNAs were significantly and differentially expressed in osteoporosis compared to healthy controls. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the host genes of significantly differentially expressed circRNAs were mainly enriched in the regulation of cell cycle process: biological process (BP), organelle part cellular components (CC), protein binding molecular function (MF), Toll-like receptor signaling pathway, tumor necrosis factor (TNF) signaling pathway, and thyroid hormone signaling pathway. circRNA-miRNA-mRNA regulatory network was constructed using the differentially expressed RNAs. Moreover, key circRNAs (hsa_circ_0042409) in osteoporosis were discovered and validated by qPCR.

Conclusions

The key cicrRNAs plays a major role in the pathogenesis of osteoporosis and could be used as potential biomarkers or targets in the diagnosis and treatment of osteoporosis.

Noncoding RNA-associated competing endogenous RNA networks in trastuzumab-induced cardiotoxicity

Trastuzumab-induced cardiotoxicity (TIC) is a common and serious disease with abnormal cardiac function. Accumulating evidence has indicated certain non-coding RNAs (ncRNAs), functioning as competing endogenous RNAs (ceRNAs), impacting the progression of cardiovascular diseases. Nonetheless, the specific involvement of ncRNA-mediated ceRNA regulatory mechanisms in TIC remains elusive. The present research aims to comprehensively investigate changes in the expressions of all ncRNA using whole-transcriptome RNA sequencing. The sequencing analysis unveiled significant dysregulation, identifying a total of 43 circular RNAs (circRNAs), 270 long noncoding RNAs (lncRNAs), 12 microRNAs (miRNAs), and 4131 mRNAs in trastuzumab-treated mouse hearts. Subsequently, circRNA-based ceRNA networks consisting of 82 nodes and 91 edges, as well as lncRNA-based ceRNA networks comprising 111 nodes and 112 edges, were constructed. Using the CytoNCA plugin, pivotal genes-miR-31-5p and miR-644-5p-were identified within these networks, exhibiting potential relevance in TIC treatment. Additionally, KEGG and GO analyses were conducted to explore the functional pathways associated with the genes within the ceRNA networks. The outcomes of the predicted ceRNAs and bioinformatics analyses elucidated the plausible involvement of ncRNAs in TIC pathogenesis. This insight contributes to a better understanding of underlying mechanisms and aids in identifying promising targets for effective prevention and treatment strategies.

Whole transcriptome analysis of RNA expression profiles reveals the potential regulating action of long noncoding RNA in lactating cows fed a high concentrate diet

Subacute ruminal acidosis (SARA) is a common metabolic disease in the dairy farming industry which is usually caused by an excessive amount of high concentrate diet. SARA not only threatens animal welfare but also leads to economic losses in the farming industry. The liver plays an important role in the distribution of nutritional substances and metabolism; however, a high concentrate diet can cause hepatic metabolic disorders and liver injury. Recently, noncoding RNA has been considered as a critical regulator of hepatic disease, however, its role in the bovine liver is limited. In this study, 12 mid-lactating dairy cows were randomly assigned to a control (CON) group (40% concentrate of dry matter, n = 6) and a SARA group (60% concentrate of dry matter, n = 6). After 21 d of treatment, all cows were sacrificed, and liver tissue samples were collected. Three dairy cows were randomly selected from the CON and SARA groups respectively to perform whole transcriptome analysis. More than 20,000 messenger RNA (mRNA), 10,000 long noncoding RNA (lncRNA), 3,500 circular RNA (circRNA) and 1,000 micro RNA (miRNA) were identified. Furthermore, 43 mRNA, 121 lncRNA and 3 miRNA were differentially expressed, whereas no obvious differentially expressed circRNA were detected between the 2 groups. Gene Ontology (GO) annotation revealed that the differentially expressed genes were mainly enriched in oxidoreductase activity, stress, metabolism, the immune response, cell apoptosis, and cell proliferation. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the deferentially expressed genes were highly enriched in the phosphatidylinositol 3 kinase (PI3K)-serine/threonine kinase (AKT) signaling pathway (P < 0.05). According to KEGG pathway analysis, the differentially expressed lncRNA (DElncRNA) target genes were mainly related to proteasomes, peroxisomes, and the hypoxia-inducible factor-1 signaling pathway (P < 0.005). Further bioinformatics and integrative analyses revealed that the lncRNA were strongly correlated with mRNA; therefore, it is reasonable to speculate that lncRNA potentially play important roles in the liver dysfunction induced by SARA. Our study provides a valuable resource for future investigations on the mechanisms of SARA to facilitate an understanding of the importance of lncRNA, and offer functional RNA information.

High throughput sequencing of whole transcriptome and construct of ceRNA regulatory network in RD cells infected with enterovirus D68

Background

With the advancement of sequencing technologies, a plethora of noncoding RNA (ncRNA) species have been widely discovered, including microRNAs (miRNAs), circular RNAs (circRNAs), and long ncRNAs (lncRNAs). However, the mechanism of these non-coding RNAs in diseases caused by enterovirus d68 (EV-D68) remains unclear. The goal of this research was to identify significantly altered circRNAs, lncRNAs, miRNAs, and mRNAs pathways in RD cells infected with EV-D68, analyze their target relationships, demonstrate the competing endogenous RNA (ceRNA) regulatory network, and evaluate their biological functions.

Methods

The total RNAs were sequenced by high-throughput sequencing technology, and differentially expressed genes between control and infection groups were screened using bioinformatics method. We discovered the targeting relationship between three ncRNAs and mRNA using bioinformatics methods, and then built a ceRNA regulatory network centered on miRNA. The biological functions of differentially expressed mRNAs (DEmRNAs) were discovered through GO and KEGG enrichment analysis. Create a protein interaction network (PPI) to seek for hub mRNAs and learn more about protein–protein interactions. The relative expression was verified using RT-qPCR. The effects of Fos and ARRDC3 on virus replication were confirmed using RT-qPCR, virus titer (TCID₅₀/ml), Western blotting.

Results

375 lncRNAs (154 upregulated and 221 downregulated), 33 circRNAs (32 upregulated and 1 downregulated), 96 miRNAs (49 upregulated and 47 downregulated), and 239 mRNAs (135 upregulated and 104 downregulated) were identified as differently in infected group compare to no-infected group. A single lncRNA or circRNA can be connected with numerous miRNAs, which subsequently coregulate additional mRNAs, according to the ceRNA regulatory network. The majority of DEmRNAs were shown to be connected to DNA binding, transcription regulation by RNA polymerase II, transcription factor, MAPK signaling pathways, Hippo signal pathway, and apoptosis pathway, according to GO and KEGG pathway enrichment analysis. The hub mRNAs with EGR1, Fos and Jun as the core were screened through PPI interaction network. We preliminarily demonstrated that the Fos and ARRDC3 genes can suppress EV-D68 viral replication in order to further verify the results of full transcriptome sequencing.

Conclusion

The results of whole transcriptome analysis after EV-D68 infection of RD cells were first reported in this study, and for the first time, a ceRNA regulation network containing miRNA at its center was established for the first time. The Fos and ARRDC3 genes were found to hinder viral in RD cells. This study establishes a novel insight host response during EV-D68 infection and further investigated potential drug targets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01686-x.

Whole transcriptome sequencing analysis of blood plasma-derived exosomes from immune-related hearing loss

BACKGROUND

Early diagnosis of immune-related hearing loss and timely treatment can prevent structural damage to the inner ear and contribute to hearing retention. Exosomal miRNAs, lncRNAs and proteins have great prospects as novel biomarkers for clinical diagnosis. Our study aimed to investigate the molecular mechanisms of exosomes or exosomal ceRNA regulatory networks in immune-related hearing loss.

METHODS

An immune-related hearing loss mice model was constructed by injection with inner ear antigen, and then the blood plasma samples of the mice were collected for exosomes isolation by ultra-centrifugation. Subsequently, the different exosomes were sent for whole transcriptome sequencing using Illumina platform. Finally, a ceRNA pair was chosen for validation by RT-qPCR and dual luciferase reporter gene assay.

RESULTS

The exosomes were successfully extracted from the blood samples of the control and the immune-related hearing loss mice. After sequencing, 94 differentially expressed (DE) lncRNAs, 612 DEmRNAs, and 100 DEmiRNAs were found in the immune-related hearing loss-associated exosomes. Afterwards, ceRNA regulatory networks consisting of 74 lncRNAs, 28 miRNAs and 256 mRNAs were proposed, and the genes in the ceRNA regulatory networks were significantly enriched in 34 GO terms of biological processes and 9 KEGG pathways. Finally, Gm9866 and Dusp7 were significantly up-regulated, while miR-185-5p level was declined in the exosomes from immune-related hearing loss, and Gm9866, miR-185-5p and Dusp7 interacted with each other.

CONCLUSIONS

Gm9866-miR-185-5p-Dusp7 was confirmed to be closely correlated with the occurrence and progression of immune-related hearing loss.

Novel Molecular Biomarkers to Guide Treatment Decision-making in Metastatic Urothelial Cancer-A Patient Cohort Analysis

The current options and recent developments in the field of systemic therapy for advanced urothelial cancer (UC) patients urge the need for selection criteria to identify the most optimal therapeutic option for individual patients. The molecular makeup of tumors, including molecular subtype, tumor microenvironment composition, and gene mutations, fusions, and amplifications, has previously been correlated with a response to immune checkpoint inhibitors, erdafitinib, or enfortumab vedotin (EV) monotherapy, and may withhold potential candidate biomarkers. In this study, we aimed to stratify metastatic UC (mUC) patients based on molecular biomarkers that might be associated with a response to EV, a fibroblast growth factor receptor inhibitor, or anti-PD-(L)1, by using whole-genome DNA-sequencing and paired RNA-sequencing data of fresh-frozen metastatic tumor biopsies of 155 mUC patients. We observed that NECTIN4 amplification, FGFR2/3 mutations, and the RNA expression-based T-cell-to-stroma enrichment (TSE) score were mutually exclusive, and may therefore reflect biologically distinct tumors and sensitivity to treatments. This finding was validated in two independent bladder cohorts: the IMvigor210 study and The Cancer Genome Atlas. Stratification of patients into subgroups based on these molecular features is possible. Our data challenge the concept of a one-treatment-fits-all paradigm and support the rationale for prospective clinical trials with biomarker-guided treatment selection of mUC patients.

Whole transcriptome analysis and construction of gene regulatory networks of granulosa cells from patients with polycystic ovary syndrome (PCOS)

OBJECTIVE

Polycystic ovary syndrome (PCOS) is a reproductive endocrine disease characterized by reproductive dysfunction and metabolic abnormalities. The purpose of this study was to explore the expression characteristics of coding and non-coding RNAs in granulosa cells of PCOS, and to provide data support for understanding the pathogenesis of PCOS.

METHODS

Three patients with PCOS (according to the 2003 Rotterdam diagnostic criteria) and three normal controls were selected. We used the standard long protocol to collect granulosa cells from two groups, who underwent assisted reproduction at the Reproductive Medicine Center of the Affiliated Hospital of Inner Mongolia Medical University, China. We performed whole-transcriptome sequencing using RNA-Seq technology to construct transcriptome patterns of messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs). These patterns were then subjected to in-depth analysis using bioinformatics tools.

RESULTS

We identified a total of 2111 mRNAs and 4328 non-coding RNAs (ncRNAs) in the PCOS group as compared with the control group. Among the ncRNAs, there were 2047 lncRNAs, 892 circRNAs, and 1389 miRNAs. Based on the condition |log2(fold_change) |≥ 1 and a P-value of ≤ 0.05, we obtained 705 differentially expressed genes (DEGs), 204 differentially expressed lncRNAs, 111 differentially expressed circRNAs, and 88 differentially expressed miRNAs. The target genes were mainly enriched in metabolic pathways such as mitogen-activated protein kinase (MAPK), Wnt, transforming growth factor-beta (TGF-β), and the cell cycle. There were three types of circRNAs, among which the number of exon-type circRNAs accounted for more than 90%. Using co-expression network analysis, we identified several important candidate gene mRNAs (VLDLR, PPP2R2B, and MYOCD), lncRNAs (FBXO30, SNHG14, and PVT1), and miRNAs (miRNA-150); these mRNAs and ncRNAs could play a regulatory role in PCOS granulosa cells.

CONCLUSION

In this study, we discovered significant alterations in mRNAs, lncRNAs, circRNAs, and miRNAs in PCOS granulosa cells, indicating dysregulation in vital pathways. Notably, genes like VLDLR, PPP2R2B, and MYOCD, along with lncRNAs FBXO30, SNHG14, and PVT1, may contribute to PCOS pathology, shedding light on potential therapeutic targets.

Characterization of the regulatory network and pathways in duodenum affecting chicken abdominal fat deposition

The excessive accumulation of abdominal fat in chickens has resulted in a reduction in both the feed conversion efficiency and the slaughter yield. To elucidate the regulatory mechanisms and metabolic pathways affecting abdominal fat deposition in the context of broiler breeding, a cohort of 400 Qingyuan partridge chickens with varying abdominal fat deposition was established. Whole transcriptome sequencing analyses were conducted on the duodenum of 20 representative chickens to ascertain the regulatory networks at this vital digestive and absorptive organ. Consequently, 116 differentially expressed genes were identified, exhibiting a trend of increasing or decreasing expression in correlation with the accumulation of abdominal fat. A total of 36 DEmRNAs, 170 DElncRNAs, 92 DEcircRNAs and 88 DEmiRNAs were identified as differentially expressed between chickens with extremely high and low abdominal fat deposition. The functional enrichment analyses demonstrated that the differentially expressed RNA in the duodenum were involved in the regulation of chicken abdominal fat deposition by mediating a series of metabolic pathways, including the Wnt signaling pathway, the PPAR signaling pathway, the Hippo signaling pathway, the FoxO signaling pathway, the MAPK signaling pathway and other signaling pathways that are involved in fatty acid metabolism and degradation. The construction of putative interaction pairs led to the suggestion of two lncRNA-miRNA-mRNA ceRNA networks comprising two mRNAs, two miRNAs, and 29 lncRNAs, as well as two circRNA-lncRNA-miRNA-mRNA ceRNA networks comprising 26 mRNAs, 12 miRNAs, 17 lncRNAs, and nine circRNAs, as core regulatory networks in the duodenum affecting chicken abdominal fat deposition. The aforementioned genes including TMEM150C, REXO1, PIK3C2G, ppp1cb, PARP12, SERPINE2, LRAT, CYP1A1, INSR and APOA4, were proposed as candidate genes, while the miRNAs, including miR-107-y, miR-22-y, miR-25-y, miR-2404-x and miR-16-x, as well as lncRNAs such as ENSGALT00000100291, TCONS_00063508, TCONS_00061201 and TCONS_00079402 were the candidate regulators associated with chicken abdominal fat deposition. The findings of this study provide a theoretical foundation for the molecular mechanisms of mRNAs and non-coding RNAs in duodenal tissues on abdominal fat deposition in chickens.