1
|
Arencibia A, Salazar LA. Microarray meta-analysis reveals IL6 and p38β/MAPK11 as potential targets of hsa-miR-124 in endothelial progenitor cells: Implications for stent re-endothelization in diabetic patients. Front Cardiovasc Med 2022; 9:964721. [PMID: 36176980 PMCID: PMC9513120 DOI: 10.3389/fcvm.2022.964721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Accepted: 08/25/2022] [Indexed: 11/13/2022] Open
Abstract
Circulating endothelial progenitor cells (EPCs) play an important role in the repair processes of damaged vessels, favoring re-endothelization of stented vessels to minimize restenosis. EPCs number and function is diminished in patients with type 2 diabetes, a known risk factor for restenosis. Considering the impact of EPCs in vascular injury repair, we conducted a meta-analysis of microarray to assess the transcriptomic profile and determine target genes during the differentiation process of EPCs into mature ECs. Five microarray datasets, including 13 EPC and 12 EC samples were analyzed, using the online tool ExpressAnalyst. Differentially expressed genes (DEGs) analysis was done by Limma method, with an | log2FC| > 1 and FDR < 0.05. Combined p-value by Fisher exact method was computed for the intersection of datasets. There were 3,267 DEGs, 1,539 up-regulated and 1,728 down-regulated in EPCs, with 407 common DEGs in at least four datasets. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed enrichment for terms related to “AGE-RAGE signaling pathway in diabetic complications.” Intersection of common DEGs, KEGG pathways genes and genes in protein-protein interaction network (PPI) identified four key genes, two up-regulated (IL1B and STAT5A) and two down-regulated (IL6 and MAPK11). MicroRNA enrichment analysis of common DEGs depicted five hub microRNA targeting 175 DEGs, including STAT5A, IL6 and MAPK11, with hsa-miR-124 as common regulator. This group of genes and microRNAs could serve as biomarkers of EPCs differentiation during coronary stenting as well as potential therapeutic targets to improve stent re-endothelization, especially in diabetic patients.
Collapse
|
2
|
Prasad B, Li X. Fused inverse-normal method for integrated differential expression analysis of RNA-seq data. BMC Bioinformatics 2022; 23:320. [PMID: 35931958 PMCID: PMC9354357 DOI: 10.1186/s12859-022-04859-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Accepted: 07/19/2022] [Indexed: 11/10/2022] Open
Abstract
Background Use of next-generation sequencing technologies to transcriptomics (RNA-seq) for gene expression profiling has found widespread application in studying different biological conditions including cancers. However, RNA-seq experiments are still small sample size experiments due to the cost. Recently, an increased focus has been on meta-analysis methods for integrated differential expression analysis for exploration of potential biomarkers. In this study, we propose a p-value combination method for meta-analysis of multiple independent but related RNA-seq studies that accounts for sample size of a study and direction of expression of genes in individual studies. Results The proposed method generalizes the inverse-normal method without an increase in statistical or computational complexity and does not pre- or post-hoc filter genes that have conflicting direction of expression in different studies. Thus, the proposed method, as compared to the inverse-normal, has better potential for the discovery of differentially expressed genes (DEGs) with potentially conflicting differential signals from multiple studies related to disease. We demonstrated the use of the proposed method in detection of biologically relevant DEGs in glioblastoma (GBM), the most aggressive brain cancer. Our approach notably enabled the identification of over-expressed tumour suppressor gene RAD51 in GBM compared to healthy controls, which has recently been shown to be a target for inhibition to enhance radiosensitivity of GBM cells during treatment. Pathway analysis identified multiple aberrant GBM related pathways as well as novel regulators such as TCF7L2 and MAPT as important upstream regulators in GBM. Conclusions The proposed meta-analysis method generalizes the existing inverse-normal method by providing a way to establish differential expression status for genes with conflicting direction of expression in individual RNA-seq studies. Hence, leading to further exploration of them as potential biomarkers for the disease. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-022-04859-9.
Collapse
Affiliation(s)
- Birbal Prasad
- National Horizons Centre, School of Health and Life Sciences, Teesside University, Darlington, DL1 1HG, UK
| | - Xinzhong Li
- National Horizons Centre, School of Health and Life Sciences, Teesside University, Darlington, DL1 1HG, UK.
| |
Collapse
|
3
|
Ye Z, Mo C, Ke H, Yan Q, Chen C, Kochunov P, Hong LE, Mitchell BD, Chen S, Ma T. Meta-Analysis of Transcriptome-Wide Association Studies across 13 Brain Tissues Identified Novel Clusters of Genes Associated with Nicotine Addiction. Genes (Basel) 2021; 13:37. [PMID: 35052378 PMCID: PMC8775257 DOI: 10.3390/genes13010037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2021] [Revised: 12/14/2021] [Accepted: 12/18/2021] [Indexed: 12/01/2022] Open
Abstract
Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci, but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as a reference panel to detect associations with the phenotype at the gene level and have been gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as a reference panel to conduct tissue-specific TWAS on cigarettes per day (CPD) over thirteen brain tissues in two large cohorts: UK Biobank (UKBB; number of participants (N) = 142,202) and the GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N = 143,210), then meta-analyzing the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues; partially homogeneous genes associated with CPD in cortex, cerebellum, and hippocampus tissues; and, lastly, the tissue-specific genes associated with CPD in only a few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.
Collapse
Affiliation(s)
- Zhenyao Ye
- Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, MD 21201, USA; (Z.Y.); (C.M.); (P.K.); (L.E.H.)
- Division of Biostatistics and Bioinformatics, Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, MD 21201, USA;
| | - Chen Mo
- Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, MD 21201, USA; (Z.Y.); (C.M.); (P.K.); (L.E.H.)
- Division of Biostatistics and Bioinformatics, Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, MD 21201, USA;
| | - Hongjie Ke
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, MD 20742, USA;
| | - Qi Yan
- Irving Medical Center, Department of Obstetrics & Gynecology, Columbia University, New York, NY 10032, USA;
| | - Chixiang Chen
- Division of Biostatistics and Bioinformatics, Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, MD 21201, USA;
| | - Peter Kochunov
- Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, MD 21201, USA; (Z.Y.); (C.M.); (P.K.); (L.E.H.)
| | - L. Elliot Hong
- Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, MD 21201, USA; (Z.Y.); (C.M.); (P.K.); (L.E.H.)
| | - Braxton D. Mitchell
- Department of Medicine, University of Maryland School of Medicine, Baltimore, MD 21201, USA;
| | - Shuo Chen
- Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, MD 21201, USA; (Z.Y.); (C.M.); (P.K.); (L.E.H.)
- Division of Biostatistics and Bioinformatics, Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, MD 21201, USA;
| | - Tianzhou Ma
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, MD 20742, USA;
| |
Collapse
|
4
|
Ye Z, Ke H, Chen S, Cruz-Cano R, He X, Zhang J, Dorgan J, Milton DK, Ma T. Biomarker Categorization in Transcriptomic Meta-Analysis by Concordant Patterns With Application to Pan-Cancer Studies. Front Genet 2021; 12:651546. [PMID: 34276766 PMCID: PMC8283696 DOI: 10.3389/fgene.2021.651546] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2021] [Accepted: 05/28/2021] [Indexed: 01/21/2023] Open
Abstract
With the increasing availability and dropping cost of high-throughput technology in recent years, many-omics datasets have accumulated in the public domain. Combining multiple transcriptomic studies on related hypothesis via meta-analysis can improve statistical power and reproducibility over single studies. For differential expression (DE) analysis, biomarker categorization by DE pattern across studies is a natural but critical task following biomarker detection to help explain between study heterogeneity and classify biomarkers into categories with potentially related functionality. In this paper, we propose a novel meta-analysis method to categorize biomarkers by simultaneously considering the concordant pattern and the biological and statistical significance across studies. Biomarkers with the same DE pattern can be analyzed together in downstream pathway enrichment analysis. In the presence of different types of transcripts (e.g., mRNA, miRNA, and lncRNA, etc.), integrative analysis including miRNA/lncRNA target enrichment analysis and miRNA-mRNA and lncRNA-mRNA causal regulatory network analysis can be conducted jointly on all the transcripts of the same category. We applied our method to two Pan-cancer transcriptomic study examples with single or multiple types of transcripts available. Targeted downstream analysis identified categories of biomarkers with unique functionality and regulatory relationships that motivate new hypothesis in Pan-cancer analysis.
Collapse
Affiliation(s)
- Zhenyao Ye
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| | - Hongjie Ke
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| | - Shuo Chen
- Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, Baltimore, MD, United States
| | - Raul Cruz-Cano
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| | - Xin He
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| | - Jing Zhang
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| | - Joanne Dorgan
- Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, Baltimore, MD, United States
| | - Donald K Milton
- Maryland Institute for Applied Environmental Health, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| | - Tianzhou Ma
- Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, College Park, MD, United States
| |
Collapse
|
5
|
Ma T, Liang F, Oesterreich S, Tseng GC. A Joint Bayesian Model for Integrating Microarray and RNA Sequencing Transcriptomic Data. J Comput Biol 2017; 24:647-662. [PMID: 28541721 PMCID: PMC5510692 DOI: 10.1089/cmb.2017.0056] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
As the sequencing cost continued to drop in the past decade, RNA sequencing (RNA-seq) has replaced microarray to become the standard high-throughput experimental tool to analyze transcriptomic profile. As more and more datasets are generated and accumulated in the public domain, meta-analysis to combine multiple transcriptomic studies to increase statistical power has received increasing popularity. In this article, we propose a Bayesian hierarchical model to jointly integrate microarray and RNA-seq studies. Since systematic fold change differences across RNA-seq and microarray for detecting differentially expressed genes have been previously reported, we replicated this finding in several real datasets and showed that incorporation of a normalization procedure to account for the bias improves the detection accuracy and power. We compared our method with the popular two-stage Fisher's method using simulations and two real applications in a histological subtype (invasive lobular carcinoma) of breast cancer comparing PR+ versus PR- and early-stage versus late-stage patients. The result showed improved detection power and more significant and interpretable pathways enriched in the detected biomarkers from the proposed Bayesian model.
Collapse
Affiliation(s)
- Tianzhou Ma
- Department of Biostatistics, University of Pittsburgh, Pittsburgh, Pennsylvania
| | - Faming Liang
- Department of Biostatistics, University of Florida, Gainesville, Florida
| | - Steffi Oesterreich
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, Pennsylvania
- Women's Cancer Research Center, Pittsburgh, Pennsylvania
| | - George C. Tseng
- Department of Biostatistics, University of Pittsburgh, Pittsburgh, Pennsylvania
- Department of Human Genetics, University of Pittsburgh, Pittsburgh, Pennsylvania
- Department of Computational Biology, University of Pittsburgh, Pittsburgh, Pennsylvania
| |
Collapse
|