scholarly journals Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

2020 ◽  
Author(s):  
Yue Zhang ◽  
Razgar Seyed Rahmani ◽  
Xingyu Yang ◽  
Jinming Chen ◽  
Tao Shi

Abstract Background Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus ( Nelumbo nucifera ), an economically important flowering plant. Results Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets, and tend to be tissue-specific. Further, our results showed that frequent interactions between miRNAs and hub gene isoforms highlight the critical roles of miRNAs in the gene network. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes. Conclusions Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

2020 ◽  
Author(s):  
Yue Zhang ◽  
Razgar Seyed Rahmani ◽  
Xingyu Yang ◽  
Jinming Chen ◽  
Tao Shi

Abstract Background: Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.Results: Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.Conclusions: Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.


2019 ◽  
Author(s):  
Yue Zhang ◽  
Razgar Seyed Rahmani ◽  
Xingyu Yang ◽  
Jinming Chen ◽  
Tao Shi

Abstract BACKGROUND. Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, coordination of interaction between miRNAs with different splicing isoforms, and the role of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus ( Nelumbo nucifera ), an economically important flowering plant. RESULTS. Through RNA-seq expression analyses among six lotus tissues, the negative regulatory roles of most miRNAs are reflected by their tissue-biased expression and the negative correlation with their targets in expression. Further, the central roles of miRNAs in the gene network was unveiled as there are more frequent interactions between miRNAs and hub isoforms than between miRNAs and non-hub isoforms. Surprisingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. The gene function enrichment analysis of miRNA targeted reveals that miRNAs are involved in regulation of lotus growth and development by regulating plant hormone-related pathway genes. CONCLUSION. Taken together, we carry out a comprehensive and deep analysis between miRNA and mRNA transcriptome to study coordination of interaction between miRNAs with different splicing isoforms. Our study on lotus highlights not only the complicate interactions between the miRNAs and transcript isoforms but also functional divergence of many transcript isoforms from the same locus in plant.


2021 ◽  
Author(s):  
Parastoo Karimi ◽  
Mohammad Reza Bakhtiarizadeh ◽  
Abdolreza Salehi ◽  
Hamid Reza Izadnia

Abstract Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs.In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.


2019 ◽  
Vol 14 (7) ◽  
pp. 591-601 ◽  
Author(s):  
Aravind K. Konda ◽  
Parasappa R. Sabale ◽  
Khela R. Soren ◽  
Shanmugavadivel P. Subramaniam ◽  
Pallavi Singh ◽  
...  

Background: Chickpea is a nutritional rich premier pulse crop but its production encounters setbacks due to various stresses and understanding of molecular mechanisms can be ascribed foremost importance. Objective: The investigation was carried out to identify the differentially expressed WRKY TFs in chickpea in response to herbicide stress and decipher their interacting partners. Methods: For this purpose, transcriptome wide identification of WRKY TFs in chickpea was done. Behavior of the differentially expressed TFs was compared between other stress conditions. Orthology based cofunctional gene networks were derived from Arabidopsis. Gene ontology and functional enrichment analysis was performed using Blast2GO and STRING software. Gene Coexpression Network (GCN) was constructed in chickpea using publicly available transcriptome data. Expression pattern of the identified gene network was studied in chickpea-Fusarium interactions. Results: A unique WRKY TF (Ca_08086) was found to be significantly (q value = 0.02) upregulated not only under herbicide stress but also in other stresses. Co-functional network of 14 genes, namely Ca_08086, Ca_19657, Ca_01317, Ca_20172, Ca_12226, Ca_15326, Ca_04218, Ca_07256, Ca_14620, Ca_12474, Ca_11595, Ca_15291, Ca_11762 and Ca_03543 were identified. GCN revealed 95 hub genes based on the significant probability scores. Functional annotation indicated role in callose deposition and response to chitin. Interestingly, contrasting expression pattern of the 14 network genes was observed in wilt resistant and susceptible chickpea genotypes, infected with Fusarium. Conclusion: This is the first report of identification of a multi-stress responsive WRKY TF and its associated GCN in chickpea.


2020 ◽  
Vol 2020 ◽  
pp. 1-8
Author(s):  
Yuntao Shi ◽  
Yingying Zhuang ◽  
Jialing Zhang ◽  
Mengxue Chen ◽  
Shangnong Wu

Objective. Although noncoding RNAs, especially the microRNAs, have been found to play key roles in CRC development in intestinal tissue, the specific mechanism of these microRNAs has not been fully understood. Methods. GEO and TCGA database were used to explore the microRNA expression profiles of normal mucosa, adenoma, and carcinoma. And the differential expression genes were selected. Computationally, we built the SVM model and multivariable Cox regression model to evaluate the performance of tumorigenic microRNAs in discriminating the adenomas from normal tissues and risk prediction. Results. In this study, we identified 20 miRNA biomarkers dysregulated in the colon adenomas. The functional enrichment analysis showed that MAPK activity and MAPK cascade were highly enriched by these tumorigenic microRNAs. We also investigated the target genes of the tumorigenic microRNAs. Eleven genes, including PIGF, TPI1, KLF4, RARS, PCBP2, EIF5A, HK2, RAVER2, HMGN1, MAPK6, and NDUFA2, were identified to be frequently targeted by the tumorigenic microRNAs. The high AUC value and distinct overall survival rates between the two risk groups suggested that these tumorigenic microRNAs had the potential of diagnostic and prognostic value in CRC. Conclusions. The present study revealed possible mechanisms and pathways that may contribute to tumorigenesis of CRC, which could not only be used as CRC early detection biomarkers, but also be useful for tumorigenesis mechanism studies.


Genes ◽  
2018 ◽  
Vol 9 (11) ◽  
pp. 545 ◽  
Author(s):  
Wei Wu ◽  
Lingxiang Wu ◽  
Mengyan Zhu ◽  
Ziyu Wang ◽  
Min Wu ◽  
...  

Somatic mutations in 3′-untranslated regions (3′UTR) do not alter amino acids and are considered to be silent in cancers. We found that such mutations can promote tumor progression by altering microRNA (miRNA) targeting efficiency and consequently affecting miRNA–mRNA interactions. We identified 67,159 somatic mutations located in the 3′UTRs of messenger RNAs (mRNAs) which can alter miRNA–mRNA interactions (functional somatic mutations, funcMutations), and 69.3% of these funcMutations (the degree of energy change > 12 kcal/mol) were identified to significantly promote loss of miRNA-mRNA binding. By integrating mRNA expression profiles of 21 cancer types, we found that the expression of target genes was positively correlated with the loss of absolute affinity level and negatively correlated with the gain of absolute affinity level. Functional enrichment analysis revealed that genes carrying funcMutations were significantly enriched in the MAPK and WNT signaling pathways, and analysis of regulatory modules identified eighteen miRNA modules involved with similar cellular functions. Our findings elucidate a complex relationship between miRNA, mRNA, and mutations, and suggest that 3′UTR mutations may play an important role in tumor development.


2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Meng Liu ◽  
Xia Li ◽  
Rui Fan ◽  
Xinhua Liu ◽  
Ju Wang

Nicotine, as the major psychoactive component of tobacco, has broad physiological effects within the central nervous system, but our understanding of the molecular mechanism underlying its neuronal effects remains incomplete. In this study, we performed a systematic analysis on a set of nicotine addiction-related genes to explore their characteristics at network levels. We found that NAGenes tended to have a more moderate degree and weaker clustering coefficient and to be less central in the network compared to alcohol addiction-related genes or cancer genes. Further, clustering of these genes resulted in six clusters with themes in synaptic transmission, signal transduction, metabolic process, and apoptosis, which provided an intuitional view on the major molecular functions of the genes. Moreover, functional enrichment analysis revealed that neurodevelopment, neurotransmission activity, and metabolism related biological processes were involved in nicotine addiction. In summary, by analyzing the overall characteristics of the nicotine addiction related genes, this study provided valuable information for understanding the molecular mechanisms underlying nicotine addiction.


2020 ◽  
Author(s):  
Yiyuan Zhang ◽  
Rongguo Yu ◽  
Jiayu Zhang ◽  
Eryou Feng ◽  
Haiyang Wang ◽  
...  

Abstract BackgroundOsteoarthritis (OA) is a common chronic disease worldwide. Subchondral bone is an important pathological change in OA and responds more rapidly to adverse loading and events compared to cartilage. However, the pathogenic genes and pathways of subchondral bone are largely unclear.ObjectiveThis study aimed to identify signature differences in genes involved in knee lateral tibial (LT) and medial tibial (MT) plateaus of subchondral bone tissue while exploring their potential molecular mechanisms via bioinformatics analysis.MethodsFirst, the gene expression data of GSE51588 was downloaded from the GEO database. Differentially expressed genes (DEGs) between knee LT and MT were identified, and functional enrichment analyses were performed. Then, a protein-protein interactive network was constructed in order to acquire the hub genes, and modules analysis was conducted using STRING and Cytoscape for further analysis. The enriched hub genes were queried in DGIdb database to find suitable drug candidates in OA.ResultsA total of 202 DEGs (112 upregulated genes and 84 downregulated genes) were determined. In the PPI network, ten hub genes were identified. Five significant modules were identified using the MCODE plugin unit. Functional enrichment analysis revealed the most important signaling pathways. Six of the ten hub genes were targetable by a total of 35 drugs, suggesting their possible therapeutic use for OA .ConclusionsThe identified hub genes and functional enrichment pathways were implicated in the development and progression of subchondral bone in OA, thus improving our understanding of OA and offering molecular targets for future therapeutic modalities.


2021 ◽  
Author(s):  
Nana Yang ◽  
Qianghua Wang ◽  
Biao Ding ◽  
Yinging Gong ◽  
Yue Wu ◽  
...  

Abstract Background: The accumulation of ROS resulting from upregulated levels of oxidative stress is commonly implicated in preeclampsia (PE). Ferroptosis is a novel form of iron-dependent cell death instigated by lipid peroxidation likely plays important role in PE pathogenesis. This study aims to investigate expression profiles and functions of the ferroptosis-related genes (FRGs) in early- and late-onset preeclampsia.Methods: The gene expression data and clinical information were downloaded from GEO database. The “limma” R package was used for screening differentially expressed genes. GO(Gene Ontology), Kyoto Encyclopedia of Genes and Genomes(KEGG) and protein protein interaction (PPI) network analyses were conducted to investigate the bioinformatics functions and molecular interactions of significantly different FRGs. Quantitative real-time reverse transcriptase PCR was used to verify the expression of hub FRGs in PE.Results: A total number of 4,215 DEGs were identified between EOPE and preterm cases and 3,356 DEGs were found between EOPE and LOPE subtypes. 20 significantly different FRGs were identified in EOPE, while only 3 in LOPE. Functional enrichment analysis revealed that the differentially expressed FRGs was mainly involved in EOPE and enriched in hypoxia- and iron-related pathways, such as response to hypoxia, iron homeostasis and iron ion binding process. The PPI network analysis and verification by RT-qPCR resulted in the identification of the following six interesting FRGs: FTH1, HIF1A, FTL, IREB2, MAPK8 and PLIN2. Conclusions: EOPE and LOPE owned distinct underlying molecular mechanisms and ferroptosis may be mainly implicated in pathogenesis of EOPE. Further studies are necessary for deeper inquiry into placental ferroptosis and its role in the pathogenesis of EOPE.


2019 ◽  
Vol 3 (s1) ◽  
pp. 111-111 ◽  
Author(s):  
Andrea Comba ◽  
Patrick Dunn ◽  
Anna E Argento ◽  
Padma Kadiyala ◽  
Sebastien Motsch ◽  
...  

OBJECTIVES/SPECIFIC AIMS: Oncostreams represent a novel growth pattern of GBM. In this study we uncovered the cellular and molecular mechanism that regulates the oncostreams function in GBM growth and invasion. METHODS/STUDY POPULATION: We studied oncostreams organization and function using genetically engineered mouse gliomas models (GEMM), mouse primary patient derived GBM model and human glioma biopsies. We evaluated the molecular landscape of oncostreams by laser capture microdissection (LCM) followed by RNA-Sequencing and bioinformatics analysis. RESULTS/ANTICIPATED RESULTS: Oncostreams are multicellular structures of 10-20 cells wide and 2-400 μm long. They are distributed throughout the tumors in mouse and human GBM. Oncostreams are heterogeneous structures positive for GFAP, Nestin, Olig2 and Iba1 cells and negative for Neurofilament. Using GEMM we found a negative correlation between oncostream density and animal survival. Moreover, examination of patient’s glioma biopsies evidenced that oncostreams are present in high grade but no in low grade gliomas. This suggests that oncostreams may play a role in tumor malignancy. Our data also indicated that oncostreams aid local invasion of normal brain. Transcriptome analysis of oncostreams revealed 43 differentially expressed (DE) genes. Functional enrichment analysis of DE genes showed that “collagen catabolic processes”, “positive regulation of cell migration”, and “extracellular matrix organization” were the most over-represented GO biological process. Network analysis indicated that Col1a1, ACTA2, MMP9 and MMP10 are primary target genes. These genes were also overexpressed in more malignant tumors (WT-IDH) compared to the less malignant (IDH1- R132H) tumors. Confocal time lapse imagining of 3D tumor slices demonstrated that oncostreams display a collective motion pattern within gliomas that has not been seen before. DISCUSSION/SIGNIFICANCE OF IMPACT: In summary, oncostreams are anatomically and molecularly distinctive, regulate glioma growth and invasion, display collective motion and are regulated by the extracellular matrix. We propose oncostreams as novel pathological markers valuable for diagnosis, prognosis and designing therapeutics for GBM patients.


Sign in / Sign up

Export Citation Format

Share Document