scholarly journals Genome-Wide Identification of Potential Biomarkers in Multiple Myeloma Using Meta-Analysis of mRNA and miRNA Expression Data

Author(s):  
Amit Katiyar ◽  
Gurvinder Kaur ◽  
Lata Rani ◽  
Harpreet Singh ◽  
Punit Kaur ◽  
...  

Abstract Multiple myeloma (MM) is a plasma cell malignancy with diverse clinical phenotypes and molecular heterogeneity that is not completely understood. Recent studies have identified differentially expressed patterns of genes (DEGs) or miRNAs (DEMs) in MM. But these signatures overlap partially, plausibly due to complexity of myeloma genome, diversity in cell lines studied, molecular technologies and analytical tools utilized in these studies. This warrants further investigations since DEGs/DEMs can impact clinical outcomes and guide personalized therapy. We conducted the genome-wide meta-analysis of expression datasets on DEGs/DEMs in MM and derive net putative signatures and potential biomarkers for MM. A set of 110 DEMs and 3,817 DEGs were identified to be differentially expressed. Among these, 86 DEMs (60 downregulated; 26 upregulated) correlated with 1,970 target DEGs (1373 downregulated; 597 upregulated). Signatures of 23 DEMs (‘Union 23’) and 196 DEGs (‘Union 196’) were deduced that shared 10 DEMs and 13 DEGs with published signatures, respectively. The study has identified five topmost nodal genes (APP, KIAA0101, CDK2, ESR1 and FN1) derived from functional modules in PPI networks and has paved the way for further studies to establish their prognostic potential and role in therapeutics for MM. The integrated bioinformatics methods and expression profiling techniques may lead to the identification of putative hub genes and expression signatures that can serve as predictive biomarkers of MM progression.

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Amit Katiyar ◽  
Gurvinder Kaur ◽  
Lata Rani ◽  
Lingaraja Jena ◽  
Harpreet Singh ◽  
...  

AbstractMultiple myeloma (MM) is a plasma cell malignancy with diverse clinical phenotypes and molecular heterogeneity not completely understood. Differentially expressed genes (DEGs) and miRNAs (DEMs) in MM may influence disease pathogenesis, clinical presentation / drug sensitivities. But these signatures overlap meagrely plausibly due to complexity of myeloma genome, diversity in primary cells studied, molecular technologies/ analytical tools utilized. This warrants further investigations since DEGs/DEMs can impact clinical outcomes and guide personalized therapy. We have conducted genome-wide meta-analysis of DEGs/DEMs in MM versus Normal Plasma Cells (NPCs) and derived unified putative signatures for MM. 100 DEMs and 1,362 DEGs were found deranged between MM and NPCs. Signatures of 37 DEMs (‘Union 37’) and 154 DEGs (‘Union 154’) were deduced that shared 17 DEMs and 22 DEGs with published prognostic signatures, respectively. Two miRs (miR-16–2-3p, 30d-2-3p) correlated with survival outcomes. PPI analysis identified 5 topmost functionally connected hub genes (UBC, ITGA4, HSP90AB1, VCAM1, VCP). Transcription factor regulatory networks were determined for five seed DEGs with ≥ 4 biomarker applications (CDKN1A, CDKN2A, MMP9, IGF1, MKI67) and three topmost up/ down regulated DEMs (miR-23b, 195, let7b/ miR-20a, 155, 92a). Further studies are warranted to establish and translate prognostic potential of these signatures for MM.


2019 ◽  
Author(s):  
Daniel Palmer ◽  
Fabio Fabris ◽  
Aoife Doherty ◽  
Alex A. Freitas ◽  
João Pedro de Magalhães

1AbstractUnderstanding the expression changes that come with age is an important step in understanding the ageing process as a whole. By combining such transcriptomic data with other sources of information, for instance protein-protein interaction (PPI) data, it is possible to make inferences about the functional changes that occur with age. To address this, we conducted a meta-analysis on 127 publicly available microarray and RNA-Seq datasets from mice, rats and humans, to identify genes that are commonly differentially expressed with age in mammals. We also conducted analyses on subsets of these datasets, to produce transcriptomic signatures for brain, heart and muscle tissues, all of which are important tissues in the pathophysiology of ageing. This approach identified the transcriptomic signatures of the ageing system, as well as brain, heart and muscle tissues. We then applied enrichment analysis and machine learning to functionally describe those signatures. This revealed a typical ageing signature including the overexpression of immune and stress response genes and the underexpression of metabolic and developmental genes. Further analysis of the ageing expression signatures revealed that genes differentially expressed with age tend to be broadly expressed across tissues, rather than be tissue-specific, and that the ageing expression signatures (particularly the overexpressed signatures) of the whole system, brain and muscle tend to include genes that are central in PPI networks. We also show that genes underexpressed in the brain are highly central in a co-expression map, suggesting that underexpression of these genes may play a part in cognitive ageing. In sum, we show numerous functional similarities between the ageing transcriptomes of these important tissues, a broad non-specific expression pattern in genes differentially expressed with age, along with altered network properties of these genes in both a PPI and co-expression network.


2019 ◽  
Vol 10 (1) ◽  
Author(s):  
Charles E. Mordaunt ◽  
Bo Y. Park ◽  
Kelly M. Bakulski ◽  
Jason I. Feinberg ◽  
Lisa A. Croen ◽  
...  

Abstract Background Autism spectrum disorder (ASD) is a neurodevelopmental disorder that affects more than 1% of children in the USA. ASD risk is thought to arise from both genetic and environmental factors, with the perinatal period as a critical window. Understanding early transcriptional changes in ASD would assist in clarifying disease pathogenesis and identifying biomarkers. However, little is known about umbilical cord blood gene expression profiles in babies later diagnosed with ASD compared to non-typically developing and non-ASD (Non-TD) or typically developing (TD) children. Methods Genome-wide transcript levels were measured by Affymetrix Human Gene 2.0 array in RNA from cord blood samples from both the Markers of Autism Risk in Babies-Learning Early Signs (MARBLES) and the Early Autism Risk Longitudinal Investigation (EARLI) high-risk pregnancy cohorts that enroll younger siblings of a child previously diagnosed with ASD. Younger siblings were diagnosed based on assessments at 36 months, and 59 ASD, 92 Non-TD, and 120 TD subjects were included. Using both differential expression analysis and weighted gene correlation network analysis, gene expression between ASD and TD, and between Non-TD and TD, was compared within each study and via meta-analysis. Results While cord blood gene expression differences comparing either ASD or Non-TD to TD did not reach genome-wide significance, 172 genes were nominally differentially expressed between ASD and TD cord blood (log2(fold change) > 0.1, p < 0.01). These genes were significantly enriched for functions in xenobiotic metabolism, chromatin regulation, and systemic lupus erythematosus (FDR q < 0.05). In contrast, 66 genes were nominally differentially expressed between Non-TD and TD, including 8 genes that were also differentially expressed in ASD. Gene coexpression modules were significantly correlated with demographic factors and cell type proportions. Limitations ASD-associated gene expression differences identified in this study are subtle, as cord blood is not the main affected tissue, it is composed of many cell types, and ASD is a heterogeneous disorder. Conclusions This is the first study to identify gene expression differences in cord blood specific to ASD through a meta-analysis across two prospective pregnancy cohorts. The enriched gene pathways support involvement of environmental, immune, and epigenetic mechanisms in ASD etiology.


Blood ◽  
2019 ◽  
Vol 134 (Supplement_1) ◽  
pp. 361-361 ◽  
Author(s):  
Anil Aktas Samur ◽  
Mehmet Kemal Samur ◽  
Michael A Lopez ◽  
Sanika Derebail ◽  
Kenneth C. Anderson ◽  
...  

Alternative splicing (AS) is a critical post-transcriptional event, which affects the number of cellular processes. Aberrant splicing of some genes has been reported in multiple myeloma (MM). However, to date, whole-transcriptome-wide AS study has not been performed. We used deep RNA-sequencing data from 16 normal plasma cells (NPC) and 360 newly-diagnosed MM patients to describe the landscape of the alternative splicing events and the molecular mechanisms driving aberrant AS in MM. Global splicing analysis showed that mutually exclusive exon (MXE) (n=510) and Skipped Exon (SE) (n=417) are the most frequent splicing events in MM compared to NPC. Among these events, ~54% were observed in genes which are not differentially expressed between MM and NPC and 46% of the AS events (SE, MXE, retained intron, alternative 3'/5' sites) were observed in differentially expressed genes targeting 203 unique genes. AS affected RNA transcription regulation genes such as IKZF1, IKZF3, and key regulatory elements in MM including, IRF3, IRF4, or key transcription factors such as MEF2C, XPB1, STAT2, and ILF3. In general, AS targetted DNA replication, cell cycle, and apoptosis pathways. MM subgroups showed a heterogeneity for AS events. Monosomy 14, t(4;14), del1p and del17p had the highest number of unique (not observed in other subgroups) AS events compared to NPC.To understanding the molecular mechanisms driving aberrant alternate splicing we next investigated115 splicing factors (SF) in MM and associated them with AS events. We observed that ~40% of SF were dysregulated (dysregulated expression and/or copy number alteration) in MM compared to NPC, including SRSF, PCBP and RBM families. To understand the key binding regions, we have performed SF binding motif enrichment analysis around AS events and found that SRFS1, SRSF9, and PCB1 motifs to be enriched among the splicing events. Importantly, SRSF1 expression was linked with survival in two independent MM datasets.We therefore explored functional role of SRSF1 in MM with perturbation studies. While upregulation of SRSF1 expression significantly increased the cell growth and survival, conversely downregulation of SRSF1 inhibited the both. To dissect the mechanisms of SRSF1-mediated MM growth induction, we utilized SRSF1 mutants lacking either of the 2 RNA-recognition motifs or the serine/argine-rich C-terminal domain involved in protein-protein interactions, and recruitment of spliceosome components. We also utilized a C-terminal fusion of SRSF1 with the nuclear-retention signal of SRSF2 (NRS1 mutant), to force SRSF1 retention in the nucleus and assess the role of its nuclear versus cytoplasmic functions. These studies suggested that SRSF1-regulated AS effects MM cell proliferation. We surprisingly also found that even NRS1 mutant failed to promote MM growth, suggesting an important role of cytoplasmic SRSF1 in promoting MM cells proliferation.We next investigated alternative splicing pattern changes induced by SRSF1 knock down.When analyzing cellular functions of SRSF1-regulated splicing events, we found that SRSF1 knock down affect's genes in the RNA processing pathway as well as genes involved in cancer-related functions such as mTOR, E2F and MYC-related pathways. Splicing analysis was corroborated with immunoprecipitation (IP) followed by mass spectrometry (MS) analysis of T7-tagged SRSF1 MM cells.Finally, using genome wide chromatin and transcription landscape mapping techniques, we have found SRSF1 to be under the transcriptional control of oncogenic E2F1 in MM cells. Consistent with these findings, we observed greater in vitro loss of viability in a large panel of MM cell lines compared with PBMCs from healthy volunteers, following exposure to the splicing modulator pladeniolide. In summary, this study for the first time reports a detailed splicing landscape in myeloma and highlights the biological and clinical importance of alternative splicing events. Moreover, these results indicate a functional role and clinical significance of a gene involved in regulation of alternate splicing in MM, highlighting the need to further understand the splicing pattern in myeloma initiation and progression. Disclosures Anderson: Takeda: Consultancy, Speakers Bureau; Celgene: Consultancy, Speakers Bureau; Amgen: Consultancy, Speakers Bureau; Janssen: Consultancy, Speakers Bureau; Oncopep: Other: Scientific Founder; Sanofi-Aventis: Other: Advisory Board; Bristol-Myers Squibb: Other: Scientific Founder. Avet-Loiseau:takeda: Consultancy, Other: travel fees, lecture fees, Research Funding; celgene: Consultancy, Other: travel fees, lecture fees, Research Funding. Munshi:Adaptive: Consultancy; Abbvie: Consultancy; Takeda: Consultancy; Janssen: Consultancy; Oncopep: Consultancy; Amgen: Consultancy; Celgene: Consultancy.


2018 ◽  
Author(s):  
Charles E. Mordaunt ◽  
Bo Y. Park ◽  
Kelly M. Bakulski ◽  
Jason I. Feinberg ◽  
Lisa A. Croen ◽  
...  

AbstractBackgroundAutism spectrum disorder (ASD) is a neurodevelopmental disorder that affects more than 1% of children in the United States. ASD risk is thought to arise from a combination of genetic and environmental factors, with the perinatal period as a critical window. Understanding early transcriptional changes in ASD would assist in clarifying disease pathogenesis and identifying biomarkers and treatments. However, little is known about umbilical cord blood gene expression profiles in babies later diagnosed with ASD compared to non-typically developing (Non-TD) or neurotypical children.MethodsGenome-wide transcript levels were measured by Affymetrix Human Gene 2.0 array in RNA from umbilical cord blood samples from both the Markers of Autism Risk in Babies--Learning Early Signs (MARBLES) and the Early Autism Risk Longitudinal Investigation (EARLI) high-risk pregnancy cohorts that enroll younger siblings of a child previously diagnosed with ASD. An algorithm-based diagnosis from 36 month assessments categorized the younger sibling as either ASD, typically developing (TD), or not ASD but non-typically developing (Non-TD). 59 ASD, 92 Non-TD, and 120 TD subjects were included and differences were identified in ASD versus TD subjects, with Non-TD versus TD as a specificity control. Meta-analysis was used to combine the results from both studies. Functional enrichments of differentially-expressed genes were examined across diagnostic groups.ResultsWhile cord blood gene expression differences comparing either ASD or Non-TD to TD did not reach genome-wide significance when adjusting for multiple comparisons, 172 genes were nominally differentially-expressed between ASD and TD cord blood (log2(fold change) > 0.1, p < 0.01). These genes were significantly enriched for toxic substance response and xenobiotic metabolism functions, and gene sets involved in chromatin regulation and systemic lupus erythematosus were significantly upregulated (FDR q < 0.05). In contrast, 66 genes were differentially-expressed between Non-TD and TD cord blood, including only 8 genes that were also differentially-expressed in ASD.ConclusionsThis is the first study to identify perinatal gene expression differences in umbilical cord blood specific to ASD. The results of this meta-analysis across two prospective ASD cohorts support involvement of environmental, immune, and epigenetic mechanisms in ASD etiology.


2020 ◽  
Author(s):  
Mariana Ferrarini ◽  
Avantika Lal ◽  
Rita Rebollo ◽  
Andreas Gruber ◽  
Andrea Guarracino ◽  
...  

Abstract The novel betacoronavirus named Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) caused a worldwide pandemic (COVID-19) after initially emerging in Wuhan, China. Here we applied a novel, comprehensive bioinformatic strategy to public RNA sequencing and viral genome sequencing data, to better understand how SARS-CoV-2 interacts with human cells. To our knowledge, this is the first meta-analysis to predict host factors that play a specific role in SARS-CoV-2 pathogenesis, distinct from other respiratory viruses. We identified differentially expressed genes, isoforms and transposable element families specifically altered in SARS-CoV-2 infected cells. Well-known immunoregulators including CSF2, IL-32, IL-6 and SERPINA3 were differentially expressed, while immunoregulatory transposable element families were overexpressed. We predicted conserved interactions between the SARS-CoV-2 genome and human RNA-binding proteins such as hnRNPA1, PABPC1 and eIF4b, which may play important roles in the viral life cycle. We also detected four viral sequence variants in the spike, polymerase, and nonstructural proteins that correlate with severity of COVID-19. The host factors we identified likely represent important mechanisms in the disease profile of this pathogen, and could be targeted by prophylactics and/or therapeutics against SARS-CoV-2.


2020 ◽  
Vol 4 (1) ◽  
pp. 181-190 ◽  
Author(s):  
Zhaohui Du ◽  
Niels Weinhold ◽  
Gregory Chi Song ◽  
Kristin A. Rand ◽  
David J. Van Den Berg ◽  
...  

Abstract Persons of African ancestry (AA) have a twofold higher risk for multiple myeloma (MM) compared with persons of European ancestry (EA). Genome-wide association studies (GWASs) support a genetic contribution to MM etiology in individuals of EA. Little is known about genetic risk factors for MM in individuals of AA. We performed a meta-analysis of 2 GWASs of MM in 1813 cases and 8871 controls and conducted an admixture mapping scan to identify risk alleles. We fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. In GWAS meta-analysis, we identified 2 suggestive novel loci located at 9p24.3 and 9p13.1 at P &lt; 1 × 10−6; however, no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. Of the 23 known EA risk variants, 20 showed directional consistency, and 9 replicated at P &lt; .05 in AA individuals. In 8 regions, we identified markers that better capture MM risk in persons with AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95% confidence interval, 1.56-2.11) increased MM risk compared with those with average risk (25%-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P = 5.1 × 10−12). Our study shows that common genetic variation contributes to MM risk in individuals with AA.


Sign in / Sign up

Export Citation Format

Share Document