Functional Enrichment Analysis of Deregulated Long Non-Coding RNAs in Cancer Based on their Genomic Neighbors

The discovery of microRNAs (miRNAs) has fundamentally transformed our understanding of gene regulation. The competing endogenous RNA (ceRNA) hypothesis postulates that not only messenger RNAs but also other RNA transcripts, such as long non-coding RNAs and pseudogenes, can act as natural miRNA sponges. These RNAs influence each other’s expression levels by competing for the same pool of miRNAs through miRNA response elements on their target transcripts, thereby modulating gene expression and protein activity. In recent years, these ceRNA regulatory networks have gained considerable attention in cancer research. Several studies have identified cancer-specific ceRNA networks. Nevertheless, prior bioinformatic analyses have focused on long non-coding RNAs-associated ceRNA networks. Here, we identify an extended-ceRNA network (including both long non-coding RNAs and pseudogenes) shared across a group of four hormone-dependent (HD) cancers, i.e., prostate, breast, colorectal, and endometrial cancers, using data from The Cancer Genome Atlas (TCGA). We performed a functional enrichment analysis for differentially expressed genes in the shared ceRNA network of HD cancers, followed by a survival analysis to determine their prognostic ability. We identified two long non-coding RNAs, nine genes, and seventy-four miRNAs in the shared ceRNA network across four HD cancers. Among them, two genes and forty-one miRNAs were associated with at least one HD cancer survival. This study is the first to investigate pseudogene associated ceRNAs across a group of related cancers and highlights the value of this approach to understanding shared molecular pathogenesis in a group of related diseases.

Download Full-text

Comparative Transcriptomic Analysis to Identify the Important Coding and Non-coding RNAs Involved in the Pathogenesis of Pterygium

Frontiers in Genetics ◽

10.3389/fgene.2021.646550 ◽

2021 ◽

Vol 12 ◽

Author(s):

Xin Liu ◽

Jing Zhang ◽

Danyao Nie ◽

Kun Zeng ◽

Huiling Hu ◽

...

Keyword(s):

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Scientific Basis ◽

Functional Enrichment ◽

Circular Rnas ◽

Messenger Rnas ◽

Protein Protein Interaction ◽

Matrix Organization ◽

Non Coding Rnas

Pterygium is a common ocular surface disease characterized by abnormal fibrovascular proliferation and invasion, similar to tumorigenesis. The formation of tumors is related to a change in the expression of various RNAs; however, whether they are involved in the formation and development of pterygium remains unclear. In this study, transcriptome analysis of messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) of paired pterygium and normal conjunctiva was performed to explore key genes regulating the development of pterygium. In total, 579 mRNAs, 275 lncRNAs, and 21 circRNAs were differentially expressed (DE) in pterygium compared with paired conjunctival tissues. Functional enrichment analysis indicated that DE RNAs were associated with extracellular matrix organization, blood vessel morphogenesis, and focal adhesion. Furthermore, through protein-protein interaction network and mRNA-lncRNA co-expression network analysis, key mRNAs including FN1, VCAM1, and MMP2, and key lncRNAs including MIR4435-2HG and LINC00968 were screened and might be involved in the pathogenesis of pterygium. In addition, several circRNAs including hsa_circ_0007482 and hsa_circ_001730 were considered to be involved in the pterygium development. This study provides a scientific basis for elucidating the pathogenesis of pterygium and will be beneficial for the development of preventive and therapeutic strategies.

Download Full-text

Transcriptome Analysis Reveals the Potential Roles of Long Non-Coding RNAs in Feed Efficiency of Chicken

10.21203/rs.3.rs-882290/v1 ◽

2021 ◽

Author(s):

Parastoo Karimi ◽

Mohammad Reza Bakhtiarizadeh ◽

Abdolreza Salehi ◽

Hamid Reza Izadnia

Keyword(s):

Feed Efficiency ◽

Liver Tissue ◽

Molecular Mechanisms ◽

Target Genes ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Production Costs ◽

Functional Enrichment ◽

Localization Analysis ◽

Non Coding Rnas

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.

Download Full-text

Identifying Complex lncRNA/Pseudogene–miRNA–mRNA Crosstalk in Hormone-Dependent Cancers

Biology ◽

10.3390/biology10101014 ◽

2021 ◽

Vol 10 (10) ◽

pp. 1014

Author(s):

Dulari K. Jayarathna ◽

Miguel E. Rentería ◽

Emilie Sauret ◽

Jyotsna Batra ◽

Neha S. Gandhi

Keyword(s):

Regulatory Networks ◽

Cancer Survival ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Messenger Rnas ◽

Protein Activity ◽

Cerna Network ◽

Non Coding Rnas

The discovery of microRNAs (miRNAs) has fundamentally transformed our understanding of gene regulation. The competing endogenous RNA (ceRNA) hypothesis postulates that messenger RNAs and other RNA transcripts, such as long non-coding RNAs and pseudogenes, can act as natural miRNA sponges. These RNAs influence each other’s expression levels by competing for the same pool of miRNAs through miRNA response elements on their target transcripts, thereby modulating gene expression and protein activity. In recent years, these ceRNA regulatory networks have gained considerable attention in cancer research. Several studies have identified cancer-specific ceRNA networks. Nevertheless, prior bioinformatic analyses have focused on long-non-coding RNA-associated ceRNA networks. Here, we identify an extended ceRNA network (including both long non-coding RNAs and pseudogenes) shared across a group of five hormone-dependent (HD) cancers, i.e., prostate, breast, colon, rectal, and endometrial cancers, using data from The Cancer Genome Atlas (TCGA). We performed a functional enrichment analysis for differentially expressed genes in the shared ceRNA network of HD cancers, followed by a survival analysis to determine their prognostic ability. We identified two long non-coding RNAs, nine genes, and seventy-four miRNAs in the shared ceRNA network across five HD cancers. Among them, two genes and forty-one miRNAs were associated with at least one HD cancer survival. This study is the first to investigate pseudogene-associated ceRNAs across a group of related cancers and highlights the value of this approach to understanding the shared molecular pathogenesis in a group of related diseases.

Download Full-text

Functional Enrichment Analysis of Deregulated Long Non-Coding RNAs in Cancer Based on their Genomic Neighbors

10.1101/2020.09.14.296921 ◽

2020 ◽

Author(s):

Gulden Olgun ◽

Oznur Tastan

Keyword(s):

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Transcriptional Regulators ◽

Functional Enrichment ◽

Biological Processes ◽

Driver Genes ◽

Cancer Driver ◽

Cancer Types ◽

Non Coding Rnas ◽

The Impact

AbstractThe dysregulation of long non-coding RNAs’ (lncRNAs) expressions has been implicated in cancer. Since most of the lncRNAs’ are not functionally characterized well, investigating the set of perturbed lncRNAs are is challenging. Existing methods that inspect lncRNAs functionally rely on the co-expressed coding genes, which are far better characterized functionally. LncRNAs can be known to act as transcriptional regulators; they may activate or repress the neighborhood’s coding genes on the genome. Based on this, in this work, we aim to analyze the deregulated lncRNAs in cancer by taking into account their ability to regulate nearby loci on the genome. We perform functional analysis on differentially expressed lncRNAs for 28 different cancers considering their adjacent coding genes. We identify that some deregulated lncRNAs are cancer-specific, but a substantial number of lncRNAs are shared across cancers. Next, we assess the similarities of the cancer types based on the functional enrichment of the deregulated lncRNA sets. We find some cancers are very similar in the functions and biological processes related to the deregulated lncRNAs. We observe that some of the cancers for which we find similarity can be linked through primary, metastatic site relations. We investigate the similarity of enriched functional terms for the deregulated lncRNAs and the mRNAs. We further assess the enriched functions’ similarity to the functions and processes that the known cancer driver genes take place. We believe that our methodology help to understand the impact of the lncRNAs in cancer functionally.

Download Full-text

Construction and analysis for dys-regulated lncRNAs and mRNAs in LPS-induced porcine PBMCs

Innate Immunity ◽

10.1177/1753425920983869 ◽

2021 ◽

Vol 27 (2) ◽

pp. 170-183

Author(s):

Jing Zhang ◽

Xin Xu ◽

Hongbo Chen ◽

Ping Kang ◽

Huiling Zhu ◽

...

Keyword(s):

Signaling Pathways ◽

Function Analysis ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Receptor Interaction ◽

Lps Challenge ◽

Tnf Α ◽

Non Coding Rnas ◽

And Function

Long non-coding RNAs (lncRNAs) are emerging as key regulators in inflammation. However, their functions and profiles in LPS-induced inflammation in pigs are largely unknown. In this study, we profiled global lncRNA and mRNA expression changes in PBMCs treated with LPS using the lncRNA-seq technique. In total 43 differentially expressed (DE) lncRNAs and 1082 DE mRNAs were identified in porcine PBMCs after LPS stimulation. Functional enrichment analysis on DE mRNAs indicated these genes were involved in inflammation-related signaling pathways, including cytokine–cytokine receptor interaction, TNF-α, NF-κB, Jak-STAT and TLR signaling pathways. In addition, co-expression network and function analysis identified the potential lncRNAs related to inflammatory response and immune response. The expressions of eight lncRNAs (ENSSSCT00000045208, ENSSSCT00000051636, ENSSSCT00000049770, ENSSSCT00000050966, ENSSSCT00000047491, ENSSSCT00000049750, ENSSSCT00000054262 and ENSSSCT00000044651) were validated in the LPS-treated PBMCs by quantitative real-time PCR (qRT-PCR). In LPS-challenged piglets, we identified that expression of three lncRNAs (ENSSSCT00000051636, ENSSSCT00000049770, and ENSSSCT00000047491) was significantly up-regulated in liver, spleen and jejunum tissues after LPS challenge, which indicated that these lncRNAs might be important regulators for inflammation. This study provides the first lncRNA and mRNA transcriptomic landscape of LPS-mediated changes in porcine PBMCs, which might provide potential insights into lncRNAs involved in regulating inflammation in pigs.

Download Full-text

Systems Biology Approaches Reveal a Multi-stress Responsive WRKY Transcription Factor and Stress Associated Gene Co-expression Networks in Chickpea

Current Bioinformatics ◽

10.2174/1574893614666190204152500 ◽

2019 ◽

Vol 14 (7) ◽

pp. 591-601 ◽

Cited By ~ 1

Author(s):

Aravind K. Konda ◽

Parasappa R. Sabale ◽

Khela R. Soren ◽

Shanmugavadivel P. Subramaniam ◽

Pallavi Singh ◽

...

Keyword(s):

Expression Pattern ◽

Gene Networks ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Wrky Transcription Factor ◽

Functional Enrichment ◽

Differentially Expressed ◽

Callose Deposition ◽

Significant Probability

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.

Download Full-text

Structural variations in papaya genomes

BMC Genomics ◽

10.1186/s12864-021-07665-4 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Zhenyang Liao ◽

Xunxiao Zhang ◽

Shengcheng Zhang ◽

Zhicong Lin ◽

Xingtan Zhang ◽

...

Keyword(s):

Agronomic Traits ◽

Enrichment Analysis ◽

Copy Number Variant ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Structural Variations ◽

Reproduction Traits ◽

Depth Study ◽

Environmental Adaptability ◽

The Impact

Abstract Background Structural variations (SVs) are a type of mutations that have not been widely detected in plant genomes and studies in animals have shown their role in the process of domestication. An in-depth study of SVs will help us to further understand the impact of SVs on the phenotype and environmental adaptability during papaya domestication and provide genomic resources for the development of molecular markers. Results We detected a total of 8083 SVs, including 5260 deletions, 552 tandem duplications and 2271 insertions with deletion being the predominant, indicating the universality of deletion in the evolution of papaya genome. The distribution of these SVs is non-random in each chromosome. A total of 1794 genes overlaps with SV, of which 1350 genes are expressed in at least one tissue. The weighted correlation network analysis (WGCNA) of these expressed genes reveals co-expression relationship between SVs-genes and different tissues, and functional enrichment analysis shows their role in biological growth and environmental responses. We also identified some domesticated SVs genes related to environmental adaptability, sexual reproduction, and important agronomic traits during the domestication of papaya. Analysis of artificially selected copy number variant genes (CNV-genes) also revealed genes associated with plant growth and environmental stress. Conclusions SVs played an indispensable role in the process of papaya domestication, especially in the reproduction traits of hermaphrodite plants. The detection of genome-wide SVs and CNV-genes between cultivated gynodioecious populations and wild dioecious populations provides a reference for further understanding of the evolution process from male to hermaphrodite in papaya.

Download Full-text

CCTs as new biomarkers for the prognosis of head and neck squamous cancer

Open Medicine ◽

10.1515/med-2020-0114 ◽

2020 ◽

Vol 15 (1) ◽

pp. 672-688

Author(s):

Yanbo Dong ◽

Siyu Lu ◽

Zhenxiao Wang ◽

Liangfa Liu

Keyword(s):

Head And Neck ◽

Survival Data ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Advanced Tumor Stage ◽

Expression Levels ◽

T Complex ◽

Advanced Tumor ◽

Squamous Cancer

AbstractThe chaperonin-containing T-complex protein 1 (CCT) subunits participate in diverse diseases. However, little is known about their expression and prognostic values in human head and neck squamous cancer (HNSC). This article aims to evaluate the effects of CCT subunits regarding their prognostic values for HNSC. We mined the transcriptional and survival data of CCTs in HNSC patients from online databases. A protein–protein interaction network was constructed and a functional enrichment analysis of target genes was performed. We observed that the mRNA expression levels of CCT1/2/3/4/5/6/7/8 were higher in HNSC tissues than in normal tissues. Survival analysis revealed that the high mRNA transcriptional levels of CCT3/4/5/6/7/8 were associated with a low overall survival. The expression levels of CCT4/7 were correlated with advanced tumor stage. And the overexpression of CCT4 was associated with higher N stage of patients. Validation of CCTs’ differential expression and prognostic values was achieved by the Human Protein Atlas and GEO datasets. Mechanistic exploration of CCT subunits by the functional enrichment analysis suggests that these genes may influence the HNSC prognosis by regulating PI3K-Akt and other pathways. This study implies that CCT3/4/6/7/8 are promising biomarkers for the prognosis of HNSC.

Download Full-text

Investigation and Functional Enrichment Analysis of the Human Host Interaction Network with Common Gram-Negative Respiratory Pathogens Predicts Possible Association with Lung Adenocarcinoma

Pathophysiology ◽

10.3390/pathophysiology28010003 ◽

2021 ◽

Vol 28 (1) ◽

pp. 20-33

Author(s):

Lydia-Eirini Giannakou ◽

Athanasios-Stefanos Giannopoulos ◽

Chrissi Hatzoglou ◽

Konstantinos I. Gourgoulianis ◽

Erasmia Rouka ◽

...

Keyword(s):

Lung Adenocarcinoma ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Cell Junctions ◽

Respiratory Pathogens ◽

Gram Negative ◽

Human Proteins ◽

Apoptotic Pathways

Haemophilus influenzae (Hi), Moraxella catarrhalis (MorCa) and Pseudomonas aeruginosa (Psa) are three of the most common gram-negative bacteria responsible for human respiratory diseases. In this study, we aimed to identify, using the functional enrichment analysis (FEA), the human gene interaction network with the aforementioned bacteria in order to elucidate the full spectrum of induced pathogenicity. The Human Pathogen Interaction Database (HPIDB 3.0) was used to identify the human proteins that interact with the three pathogens. FEA was performed via the ToppFun tool of the ToppGene Suite and the GeneCodis database so as to identify enriched gene ontologies (GO) of biological processes (BP), cellular components (CC) and diseases. In total, 11 human proteins were found to interact with the bacterial pathogens. FEA of BP GOs revealed associations with mitochondrial membrane permeability relative to apoptotic pathways. FEA of CC GOs revealed associations with focal adhesion, cell junctions and exosomes. The most significantly enriched annotations in diseases and pathways were lung adenocarcinoma and cell cycle, respectively. Our results suggest that the Hi, MorCa and Psa pathogens could be related to the pathogenesis and/or progression of lung adenocarcinoma via the targeting of the epithelial cellular junctions and the subsequent deregulation of the cell adhesion and apoptotic pathways. These hypotheses should be experimentally validated.

Download Full-text