Intraspecific transcriptome variation and sex-biased expression in Anopheles arabiensis

Abstract The magnitude and functional patterns of intraspecific transcriptional variation in the anophelines, including those of sex-biased genes underlying sex-specific traits relevant for malaria transmission, remain understudied. As a result, how changes in expression levels drive adaptation in these species is poorly understood. We sequenced the female, male, and larval transcriptomes of three populations of Anopheles arabiensis from Burkina Faso. One-third of the genes were differentially expressed between populations, often involving insecticide resistance-related genes in a sample type-specific manner, and with the females showing the largest number of differentially expressed genes. At the genomic level, the X chromosome appears depleted of differentially expressed genes compared to the autosomes, chromosomes harbouring inversions do not exhibit evidence for enrichment of such genes, and genes that are top contributors to functional enrichment patterns of population differentiation tend to be clustered in the genome. Further, the magnitude of variation for the sex expression ratio across populations did not substantially differ between male- and female-biased genes, except for some populations in which male-limited expressed genes showed more variation than their female counterparts. In fact, female-biased genes exhibited a larger level of interpopulation variation than male-biased genes, both when assayed in males and females. Beyond uncovering extensive adaptive potential of transcriptional variation in An. arabiensis, our findings suggest that the evolutionary rate of changes in expression levels on the X chromosome exceeds that on the autosomes, while pointing to female-biased genes as the most variable component of the An. arabiensis transcriptome.

Download Full-text

Differential gene expression analysis of ankylosing spondylitis shows deregulation of the HLA-DRB, HLA-DQB, ITM2A, and CTLA4 genes

Egyptian Journal of Medical Human Genetics ◽

10.1186/s43042-021-00161-0 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Rowan AlEjielat ◽

Anas Khaleel ◽

Amneh H. Tarkhan

Keyword(s):

Gene Expression ◽

Ankylosing Spondylitis ◽

Differentially Expressed Genes ◽

Gene Network ◽

Disease Diagnosis ◽

Receptor Activation ◽

Functional Enrichment ◽

Differentially Expressed ◽

Inflammatory Disorder ◽

Data Set

Abstract Background Ankylosing spondylitis (AS) is a rare inflammatory disorder affecting the spinal joints. Although we know some of the genetic factors that are associated with the disease, the molecular basis of this illness has not yet been fully elucidated, and the genes involved in AS pathogenesis have not been entirely identified. The current study aimed at constructing a gene network that may serve as an AS gene signature and biomarker, both of which will help in disease diagnosis and the identification of therapeutic targets. Previously published gene expression profiles of 16 AS patients and 16 gender- and age-matched controls that were profiled on the Illumina HumanHT-12 V3.0 Expression BeadChip platform were mined. Patients were Portuguese, 21 to 64 years old, were diagnosed based on the modified New York criteria, and had Bath Ankylosing Spondylitis Disease Activity Index scores > 4 and Bath Ankylosing Spondylitis Functional Index scores > 4. All patients were receiving only NSAIDs and/or sulphasalazine. Functional enrichment and pathway analysis were performed to create an interaction network of differentially expressed genes. Results ITM2A, ICOS, VSIG10L, CD59, TRAC, and CTLA-4 were among the significantly differentially expressed genes in AS, but the most significantly downregulated genes were the HLA-DRB6, HLA-DRB5, HLA-DRB4, HLA-DRB3, HLA-DRB1, HLA-DQB1, ITM2A, and CTLA-4 genes. The genes in this study were mostly associated with the regulation of the immune system processes, parts of cell membrane, and signaling related to T cell receptor and antigen receptor, in addition to some overlaps related to the IL2 STAT signaling, as well as the androgen response. The most significantly over-represented pathways in the data set were associated with the “RUNX1 and FOXP3 which control the development of regulatory T lymphocytes (Tregs)” and the “GABA receptor activation” pathways. Conclusions Comprehensive gene analysis of differentially expressed genes in AS reveals a significant gene network that is involved in a multitude of important immune and inflammatory pathways. These pathways and networks might serve as biomarkers for AS and can potentially help in diagnosing the disease and identifying future targets for treatment.

Download Full-text

Behavioral sensitization induced by methamphetamine causes differential alterations in gene expression and histone acetylation of the prefrontal cortex in rats

BMC Neuroscience ◽

10.1186/s12868-021-00616-5 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Hui Li ◽

Jing-An Chen ◽

Qian-Zhi Ding ◽

Guan-Yi Lu ◽

Ning Wu ◽

...

Keyword(s):

Gene Expression ◽

Prefrontal Cortex ◽

Histone Acetylation ◽

Differentially Expressed Genes ◽

Behavioral Sensitization ◽

Functional Enrichment ◽

Differentially Expressed ◽

Adult Rats ◽

Mrna Microarray ◽

H3 Acetylation

Abstract Background Methamphetamine (METH) is one of the most widely abused illicit substances worldwide; unfortunately, its addiction mechanism remains unclear. Based on accumulating evidence, changes in gene expression and chromatin modifications might be related to the persistent effects of METH on the brain. In the present study, we took advantage of METH-induced behavioral sensitization as an animal model that reflects some aspects of drug addiction and examined the changes in gene expression and histone acetylation in the prefrontal cortex (PFC) of adult rats. Methods We conducted mRNA microarray and chromatin immunoprecipitation (ChIP) coupled to DNA microarray (ChIP-chip) analyses to screen and identify changes in transcript levels and histone acetylation patterns. Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, were performed to analyze the differentially expressed genes. We then further identified alterations in ANP32A (acidic leucine-rich nuclear phosphoprotein-32A) and POU3F2 (POU domain, class 3, transcription factor 2) using qPCR and ChIP-PCR assays. Results In the rat model of METH-induced behavioral sensitization, METH challenge caused 275 differentially expressed genes and a number of hyperacetylated genes (821 genes with H3 acetylation and 10 genes with H4 acetylation). Based on mRNA microarray and GO and KEGG enrichment analyses, 24 genes may be involved in METH-induced behavioral sensitization, and 7 genes were confirmed using qPCR. We further examined the alterations in the levels of the ANP32A and POU3F2 transcripts and histone acetylation at different periods of METH-induced behavioral sensitization. H4 hyperacetylation contributed to the increased levels of ANP32A mRNA and H3/H4 hyperacetylation contributed to the increased levels of POU3F2 mRNA induced by METH challenge-induced behavioral sensitization, but not by acute METH exposure. Conclusions The present results revealed alterations in transcription and histone acetylation in the rat PFC by METH exposure and provided evidence that modifications of histone acetylation contributed to the alterations in gene expression caused by METH-induced behavioral sensitization.

Download Full-text

Genomic Expression Profiling of Colorectal Cancer in Silico

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

2021 ◽

Author(s):

Feifei Liu ◽

Yu Wang ◽

Wenxue Li ◽

Diancheng Li ◽

Yuwei Xin ◽

...

Keyword(s):

Colorectal Cancer ◽

Survival Analysis ◽

Mrna Expression ◽

Differentially Expressed Genes ◽

Expression Profiles ◽

Functional Enrichment ◽

Differentially Expressed ◽

Analysis Tool ◽

Hub Genes ◽

Normal Tissues

Abstract Background: Colorectal cancer (CRC) is one of the most common malignancies of the digestive system; the progression and prognosis of which are affected by a complicated network of genes and pathways. The aim of this study was to identify potential hub genes associated with the progression and prognosis of colorectal cancer (CRC).Methods: We obtained gene expression profiles from GEO database to search differentially expressed genes (DEGs) between CRC tissues and normal tissue. Subsequently, we conducted a functional enrichment analysis, generated a protein–protein interaction (PPI) network to identify the hub genes, and analyzed the expression validation of the hub genes. Kaplan–Meier plotter survival analysis tool was performed to evaluate the prognostic value of hub genes expression in CRC patients.Results: A total of 370 samples, involving CRC and normal tissues were enrolled in this article. 283 differentially expressed genes (DEGs), including 62 upregulated genes and 221 downregulated genes between CRC and normal tissues were selected. We finally filtered out 6 hub genes, including INSL5, MTIM, GCG, SPP1, HSD11B2, and MAOB. In the database of TCGA-COAD, the mRNA expression of INSL5, MT1M, HSD11B2, MAOB in tumor is lower than that in normal; the mRNA expression of SPP1 in tumor is higher than that in normal. In the HPA database, the expression of INSL5, GCG, HSD11B2, MAOB in tumor is lower than that in normal tissues; the expression of SPP1 in the tumor is higher than that in normal tissues. Survival analysis revealed that INSL5, GCG, SPP1 and MT1M may serve as prognostic biomarkers in CRC. Conclusions: We screened out six hub genes to predict the occurrence and prognosis of patients with CRC using bioinformatics methods, which may provide new targets and ideas for diagnosis, prognosis and individualized treatment for CRC.

Download Full-text

Diagnostic model of combined ceRNA and DNA methylation related genes in esophageal carcinoma

PeerJ ◽

10.7717/peerj.8831 ◽

2020 ◽

Vol 8 ◽

pp. e8831 ◽

Cited By ~ 1

Author(s):

Xiaojiao Guan ◽

Yao Yao ◽

Guangyao Bao ◽

Yue Wang ◽

Aimeng Zhang ◽

...

Keyword(s):

Gene Expression ◽

Esophageal Cancer ◽

Esophageal Carcinoma ◽

Differentially Expressed Genes ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Diagnostic Model ◽

Link Type ◽

Key Genes

Esophageal cancer is a common malignant tumor in the world, and the aim of this study was to screen key genes related to the development of esophageal cancer using a variety of bioinformatics analysis tools and analyze their biological functions. The data of esophageal squamous cell carcinoma from the Gene Expression Omnibus (GEO) were selected as the research object, processed and analyzed to screen differentially expressed microRNAs (miRNAs) and differential methylation genes. The competing endogenous RNAs (ceRNAs) interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, String, and Cytoscape. Biofunctional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of the screened genes and the survival of the patients were verified. By analyzing GSE59973 and GSE114110, we found three down-regulated and nine up-regulated miRNAs. The gene expression matrix of GSE120356 was calculated by Pearson correlation coefficient, and the 11696 pairs of ceRNA relation were determined. In the ceRNA network, 643 lncRNAs and 147 mRNAs showed methylation difference. Functional enrichment analysis showed that these differentially expressed genes were mainly concentrated in the FoxO signaling pathway and were involved in the corresponding cascade of calcineurin. By analyzing the clinical data in The Cancer Genome Atlas (TCGA) database, it was found that four lncRNAs had an important impact on the survival and prognosis of esophageal carcinoma patients. QRT-PCR was also conducted to identify the expression of the key lncRNAs (RNF217-AS1, HCP5, ZFPM2-AS1 and HCG22) in ESCC samples. The selected key genes can provide theoretical guidance for further research on the molecular mechanism of esophageal carcinoma and the screening of molecular markers.

Download Full-text

Identification of differentially expressed genes between the colon and ileum of patients with inflammatory bowel disease by gene co-expression analysis

Journal of International Medical Research ◽

10.1177/0300060519887268 ◽

2019 ◽

Vol 48 (5) ◽

pp. 030006051988726

Author(s):

Yuting Zhang ◽

Bo Shen ◽

Liya Zhuge ◽

Yong Xie

Keyword(s):

Inflammatory Bowel Disease ◽

Differentially Expressed Genes ◽

Bowel Disease ◽

Gene Expression Pattern ◽

Functional Enrichment ◽

Wilcoxon Test ◽

Differentially Expressed ◽

Ppi Network ◽

Hub Genes ◽

Inflammatory Bowel

Objective We aimed to identify differentially expressed genes (DEG) in patients with inflammatory bowel disease (IBD). Methods RNA-seq data were obtained from the Array Express database. DEG were identified using the edgeR package. A co-expression network was constructed and key modules with the highest correlation with IBD inflammatory sites were identified for analysis. The Cytoscape MCODE plugin was used to identify key sub-modules of the protein–protein interaction (PPI) network. The genes in the sub-modules were considered hub genes, and functional enrichment analysis was performed. Furthermore, we constructed a drug–gene interaction network. Finally, we visualized the hub gene expression pattern between the colon and ileum of IBD using the ggpubr package and analyzed it using the Wilcoxon test. Results DEG were identified between the colon and ileum of IBD patients. Based on the co-expression network, the green module had the highest correlation with IBD inflammatory sites. In total, 379 DEG in the green module were identified for the PPI network. Nineteen hub genes were differentially expressed between the colon and ileum. The drug–gene network identified these hub genes as potential drug targets. Conclusion Nineteen DEG were identified between the colon and ileum of IBD patients.

Download Full-text

Association of differentially expressed genes and autoantibody type in patients with systemic sclerosis

Rheumatology ◽

10.1093/rheumatology/keaa447 ◽

2020 ◽

Author(s):

Jun Inamo

Keyword(s):

Gene Expression ◽

Differentially Expressed Genes ◽

Cellular Response ◽

Rna Polymerase Iii ◽

Skin Lesions ◽

Functional Enrichment ◽

Differentially Expressed ◽

Study Cohort ◽

Healthy Controls ◽

Specific Autoantibody

Abstract Objectives The aims of this study were to investigate the relationship between the type of autoantibody and gene expression profile in skin lesions from patients with SSc, and to identify specific dysregulated pathways in SSc patients compared with healthy controls. Methods Sixty-one patients with SSc from the Genetics vs Environment in Scleroderma Outcome Study cohort and 36 healthy controls were included in this study. Differentially expressed genes were extracted and functional enrichment and pathway analysis were conducted. Results Compared with healthy controls, lists containing 2, 71, 10, 144 and 78 differentially expressed genes were created for patients without specific autoantibody, ACA, anti-U1 RNP antibody (RNP), anti-RNA polymerase III antibody (RNAP) and anti-topoisomerase I antibody (ATA), respectively. While part of the enriched pathways overlapped, distinct pathways were identified except in those patients lacking specific autoantibody. The distinct enriched pathways included ‘keratinocyte differentiation’ for ACA, ‘nuclear factor κB signalling’ and ‘cellular response to TGF-β stimulus’ for RNAP, ‘interferon α/β signalling’ for RNP, and ‘cellular response to stress’ for ATA. Cell type signature score analysis revealed that macrophages/monocytes, endothelial cells and fibroblasts were associated with ACA, RNAP, ATA and the severity of the SSc skin lesions. Conclusion Pathogenic pathways were identified according to the type of autoantibody by leveraging gene expression data of patients and controls from a multicentre cohort. The current study may promote the search for new therapeutic targets for SSc.

Download Full-text

Deep phenotyping of synovial molecular signatures by integrative systems analysis in rheumatoid arthritis

Rheumatology ◽

10.1093/rheumatology/keaa751 ◽

2020 ◽

Author(s):

Seung Min Jung ◽

Kyung-Su Park ◽

Ki-Jo Kim

Keyword(s):

Differentially Expressed Genes ◽

Systems Approach ◽

Synovial Fibroblasts ◽

Functional Enrichment ◽

Differentially Expressed ◽

Driver Genes ◽

Key Driver ◽

Driver Analysis ◽

Synovial Tissues ◽

Deep Phenotyping

Abstract Objective RA encompasses a complex, heterogeneous and dynamic group of diseases arising from molecular and cellular perturbations of synovial tissues. The aim of this study was to decipher this complexity using an integrative systems approach and provide novel insights for designing stratified treatments. Methods An RNA sequencing dataset of synovial tissues from 152 RA patients and 28 normal controls was imported and subjected to filtration of differentially expressed genes, functional enrichment and network analysis, non-negative matrix factorization, and key driver analysis. A naïve Bayes classifier was applied to the independent datasets to investigate the factors associated with treatment outcome. Results A matrix of 1241 upregulated differentially expressed genes from RA samples was classified into three subtypes (C1–C3) with distinct molecular and cellular signatures. C3 with prominent immune cells and proinflammatory signatures had a stronger association with the presence of ACPA and showed a better therapeutic response than C1 and C2, which were enriched with neutrophil and fibroblast signatures, respectively. C2 was more occupied by synovial fibroblasts of destructive phenotype and carried highly expressed key effector molecules of invasion and osteoclastogenesis. CXCR2, JAK3, FYN and LYN were identified as key driver genes in C1 and C3. HDAC, JUN, NFKB1, TNF and TP53 were key regulators modulating fibroblast aggressiveness in C2. Conclusions Deep phenotyping of synovial heterogeneity captured comprehensive and discrete pathophysiological attributes of RA regarding clinical features and treatment response. This result could serve as a template for future studies to design stratified approaches for RA patients.

Download Full-text

Identification of key genes, pathways and potential therapeutic agents for IgA nephropathy using an integrated bioinformatics analysis

Journal of the Renin-Angiotensin-Aldosterone System ◽

10.1177/1470320320919635 ◽

2020 ◽

Vol 21 (2) ◽

pp. 147032032091963

Author(s):

Xiaoxue Chen ◽

Mindan Sun

Keyword(s):

Differentially Expressed Genes ◽

Immunoglobulin A ◽

Therapeutic Targets ◽

Gene Expression Omnibus ◽

Functional Enrichment ◽

Differentially Expressed ◽

Signal Pathways ◽

Immunoglobulin A Nephropathy ◽

Hub Genes ◽

Limma Package

Purpose: This study aims to identify immunoglobulin-A-nephropathy-related genes based on microarray data and to investigate novel potential gene targets for immunoglobulin-A-nephropathy treatment. Methods: Immunoglobulin-A-nephropathy chip data was obtained from the Gene Expression Omnibus database, which included 10 immunoglobulin-A-nephropathy and 22 normal samples. We used the limma package of R software to screen differentially expressed genes in immunoglobulin-A-nephropathy and normal glomerular compartment tissues. Functional enrichment (including cellular components, molecular functions, biological processes) and signal pathways were performed for the differentially expressed genes. The online analysis database (STRING) was used to construct the protein-protein interaction networks of differentially expressed genes, and Cytoscape software was used to identify the hub genes of the signal pathway. In addition, we used the Connectivity Map database to predict possible drugs for the treatment of immunoglobulin-A-nephropathy. Results: A total of 348 differentially expressed genes were screened including 107 up-regulated and 241 down-regulated genes. Functional analysis showed that up-regulated differentially expressed genes were mainly concentrated on leukocyte migration, and the down-regulated differentially expressed genes were significantly enriched in alpha-amino acid metabolic process. A total of six hub genes were obtained: JUN, C3AR1, FN1, AGT, FOS, and SUCNR1. The small-molecule drugs thapsigargin, ciclopirox and ikarugamycin were predicted therapeutic targets against immunoglobulin-A-nephropathy. Conclusion: Differentially expressed genes and hub genes can contribute to understanding the molecular mechanism of immunoglobulin-A-nephropathy and providing potential therapeutic targets and drugs for the diagnosis and treatment of immunoglobulin-A-nephropathy.

Download Full-text

Identification of nine key genes by bioinformatics analysis for predicting poor prognosis in smoking-induced lung adenocarcinoma

Lung Cancer Management ◽

10.2217/lmt-2020-0009 ◽

2020 ◽

Vol 9 (2) ◽

pp. LMT30

Author(s):

Chuanli Ren ◽

Weixiu Sun ◽

Xu Lian ◽

Chongxu Han

Keyword(s):

Lung Adenocarcinoma ◽

Differentially Expressed Genes ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Kaplan Meier ◽

Key Genes ◽

Core Genes

Aim: To screen and identify key genes related to the development of smoking-induced lung adenocarcinoma (LUAD). Materials & methods: We obtained data from the GEO chip dataset GSE31210. The differentially expressed genes were screened by GEO2R. The protein interaction network of differentially expressed genes was constructed by STRING and Cytoscape. Finally, core genes were screened. The overall survival time of patients with the core genes was analyzed by Kaplan–Meier method. Gene ontology and Kyoto encyclopedia of genes and genomes bioaccumulation was calculated by DAVID. Results: Functional enrichment analysis indicated that nine key genes were actively involved in the biological process of smoking-related LUAD. Conclusion: 23 core genes and nine key genes among them were correlated with adverse prognosis of LUAD induced by smoking.

Download Full-text

Significance of tumor protein p53 mutation in cellular process and drug selection in brain lower grade (WHO grades II and III) glioma

Biomarkers in Medicine ◽

10.2217/bmm-2020-0331 ◽

2020 ◽

Vol 14 (12) ◽

pp. 1139-1150

Author(s):

Chang-feng Guo ◽

Yugang Zhuang ◽

Yuanzhuo Chen ◽

Sheng Chen ◽

Hu Peng ◽

...

Keyword(s):

Differentially Expressed Genes ◽

Tp53 Mutation ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Differentially Expressed ◽

World Health ◽

Outcome Analysis ◽

Lower Grade ◽

Protein P53 ◽

Tumor Protein P53

Aim: Tumor protein p53 ( TP53) mutant is one of the most frequently mutated genes in glioma. Results: The Cancer Genome Atlas data has shown that TP53 mutation is present in 49% of lower grade (World Health Organization [WHO] grades II and III) glioma patients. Data from The Genomics of Drug Sensitivity in Cancer database showed that three drugs: (5Z)-7-oxozeaenol, dabrafenib and nutlin-3a (−), have shown more resistance in patients with TP53 mutation. We identified 1100 differentially expressed genes. Functional enrichment analysis showed that the differentially expressed genes are mainly concentrated in the transport of ionic and cancer-related pathways. The top ten hub genes were identified and an outcome analysis revealed the most critical genes related to prognosis. Conclusion: Our results identified the key genes and pathways that might provide the basic proof to improve individualized treatment in patients with glioma.

Download Full-text