Next generation sequencing (NGS) has been increasingly used to generate mutation, transcriptome and epigenomic profiles, as well demonstrated by The Cancer Genome Atlas (TCGA) () and the International Cancer Genome Consortium (ICGC) in major cancer types (). It is evident that utilizing NGS-based omics data, individually or in combination, along with clinical metadata, can foster the development of robust biomarkers, such as tumor mutational burden, gene mutation and expression signature, and the classification of disease subtypes, thus benefiting patients in diagnosis, risk evaluation and potentially individualized therapy. In practice, however, prioritization on causal variants and genes still faces key challenges in data processing, harmonization, and clinical interpretation. Misinterpretation of genetic testing results remains a major bottleneck in cases of challenges (). This topic covers research articles that, as we described below, aimed to identify potentially functional variants and genes, or to build models for risk prediction.
Nomogram is a predictive model that is widely used to predict individual’s risk of recurrence, metastases and overall survival (). To build a nomogram for early-stage hepatocellular carcinoma (HCC), Huang et al. downloaded transcriptome, mutation and clinical data for patients from a single cohort in TCGA and another four in ICGC. Cox regression analysis identified seven significant variables, including mutation status of TP53, MACF1, EYS and DOCK2, that were used to build the nomogram. The patients were then divided into low-versus high-risk group, with the former being associated with a better overall survival. Focused analysis of the cohort from TCGA revealed clear differences between the two risk groups in the abundance for seven of the 22 tumor-infiltrating hematopoietic cell subpopulations (); also, the low-risk group had significantly lower Tumor Immune Dysfunction and Exclusion (TIDE) scores (), suggestive of a better immunotherapy response. This study demonstrated a risk stratification nomogram that is potentially linked to the infiltrating immune cell composition in HCC.
Starting with a public RNA-seq data of 117 Ewing sarcoma (ES) patients, Zhou et al. first calculated, for each sample, an immune enrichment score across each of the 28 infiltrating immune cell subpopulations (), followed by unsupervised sample clustering. Two clusters with the highest and lowest overall score were retained. Of the differentially expressed genes (DEGs) between the two clusters, 862 formed a distinct immune-related module that showed the strongest negative correlation with immune score (estimated via the ESTIMATE package). About 10% (85 genes) were DEGs between normal skeletal muscle tissue and ES. They focused on NPM1 (nucleophosmin 1) involved in DNA repair and cell proliferation, showing that its mRNA and protein expression levels were markedly higher in ES cell lines compared to mesenchymal stem cells. The higher mRNA expression correlated with lower immune score, TIDE score and PD-L1 expression, as well as worse prognosis in ES. Importantly, NSC348884, a nucleophosmin inhibitor (), can induce apoptosis in treated ES cells. This work recapitulates the previous finding that NPM1, a drug-targetable gene, is a prognostic biomarker in ES ().
Through total RNA and miRNA sequencing, Wang et al. identified mRNAs, IncRNAs, and miRNAs differentially expressed between acute myeloid leukemia (AML) patients and healthy subjects. They used RAID, a comprehensive RNA-associated interaction database (), to predict mRNAs and lncRNAs targeted by the differential miRNAs. The analysis revealed a potential network of the top 25 hub mRNAs with 15 miRNAs and 12 lncRNAs, including at least four mRNAs and two lncRNAs that are associated with overall survival. Notably, the expression of CCL5 and lncRNA UCA1, known to play key roles in the proliferation of AML, correlated with the fraction of infiltrating immune and stromal cells (). The analysis also revealed a novel interaction between UCA1 and miR-16-5p, expanding the known UCA1-miRNA crosstalk in AML (). Together, this study supports CCL5 and UCA1 as potential diagnostic biomarkers in AML.
Biomarker discovery often relies on the integration of different datasets. In ulcerative colitis (UC), Chen et al. selected six microarray gene expression data from GEO, including 22–162 patients and 11–21 controls. After batch effect correction, 231–436 DEGs were identified from each dataset, with only 79 DEGs in common by a simple intersection approach. To effectively integrate the results, the authors applied the robust rank aggregation (RRA) method, which is robust to outliers and noises (), on the ranked DEG lists. Of the 208 RRA-identified DEGs, six hub genes were selected and confirmed to be upregulated in a UC mouse model. Indeed, these six genes are known to be associated with UC. Thus, to extract biological signatures shared across multiple datasets, one should consider robust meta-analysis approaches for high reproducibility.
Finally, Shestak et al. reported the genetic test of a 14-year-old female athlete, who was suspected to have long QT syndrome (LQTS). WES identified a rare mutation (c.647C > T, p. S216L, chr3:38655522-38655522) in the non-canonical exon 6 of SCN5A. SCN5A is a cardiac ion channel gene implicated in multiple cardiac diseases, with conclusive evidence for its causation in congenital LQTS (). The clinical report, however, mistook this variant for the one previously reported in the canonical exon 6 (c.647C > T, p. S216L, chr3:38655290-38655290) (), leading to misinterpretation. Subsequent Sanger sequencing confirmed a lack of mutation in canonical exon 6. Two more tests were ordered, and both identified the mutation only in the non-canonical exon 6. First, DNA was sequenced in a targeted panel of 11 genes including SCN5A, followed by Sanger sequencing validation. Second, Sanger sequencing revealed the mutation in the mother, but not in the father. The variant was classified as benign, suggesting negative result of the genetic testing. This study highlights the importance of variant validation. Obviously, the collaboration between clinicians and bioinformaticians is vital for genetic counseling. With the ongoing efforts, we are expecting the development of systems for accurately prioritizing causal variants and genes in accelerating biomarker discovery.
Statements
Author contributions
All authors listed have made a substantial contribution to the work and approved it for publication.
Funding
This work is supported by the Mayo Clinic Center for Individualized Medicine.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AdlerA.NovelliV.AminA. S.AbiusiE.CareM.NannenbergE. A.et al (2020). An International, Multicentered, Evidence-Based Reappraisal of Genes Reported to Cause Congenital Long QT Syndrome. Circulation141, 418–428. 10.1161/circulationaha.119.043132
2
BalachandranV. P.GonenM.SmithJ. J.DematteoR. P. (2015). Nomograms in Oncology: More Than Meets the Eye. Lancet Oncol.16, e173–e180. 10.1016/s1470-2045(14)71116-7
3
FarmerM. B.BonadiesD. C.PedersonH. J.MrazK. A.WhatleyJ. W.DarnesD. R.et al (2021). Challenges and Errors in Genetic Testing. Cancer J.27, 417–422. 10.1097/ppo.0000000000000553
4
JiaQ.WuW.WangY.AlexanderP. B.SunC.GongZ.et al (2018). Local Mutational Diversity Drives Intratumoral Immune Heterogeneity in Non-small Cell Lung Cancer. Nat. Commun.9, 5361. 10.1038/s41467-018-07767-w
5
JiangP.GuS.PanD.FuJ.SahuA.HuX.et al (2018). Signatures of T Cell Dysfunction and Exclusion Predict Cancer Immunotherapy Response. Nat. Med.24, 1550–1558. 10.1038/s41591-018-0136-1
6
KikutaK.TochigiN.ShimodaT.YabeH.MoriokaH.ToyamaY.et al (2009). Nucleophosmin as a Candidate Prognostic Biomarker of Ewing's Sarcoma Revealed by Proteomics. Clin. Cancer Res.15, 2885–2894. 10.1158/1078-0432.ccr-08-1913
7
KoldeR.LaurS.AdlerP.ViloJ. (2012). Robust Rank Aggregation for Gene List Integration and Meta-Analysis. Bioinformatics28, 573–580. 10.1093/bioinformatics/btr709
8
MarangoniS.Di RestaC.RocchettiM.BarileL.RizzettoR.SummaA.et al (2011). A Brugada Syndrome Mutation (p.S216L) and its Modulation by p.H558R Polymorphism: Standard and Dynamic Characterization. Cardiovasc. Res.91, 606–616. 10.1093/cvr/cvr142
9
MiliusD.DoveE. S.ChalmersD.DykeS. O. M.KatoK.NicolásP.et al (2014). The International Cancer Genome Consortium's Evolving Data-protection Policies. Nat. Biotechnol.32, 519–523. 10.1038/nbt.2926
10
NewmanA. M.LiuC. L.GreenM. R.GentlesA. J.FengW.XuY.et al (2015). Robust Enumeration of Cell Subsets from Tissue Expression Profiles. Nat. Methods12, 453–457. 10.1038/nmeth.3337
11
QiW.ShakalyaK.StejskalA.GoldmanA.BeeckS.CookeL.et al (2008). NSC348884, a Nucleophosmin Inhibitor Disrupts Oligomer Formation and Induces Apoptosis in Human Cancer Cells. Oncogene27, 4210–4220. 10.1038/onc.2008.54
12
SunM. D.ZhengY. Q.WangL. P.ZhaoH. T.YangS. (2018). Long Noncoding RNA UCA1 Promotes Cell Proliferation, Migration and Invasion of Human Leukemia Cells via Sponging miR-126. Eur. Rev. Med. Pharmacol. Sci.22, 2233–2245. 10.26355/eurrev_201804_14809
13
TomczakK.CzerwińskaP.WiznerowiczM. (2015). Review the Cancer Genome Atlas (TCGA): an Immeasurable Source of Knowledge. Wspólczesna Onkologia1A, 68–77. 10.5114/wo.2014.47136
14
YiY.ZhaoY.LiC.ZhangL.HuangH.LiY.et al (2017). RAID v2.0: an Updated Resource of RNA-Associated Interactions across Organisms. Nucleic Acids Res.45, D115–d118. 10.1093/nar/gkw1052
15
YoshiharaK.ShahmoradgoliM.MartínezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun.4, 2612. 10.1038/ncomms3612
Summary
Keywords
bioinformatics, biomarker, next-generation sequencing, nomogram, RNA sequencing, microarray, whole-exome sequencing
Citation
Tian S, Tu ZJ, Yan H and Klee EW (2022) Editorial: Clinical Genome Sequencing: Bioinformatics Challenges and Key Considerations. Front. Genet. 13:896032. doi: 10.3389/fgene.2022.896032
Received
14 March 2022
Accepted
16 March 2022
Published
31 March 2022
Volume
13 - 2022
Edited and reviewed by
Stephen J. Bush, University of Oxford, United Kingdom
Updates
Copyright
© 2022 Tian, Tu, Yan and Klee.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shulan Tian, tian.shulan@mayo.edu; Eric W. Klee, Klee.Eric@mayo.edu
This article was submitted to Human and Medical Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.