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

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

Identifcation of The Ferroptosis-Related Gene Signature In Placenta of Patients With Early-Onset Preeclampsia

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

2021 ◽

Author(s):

Nana Yang ◽

Qianghua Wang ◽

Biao Ding ◽

Yinging Gong ◽

Yue Wu ◽

...

Keyword(s):

Iron Homeostasis ◽

Molecular Mechanisms ◽

Late Onset ◽

Expression Profiles ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Ppi Network ◽

Network Analyses

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.

Download Full-text

Identification and Interaction Analysis of Significant Genes and MicroRNAs in Pterygium

BioMed Research International ◽

10.1155/2019/2767512 ◽

2019 ◽

Vol 2019 ◽

pp. 1-12 ◽

Cited By ~ 2

Author(s):

Siying He ◽

Hui Sun ◽

Yifang Huang ◽

Shiqi Dong ◽

Chen Qiao ◽

...

Keyword(s):

Differentially Expressed Genes ◽

Molecular Mechanisms ◽

Interaction Analysis ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Gene Expression Omnibus ◽

Functional Enrichment ◽

Differentially Expressed ◽

Differentially Expressed Mirnas

Purpose. MiRNAs have been widely analyzed in the occurrence and development of many diseases, including pterygium. This study aimed to identify the key genes and miRNAs in pterygium and to explore the underlying molecular mechanisms. Methods. MiRNA expression was initially extracted and pooled by published literature. Microarray data about differentially expressed genes was downloaded from Gene Expression Omnibus (GEO) database and analyzed with the R programming language. Functional and pathway enrichment analyses were performed using the database for Annotation, Visualization and Integrated Discovery (DAVID). The protein-protein interaction network was constructed with the STRING database. The associations between chemicals, differentially expressed miRNAs, and differentially expressed genes were predicted using the online resource. All the networks were constructed using Cytoscape. Results. We found that 35 miRNAs and 301 genes were significantly differentially expressed. Functional enrichment analysis showed that upregulated genes were significantly enriched in extracellular matrix (ECM) organization, while downregulated genes were mainly involved in cell death and apoptotic process. Finally, we concluded the chemical-gene affected network, miRNA-mRNA interacted networks, and significant pathway network. Conclusion. We identified lists of differentially expressed miRNAs and genes and their possible interaction in pterygium. The networks indicated that ECM breakdown and EMT might be two major pathophysiological mechanisms and showed the potential significance of PI3K-Akt signalling pathway. MiR-29b-3p and collagen family (COL4A1 and COL3A1) might be new treatment target in pterygium.

Download Full-text

Integrated analysis of lncRNA-miRNA-mRNA ceRNA network in human aortic dissection

BMC Genomics ◽

10.1186/s12864-021-08012-3 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Hao Zhang ◽

Ce Bian ◽

Simei Tu ◽

Fanxing Yin ◽

Panpan Guo ◽

...

Keyword(s):

Aortic Dissection ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Integrated Analysis ◽

Down Regulation ◽

Tissue Samples ◽

Cerna Network

Abstract Background Many studies on long chain non-coding RNAs (lncRNAs) are published in recent years. But the roles of lncRNAs in aortic dissection (AD) are still unclear and should be further examined. The present work focused on determining the molecular mechanisms underlying lncRNAs regulation in aortic dissection on the basis of the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network. Methods This study collected the lncRNAs (GSE52093), mRNAs (GSE52093) and miRNAs (GSE92427) expression data within human tissue samples with aortic dissection group and normal group based on Gene Expression Omnibus (GEO) database. Results This study identified three differentially expressed lncRNAs (DELs), 19 differentially expressed miRNAs (DEmiRs) and 1046 differentially expressed mRNAs (DEGs) identified regarding aortic dissection. Furthermore, we constructed a lncRNA-miRNA-mRNA network through three lncRNAs (including two with up-regulation and one with down-regulation), five miRNAs (five with up-regulation), as well as 211 mRNAs (including 103 with up-regulation and 108 with down-regulation). Simultaneously, we conducted functional enrichment and pathway analyses on genes within the as-constructed ceRNA network. According to our PPI/ceRNA network and functional enrichment analysis results, four critical genes were found (E2F2, IGF1R, BDNF and PPP2R1B). In addition, E2F2 level was possibly modulated via lncRNA FAM87A-hsa-miR-31-5p/hsa-miR-7-5p or lncRNA C9orf106-hsa-miR-7-5p. The expression of IGF1R may be regulated by lncRNA FAM87A-hsa-miR-16-5p/hsa-miR-7-5p or lncRNA C9orf106-hsa-miR-7-5p. Conclusion In conclusion, the ceRNA interaction axis we identified is a potentially critical target for treating AD. Our results shed more lights on the possible pathogenic mechanism in AD using a lncRNA-associated ceRNA network.

Download Full-text

Identification of Biomarkers and Pathogenesis in Severe Asthma by Co-expression Network Analysis

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

2020 ◽

Author(s):

Zeyi Zhang ◽

Ou Chen ◽

Jingjing Wang

Keyword(s):

Network Analysis ◽

Severe Asthma ◽

Gene Networks ◽

Molecular Mechanisms ◽

Area Under The Curve ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Asthma Severity ◽

Functional Enrichment ◽

Hub Genes

Abstract BackgroundSevere asthma is a heterogeneous inflammatory disease. The rise of precise immunotherapy for severe asthmatics underlines more understanding of molecular mechanisms and biomarkers. In this study, we aim to identify underlying mechanisms and hub genes that define asthma severity.MethodsDifferentially expressed genes were screened out based on bronchial epithelial brushings from mild and severe asthmatics. Then, the weighted gene co-expression network analysis was adopted to identify gene networks and the most significant module associated with asthma severity. Meanwhile, hub genes screening and functional enrichment analysis was performed. Receiver operating characteristic was conducted to validate the hub genes.ResultsWeighted gene co-expression network analysis identified 6 modules associated with asthma severity. Three modules were positively correlated (P < 0.001) with asthma severity, containing genes upregulated in severe asthmatics. Functional enrichment analysis found genes in the highlighted module mainly enriched in neutrophil degranulation and activation, leukocyte migration and chemotaxis. Hub genes identified in the module were CXCR1, CXCR2, CCR1, CCR7, TLR2, FPR1, FCGR3B, FCGR2A, ITGAM, and PLEK. Combining these hub genes possessed a moderate ability for discriminating between severe asthmatics and mild-moderate asthmatics with an area under the curve of 0.75.ConclusionOur results identified biomarkers and potential pathogenesis of severe asthma, which provides sight into treatment targets and prognostic markers.

Download Full-text

A Systematic Analysis of Candidate Genes Associated with Nicotine Addiction

BioMed Research International ◽

10.1155/2015/313709 ◽

2015 ◽

Vol 2015 ◽

pp. 1-9 ◽

Cited By ~ 2

Author(s):

Meng Liu ◽

Xia Li ◽

Rui Fan ◽

Xinhua Liu ◽

Ju Wang

Keyword(s):

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Nicotine Addiction ◽

Clustering Coefficient ◽

Functional Enrichment ◽

Moderate Degree ◽

Cancer Genes ◽

Systematic Analysis ◽

The Central Nervous System

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.

Download Full-text

Identification of Key Genes and Pathways Associated with Differences of Subchondral Bone in Osteoarthritis

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

2020 ◽

Author(s):

Yiyuan Zhang ◽

Rongguo Yu ◽

Jiayu Zhang ◽

Eryou Feng ◽

Haiyang Wang ◽

...

Keyword(s):

Subchondral Bone ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Hub Genes ◽

Common Chronic Disease ◽

Drug Candidates ◽

Therapeutic Modalities ◽

Pathogenic Genes

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.

Download Full-text

Identification of a prognostic long non-coding RNA signature in breast cancer

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

2020 ◽

Author(s):

Gaochen Lan ◽

Xiaoling Yu ◽

Yanna Zhao ◽

Jinjian Lan ◽

Wan Li ◽

...

Keyword(s):

Breast Cancer ◽

Cox Regression ◽

Prognostic Biomarker ◽

Expression Profiles ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Differentially Expressed ◽

Cox Regression Analysis

Abstract Background: Breast cancer is the most common malignant disease among women. At present, more and more attention has been paid to long non-coding RNAs (lncRNAs) in the field of breast cancer research. We aimed to investigate the expression profiles of lncRNAs and construct a prognostic lncRNA for predicting the overall survival (OS) of breast cancer.Methods: The expression profiles of lncRNAs and clinical data with breast cancer were obtained from The Cancer Genome Atlas (TCGA). Differentially expressed lncRNAs were screened out by R package (limma). The survival probability was estimated by the Kaplan‑Meier Test. The Cox Regression Model was performed for univariate and multivariate analysis. The risk score (RS) was established on the basis of the lncRNAs’ expression level (exp) multiplied regression coefficient (β) from the multivariate cox regression analysis with the following formula: RS=exp a1 * β a1 + exp a2 * β a2 +……+ exp an * β an. Functional enrichment analysis was performed by Metascape.Results: A total of 3404 differentially expressed lncRNAs were identified. Among them, CYTOR, MIR4458HG and MAPT-AS1 were significantly associated with the survival of breast cancer. Finally, The RS could predict OS of breast cancer (RS=exp CYTOR * β CYTOR + exp MIR4458HG * β MIR4458HG + exp MAPT-AS1 * β MAPT-AS1). Moreover, it was confirmed that the three-lncRNA signature could be an independent prognostic biomarker for breast cancer (HR=3.040, P=0.000).Conclusions: This study established a three-lncRNA signature, which might be a novel prognostic biomarker for breast cancer.

Download Full-text

Immunological pathways of macrophage response to Brucella ovis infection

Innate Immunity ◽

10.1177/1753425920958179 ◽

2020 ◽

Vol 26 (7) ◽

pp. 635-648

Author(s):

Zhixiong Zhou ◽

Guojing Gu ◽

Yichen Luo ◽

Wenjie Li ◽

Bowen Li ◽

...

Keyword(s):

Immune Response ◽

Molecular Mechanisms ◽

Group Versus ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Sequencing Data ◽

Brucella Ovis ◽

Qrt Pcr ◽

Raw264.7 Macrophage

As the molecular mechanisms of Brucella ovis pathogenicity are not completely clear, we have applied a transcriptome approach to identify the differentially expressed genes (DEGs) in RAW264.7 macrophage infected with B. ovis. The DEGs related to immune pathway were identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) functional enrichment analysis. Quantitative real-time PCR (qRT-PCR) was performed to validate the transcriptome sequencing data. In total, we identified 337 up-regulated and 264 down-regulated DEGs in B. ovis-infected group versus mock group. Top 20 pathways were enriched by KEGG analysis and 20 GO by functional enrichment analysis in DEGs involved in the molecular function, cellular component, and biological process and so on, which revealed multiple immunological pathways in RAW264.7 macrophage cells in response to B. ovis infection, including inflammatory response, immune system process, immune response, cytokine activity, chemotaxis, chemokine-mediated signaling pathway, chemokine activity, and CCR chemokine receptor binding. qRT-PCR results showed Ccl2 (ENSMUST00000000193), Ccl2 (ENSMUST00000124479), Ccl3 (ENSMUST00000001008), Hmox1 (ENSMUST00000005548), Hmox1 (ENSMUST00000159631), Cxcl2 (ENSMUST00000075433), Cxcl2 (ENSMUST00000200681), Cxcl2 (ENSMUST00000200919), and Cxcl2 (ENSMUST00000202317). Our findings firstly elucidate the pathways involved in B. ovis-induced host immune response, which may lay the foundation for revealing the bacteria–host interaction and demonstrating the pathogenic mechanism of B. ovis.

Download Full-text

Comparative Transcriptome Analysis Suggests Key Roles for 5-Hydroxytryptamlne Receptors in Control of Goose Egg Production

Genes ◽

10.3390/genes11040455 ◽

2020 ◽

Vol 11 (4) ◽

pp. 455 ◽

Cited By ~ 1

Author(s):

Qingyuan Ouyang ◽

Shenqiang Hu ◽

Guosong Wang ◽

Jiwei Hu ◽

Jiaman Zhang ◽

...

Keyword(s):

Egg Production ◽

Expression Profiles ◽

Expression Patterns ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Production Performance ◽

Functional Enrichment ◽

Differentially Expressed ◽

Receptor Interaction ◽

Anser Anser

To date, research on poultry egg production performance has only been conducted within inter or intra-breed groups, while those combining both inter- and intra-breed groups are lacking. Egg production performance is known to differ markedly between Sichuan white goose (Anser cygnoides) and Landes goose (Anser anser). In order to understand the mechanism of egg production performance in geese, we undertook this study. Here, 18 ovarian stromal samples from both Sichuan white goose and Landes goose at the age of 145 days (3 individuals before egg production initiation for each breed) and 730 days (3 high- and low egg production individuals during non-laying periods for each breed) were collected to reveal the genome-wide expression profiles of ovarian mRNAs and lncRNAs between these two geese breeds at different physiological stages. Briefly, 58, 347, 797, 777, and 881 differentially expressed genes (DEGs) and 56, 24, 154, 105, and 224 differentially expressed long non-coding RNAs (DElncRNAs) were found in LLD vs. HLD (low egg production Landes goose vs. high egg production Landes goose), LSC vs. HSC (low egg production Sichuan White goose vs. high egg production Sichuan white goose), YLD vs. YSC (young Landes goose vs. young Sichuan white goose), HLD vs. HSC (high egg production Landes goose vs. high egg production Sichuan white goose), and LLD vs. LSC (low egg production Landes goose vs. low egg production Sichuan white goose) groups, respectively. Functional enrichment analysis of these DEGs and DElncRNAs suggest that the “neuroactive ligand–receptor interaction pathway” is crucial for egg production, and particularly, members of the 5-hydroxytryptamine receptor (HTR) family affect egg production by regulating ovarian metabolic function. Furthermore, the big differences in the secondary structures among HTR1F and HTR1B, HTR2B, and HTR7 may lead to their different expression patterns in goose ovaries of both inter- and intra-breed groups. These results provide novel insights into the mechanisms regulating poultry egg production performance.

Download Full-text

Integrative Analysis Reveals Common and Unique Roles of Tetraspanins in Fibrosis and Emphysema

Frontiers in Genetics ◽

10.3389/fgene.2020.585998 ◽

2020 ◽

Vol 11 ◽

Author(s):

Lokesh P. Tripathi ◽

Mari N. Itoh ◽

Yoshito Takeda ◽

Kazuyuki Tsujino ◽

Yasushi Kondo ◽

...

Keyword(s):

Pulmonary Fibrosis ◽

Molecular Mechanisms ◽

Inflammatory Responses ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Integrative Approach ◽

Chronic Obstructive ◽

Lung Pathology ◽

Obstructive Pulmonary Disease

While both chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are multifactorial disorders characterized by distinct clinical and pathological features, their commonalities and differences have not been fully elucidated. We sought to investigate the preventive roles of tetraspanins Cd151 and Cd9 -that are involved in diverse cellular processes in lung pathophysiology- in pulmonary fibrosis and emphysema, respectively, and to obtain a deeper understanding of their underlying molecular mechanisms toward facilitating improved therapeutic outcomes. Using an integrative approach, we examined the transcriptomic changes in the lungs of Cd151- and Cd9-deficient mice using functional-enrichment-analysis, pathway-perturbation-analysis and protein-protein-interaction (PPI) network analysis. Circadian-rhythm, extracellular-matrix (ECM), cell-adhesion and inflammatory responses and associated factors were prominently influenced by Cd151-deletion. Conversely, cellular-junctions, focal-adhesion, vascular-remodeling, and TNF-signaling were deeply impacted by Cd9-deletion. We also highlighted a “common core” of factors and signaling cascades that underlie the functions of both Cd151 and Cd9 in lung pathology. Circadian dysregulation following Cd151-deletion seemingly facilitated progressive fibrotic lung phenotype. Conversely, TGF-β signaling attenuation and TNF-signaling activation emerged as potentially novel functionaries of Cd9-deletion-induced emphysema. Our findings offer promising avenues for developing novel therapeutic treatments for pulmonary fibrosis and emphysema.

Download Full-text