Restricted accessResearch articleFirst published online 2018-06
Harnessing Omics Big Data in Nine Vertebrate Species by Genome-Wide Prioritization of Sequence Variants with the Highest Predicted Deleterious Effect on Protein Function
Harnessing the genomics big data requires innovation in how we extract and interpret biologically relevant variants. Currently, there is no established catalog of prioritized missense variants associated with deleterious protein function phenotypes. We report in this study, to the best of our knowledge, the first genome-wide prioritization of sequence variants with the most deleterious effect on protein function (potentially deleterious variants [pDelVars]) in nine vertebrate species: human, cattle, horse, sheep, pig, dog, rat, mouse, and zebrafish. The analysis was conducted using the Ensembl/BioMart tool. Genes comprising pDelVars in the highest number of examined species were identified using a Python script. Multiple genomic alignments of the selected genes were built to identify interspecies orthologous potentially deleterious variants, which we defined as the “ortho-pDelVars.” Genome-wide prioritization revealed that in humans, 0.12% of the known variants are predicted to be deleterious. In seven out of nine examined vertebrate species, the genes encoding the multiple PDZ domain crumbs cell polarity complex component (MPDZ) and the transforming acidic coiled-coil containing protein 2 (TACC2) comprise pDelVars. Five interspecies ortho-pDelVars were identified in three genes. These findings offer new ways to harness genomics big data by facilitating the identification of functional polymorphisms in humans and animal models and thus provide a future basis for optimization of protocols for whole genome prioritization of pDelVars and screening of orthologous sequence variants. The approach presented here can inform various postgenomic applications such as personalized medicine and multiomics study of health interventions (iatromics).
Get full access to this article
View all access options for this article.
References
1.
Genomes ProjectConsortium, AbecasisGR, AltshulerD, et al. (2010). A map of human genome variation from population-scale sequencing. Nature, 467, 1061–1073.
2.
AdzhubeiI, JordanDM, and SunyaevSR. (2013). Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet, 7, Unit7.20.
3.
AdzhubeiIA, SchmidtS, PeshkinL, et al. (2010). A method and server for predicting damaging missense mutations. Nat Methods, 7, 248–249.
4.
BalbiKJ, and FeilEJ. (2007). The rise and fall of deleterious mutation. Res Microbiol, 158, 779–786.
5.
BarchetW, KrugA, CellaM, et al. (2005). Dendritic cells respond to influenza virus through TLR7- and PKR-independent pathways. Eur J Immunol, 35, 236–242.
6.
BrandomBW, and CallahanPM. (2015). Malignant hyperthermia: An update. Adv Anesth, 33, 113–128.
7.
CareMA, BradfordJR, NeedhamCJ, BulpittAJ, and WestheadDR. (2009). Combining the interactome and deleterious SNP predictions to improve disease gene identification. Hum Mutat, 30, 485–492.
8.
ChenCY, HungLY, WuCS, and ChuangTJ. (2016). Purifying selection shapes the coincident SNP distribution of primate coding sequences. Sci Rep 6,27272.
9.
ChengS, Douglas-JonesA, YangX, ManselRE, and JiangWG. (2010). Transforming acidic coiled-coil-containing protein 2 (TACC2) in human breast cancer, expression pattern and clinical/prognostic relevance. Cancer Genomics Proteomics, 7, 67–73.
10.
ChoiY, and ChanAP. (2015). PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics, 31, 2745–2747.
11.
ChuM, GregorioCC, and PappasCT. (2016). Nebulin, a multi-functional giant. J Exp Biol, 219, 146–152.
12.
ColakE, LeslieA, ZausmerK, et al. (2014). RNA and imidazoquinolines are sensed by distinct TLR7/8 ectodomain sites resulting in functionally disparate signaling events. J Immunol, 192, 5963–5973.
13.
CorreiaACC, SilvaPCB, and da SilvaBA. (2012). Malignant hyperthermia: Clinical and molecular aspects. Rev Bras Anestesiol, 62, 820–837.
FagerlundT, OrdingH, BendixenD, and BergK. (1994). Search for three known mutations in the RYR1 gene in 48 Danish families with malignant hyperthermia. Clin Genet, 46, 401–404.
16.
FehrC, ShirleyRL, BelknapJK, CrabbeJC, and BuckKJ. (2002). Congenic mapping of alcohol and pentobarbital withdrawal liability loci to a <1 centimorgan interval of murine chromosome 4: Identification of Mpdz as a candidate gene. J Neurosci, 22, 3730–3738.
17.
FeldnerA, AdamMG, TetzlaffF, et al. (2017). Loss of Mpdz impairs ependymal cell integrity leading to perinatal-onset hydrocephalus in mice. EMBO Mol Med, 9, 890–905.
18.
GajdaA, HorváthE, HortobágyiT, et al. (2015). Nemaline myopathy type 2 (NEM2): Two novel mutations in the nebulin (NEB) gene. J Child Neurol, 30, 627–630.
19.
GantierMP, TongS, BehlkeMA, et al. (2008). TLR7 is involved in sequence-specific sensing of single-stranded RNAs in human macrophages. J Immunol, 180, 2117–2124.
20.
GillardEF, OtsuK, FujiiJ, et al. (1991). A substitution of cysteine for arginine 614 in the ryanodine receptor is potentially causative of human malignant hyperthermia. Genomics, 11, 751–755.
21.
HaGH, KimJL, and BreuerEKY. (2013). Transforming acidic coiled-coil proteins (TACCs) in human cancer. Cancer Lett, 336, 24–33.
22.
HekimN, and ÖzdemirV. (2017). A general theory for ‘post’ systems biology: Iatromics and the environtome. OMICS, 21, 359–360.
23.
Hernández-OchoaEO, PrattSJP, LoveringRM, and SchneiderMF. (2016). Critical role of intracellular RyR1 calcium release channels in skeletal muscle function and disease. Front Physiol 6,420.
24.
HodgkinsonA, and Eyre-WalkerA. (2010). The genomic distribution and local context of coincident SNPs in human and chimpanzee. Genome Biol Evol, 2, 547–557.
25.
HodgkinsonA, LadoukakisE, and Eyre-WalkerA. (2009). Cryptic variation in the human mutation rate. PLoS Biol 7,e1000027.
26.
HoganK, CouchF, PowersPA, and GreggRG. (1992). A cysteine-for-arginine substitution (R614C) in the human skeletal muscle calcium release channel cosegregates with malignant hyperthermia. Anesth Analg, 75, 441–448.
27.
JagadeeshKA, WengerAM, BergerMJ, et al. (2016). M-CAP eliminates a majority of variants of uncertain significance in clinical exomes at high sensitivity. Nat Genet, 48, 1581–1586.
28.
JiaoX, ShermanBT, Huang daW, et al. (2012). DAVID-WS: A stateful web service to facilitate gene/protein list analysis. Bioinformatics, 28, 1805–1806.
29.
KarpyakVM, KimJH, BiernackaJM, et al. (2009). Sequence variations of the human MPDZ gene and association with alcoholism in subjects with European Ancestry. Alcohol Clin Exp Res, 33, 712–721.
30.
KruseLC, WalterNAR, and BuckKJ. (2014). MPDZ expression in the caudolateral substantia nigra pars reticulata is crucially involved in alcohol withdrawal. Genes Brain Behav, 13, 769–776.
31.
KumarP, HenikoffS, and NgPC. (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc, 4, 1073–1082.
32.
LauffartB, GangisettyO, and StillIH. (2003). Molecular cloning, genomic structure and interactions of the putative breast tumor suppressor TACC2. Genomics, 81, 192–201.
33.
LekM, KarczewskiKJ, MinikelEV, et al. (2016). Analysis of protein-coding genetic variation in 60,706 humans. Nature, 536, 285–291.
34.
MakinoT, RubinCJ, CarneiroM, AxelssonE, AnderssonL, and WebsterMT. (2018). Elevated proportions of deleterious genetic variation in domestic animals and plants. Genome Biol Evol, 10, 276–290.
35.
MalfattiE, MongesS, LehtokariVL, et al. (2015). Bilateral foot-drop as predominant symptom in nebulin (NEB) gene related ‘core-Rod’ congenital myopathy. Eur J Med Genet, 58, 556–561.
36.
MilnerLC, ShirleyRL, KozellLB, et al. (2015). Novel MPDZ/MUPP1 transgenic and knockdown models confirm Mpdz's role in ethanol withdrawal and support its role in voluntary ethanol consumption. Addict Biol, 20, 143–147.
37.
MonnierN, MartyI, FaureJ, et al. (2008). Null Mutations causing depletion of the type 1 ryanodine receptor (RYR1) are commonly associated with recessive structural congenital myopathies with cores. Hum Mutat, 29, 670–678.
OnoderaY, TakagiK, MikiY, et al. (2016). TACC2 (transforming acidic coiled-coil protein 2) in breast carcinoma as a potent prognostic predictor associated with cell proliferation. Cancer Med, 5, 1973–1982.
40.
OtsuK, PhillipsSM, KhannaKV, de LeonS, and MacLennanHD. (1992). Refinement of diagnostic assays for a probable causal mutation for porcine and human malignant hyperthermia. Genomics, 13, 835–837.
41.
OzdemirV, MotulskyAG, KolkerE, and GodardB. (2009). Genome–environment interactions and prospective technology assessment: Evolution from pharmacogenomics to nutrigenomics and ecogenomics. OMICS, 13, 1–6.
42.
PiletičK, and KunejT. (2017). Minimal standards for reporting microRNA:target interactions. OMICS, 21, 197–206.
43.
PirihN, and KunejT. (2017). Toward a taxonomy for multi-omics science? Terminology development for whole genome study approaches by omics technology and hierarchy. OMICS, 21, 1–16.
44.
RedenšekS, DolžanV, and KunejT. (2018). From genomics to omics landscapes of Parkinson's disease: Revealing the molecular mechanisms. OMICS, 21, 1–16.
45.
SachidanandamR, WeissmanD, SchmidtSC, et al. (2001). A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms. Nature, 409, 928–933.
46.
Saugier-VeberP, MarguetF, LecoquierreF, et al. (2017). Hydrocephalus due to multiple ependymal malformations is caused by mutations in the MPDZ gene. Acta Neuropathol Commun 5,36.
47.
SkrljB, KoncJ, and KunejT. (2017). Identification of sequence variants within experimentally validated protein interaction sites provides new insights into molecular mechanisms of disease development. Mol Inform, 37, 1–11.
48.
SkrljB, and KunejT. (2016). Computational identification of non-synonymous polymorphisms within regions corresponding to protein interaction sites. Comput Biol Med, 79, 30–35.
49.
TakayamaK, Horie-InoueK, SuzukiT, et al. (2012). TACC2 is an androgen-responsive cell cycle regulator promoting androgen-mediated and castration-resistant growth of prostate cancer. Mol Endocrinol, 26, 748–761.
50.
WongKC, and ZhangZ. (2014). SNPdryad: Predicting deleterious non-synonymous human SNPs using only orthologous protein sequences. Bioinformatics, 30, 1112–1119.
51.
XueY, ChenY, AyubQ, et al. (2012). Deleterious- and disease-allele prevalence in healthy individuals: Insights from current predictions, mutation databases, and population-scale resequencing. Am J Hum Genet, 91, 1022–1032.
52.
ZhouH, JungbluthH, SewryCA, et al. (2007). Molecular mechanisms and phenotypic variation in RYR1-related congenital myopathies. Brain, 130, 2024–2036.
53.
ZorcM, and KunejT. (2016). In silico screening of the chicken genome for overlaps between genomic regions: microRNA genes, coding and non-coding transcriptional units, QTL, and genetic variations. Chromosome Res, 24, 225–230.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
