Comparison study of microarray meta-analysis methods

Microarray meta-analysis reveals comprehensive effects of 3,4,5-tricaffeolyquinic acid in cell differentiation and signaling

Caffeoylquinic acids (CQA) are polyphenolic compounds found in fruits, vegetables, coffee, and spices that have exhibited several beneficial activities, including antioxidant, antibacterial, neuroprotective, anti-inflammatory, anticancer, antiviral, antidiabetic, and cardiovascular effects. A derivative, TCQA (3,4,5-Tri-O-caffeoylquinic acid), has also shown both neurogenic and pigment differentiation potential. A transcriptomic-based meta-analysis was conducted to explore potential biochemical processes and molecular targets of TCQA. This approach involved integrating data from various cell and tissue types, including human amniotic stem cells, human neural stem cells, human dermal papilla cells, and the brain cortex of aging model mice. It offered a comprehensive perspective on the significant gene regulations in response to TCQA treatment. The objective was to uncover the mechanism and novel targets of TCQA, facilitating a further understanding of its functions. New areas of interest found were TCQA's effect on adipogenesis, heart, and muscle tissue development. In addition, significantly enhanced biological activities found through meta-analysis included cell cycle, VEGFA-VEGFR2 pathway, and BMP signaling. Overall, a comprehensive functional and visual analysis using available biological databases uncovered the multi-target potential of this natural compound.

CoSIA: an R Bioconductor package for CrOss Species Investigation and Analysis

SUMMARY

High-throughput sequencing technologies have enabled cross-species comparative transcriptomic studies; however, there are numerous challenges for these studies due to biological and technical factors. We developed CoSIA (Cross-Species Investigation and Analysis), a Bioconductor R package and Shiny app that provides an alternative framework for cross-species transcriptomic comparison of non-diseased wild-type RNA sequencing gene expression data from Bgee across tissues and species (human, mouse, rat, zebrafish, fly, and nematode) through visualization of variability, diversity, and specificity metrics.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lasseignelab/CoSIA.

CoSIA: an R Bioconductor package for CrOss Species Investigation and Analysis

High throughput sequencing technologies have enabled cross-species comparative transcriptomic studies; however, there are numerous challenges for these studies due to biological and technical factors. We developed CoSIA (Cross-Species Investigation and Analysis), an Bioconductor R package and Shiny app that provides an alternative framework for cross-species transcriptomic comparison of non-diseased wild-type RNA sequencing gene expression data from Bgee across tissues and species (human, mouse, rat, zebrafish, fly, and nematode) through visualization of variability, diversity, and specificity metrics.

Transcriptomics of the Carotid Body

The carotid body (CB) has emerged as a potential therapeutic target for treating sympathetically mediated cardiovascular, respiratory, and metabolic diseases. In adjunct to its classical role as an arterial O₂ sensor, the CB is a multimodal sensor activated by a range of stimuli in the circulation. However, consensus on how CB multimodality is achieved is lacking; even the best studied O₂-sensing appears to involve multiple convergent mechanisms. A strategy to understand multimodal sensing is to adopt a hypothesis-free, high-throughput transcriptomic approach. This has proven instrumental for understanding fundamental mechanisms of CB response to hypoxia and other stimulants, its developmental niche, cellular heterogeneity, laterality, and pathophysiological remodeling in disease states. Herein, we review this published work that reveals novel molecular mechanisms underpinning multimodal sensing and reveals numerous gaps in knowledge that require experimental testing.

Meta-analysis of gene expression data in adipose tissue reveals new obesity associated genes

High-throughput transcriptomic and proteomic data like microarray data are deposited in public databases such as Gene Expression Omnibus (GEO). Omics data integration and processing from different and independent studies is achieved by using efficient and effective computational tools through meta-analysis. Meta-analysis is a statistical powerful tool combining data from numerous studies, minimizes bias and increases statistical power by increasing sample size compared to individual studies. Therefore, we performed a meta-analysis of gene expression data in adipose tissue to identify genes that are differentially expressed between obese and non-obese subjects as well as to detect gene expression signatures, pathways and networks associated with obesity. We identified 821 differentially expressed genes (DEGs) in adipose tissue of obese subjects compared to non-obese. A protein-protein interactions (PPIs) network was reconstructed consisting of 168 proteins. Functional enrichment analysis in the network revealed proteins involved in RNA and energy metabolism. The KEGG pathway analysis revealed 15 enriched pathway terms. Furthermore, multiple testing correction methods identified five statistically significant obesity associated genes (NDUFA12, SFI1, SSB, FAR2 and LACE1) that require further investigation.

A systems biology study unveils the association between a melatonin biosynthesis gene, O-methyl transferase 1 (OMT1) and wheat (Triticum aestivum L.) combined drought and salinity stress tolerance

Enhanced levels of endogenous melatonin in the root of wheat, mainly through the OMT1 gene, augment the antioxidant system, reestablish redox homeostasis and are associated with combined stress tolerance. A systems biology approach, including a collection of computational analyses and experimental assays, led us to uncover some aspects of a poorly understood phenomenon, namely wheat (Triticum aestivum L.) combined drought and salinity stress tolerance. Accordingly, a cross-study comparison of stress experiments was performed via a meta-analysis of Expressed Sequence Tags (ESTs) data from wheat roots to uncover the overlapping gene network of drought and salinity stresses. Identified differentially expressed genes were functionally annotated by gene ontology enrichment analysis and gene network analysis. Among those genes, O-methyl transferase 1 (OMT1) was highlighted as a more important (hub) gene in the dual-stress response gene network. Afterwards, the potential roles of OMT1 in mediating physiochemical indicators of stress tolerance were investigated in two wheat genotypes differing in abiotic stress tolerance. Regression analysis and correspondence analysis (CA) confirmed that the expression profiles of the OMT1 gene and variations in melatonin content, antioxidant enzyme activities, proline accumulation, H₂O₂ and malondialdehyde (MDA) contents are significantly associated with combined stress tolerance. These results reveal that the OMT1 gene may contribute to wheat combined drought and salinity stress tolerance through augmenting the antioxidant system and re-establishing redox homeostasis, probably via the regulation of melatonin biosynthesis as a master regulator molecule. Our findings provide new insights into the roles of melatonin in wheat combined drought and salinity stress tolerance and suggest a novel plausible regulatory node through the OMT1 gene to improve multiple-stress tolerant crops.

Reviewing and assessing existing meta-analysis models and tools

Over the past few years, meta-analysis has become popular among biomedical researchers for detecting biomarkers across multiple cohort studies with increased predictive power. Combining datasets from different sources increases sample size, thus overcoming the issue related to limited sample size from each individual study and boosting the predictive power. This leads to an increased likelihood of more accurately predicting differentially expressed genes/proteins or significant biomarkers underlying the biological condition of interest. Currently, several meta-analysis methods and tools exist, each having its own strengths and limitations. In this paper, we survey existing meta-analysis methods, and assess the performance of different methods based on results from different datasets as well as assessment from prior knowledge of each method. This provides a reference summary of meta-analysis models and tools, which helps to guide end-users on the choice of appropriate models or tools for given types of datasets and enables developers to consider current advances when planning the development of new meta-analysis models and more practical integrative tools.

Comprehensive Analysis of Differentially Expressed Genes in Clinically Diagnosed Irreversible Pulpitis by Multiplatform Data Integration Using a Robust Rank Aggregation Approach

INTRODUCTION

Molecular diagnosis may overcome the limitations of clinical and histologic diagnosis in pulpitis, thereby benefiting many treatment techniques, such as vital pulp therapies. In this study, integrated microarray data on pulpitis were used to obtain a list of normalized differentially expressed (DE) genes for analyzing the molecular mechanisms underlying pulpitis and identifying potential diagnostic biomarkers.

METHODS

A systematic search of public microarray and sequencing databases was performed to obtain expression data of pulpitis. Robust rank aggregation (RRA) was used to obtain DE gene lists (RRA_DEmRNAs and RRA_DElncRNAs) between inflamed pulp and normal samples. DE genes were evaluated by functional enrichment analyses, correlation analyses for inflammation-related RRA_DEmRNAs, and protein-protein interaction and competing endogenous RNA network construction. Quantitative real-time polymerase chain reaction validation was applied in snap-frozen pulp tissues.

RESULTS

Using the GSE77459 and GSE92681 data sets, 280 RRA_DEmRNAs and 90 RRA_DElncRNAs were identified. RRA_DEmRNAs were significantly enriched in inflammation-related biological processes and osteoclast differentiation and tumor necrosis factor, chemokine, and B-cell receptor signaling pathways. The molecular complex detection and cytoHubba methods identified 2 clusters and 10 hub genes in the protein-protein interaction network. The competing endogenous RNA network was composed of 2 long noncoding RNAs (ADAMTS9-AS2 and LINC00290), 2 microRNAs (hsa-miR-30a-5p and hsa-miR-128-3p), and 3 messenger RNAs (ABCA1, FBLN5, and SOCS3). The expression between most top inflammation-related RRA_DEmRNAs in pulpitis showed positive correlations. Quantitative real-time polymerase chain reacation validated the expression trends of selected genes, including ITGAX, TREM1, CD86, FCGR2A, ADAMTS9-AS2, LINC00290, hsa-miR-30a-5p, hsa-miR-128-3p, RASGRP3, IL3RA, CCDC178, CRISPLD1, LINC01857, AC007991.2, ARHGEF26-AS1, and AL021408.1.

CONCLUSIONS

The identified biomarkers provide insight into the pathology and may aid in the molecular diagnosis of pulpitis.

Identification of genes and pathways related to breast cancer metastasis in an integrated cohort

BACKGROUND

Breast cancer is the most common malignant disease in women. Metastasis is the most common cause of death from this cancer. Screening genes related to breast cancer metastasis may help elucidate the mechanisms governing metastasis and identify molecular targets for antimetastatic therapy. The development of advanced algorithms enables us to perform cross-study analysis to improve the robustness of the results.

MATERIALS AND METHODS

Ten data sets meeting our criteria for differential expression analyses were obtained from the Gene Expression Omnibus (GEO) database. Among these data sets, five based on the same platform were formed into a large cohort using the XPN algorithm. Differentially expressed genes (DEGs) associated with breast cancer metastasis were identified using the differential expression via distance synthesis (DEDS) algorithm. A cross-platform method was employed to verify these DEGs in all ten selected data sets. The top 50 validated DEGs are represented with heat maps. Based on the validated DEGs, Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Protein interaction (PPI) networks were constructed to further illustrate the direct and indirect associations among the DEGs. Survival analysis was performed to explore whether these genes can affect breast cancer patient prognosis.

RESULTS

A total of 817 DEGs were identified using the DEDS algorithm. Of these DEGs, 450 genes were validated by the second algorithm. Enriched KEGG pathway terms demonstrated that these 450 DEGs may be involved in the cell cycle and oocyte meiosis in addition to their functions in ECM-receptor interaction and protein digestion and absorption. PPI network analysis for the proteins encoded by the DEGs indicated that these genes may be primarily involved in the cell cycle and extracellular matrix. In particular, several genes played roles in multiple signalling pathways and were related to patient survival. These genes were also observed to be targetable in the CTD2 database.

CONCLUSIONS

Our study analysed multiple cross-platform data sets using two different algorithms, helping elucidate the molecular mechanisms and identify several potential therapeutic targets of metastatic breast cancer. In addition, several genes exhibited promise for applications in targeted therapy against metastasis in future research.

Rubisco activase A (RcaA) is a central node in overlapping gene network of drought and salinity in Barley (Hordeum vulgare L.) and may contribute to combined stress tolerance

Co-occurrence of abiotic stresses, especially drought and salinity, is a natural phenomenon in field conditions and is worse for crop production than any single stress. Nowadays, rigorous methods of meta-analysis and systems biology have made it possible to perform cross-study comparisons of single stress experiments, which can uncover main overlapping mechanisms underlying tolerance to combined stress. In this study, a meta-analysis of RNA-Seq data was conducted to obtain the overlapping gene network of drought and salinity stresses in barley (Hordeum vulgare L.), which identified Rubisco activase A (RcaA) as a hub gene in the dual-stress response. Thereafter, a greenhouse experiment was carried out using two barley genotypes with different abiotic stress tolerance and evaluated several physiochemical properties as well as the expression profile and protein activity of RcaA. Finally, machine learning analysis was applied to uncover relationships among combined stress tolerance and evaluated properties. We identified 441 genes which were differentially expressed under both drought and salinity stress. Results revealed that the photosynthesis pathway and, in particular, the RcaA gene are major components of the dual-stress responsive transcriptome. Comparative physiochemical and molecular evaluations further confirmed that enhanced photosynthesis capability, mainly through regulation of RcaA expression and activity as well as accumulation of proline content, have a significant association with combined drought and salinity stress tolerance in barley. Overall, our results clarify the importance of RcaA in combined stress tolerance and may provide new insights for future investigations.

Comparative transcriptome analysis of Parkinson's disease and Hutchinson-Gilford progeria syndrome reveals shared susceptible cellular network processes

BACKGROUND

Parkinson's Disease (PD) and Hutchinson-Gilford Progeria Syndrome (HGPS) are two heterogeneous disorders, which both display molecular and clinical alterations associated with the aging process. However, similarities and differences between molecular changes in these two disorders have not yet been investigated systematically at the level of individual biomolecules and shared molecular network alterations.

METHODS

Here, we perform a comparative meta-analysis and network analysis of human transcriptomics data from case-control studies for both diseases to investigate common susceptibility genes and sub-networks in PD and HGPS. Alzheimer's disease (AD) and primary melanoma (PM) were included as controls to confirm that the identified overlapping susceptibility genes for PD and HGPS are non-generic.

RESULTS

We find statistically significant, overlapping genes and cellular processes with significant alterations in both diseases. Interestingly, the majority of these shared affected genes display changes with opposite directionality, indicating that shared susceptible cellular processes undergo different mechanistic changes in PD and HGPS. A complementary regulatory network analysis also reveals that the altered genes in PD and HGPS both contain targets controlled by the upstream regulator CDC5L.

CONCLUSIONS

Overall, our analyses reveal a significant overlap of affected cellular processes and molecular sub-networks in PD and HGPS, including changes in aging-related processes that may reflect key susceptibility factors associated with age-related risk for PD.

Dysregulation of microRNA in cholangiocarcinoma identified through a meta-analysis of microRNA profiling

BACKGROUND

In the past decades, the potential of microRNA (miRNA) in cancer diagnostics and prognostics has gained a lot of interests. In this study, a meta-analysis was conducted upon the pooled miRNA microarray data of cholangiocarcinoma (CCA).

AIM

To identify differentially expressed (DE) miRNAs and perform functional analyses in order to gain insights to understanding miRNA-target interactions involved in tumorigenesis pathways of CCA.

METHODS

Raw data from 8 CCA miRNA microarray datasets, consisting of 443 samples in total, were integrated and statistically analyzed to identify DE miRNAs via comparison of levels of miRNA expression between CCA and normal bile duct samples using t-tests (P < 0.001). The 10-fold cross validation was performed in order to increase the robustness of the t-test results.

RESULTS

Our data showed 70 up-regulated and 48 down-regulated miRNAs in CCA. Gene Ontology and pathway enrichment analyses revealed that mRNA targets of DE miRNAs were significantly involved in several biological processes. The most prominent dysregulated pathways included phosphatidylinositol-3 kinases/Akt, mitogen-activated protein kinase and Ras signaling pathways.

CONCLUSION

DE miRNAs found in our meta-analysis revealed dysregulation in major cancer pathways involved in the development of CCA. These results indicated the necessity of understanding the miRNA-target interactions and the significance of dysregulated miRNAs in terms of diagnostics and prognostics of cancers.

Cross-Species Annotation of Expressed Genes and Detection of Different Functional Gene Modules Between 10 Cold- and 10 Hot-Propertied Chinese Herbal Medicines

According to the traditional Chinese medicine (TCM) system, Chinese herbal medicines (HMs) can be divided into four categories: hot, warm, cold, and cool. A cool nature usually is categorized as a cold nature, and a warm nature is classified as a hot nature. However, the detectable characteristics of the gene expression profile associated with the cold and hot properties have not been studied. To address this question, a strategy for the cross-species annotation of conserved genes was established in the present study by using transcriptome data of 20 HMs with cold and hot properties. Functional enrichment analysis was performed on group-specific expressed genes inferred from the functional genome of the reference species (i.e., Arabidopsis). Results showed that metabolic pathways relevant to chrysoeriol, luteolin, paniculatin, and wogonin were enriched for cold-specific genes, and pathways of inositol, heptadecane, lauric acid, octanoic acid, hexadecanoic acid, and pentadecanoic acid were enriched for hot-specific genes. Six functional modules were identified in the HMs with the cold property: nucleotide biosynthetic process, peptidy-L-cysteine S-palmitoylation, lipid modification, base-excision repair, dipeptide transport, and response to endoplasmic reticulum stress. For the hot HMs, another six functional modules were identified: embryonic meristem development, embryonic pattern specification, axis specification, regulation of RNA polymerase II transcriptional preinitiation complex assembly, mitochondrial RNA modification, and cell redox homeostasis. The research provided a new insight into HMs’ cold and hot properties from the perspective of the gene expression profile of plants.

Twist-1 is upregulated by NSD2 and contributes to tumour dissemination and an epithelial-mesenchymal transition-like gene expression signature in t(4;14)-positive multiple myeloma

Approximately 15% of patients with multiple myeloma (MM) harbour the t(4;14) chromosomal translocation, leading to the overexpression of the histone methyltransferase NSD2. Patients with this translocation display increased tumour dissemination, accelerated disease progression and rapid relapse. Using publicly available gene expression profile data from NSD2^high (n = 135) and NSD2^low (n = 878) MM patients, we identified 39 epithelial-mesenchymal transition (EMT)-associated genes which are overexpressed in NSD2^high MM plasma cells. In addition, our analyses identified Twist-1 as a key transcription factor upregulated in NSD2^high MM patients and t(4;14)-positive cell lines. Overexpression and knockdown studies confirmed that Twist-1 is involved in driving the expression of EMT-associated genes in the human MM cell line KMS11 and promoted the migration of myeloma cell lines in vitro. Notably, Twist-1 overexpression in the mouse MM cell line 5TGM1 significantly increased tumour dissemination in an intratibial tumour model. These findings demonstrate that Twist-1, downstream of NSD2, contributes to the induction of an EMT-like signature in t(4;14)-positive MM and enhances the dissemination of MM plasma cells in vivo, which may, in part, explain the aggressive disease features associated with t(4;14)-positive MM.

Integrated PPI- and WGCNA-Retrieval of Hub Gene Signatures Shared Between Barrett's Esophagus and Esophageal Adenocarcinoma

Esophageal adenocarcinoma (EAC) is a deadly cancer with high mortality rate, especially in economically advanced countries, while Barrett's esophagus (BE) is reported to be a precursor that strongly increases the risk of EAC. Due to the complexity of these diseases, their molecular mechanisms have not been revealed clearly. This study aims to explore the gene signatures shared between BE and EAC based on integrated network analysis. We obtained EAC- and BE-associated microarray datasets GSE26886, GSE1420, GSE37200, and GSE37203 from the Gene Expression Omnibus and ArrayExpress using systematic meta-analysis. These data were accompanied by clinical data and RNAseq data from The Cancer Genome Atlas (TCGA). Weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) analysis were conducted to explore the relationship between gene sets and clinical traits as well as to discover the key relationships behind the co-expression modules. A differentially expressed gene-based protein-protein interaction (PPI) complex was used to extract hub genes through Cytoscape plugins. As a result, 403 DEGs were excavated, comprising 236 upregulated and 167 downregulated genes, which are involved in the cell cycle and replication pathways. Forty key genes were identified using modules of MCODE, CytoHubba, and CytoNCA with different algorithms. A dark-gray module with 207 genes was identified which having a high correlation with phenotype (gender) in the WGCNA. Furthermore, five shared hub gene signatures (SHGS), namely, pre-mRNA processing factor 4 (PRPF4), serine and arginine-rich splicing factor 1 (SRSF1), heterogeneous nuclear ribonucleoprotein M (HNRNPM), DExH-Box Helicase 9 (DHX9), and origin recognition complex subunit 2 (ORC2), were identified between BE and EAC. SHGS enrichment denotes that RNA metabolism and splicosomes play a key role in esophageal cancer development and progress. We conclude that the PPI complex and WGCNA co-expression network highlight the importance of phenotypic identifying hub gene signatures for BE and EAC.

Induction of PrMADS10 on the lower side of bent pine tree stems: potential role in modifying plant cell wall properties and wood anatomy

The molecular mechanisms underlying inclination responses in trees are unclear. In this study, we identified a MADS-box transcription factor differentially expressed early after inclination in the stems of Pinus radiata D. Don. PrMADS10 has a CDS of 582 bp and encodes a group II MADS-box transcription factor. We measured highest accumulation of this transcript on the lower side of inclined pine stems. In an effort to identify putative targets, we stably transformed Arabidopsis thaliana with a 35S::PrMADS10 construct. Transcriptome analysis revealed 1,219 genes differentially-expressed, with 690 and 529 genes up- and down-regulated respectively, when comparing the transgenic and wild-type. Differentially-expressed genes belong to different biological processes, but were enriched in cell wall remodeling and phenylpropanoid metabolic functions. Interestingly, lignin content was 30% higher in transgenic as compared to wild-type plants consistent with observed changes in gene expression. Differentially expressed transcription factors and phenylpropanoid genes were analyzed using STRING. Several MYB and NAC transcription factors showed interactions with genes of the phenylpropanoid pathway. Together, these results implicate PrMADS10 as a regulatory factor, triggering the expression of other transcription factors and genes involved in the synthesis of lignin.

Leveraging Big Data to Transform Drug Discovery

The surge of public disease and drug-related data availability has facilitated the application of computational methodologies to transform drug discovery. In the current chapter, we outline and detail the various resources and tools one can leverage in order to perform such analyses. We further describe in depth the in silico workflows of two recent studies that have identified possible novel indications of existing drugs. Lastly, we delve into the caveats and considerations of this process to enable other researchers to perform rigorous computational drug discovery experiments of their own.

Identification of gene expression profiles in myocardial infarction: a systematic review and meta-analysis

BACKGROUND

Myocardial infarction (MI) is a multifactorial disease with complex pathogenesis, mainly the result of the interplay of genetic and environmental risk factors. The regulation of thrombosis, inflammation and cholesterol and lipid metabolism are the main factors that have been proposed thus far to be involved in the pathogenesis of MI. Traditional risk-estimation tools depend largely on conventional risk factors but there is a need for identification of novel biochemical and genetic markers. The aim of the study is to identify differentially expressed genes that are consistently associated with the incidence myocardial infarction (MI), which could be potentially incorporated into the traditional cardiovascular diseases risk factors models.

METHODS

The biomedical literature and gene expression databases, PubMed and GEO, respectively, were searched following the PRISMA guidelines. The key inclusion criteria were gene expression data derived from case-control studies on MI patients from blood samples. Gene expression datasets regarding the effect of medicinal drugs on MI were excluded. The t-test was applied to gene expression data from case-control studies in MI patients.

RESULTS

A total of 162 articles and 174 gene expression datasets were retrieved. Of those a total of 4 gene expression datasets met the inclusion criteria, which contained data on 31,180 loci in 93 MI patients and 89 healthy individuals. Collectively, 626 differentially expressed genes were detected in MI patients as compared to non-affected individuals at an FDR q-value = 0.01. Of those, 88 genes/gene products were interconnected in an interaction network. Totally, 15 genes were identified as hubs of the network.

CONCLUSIONS

Functional enrichment analyses revealed that the DEGs and that they are mainly involved in inflammatory/wound healing, RNA processing/transport mechanisms and a yet not fully characterized pathway implicated in RNA transport and nuclear pore proteins. The overlap between the DEGs identified in this study and the genes identified through genetic-association studies is minimal. These data could be useful in future studies on the molecular mechanisms of MI as well as diagnostic and prognostic markers.

Unified Transcriptomic Signature of Arbuscular Mycorrhiza Colonization in Roots of Medicago truncatula by Integration of Machine Learning, Promoter Analysis, and Direct Merging Meta-Analysis

Plant root symbiosis with Arbuscular mycorrhizal (AM) fungi improves uptake of water and mineral nutrients, improving plant development under stressful conditions. Unraveling the unified transcriptomic signature of a successful colonization provides a better understanding of symbiosis. We developed a framework for finding the transcriptomic signature of Arbuscular mycorrhiza colonization and its regulating transcription factors in roots of Medicago truncatula. Expression profiles of roots in response to AM species were collected from four separate studies and were combined by direct merging meta-analysis. Batch effect, the major concern in expression meta-analysis, was reduced by three normalization steps: Robust Multi-array Average algorithm, Z-standardization, and quartiling normalization. Then, expression profile of 33685 genes in 18 root samples of Medicago as numerical features, as well as study ID and Arbuscular mycorrhiza type as categorical features, were mined by seven models: RELIEF, UNCERTAINTY, GINI INDEX, Chi Squared, RULE, INFO GAIN, and INFO GAIN RATIO. In total, 73 genes selected by machine learning models were up-regulated in response to AM (Z-value difference > 0.5). Feature weighting models also documented that this signature is independent from study (batch) effect. The AM inoculation signature obtained was able to differentiate efficiently between AM inoculated and non-inoculated samples. The AP2 domain class transcription factor, GRAS family transcription factors, and cyclin-dependent kinase were among the highly expressed meta-genes identified in the signature. We found high correspondence between the AM colonization signature obtained in this study and independent RNA-seq experiments on AM colonization, validating the repeatability of the colonization signature. Promoter analysis of upregulated genes in the transcriptomic signature led to the key regulators of AM colonization, including the essential transcription factors for endosymbiosis establishment and development such as NF-YA factors. The approach developed in this study offers three distinct novel features: (I) it improves direct merging meta-analysis by integrating supervised machine learning models and normalization steps to reduce study-specific batch effects; (II) seven attribute weighting models assessed the suitability of each gene for the transcriptomic signature which contributes to robustness of the signature (III) the approach is justifiable, easy to apply, and useful in practice. Our integrative framework of meta-analysis, promoter analysis, and machine learning provides a foundation to reveal the transcriptomic signature and regulatory circuits governing Arbuscular mycorrhizal symbiosis and is transferable to the other biological settings.

Identification of target gene and prognostic evaluation for lung adenocarcinoma using gene expression meta-analysis, network analysis and neural network algorithms

BACKGROUND

Lung adenocarcinoma (LUAD) is a heterogeneous disease with poor survival in the advanced stage and a high incidence rate in the world. Novel drug targets are urgently required to improve patient treatment. Therefore, we aimed to identify therapeutic targets for LUAD based on protein-protein and protein-drug interaction network analysis with neural network algorithms using mRNA expression profiles.

RESULTS

A comprehensive meta-analysis of selective non-small cell lung cancer (NSCLC) mRNA expression profile datasets from Gene Expression Omnibus were used to identify potential biomarkers and the molecular mechanisms related to the prognosis of NSCLC patients. Using the Network Analyst tool, based on combined effect size (ES) methods, we recognized 6566 differentially expressed genes (DEGs), which included 3036 downregulated and 3530 upregulated genes linked to NSCLC patient survival. ClueGO, a Cytoscape plugin, was exploited to complete the function and pathway enrichment analysis, which disclosed "regulated exocytosis", "purine nucleotide binding", "pathways in cancer", and "cell cycle" between exceptionally supplemented terms. Enrichr, a web tool examination, demonstrated "early growth response protein 1 (EGR-1)", "hepatocyte nuclear factor 4α (HNF4A)", "mitogen-activated protein kinase 14 (MAP3K14)", and "cyclin-dependent kinase 1 (CDK1)" to be among the most prevalent TFs and kinases associated with NSCLC. Our meta-analysis identified that MAPK1 and aurora kinase (AURKA) are the most obvious class of hub nodes. Furthermore, protein-drug interaction network and neural network algorithms identified candidate drugs such as phosphothreonine and 4-(4-methylpiperazin-1-yl)-n-[5-(2-thienylacetyl)-1,5-dihydropyrrolo[3,4-c]pyrazol-3-yl] benzamide and for the targets MAPK1 and AURKA, respectively.

CONCLUSION

Our study has identified novel candidate biomarkers, pathways, transcription factors (TFs), and kinases associated with NSCLC prognosis, as well as drug candidates, which may assist treatment strategy for NSCLC patients.

Meta-analysis of transcriptomic responses to biotic and abiotic stress in tomato

A wide range of biotic stresses (BS) and abiotic stresses (AS) adversely affect plant growth and productivity worldwide. The study of individual genes cannot be considered as an effective approach for the understanding of tolerance mechanisms, since these stresses are frequent and often in combination with each other, and a large number of genes are involved in these mechanisms. The availability of high-throughput genomic data has enabled the discovery of the role of transcription factors (TFs) in regulatory networks. A meta-analysis of BS and AS responses was performed by analyzing a total of 391 microarray samples from 23 different experiments and 2,336 differentially expressed genes (DEGs) involved in multiple stresses were identified. We identified 1,862 genes differentially regulated in response to BS was much greater than that regulated by AS, 835 genes, and found 15.4% or 361 DEGs with the conserved expression between AS and BS. The greatest percent of genes related to the cellular process (>76% genes), metabolic process (>76% genes) and response to stimulus (>50%). About 4.2% of genes involved in BS and AS responses belonged to the TF families. We identified several genes, which encode TFs that play an important role in AS and BS responses. These proteins included Jasmonate Ethylene Response Factor 1 (JERF1), SlGRAS6, MYB48, SlERF4, EIL2, protein LATE ELONGATED HYPOCOTYL (LHY), SlERF1, WRKY 26, basic leucine zipper TF, inducer of CBF expression 1-like, pti6, EIL3 and WRKY 11. Six of these proteins, JERF1, MYB48, protein LHY, EIL3, EIL2 and SlGRAS6, play central roles in these mechanisms. This research promoted a new approach to clarify the expression profiles of various genes under different conditions in plants, detected common genes from differentially regulated in response to these conditions and introduced them as candidate genes for improving plant tolerance through genetic engineering approach.

Genomic Approaches to Posttraumatic Stress Disorder: The Psychiatric Genomic Consortium Initiative

Posttraumatic stress disorder (PTSD) after exposure to a traumatic event is a highly prevalent psychiatric disorder. Heritability estimates from twin studies as well as from recent molecular data (single nucleotide polymorphism-based heritability) indicate moderate to high heritability, yet robust genetic variants for PTSD have not yet been identified and the genetic architecture of this polygenic disorder remains largely unknown. To date, fewer than 10 large-scale genome-wide association studies of PTSD have been published, with findings that highlight the unique challenges for PTSD genomics, including a complex diagnostic entity with contingency of PTSD diagnosis on trauma exposure and the large genetic diversity of the study populations. The Psychiatric Genomics Consortium PTSD group has brought together more than 200 scientists with the goal to increase sample size for genome-wide association studies and other genomic analyses to sufficient numbers where robust discoveries of molecular signatures can be achieved. The sample currently includes more than 32,000 PTSD cases and 100,000 trauma-exposed control subjects, and collection is ongoing. The first results found a significant shared genetic risk of PTSD with other psychiatric disorders and sex-biased heritability estimates with higher heritability in female individuals compared with male individuals. This review describes the scope and current focus of the Psychiatric Genomics Consortium PTSD group and its expansion from the initial genome-wide association study group to nine working groups, including epigenetics, gene expression, imaging, and integrative systems biology. We further briefly outline recent findings and future directions of "omics"-based studies of PTSD, with the ultimate goal of elucidating the molecular architecture of this complex disorder to improve prevention and intervention strategies.

Multi-cohort analysis of colorectal cancer metagenome identified altered bacteria across populations and universal bacterial markers

BACKGROUND

Alterations of gut microbiota are associated with colorectal cancer (CRC) in different populations and several bacterial species were found to contribute to the tumorigenesis. The potential use of gut microbes as markers for early diagnosis has also been reported. However, cohort specific noises may distort the structure of microbial dysbiosis in CRC and lead to inconsistent results among studies. In this regard, our study targeted at exploring changes in gut microbiota that are universal across populations at species level.

RESULTS

Based on the combined analysis of 526 metagenomic samples from Chinese, Austrian, American, and German and French cohorts, seven CRC-enriched bacteria (Bacteroides fragilis, Fusobacterium nucleatum, Porphyromonas asaccharolytica, Parvimonas micra, Prevotella intermedia, Alistipes finegoldii, and Thermanaerovibrio acidaminovorans) have been identified across populations. The seven enriched bacterial markers classified cases from controls with an area under the receiver-operating characteristics curve (AUC) of 0.80 across the different populations. Abundance correlation analysis demonstrated that CRC-enriched and CRC-depleted bacteria respectively formed their own mutualistic networks, in which the latter was disjointed in CRC. The CRC-enriched bacteria have been found to be correlated with lipopolysaccharide and energy biosynthetic pathways.

CONCLUSIONS

Our study identified potential diagnostic bacterial markers that are robust across populations, indicating their potential universal use for non-invasive CRC diagnosis. We also elucidated the ecological networks and functional capacities of CRC-associated microbiota.

Methods of Analysis and Meta-Analysis for Identifying Differentially Expressed Genes

Microarray approaches are widely used high-throughput techniques to assess simultaneously the expression of thousands of genes under certain conditions and study the effects of certain treatments, diseases, and developmental stages. The traditional way to perform such experiments is to design oligonucleotide hybridization probes that correspond to specific genes and then measure the expression of the genes in order to determine which of them are up- or down-regulated compared to a condition that is used as a control. Hitherto, individual experiments cannot capture the bigger picture of how a biological system works and, therefore, data integration from multiple experimental studies and external data repositories is necessary to understand the function of genes and their expression patterns under certain conditions. Therefore, the development of methods for handling, integrating, comparing, interpreting and visualizing microarray data is necessary. The selection of an appropriate method for analysing microarray datasets is not an easy task. In this chapter, we provide an overview of the various methods developed for microarray data analysis, as well as suggestions for choosing the appropriate method for microarray meta-analysis.

Comparison of multiple transcriptomes exposes unified and divergent features of quiescent and activated skeletal muscle stem cells

BACKGROUND

Skeletal muscle satellite (stem) cells are quiescent in adult mice and can undergo multiple rounds of proliferation and self-renewal following muscle injury. Several labs have profiled transcripts of myogenic cells during the developmental and adult myogenesis with the aim of identifying quiescent markers. Here, we focused on the quiescent cell state and generated new transcriptome profiles that include subfractionations of adult satellite cell populations, and an artificially induced prenatal quiescent state, to identify core signatures for quiescent and proliferating.

METHODS

Comparison of available data offered challenges related to the inherent diversity of datasets and biological conditions. We developed a standardized workflow to homogenize the normalization, filtering, and quality control steps for the analysis of gene expression profiles allowing the identification up- and down-regulated genes and the subsequent gene set enrichment analysis. To share the analytical pipeline of this work, we developed Sherpa, an interactive Shiny server that allows multi-scale comparisons for extraction of desired gene sets from the analyzed datasets. This tool is adaptable to cell populations in other contexts and tissues.

RESULTS

A multi-scale analysis comprising eight datasets of quiescent satellite cells had 207 and 542 genes commonly up- and down-regulated, respectively. Shared up-regulated gene sets include an over-representation of the TNFα pathway via NFKβ signaling, Il6-Jak-Stat3 signaling, and the apical surface processes, while shared down-regulated gene sets exhibited an over-representation of Myc and E2F targets and genes associated to the G2M checkpoint and oxidative phosphorylation. However, virtually all datasets contained genes that are associated with activation or cell cycle entry, such as the immediate early stress response genes Fos and Jun. An empirical examination of fixed and isolated satellite cells showed that these and other genes were absent in vivo, but activated during procedural isolation of cells.

CONCLUSIONS

Through the systematic comparison and individual analysis of diverse transcriptomic profiles, we identified genes that were consistently differentially expressed among the different datasets and shared underlying biological processes key to the quiescent cell state. Our findings provide impetus to define and distinguish transcripts associated with true in vivo quiescence from those that are first responding genes due to disruption of the stem cell niche.

Microarray meta-analysis to explore abiotic stress-specific gene expression patterns in Arabidopsis

BACKGROUND

Abiotic stresses are the major limiting factors that affect plant growth, development, yield and final quality. Deciphering the underlying mechanisms of plants' adaptations to stresses using few datasets might overlook the different aspects of stress tolerance in plants, which might be simultaneously and consequently operated in the system. Fortunately, the accumulated microarray expression data offer an opportunity to infer abiotic stress-specific gene expression patterns through meta-analysis. In this study, we propose to combine microarray gene expression data under control, cold, drought, heat, and salt conditions and determined modules (gene sets) of genes highly associated with each other according to the observed expression data.

RESULTS

By analyzing the expression variations of the Eigen genes from different conditions, we had identified two, three, and five gene modules as cold-, heat-, and salt-specific modules, respectively. Most of the cold- or heat-specific modules were differentially expressed to a particular degree in shoot samples, while most of the salt-specific modules were differentially expressed to a particular degree in root samples. A gene ontology (GO) analysis on the stress-specific modules suggested that the gene modules exclusively enriched stress-related GO terms and that different genes under the same GO terms may be alternatively disturbed in different conditions. The gene regulatory events for two genes, DREB1A and DEAR1, in the cold-specific gene module had also been validated, as evidenced through the literature search.

CONCLUSIONS

Our protocols study the specificity of the gene modules that were specifically activated under a particular type of abiotic stress. The biplot can also assist to visualize the stress-specific gene modules. In conclusion, our approach has the potential to further elucidate mechanisms in plants and beneficial for future experiments design under different abiotic stresses.

Meta-analysis approach as a gene selection method in class prediction: does it improve model performance? A case study in acute myeloid leukemia

BACKGROUND

Aggregating gene expression data across experiments via meta-analysis is expected to increase the precision of the effect estimates and to increase the statistical power to detect a certain fold change. This study evaluates the potential benefit of using a meta-analysis approach as a gene selection method prior to predictive modeling in gene expression data.

RESULTS

Six raw datasets from different gene expression experiments in acute myeloid leukemia (AML) and 11 different classification methods were used to build classification models to classify samples as either AML or healthy control. First, the classification models were trained on gene expression data from single experiments using conventional supervised variable selection and externally validated with the other five gene expression datasets (referred to as the individual-classification approach). Next, gene selection was performed through meta-analysis on four datasets, and predictive models were trained with the selected genes on the fifth dataset and validated on the sixth dataset. For some datasets, gene selection through meta-analysis helped classification models to achieve higher performance as compared to predictive modeling based on a single dataset; but for others, there was no major improvement. Synthetic datasets were generated from nine simulation scenarios. The effect of sample size, fold change and pairwise correlation between differentially expressed (DE) genes on the difference between MA- and individual-classification model was evaluated. The fold change and pairwise correlation significantly contributed to the difference in performance between the two methods. The gene selection via meta-analysis approach was more effective when it was conducted using a set of data with low fold change and high pairwise correlation on the DE genes.

CONCLUSION

Gene selection through meta-analysis on previously published studies potentially improves the performance of a predictive model on a given gene expression data.

E2F transcription factors associated with up-regulated genes in glioblastoma

BACKGROUND

Glioblastoma is considered to the most common and malignant brain tumor in adults. Patients have a median survival of approximately one year from diagnosis due to poor response to therapy.

OBJECTIVE

We applied bioinformatics approaches to predict transcription factors (TF) that are deregulated in glioblastoma in an attempt to point out molecular targets for therapy.

METHODS

Up-regulated genes in glioblastoma selected from public microarray data were submitted to two TF association analyses. Thereafter, the expression levels of TF obtained in the overlap of analyses were assessed by RT-qPCR carried out in seven glioblastoma cell lines (T98, U251, U138, U87, U343, M059J, and M059K).

RESULTS

E2F1 and E2F4 were highlighted in both TF analyses. However, only E2F1 was confirmed as significantly up-regulated in all glioblastoma cell lines in vitro.

CONCLUSION

E2F1 is a potential common regulator of differentially expressed genes in glioblastoma, despite the genetic heterogeneity of tumor cells.

Searching for a Lifeline: Transcriptome Profiling Studies of Influenza Susceptibility and Resistance

Excess or dysregulated host inflammatory responses cause much of the morbidity and mortality caused by severe influenza. Given the limitations of vaccines and antiviral drugs, novel therapeutics to modulate host responses and improve outcomes in severe influenza are needed. One strategy is to learn from the direct comparison of high-survivor versus high-mortality animal models. This review surveys the results of lung transcriptome profiling studies in murine models that directly compare susceptible versus resistant hosts challenged with identical influenza infections. The potential contributions and limitations of these studies are discussed. To amplify their power, the studies are subjected to a meta-analysis, which helps identify frequently dysregulated pathways and potentially novel areas for investigation. Using connectivity map-based tools (LINCS), transcriptome signatures linked to susceptibility can identify candidate drugs that merit testing for in vivo efficacy.

Analyses of caspase-1-regulated transcriptomes in various tissues lead to identification of novel IL-1β-, IL-18- and sirtuin-1-independent pathways

BACKGROUND

It is well established that caspase-1 exerts its biological activities through its downstream targets such as IL-1β, IL-18, and Sirt-1. The microarray datasets derived from various caspase-1 knockout tissues indicated that caspase-1 can significantly impact the transcriptome. However, it is not known whether all the effects exerted by caspase-1 on transcriptome are mediated only by its well-known substrates. Therefore, we hypothesized that the effects of caspase-1 on transcriptome may be partially independent from IL-1β, IL-18, and Sirt-1.

METHODS

To determine new global and tissue-specific gene regulatory effects of caspase-1, we took novel microarray data analysis approaches including Venn analysis, cooperation analysis, and meta-analysis methods. We used these statistical methods to integrate different microarray datasets conducted on different caspase-1 knockout tissues and datasets where caspase-1 downstream targets were manipulated.

RESULTS

We made the following important findings: (1) Caspase-1 exerts its regulatory effects on the majority of genes in a tissue-specific manner; (2) Caspase-1 regulatory genes partially cooperates with genes regulated by sirtuin-1 during organ injury and inflammation in adipose tissue but not in the liver; (3) Caspase-1 cooperates with IL-1β in regulating less than half of the genes involved in cardiovascular disease, organismal injury, and cancer in mouse liver; (4) The meta-analysis identifies 40 caspase-1 globally regulated genes across tissues, suggesting that caspase-1 globally regulates many novel pathways; and (5) The meta-analysis identified new cooperatively and non-cooperatively regulated genes in caspase-1, IL-1β, IL-18, and Sirt-1 pathways.

CONCLUSIONS

Our findings suggest that caspase-1 regulates many new signaling pathways potentially via its known substrates and also via transcription factors and other proteins that are yet to be identified.

Identification of Polycystic Ovary Syndrome (PCOS) Specific Genes in Cumulus and Mural Granulosa Cells

Polycystic ovary syndrome (PCOS) is a metabolic and endocrine disorder which affects women of reproductive age with prevalence of 8–18%. The oocyte within the follicle is surrounded by cumulus cells (CCs), which connect with mural granulosa cells (MGCs) that are responsible for secreting steroid hormones. The main aim of this study is comparing gene expression profiles of MGCs and CCs in PCOS and control samples to identify PCOS-specific differentially expressed genes (DEGs). In this study, two microarray databases were searched for mRNA expression microarray studies performed with CCs and MGCs obtained from PCOS patients and control samples. Three independent studies were selected to be integrated with naive meta-analysis since raw meta-data from these studies were found to be highly correlated. DEGs in these somatic cells were identified for PCOS and control groups. This study enabled us to reveal dysregulation in MAPK (mitogen activated protein kinase), insulin and Wnt signaling pathways between CCs and MGCs in PCOS. The meta-analysis results together with qRT-PCR validations provide evidence that molecular signaling is dysregulated through MGCs and CCs in PCOS, which is important for follicle and oocyte maturation and may contribute to the pathogenesis of the syndrome.

Transcriptome of Cultured Lung Fibroblasts in Idiopathic Pulmonary Fibrosis: Meta-Analysis of Publically Available Microarray Datasets Reveals Repression of Inflammation and Immunity Pathways

Heritable profibrotic differentiation of lung fibroblasts is a key mechanism of idiopathic pulmonary fibrosis (IPF). Its mechanisms are yet to be fully understood. In this study, individual data from four independent microarray studies comparing the transcriptome of fibroblasts cultured in vitro from normal (total n = 20) and IPF (total n = 20) human lung were compiled for meta-analysis following normalization to z-scores. One hundred and thirteen transcripts were upregulated and 115 were downregulated in IPF fibroblasts using the Significance Analysis of Microrrays algorithm with a false discovery rate of 5%. Downregulated genes were highly enriched for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional classes related to inflammation and immunity such as Defense response to virus, Influenza A, tumor necrosis factor (TNF) mediated signaling pathway, interferon-inducible absent in melanoma2 (AIM2) inflammasome as well as Apoptosis. Although upregulated genes were not enriched for any functional class, select factors known to play key roles in lung fibrogenesis were overexpressed in IPF fibroblasts, most notably connective tissue growth factor (CTGF) and serum response factor (SRF), supporting their role as drivers of IPF. The full data table is available as a supplement.

Comprehensive Gene expression meta-analysis and integrated bioinformatic approaches reveal shared signatures between thrombosis and myeloproliferative disorders

Thrombosis is a leading cause of morbidity and mortality in patients with myeloproliferative disorders (MPDs), particularly polycythemia vera (PV) and essential thrombocythemia (ET). Despite the attempts to establish a link between them, the shared biological mechanisms are yet to be characterized. An integrated gene expression meta-analysis of five independent publicly available microarray data of the three diseases was conducted to identify shared gene expression signatures and overlapping biological processes. Using INMEX bioinformatic tool, based on combined Effect Size (ES) approaches, we identified a total of 1,157 differentially expressed genes (DEGs) (697 overexpressed and 460 underexpressed genes) shared between the three diseases. EnrichR tool's rich library was used for comprehensive functional enrichment and pathway analysis which revealed "mRNA Splicing" and "SUMO E3 ligases SUMOylate target proteins" among the most enriched terms. Network based meta-analysis identified MYC and FN1 to be the most highly ranked hub genes. Our results reveal that the alterations in biomarkers of the coagulation cascade like F2R, PROS1, SELPLG and ITGB2 were common between the three diseases. Interestingly, the study has generated a novel database of candidate genetic markers, pathways and transcription factors shared between thrombosis and MPDs, which might aid in the development of prognostic therapeutic biomarkers.

Methods to increase reproducibility in differential gene expression via meta-analysis

Findings from clinical and biological studies are often not reproducible when tested in independent cohorts. Due to the testing of a large number of hypotheses and relatively small sample sizes, results from whole-genome expression studies in particular are often not reproducible. Compared to single-study analysis, gene expression meta-analysis can improve reproducibility by integrating data from multiple studies. However, there are multiple choices in designing and carrying out a meta-analysis. Yet, clear guidelines on best practices are scarce. Here, we hypothesized that studying subsets of very large meta-analyses would allow for systematic identification of best practices to improve reproducibility. We therefore constructed three very large gene expression meta-analyses from clinical samples, and then examined meta-analyses of subsets of the datasets (all combinations of datasets with up to N/2 samples and K/2 datasets) compared to a ‘silver standard’ of differentially expressed genes found in the entire cohort. We tested three random-effects meta-analysis models using this procedure. We showed relatively greater reproducibility with more-stringent effect size thresholds with relaxed significance thresholds; relatively lower reproducibility when imposing extraneous constraints on residual heterogeneity; and an underestimation of actual false positive rate by Benjamini–Hochberg correction. In addition, multivariate regression showed that the accuracy of a meta-analysis increased significantly with more included datasets even when controlling for sample size.

Drug repositioning in SLE: crowd-sourcing, literature-mining and Big Data analysis

Lupus patients are in need of modern drugs to treat specific manifestations of their disease effectively and safely. In the past half century, only one new treatment has been approved by the US Food and Drug Administration (FDA) for systemic lupus erythematosus (SLE). In 2014–2015, the FDA approved 71 new drugs, only one of which targeted a rheumatic disease and none of which was approved for use in SLE. Repositioning/repurposing drugs approved for other diseases using multiple approaches is one possible means to find new treatment options for lupus patients. “Big Data” analysis approaches this challenge from an unbiased standpoint whereas literature mining and crowd sourcing for candidates assessed by the CoLTs (Combined Lupus Treatment Scoring) system provide a hypothesis-based approach to rank potential therapeutic candidates for possible clinical application. Both approaches mitigate risk since the candidates assessed have largely been extensively tested in clinical trials for other indications. The usefulness of a multi-pronged approach to drug repositioning in lupus is highlighted by orthogonal confirmation of hypothesis-based drug repositioning predictions by “Big Data” analysis of differentially expressed genes from lupus patient samples. The goal is to identify novel therapies that have the potential to affect disease processes specifically. Involvement of SLE patients and the scientists that study this disease in thinking about new drugs that may be effective in lupus though crowd-sourcing sites such as LRxL-STAT ( www.linkedin.com/in/lrxlstat ) is important in stimulating the momentum needed to test these novel drug targets for efficacy in lupus rapidly in small, proof-of-concept trials conducted by LuCIN, the Lupus Clinical Investigators Network ( www.linkedin.com/in/lucinstat ).

A comparative evaluation of data-merging and meta-analysis methods for reconstructing gene-gene interactions

BACKGROUND

We address the problem of integratively analyzing multiple gene expression, microarray datasets in order to reconstruct gene-gene interaction networks. Integrating multiple datasets is generally believed to provide increased statistical power and to lead to a better characterization of the system under study. However, the presence of systematic variation across different studies makes network reverse-engineering tasks particularly challenging. We contrast two approaches that have been frequently used in the literature for addressing systematic biases: meta-analysis methods, which first calculate opportune statistics on single datasets and successively summarize them, and data-merging methods, which directly analyze the pooled data after removing eventual biases. This comparative evaluation is performed on both synthetic and real data, the latter consisting of two manually curated microarray compendia comprising several E. coli and Yeast studies, respectively. Furthermore, the reconstruction of the regulatory network of the transcription factor Ikaros in human Peripheral Blood Mononuclear Cells (PBMCs) is presented as a case-study.

RESULTS

The meta-analysis and data-merging methods included in our experimentations provided comparable performances on both synthetic and real data. Furthermore, both approaches outperformed (a) the naïve solution of merging data together ignoring possible biases, and (b) the results that are expected when only one dataset out of the available ones is analyzed in isolation. Using correlation statistics proved to be more effective than using p-values for correctly ranking candidate interactions. The results from the PBMC case-study indicate that the findings of the present study generalize to different types of network reconstruction algorithms.

CONCLUSIONS

Ignoring the systematic variations that differentiate heterogeneous studies can produce results that are statistically indistinguishable from random guessing. Meta-analysis and data merging methods have proved equally effective in addressing this issue, and thus researchers may safely select the approach that best suit their specific application.

Gene expression analysis supports tumor threshold over 2.0 cm for T-category breast cancer

Tumor size, as indicated by the T-category, is known as a strong prognostic indicator for breast cancer. It is common practice to distinguish the T1 and T2 groups at a tumor size of 2.0 cm. We investigated the 2.0-cm rule from a new point of view. Here, we try to find the optimal threshold based on the differences between the gene expression profiles of the T1 and T2 groups (as defined by the threshold). We developed a numerical algorithm to measure the overall differential gene expression between patients with smaller tumors and those with larger tumors among multiple expression datasets from different studies. We confirmed the performance of the proposed algorithm by a simulation study and then applied it to three different studies conducted at two Norwegian hospitals. We found that the maximum difference in gene expression is obtained at a threshold of 2.2–2.4 cm, and we confirmed that the optimum threshold was over 2.0 cm, as indicated by a validation study using five publicly available expression datasets. Furthermore, we observed a significant differentiation between the two threshold groups in terms of time to local recurrence for the Norwegian datasets. In addition, we performed an associated network and canonical pathway analyses for the genes differentially expressed between tumors below and above the given thresholds, 2.0 and 2.4 cm, using the Norwegian datasets. The associated network function illustrated a cellular assembly of the genes for the 2.0-cm threshold: an energy production for the 2.4-cm threshold and an enrichment in lipid metabolism based on the genes in the intersection for the 2.0- and 2.4-cm thresholds.

Comprehensive identification of sexually dimorphic genes in diverse cattle tissues using RNA-seq

BACKGROUND

Molecular mechanisms associated with sexual dimorphism in cattle have not been well elucidated. Furthermore, as recent studies have implied that gene expression patterns are highly tissue specific, it is essential to investigate gene expression in a variety of tissues using RNA-seq. Here, we employed and compared two statistical methods, a simple two group test and Analysis of deviance (ANODEV), in order to investigate bovine sexually dimorphic genes in 40 RNA-seq samples distributed across two factors: sex and tissue.

RESULTS

As a result, we detected 752 sexually dimorphic genes across tissues from two statistical approaches and identified strong tissue-specific patterns of gene expression. Additionally, significantly detected sex-related genes shared between two mammal species (cattle and rat) were identified using qRT-PCR.

CONCLUSIONS

Results of our analyses reveal that sexual dimorphism of metabolic tissues and pituitary gland in cattle involves various biological processes. Several differentially expressed genes between sexes in cattle and rat species are shared, but show tissue-specific patterns. Finally, we concluded that two distinct statistical approaches have their advantages and disadvantages in RNA-seq studies investigating multiple tissues.

Analyzing Illumina Gene Expression Microarray Data Obtained From Human Whole Blood Cell and Blood Monocyte Samples

Microarray profiling of gene expression is widely applied to studies in molecular biology and functional genomics. Experimental and technical variations make not only the statistical analysis of single studies but also meta-analyses of different studies very challenging. Here, we describe the analytical steps required to substantially reduce the variations of gene expression data without affecting true effect sizes. A software pipeline has been established using gene expression data from a total of 3358 whole blood cell and blood monocyte samples, all from three German population-based cohorts, measured on the Illumina HumanHT-12 v3 BeadChip array. In summary, adjustment for a few selected technical factors greatly improved reliability of gene expression analyses. Such adjustments are particularly required for meta-analyses of different studies.

Type I interferon related genes are common genes on the early stage after vaccination by meta-analysis of microarray data

The objective of this study was to find common immune mechanism across different kinds of vaccines. A meta-analysis of microarray datasets was performed using publicly available microarray Gene Expression Omnibus (GEO) and Array Express data sets of vaccination records. Seven studies (out of 35) were selected for this meta-analysis. A total of 447 chips (145 pre-vaccination and 302 post-vaccination) were included. Significance analysis of microarrays (SAM) program was used for screening differentially expressed genes (DEGs). Functional pathway enrichment for the DEGs was conducted in DAVID Gene Ontology (GO) database. Twenty DEGs were identified, of which 10 up-regulated genes involved immune response. Six of which were type I interferon (IFN) related genes, including LY6E, MX1, OAS3, IFI44L, IFI6 and IFITM3. Ten down-regulated genes mainly mediated negative regulation of cell proliferation and cell motion. Results of a subgroup analysis showed that although the kinds of genes varied widely between days 3 and 7 post vaccination, the pathways between them are basically the same, such as immune response and response to viruses, etc. For an independent verification of these 6 type I IFN related genes, peripheral blood mononuclear cells (PBMCs) were collected at baseline and day 3 after the vaccination from 8 Enterovirus 71(EV71) vaccinees and were assayed by RT-PCR. Results showed that the 6 DEGs were also upregulated in EV71 vaccinees. In summary, meta-analysis methods were used to explore the immune mechanism of vaccines and results indicated that the type I IFN related genes and corresponding pathways were common in early immune responses for different kinds of vaccines.

Microarray Meta-Analysis and Cross-Platform Normalization: Integrative Genomics for Robust Biomarker Discovery

The diagnostic and prognostic potential of the vast quantity of publicly-available microarray data has driven the development of methods for integrating the data from different microarray platforms. Cross-platform integration, when appropriately implemented, has been shown to improve reproducibility and robustness of gene signature biomarkers. Microarray platform integration can be conceptually divided into approaches that perform early stage integration (cross-platform normalization) versus late stage data integration (meta-analysis). A growing number of statistical methods and associated software for platform integration are available to the user, however an understanding of their comparative performance and potential pitfalls is critical for best implementation. In this review we provide evidence-based, practical guidance to researchers performing cross-platform integration, particularly with an objective to discover biomarkers.

A Simple Rank Product Approach for Analyzing Two Classes

The rank product statistic has been widely used to detect differentially expressed genes in replicated microarrays and a one-class setting. The objective of this article is to apply a rank product statistic to approximate the P-value of differential expression in a two-class setting, such as in normal and cancer cells. For this purpose, we introduce a simple statistic that compares the P-values of each class’s rank product statistic. Its null distribution is straightforwardly derived using the change-of-variable technique.

A ranking-based meta-analysis reveals let-7 family as a meta-signature for grade classification in breast cancer

Breast cancer is one of the most important causes of cancer-related deaths worldwide in women. In addition to gene expression studies, the progressing work in the miRNA area including miRNA microarray studies, brings new aspects to the research on the cancer development and progression. Microarray technology has been widely used to find new biomarkers in research and many transcriptomic microarray studies are available in public databases. In this study, the breast cancer miRNA and mRNA microarray studies were collected according to the availability of their data and clinical information, and combined by a newly developed ranking-based meta-analysis approach to find out candidate miRNA biomarkers (meta-miRNAs) that classify breast cancers according to their grades and explain the relation between miRNAs and mRNAs. This approach provided meta-miRNAs specific to breast cancer grades, pointing out let-7 family members as grade classifiers. The qRT-PCR studies performed with independent breast tumors confirmed the potential biomarker role of let-7 family members (meta-miRNAs). The concordance between the meta-mRNAs and miRNA target genes specific to tumor grade (common genes) supported the idea of mRNAs as miRNA targets. The pathway analysis results showed that most of the let-7 family miRNA targets, and also common genes, were significantly taking part in cancer-related pathways. The qRT-PCR studies, together with bioinformatic analyses, confirmed the results of meta-analysis approach, which is dynamic and allows combining datasets from different platforms.

Identification of suitable genes contributes to lung adenocarcinoma clustering by multiple meta-analysis methods

BACKGROUND

With the widespread application of high-throughput technology, numerous meta-analysis methods have been proposed for differential expression profiling across multiple studies.

OBJECTIVES

We identified the suitable differentially expressed (DE) genes that contributed to lung adenocarcinoma (ADC) clustering based on seven popular multiple meta-analysis methods.

METHODS

Seven microarray expression profiles of ADC and normal controls were extracted from the ArrayExpress database. The Bioconductor was used to perform the data preliminary preprocessing. Then, DE genes across multiple studies were identified. Hierarchical clustering was applied to compare the classification performance for microarray data samples. The classification efficiency was compared based on accuracy, sensitivity and specificity.

RESULTS

Across seven datasets, 573 ADC cases and 222 normal controls were collected. After filtering out unexpressed and noninformative genes, 3688 genes were remained for further analysis. The classification efficiency analysis showed that DE genes identified by sum of ranks method separated ADC from normal controls with the best accuracy, sensitivity and specificity of 0.953, 0.969 and 0.932, respectively. The gene set with the highest classification accuracy mainly participated in the regulation of response to external stimulus (P = 7.97E-04), cyclic nucleotide-mediated signaling (P = 0.01), regulation of cell morphogenesis (P = 0.01) and regulation of cell proliferation (P = 0.01).

CONCLUSIONS

Evaluation of DE genes identified by different meta-analysis methods in classification efficiency provided a new perspective to the choice of the suitable method in a given application. Varying meta-analysis methods always present varying abilities, so synthetic consideration should be taken when providing meta-analysis methods for particular research.

A common gene signature across multiple studies relate biomarkers and functional regulation in tolerance to renal allograft

Patients tolerant to a kidney graft display a specific blood cell transcriptional pattern but results from five different studies were inconsistent, raising the question of relevance for future clinical application. To resolve this, we sought to identify a common gene signature, specific functional and cellular components, and discriminating biomarkers for tolerance following kidney transplantation. A meta-analysis of studies identified a robust gene signature involving proliferation of B and CD4 T cells, and inhibition of CD14 monocyte related functions among 96 tolerant samples. This signature was further supported through a cross-validation approach, yielding 92.5% accuracy independent of the study of origin. Experimental validation, performed on new tolerant samples and using a selection of the top-20 biomarkers, returned 91.7% of good classification. Beyond the confirmation of B-cell involvement, our data also indicated participation of other cell subsets in tolerance. Thus, the use of the top 20 biomarkers, mostly centered on B cells, may provide a common and standardized tool towards personalized medicine for the monitoring of tolerant or low-risk patients among kidney allotransplant recipients. These data point to a global preservation of genes favoring the maintenance of a homeostatic and ‘healthy' environment in tolerant patients and may contribute to a better understanding of tolerance maintenance mechanisms.

Integrative correlation: Properties and relation to canonical correlations

The integrative correlation coefficient was developed to facilitate the validation of expression microarray results in public datasets, by identifying genes that are reproducibly measured across studies and even across microarray platforms. In the current study, we develop a number of interesting and important mathematical and statistical properties of the integrative correlation coefficient, including a unique permutation-based null distribution with the unusual property that the variance does not shrink as the sample size increases, discussing how these findings impact its use and interpretation, and what they have to say about any method for identifying reproducible genes in a meta-analysis.

Meta-analysis methods for combining multiple expression profiles: comparisons, statistical characterization and an application guideline

BACKGROUND

As high-throughput genomic technologies become accurate and affordable, an increasing number of data sets have been accumulated in the public domain and genomic information integration and meta-analysis have become routine in biomedical research. In this paper, we focus on microarray meta-analysis, where multiple microarray studies with relevant biological hypotheses are combined in order to improve candidate marker detection. Many methods have been developed and applied in the literature, but their performance and properties have only been minimally investigated. There is currently no clear conclusion or guideline as to the proper choice of a meta-analysis method given an application; the decision essentially requires both statistical and biological considerations.

RESULTS

We performed 12 microarray meta-analysis methods for combining multiple simulated expression profiles, and such methods can be categorized for different hypothesis setting purposes: (1) HS(A): DE genes with non-zero effect sizes in all studies, (2) HS(B): DE genes with non-zero effect sizes in one or more studies and (3) HS(r): DE gene with non-zero effect in "majority" of studies. We then performed a comprehensive comparative analysis through six large-scale real applications using four quantitative statistical evaluation criteria: detection capability, biological association, stability and robustness. We elucidated hypothesis settings behind the methods and further apply multi-dimensional scaling (MDS) and an entropy measure to characterize the meta-analysis methods and data structure, respectively.

CONCLUSIONS

The aggregated results from the simulation study categorized the 12 methods into three hypothesis settings (HS(A), HS(B), and HS(r)). Evaluation in real data and results from MDS and entropy analyses provided an insightful and practical guideline to the choice of the most suitable method in a given application. All source files for simulation and real data are available on the author's publication website.