A global meta-analysis of microarray expression data to predict unknown gene functions and estimate the literature-data divide

Gene Ontology (GO) analysis can provide a comprehensive function analysis for investigating genes, allowing us to identify the potential biological roles of genes. The present study conducted GO analysis to explore the biological function of IRAK2 and performed a case analysis to define its clinical role in disease progression and mediating tumor response to RT. Methods: We performed a GO enrichment analysis on the RNA-seq data to validate radiation-induced gene expression. A total of 172 I-IVB specimens from oral squamous cell carcinoma patients were collected for clinical analysis, from which IRAK2 expression was analyzed by immunohistochemistry. This was a retrospective study conducted between IRAK2 expression and the outcomes of oral squamous cell carcinoma patients after radiotherapy treatment. We conducted Gene Ontology (GO) analysis to explore the biological function of IRAK2 and performed a case analysis to define its clinical role in mediating tumor response to radiotherapy. GO enrichment analysis to validate radiation-induced gene expression was performed. Clinically, 172 stage I-IVB resected oral cancer patients were used to validate IRAK2 expression in predicting clinical outcomes. GO enrichment analysis showed that IRAK2 is involved in 10 of the 14 most enriched GO categories for post-irradiation biological processes, focusing on stress response and immune modulation. Clinically, high IRAK2 expression was correlated with adverse disease features, including pT3-4 status (p = 0.01), advanced overall stage (p = 0.02), and positive bone invasion (p = 0.01). In patients who underwent radiotherapy, the IRAK2-high group was associated with reduced post-irradiation local recurrence (p = 0.025) compared to the IRAK2-low group. IRAK2 plays a crucial role in the radiation-induced response. Patients with high IRAK2 expression demonstrated more advanced disease features but predicted higher post-irradiation local control in a clinical setting. These findings support IRAK2 as a potential predictive biomarker for radiotherapy response in non-metastatic and resected oral cancer patients.

Study question What are the molecular pathways overactivated in arrested cleavage stage human embryos? Summary answer There is an upregulation of mitochondrial activity and cellular stress-related pathways in arrested cleavage embryos, which is in agreement with the “quiet hypothesis”. What is known already mtDNA content decreases during embryo development, however there is a high increase in mtDNA content in arrested cleavage embryos that may correspond to a response to intrinsic or extrinsic factor creating stress. This reasoning would be in agreement with stablished hypothesis showing that a basal metabolism with a moderate-low energy consumption is actually a sign of embryo health, so a viable embryo does not need to use an extra energy to accommodate to the environment (Leese, 2012). The study of the transcriptional profile during human embryo development will give further information about key molecular process involved in in-vitro embryonic competence. Study design, size, duration A prospective cohort study was performed with 11 MII oocytes (average age= 22,9 years), 10 non-arrested cleavage embryos (average age= 29,9 years), 5 arrested cleavage embryos (average age= 38,8 years) and 8 blastocysts (average age= 39,1 years). All specimens were warmed and sampled in PCR tubes with 2 μl of suitable reaction buffer for the RNA sequencing protocol. Participants/materials, setting, methods Specimens were analyzed by single-cell RNA sequencing (scRNA-seq). Correlation studies, principal component and differential expression analysis were performed with DESeq2 package. Differential gene expression analyses were done using the parametric Wald test, with Benjamini-Hochberg multiple test correction (padj). Finally, Fgsea algorithm was used for enrichment analysis on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene ontology (GO) terms. Main results and the role of chance We do not observe mitochondria-related activity pathways significantly (P>.05) deregulated between MII oocytes and non-arrested cleavage embryos considering GO and KEGG categories. When comparing non-arrested cleavage embryos versus blastocysts, we observe several ATP production/consumption and cristae formation-related pathways significantly (P<.05) upregulated in blastocysts compared to non-arrested cleavage embryos considering GO and KEGG categories. This change in activity coincides with the metabolic activation event that occurs in the blastocyst stage. However, when we analyze arrested cleavage embryos versus non-arrested cleavage embryos, we observe several ATP production related-pathways and mitochondria-related apoptosis pathways significantly (P<.05) upregulated in arrested cleavage embryos compared to non-arrested cleavage embryos considering GO categories. With KEGG categories, we notice a significant (P<.05) upregulation of oxidative phosphorylation in arrested cleavage embryos. On the other hand, when we analyze the differences between arrested cleavage embryos and blastocysts, taking into account the differences related to the change of embryo stage, we do not observe ATP production or consumption-related pathways significantly (P>.05) deregulated considering GO and KEGG categories. Then, human arrested embryos in parallel with the increase in mtDNA content, display an upregulation of mitochondrial activity and cellular stress which is in line with the expected overactive metabolism of non-viable embryos. Limitations, reasons for caution All analyzed blastocyst were aneuploid, so we are unable to determine what the results would be with euploid blastocysts. Also, although age can it be ruled out, no differences were observed between mean age from cleavage stage embryos (arrested and non-arrested ones) and blastocyst. Wider implications of the findings: This study evidences the relation between extreme mtDNA content increase and the identification of the pathways involved in active metabolism and apoptosis in arrested cleavage stage human embryos. Trial registration number Not applicable

Previous studies have shown that asthma is a risk factor for lung cancer, while the mechanisms involved remain unclear. We attempted to further explore the association between asthma and non-small cell lung cancer (NSCLC) via bioinformatics analysis. We obtained GSE143303 and GSE18842 from the GEO database. Lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) groups were downloaded from the TCGA database. Based on the results of differentially expressed genes (DEGs) between asthma and NSCLC, we determined common DEGs by constructing a Venn diagram. Enrichment analysis was used to explore the common pathways of asthma and NSCLC. A protein-protein interaction (PPI) network was constructed to screen hub genes. KM survival analysis was performed to screen prognostic genes in the LUAD and LUSC groups. A Cox model was constructed based on hub genes and validated internally and externally. Tumor Immune Estimation Resource (TIMER) was used to evaluate the association of prognostic gene models with the tumor microenvironment (TME) and immune cell infiltration. Nomogram model was constructed by combining prognostic genes and clinical features. 114 common DEGs were obtained based on asthma and NSCLC data, and enrichment analysis showed that significant enrichment pathways mainly focused on inflammatory pathways. Screening of 5 hub genes as a key prognostic gene model for asthma progression to LUAD, and internal and external validation led to consistent conclusions. In addition, the risk score of the 5 hub genes could be used as a tool to assess the TME and immune cell infiltration. The nomogram model constructed by combining the 5 hub genes with clinical features was accurate for LUAD. Five-hub genes enrich our understanding of the potential mechanisms by which asthma contributes to the increased risk of lung cancer.

BackgroundGrapevine is a major food crop that is affected by global climate change. Consistent with field studies, dehydration assays of grapevine leaves can reveal valuable information of the plant’s response at physiological, transcript, and protein levels. There are well-known differences in grapevine rootstocks responses to dehydration. We used time-series transcriptomic approaches combined with network analyses to elucidate and identify important physiological processes and network hubs that responded to dehydration in three different grapevine species differing in their drought tolerance.ResultsTranscriptomic analyses of the leaves of Cabernet Sauvignon, Riparia Gloire, and Ramsey were evaluated at different times during a 24-h controlled dehydration. Analysis of variance (ANOVA) revealed that approximately 11,000 transcripts changed significantly with respect to the genotype x treatment interaction term and approximately 6000 transcripts changed significantly according to the genotype x treatment x time interaction term indicating massive differential changes in gene expression over time. Standard analyses determined substantial effects on the transcript abundance of genes involved in the metabolism and signaling of two known plant stress hormones, abscisic acid (ABA) and ethylene. ABA and ethylene signaling maps were constructed and revealed specific changes in transcript abundance that were associated with the known drought tolerance of the genotypes including genes such as VviABI5, VviABF2, VviACS2, and VviWRKY22. Weighted-gene coexpression network analysis (WGCNA) confirmed these results. In particular, WGCNA identified 30 different modules, some of which had highly enriched gene ontology (GO) categories for photosynthesis, phenylpropanoid metabolism, ABA and ethylene signaling. The ABA signaling transcription factors, VviABI5 and VviABF2, were highly connected hubs in two modules, one being enriched in gaseous transport and the other in ethylene signaling. VviABI5 was distinctly correlated with an early response and high expression for the drought tolerant Ramsey and with little response from the drought sensitive Riparia Gloire. These ABA signaling transcription factors were highly connected to VviSnRK1 and other gene hubs associated with sugar, ethylene and ABA signaling.ConclusionA leaf dehydration assay provided transcriptomic evidence for differential leaf responses to dehydration between genotypes differing in their drought tolerance. WGCNA proved to be a powerful network analysis approach; it identified 30 distinct modules (networks) with highly enriched GO categories and enabled the identification of gene hubs in these modules. Some of these genes were highly connected hubs in both the ABA and ethylene signaling pathways, supporting the hypothesis that there is substantial crosstalk between the two hormone pathways. This study identifies solid gene candidates for future investigations of drought tolerance in grapevine.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-016-0804-6) contains supplementary material, which is available to authorized users.

Mild traumatic brain injuries (mTBIs) are one of the most prevalent neurological disorders, and humans are severely limited in their ability to repair and regenerate central nervous system (CNS) tissue postinjury. However, zebrafish (Danio rerio) maintain the remarkable ability to undergo complete and functional neuroregeneration as an adult. We wish to extend knowledge of the known mechanisms of neuroregeneration by analyzing the differentially expressed genes (DEGs) in a novel adult zebrafish model of mTBI. In this study, a rodent weight drop model of mTBI was adapted to the adult zebrafish. A memory test showed significant deficits in spatial memory in the mTBI group. We identified DEGs at 3 and 21 days postinjury (dpi) through RNA-sequencing analysis. The resulting DEGs were categorized according to gene ontology (GO) categories. At 3 dpi, GO categories consisted of peak injury response pathways. Significantly, at 21 dpi, GO categories consisted of neuroregeneration pathways. Ultimately, these results validate a novel zebrafish model of mTBI and elucidate significant DEGs of interest in CNS injury and neuroregeneration.

Molecular fingerprinting of the podocyte reveals novel gene and protein regulatory networks

Motivation: Several functional gene annotation databases have been developed in the recent years, and are widely used to infer the biological function of gene sets, by scrutinizing the attributes that appear over- and underrepresented. However, this strategy is not directly applicable to the study of non-coding DNA, as the non-coding sequence span varies greatly among different gene loci in the human genome and longer loci have a higher likelihood of being selected purely by chance. Therefore, conclusions involving the function of non-coding elements that are drawn based on the annotation of neighboring genes are often biased. We assessed the systematic bias in several particular Gene Ontology (GO) categories using the standard hypergeometric test, by randomly sampling non-coding elements from the human genome and inferring their function based on the functional annotation of the closest genes. While no category is expected to occur significantly over- or underrepresented for a random selection of elements, categories such as ‘cell adhesion’, ‘nervous system development’ and ‘transcription factor activities’ appeared to be systematically overrepresented, while others such as ‘olfactory receptor activity’—underrepresented. Results: Our results suggest that functional inference for non-coding elements using gene annotation databases requires a special correction. We introduce a set of correction coefficients for the probabilities of the GO categories that accounts for the variability in the length of the non-coding DNA across different loci and effectively eliminates the ascertainment bias from the functional characterization of non-coding elements. Our approach can be easily generalized to any other gene annotation database. Contact: ovcharei@ncbi.nlm.nih.gov Supplementary information: Supplementary data are available at Bioinformatics Online.

Soil salinity affects the plant growth and productivity detrimentally, but Solanum chilense, a wild relative of cultivated tomato (Solanum lycopersicum L.), is known to have exceptional salt tolerance. It has precise adaptations against direct exposure to salt stress conditions. Hence, a better understanding of the mechanism to salinity stress tolerance by S. chilense can be accomplished by comprehensive gene expression studies. In this study 1-month-old seedlings of S. chilense and S. lycopersicum were subjected to salinity stress through application of sodium chloride (NaCl) solution. Through RNA-sequencing here we have studied the differences in the gene expression patterns. A total of 386 million clean reads were obtained through RNAseq analysis using the Illumina HiSeq 2000 platform. Clean reads were further assembled de novo into a transcriptome dataset comprising of 514,747 unigenes with N50 length of 578 bp and were further aligned to the public databases. Genebank non-redundant (Nr), Viridiplantae, Gene Ontology (GO), KOG, and KEGG databases classification suggested enrichment of these unigenes in 30 GO categories, 26 KOG, and 127 pathways, respectively. Out of 265,158 genes that were differentially expressed in response to salt treatment, 134,566 and 130,592 genes were significantly up and down-regulated, respectively. Upon placing all the differentially expressed genes (DEG) in known signaling pathways, it was evident that most of the DEGs involved in cytokinin, ethylene, auxin, abscisic acid, gibberellin, and Ca2+ mediated signaling pathways were up-regulated. Furthermore, GO enrichment analysis was performed using REVIGO and up-regulation of multiple genes involved in various biological processes in chilense under salinity were identified. Through pathway analysis of DEGs, “Wnt signaling pathway” was identified as a novel pathway for the response to the salinity stress. Moreover, key genes for salinity tolerance, such as genes encoding proline and arginine metabolism, ROS scavenging system, transporters, osmotic regulation, defense and stress response, homeostasis and transcription factors were not only salt-induced but also showed higher expression in S. chilense as compared to S. lycopersicum. Thus indicating that these genes may have an important role in salinity tolerance in S. chilense. Overall, the results of this study improve our understanding on possible molecular mechanisms underlying salt tolerance in plants in general and tomato in particular.

Synaptic loss and dysfunction are one of the earliest signs of neurodegeneration associated with cognitive decline in Alzheimer’s disease (AD) and other neurodegenerative diseases. This study aimed to assess the relationships between biological processes of the synaptic pathology underlying AD, molecular functions, and dynamics of the change concentrations of selected proteins reflecting synaptic and axonal pathology in dementia stages. Neurogranin (Ng), neuronal pentraxin receptor (NPTXR), and Visinin-like protein 1 (VILIP1) concentrations were measured in the cerebrospinal fluid (CSF) of MCI, AD, and non-demented controls (CTRL) using quantitative immunological methods. Gene ontology (GO) enrichment analysis was used for the functional analysis of tested proteins. The CSF Aβ42/Ng ratio was significantly different between all the compared groups. The CSF NPTXR/Ng ratio was significantly different between MCI compared to CTRL and AD compared to CTRL. The GO enrichment analysis revealed that two terms (the Biological Process (BP) and Cellular Component (CC) levels) are significantly enriched for NPTXR and Ng but not for VILIP1. Both Ng and NPTXR concentrations in CSF are promising synaptic dysfunction biomarkers for the early diagnosis of the disease. Moreover, both proteins are biochemically associated with classical biomarkers and VILIP-1. Mapping shared molecular and biological functions for the tested proteins by GO enrichment analysis may be beneficial in screening and setting new research targets.

Non-negative matrix factorization and differential expression analyses identify hub genes linked to progression and prognosis of glioblastoma multiforme

BackgroundDifferential gene expression analysis using RNA-seq data is a popular approach for discovering specific regulation mechanisms under certain environmental settings. Both gene ontology (GO) and KEGG pathway enrichment analysis are major processes for investigating gene groups that participate in common biological responses or possess related functions. However, traditional approaches based on differentially expressed genes only detect a few significant GO terms and pathways, which are frequently insufficient to explain all-inclusive gene regulation mechanisms.MethodsTranscriptomes of survivin (birc5) gene knock-down experimental and wild-type control zebrafish embryos were sequenced and assembled, and a differential expression (DE) gene list was obtained for traditional functional enrichment analysis. In addition to including DE genes with significant fold-change levels, we considered additional associated genes near or overlapped with differentially expressed long noncoding RNAs (DE lncRNAs), which may directly or indirectly activate or inhibit target genes and play important roles in regulation networks. Both the original DE gene list and the additional DE lncRNA-associated genes were combined to perform a comprehensive overrepresentation analysis.ResultsIn this study, a total of 638 DE genes and 616 DE lncRNA-associated genes (lncGenes) were leveraged simultaneously in searching for significant GO terms and KEGG pathways. Compared to the traditional approach of only using a differential expression gene list, the proposed method of employing DE lncRNA-associated genes identified several additional important GO terms and KEGG pathways. In GO enrichment analysis, 60% more GO terms were obtained, and several neuron development functional terms were retrieved as complete annotations. We also observed that additional important pathways such as the FoxO and MAPK signaling pathways were retrieved, which were shown in previous reports to play important roles in apoptosis and neuron development functions regulated by the survivin gene.ConclusionsWe demonstrated that incorporating genes near or overlapped with DE lncRNAs into the DE gene list outperformed the traditional enrichment analysis method for effective biological functional interpretations. These hidden interactions between lncRNAs and target genes could facilitate more comprehensive analyses.

Advances in mass spectrometry-based proteomics have yielded a substantial mapping of the tyrosine phosphoproteome and thus provided an important step toward a systematic analysis of intracellular signaling networks in higher eukaryotes. In this study we decomposed an uncharacterized proteomics data set of 481 unique phosphotyrosine (Tyr(P)) peptides by sequence similarity to known ligands of the Src homology 2 (SH2) and the phosphotyrosine binding (PTB) domains. From 20 clusters we extracted 16 known and four new interaction motifs. Using quantitative mass spectrometry we pulled down Tyr(P)-specific binding partners for peptides corresponding to the extracted motifs. We confirmed numerous previously known interaction motifs and found 15 new interactions mediated by phosphosites not previously known to bind SH2 or PTB. Remarkably, a novel hydrophobic N-terminal motif ((L/V/I)(L/V/I)pY) was identified and validated as a binding motif for the SH2 domain-containing inositol phosphatase SHIP2. Our decomposition of the in vivo Tyr(P) proteome furthermore suggests that two-thirds of the Tyr(P) sites mediate interaction, whereas the remaining third govern processes such as enzyme activation and nucleic acid binding.

Identification of Hub genes associated with infection of three lung cell lines by SARS-CoV-2 with integrated bioinformatics analysis.

The aim of the present study was to identify differentially expressed (DE) genes in patients with osteoarthritis (OA), and biological processes associated with changes in gene expression that occur in this disease. Using the INMEX (integrative meta-analysis of expression data) software tool, a meta-analysis of publicly available microarray Gene Expression Omnibus (GEO) datasets of OA was performed. Gene ontology (GO) enrichment analysis was performed in order to detect enriched functional attributes based on gene-associated GO terms. Three GEO datasets, containing 137 patients with OA and 52 healthy controls, were included in the meta-analysis. The analysis identified 85 genes that were consistently differentially expressed in OA (30 genes were upregulated and 55 genes were downregulated). The upregulated gene with the lowest P-value (P=5.36E-07) was S-phase kinase-associated protein 2, E3 ubiquitin protein ligase (SKP2). The downregulated gene with the lowest P-value (P=4.42E-09) was Proline rich 5 like (PRR5L). Among the 210 GO terms that were associated with the set of DE genes, the most significant two enrichments were observed in the GO categories of 'Immune response', with a P-value of 0.000129438, and 'Immune effectors process', with a P-value of 0.000288619. The current meta-analysis identified genes that were consistently DE in OA, in addition to biological pathways associated with changes in gene expression that occur during OA, which may provide insight into the molecular mechanisms underlying the pathogenesis of this disease.

Tongue squamous cell carcinoma (TSCC) is significantly more malignant than other type of oral squamous cell carcinoma (OSCC). In this study, we aimed to identify specific global gene expression signatures of TSCC to investigate the more invasive behavior of the deeply infiltrating cancer. Using RNA-seq technology, we detected gene expression of 20 TSCCs, 20 matched paratumor tissues, and 10 healthy normal mucosa tissues. Enrichment analysis of gene ontology (GO) and pathway was conducted using online tools DAVID for the dysregulated genes. Additionally, we performed the quantitative real-time RT-PCR (qRT-PCR) to validate the findings of RNA-Seq in 10 samples of TSCC, matched paratumor, and normal mucosa, respectively. We detected 252 differentially expressed genes (DEGs) between TSCC and matched paratumor tissue, including 117 up-regulated and 135 down-regulated genes. For comparison between TSCC and normal mucosa, 234 DEGS were identified, consisting of 67 up-regulated and 167 down-regulated genes. For both two comparisons, GO categories of muscle contraction (GO: 0006936), epidermis development (GO: 0008544), epithelial cell differentiation (GO: 0030855), and keratinization (GO: 0031424) were commonly enriched. Altered gene expression affected some cancer-related pathways, such as tight junction. The qRT-PCR validation showed that gene expression patterns of FOLR1, NKX3-1, TFF3, PIGR, NEFL, MMP13, and HMGA2 were fully in concordance with RNA-Seq results. Findings in this study demonstrated the genetic and molecular alterations associated with TSCC, providing new clues for understanding the molecular mechanisms of TSCC pathogenesis.