Novel genes and variants associated with longevity in Bulgarian centenarians revealed by whole exome sequencing DNA pools: a pilot study

Abstract

Aim: To determine specific genetic loci that might be associated with longevity in Bulgarian population by analyzing exome pool-seq data from centenarians and a control group.

Methods: We performed whole-exome sequencing of two DNA pools, set up of 32 Bulgarian centenarians and 61 young healthy controls, respectively, and 59935 quality filtered variants were concurrently detected in both pools. Fisher’s exact test was employed to establish the significance of allele frequency difference between the pools.

Results: Forty seven variants showed significantly higher allele frequency in the centenarian compared to the control pool, and these can be considered to be positively associated with longevity in Bulgarian populaton. Based on their assigned functional role, three genes containing three of these variants were further investigated. These genes, RNF43, WNK1 and NADSYN1, are involved in evolutionary conserved processes with well ascertained association with longevity, i.e., Wnt signaling pathway, insulin/IGF-1 signal pathway and redox balancing processes, respectively. For the remaining genes exhibiting variants with significantly higher allele frequency in the Bulgarian centenarian pool there is not enough evidence about their functional role in determining longevity and further research is needed.

Conclusion: The results confirm the importance of studying centenarians in different populations to discover those combinations of variants that associate with longer health span.

Keywords

Centenarian / exome / Bulgarian

References

[1] Eatock D. Demographic outlook for the European Union 2019. European Union; 2019. Available from: https://www.europarl.europa.eu/thinktank/en/document.html?reference=EPRS_IDA(2019)637955. [Last accessed on 27 Sep 2020]

[2] NSI. Statistical Reference Book. Bulgaria: National Statistical Institute; 2020. Available from: https://www.nsi.bg/en/content/18296/%D0%BF%D1%83%D0%B1%D0%BB%D0%B8%D0%BA%D0%B0%D1%86%D0%B8%D1%8F/statistical-reference-book-2020-bulgarian-version. [Last accessed on 27 Sep 2020]

[3] HerskindAM,HolmNV,HarvaldB.The heritability of human longevity: a population-based study of 2872 Danish twin pairs born 1870-1900..Hum Genet1996;97:319-23

[4] DeelenJ,UhHW,KuningasM.Genome-wide association study identifies a single major locus contributing to survival into old age; the APOE locus revisited..Aging Cell2011;10:686-98 PMCID:PMC3193372

[5] SebastianiP,DewanAT,PucaA.Genetic signatures of exceptional longevity in humans..PLoS One2012;7:e29848 PMCID:PMC3261167

[6] GkikasI,TavernarakisN.Longevity pathways and memory aging..Front Genet2014;5:155 PMCID:PMC4044971

[7] SerbezovD,HadjidekovaS.Genomics of longevity: recent insights from research on centenarians..Biotechnol Biotechnol Equip2018;32:1359-66

[8] Santos-LozanoA,Pareja-GaleanoH,EmanueleE.Where are supercentenarians located? A worldwide demographic study..Rejuvenation Res2015;18:14-9

[9] SchlöttererC,KoflerR.Sequencing pools of individuals - mining genome-wide polymorphism data without big funding..Nat Rev Genet2014;15:749-63

[10] FracassettiM,WilliY.Validation of pooled whole-genome re-sequencing in arabidopsis lyrata..PLoS One2015;10:e0140462 PMCID:PMC4604096

[11] YangH.Genomic variant annotation and prioritization with ANNOVAR and wANNOVAR..Nat Protoc2015;10:1556-66 PMCID:PMC4718734

[12] KarczewskiKJ,TiaoG,AlföldiJ.Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes..bioRxiv2019;531210

[13] BenjaminiY.Controlling the false discovery rate: a practical and powerful approach to multiple testing..J R Stat Soc B1995;57:289-300

[14] Team RC. R: a language and environment for statistical computing Vienna 2018 Available from: https://www.gbif.org/zh/tool/81287/r-a-language-and-environment-for-statistical-computing. [Last accessed on 27 Sep 2020]

[15] NSI. European health interview. In: Institute NS, editor. 2008. Available from: https://ec.europa.eu/eurostat/web/microdata/european-health-interview-survey. [Last accessed on 27 Sep 2020]

[16] BudovskyA,WangJ,CsordasA.LongevityMap: a database of human genetic variants associated with longevity..Trends Genet2013;29:559-60

[17] ChenJ,AronowBJ.ToppGene suite for gene list enrichment analysis and candidate gene prioritization..Nucleic Acids Res2009;37:W305-11 PMCID:PMC2703978

[18] FabregatA,MatthewsL,GillespieM.The reactome pathway knowledgebase..Nucleic Acids Res2018;46:D649-55 PMCID:PMC5753187

[19] ChenY,LinAC,WeiQ.Beta-catenin signaling plays a disparate role in different phases of fracture repair: implications for therapy to improve bone healing..PLoS Med2007;4:e249 PMCID:PMC1950214

[20] ItoM,AndlT,KimN.Wnt-dependent de novo hair follicle regeneration in adult mouse skin after wounding..Nature2007;447:316-20

[21] ReyaT.Wnt signalling in stem cells and cancer..Nature2005;434:843-50

[22] GiannakisM,JasmineMu X,RosenbluhJ.RNF43 is frequently mutated in colorectal and endometrial cancers..Nat Genet2014;46:1264-6 PMCID:PMC4283570

[23] MandaiS,NomuraN,AraiY.WNK1 regulates skeletal muscle cell hypertrophy by modulating the nuclear localization and transcriptional activity of FOXO4..Sci Rep2018;8:9101 PMCID:PMC6002401

[24] KlotzLO,Prieto-ArroyoI,SteinbrennerH.Redox regulation of FoxO transcription factors..Redox Biol2015;6:51-72 PMCID:PMC4511623

[25] NewhouseS,WallaceC,BurkeB.Polymorphisms in the WNK1 gene are associated with blood pressure variation and urinary potassium excretion..PLoS One2009;4:e5003 PMCID:PMC2661139

[26] PutkuM,OrgE,ComasD.HYPertension in ESTonia (HYPEST), BRItish Genetics of HyperTension (BRIGHT). Novel polymorphic AluYb8 insertion in the WNK1 gene is associated with blood pressure variation in Europeans..Hum Mutat2011;32:806-14 PMCID:PMC3298642

[27] SchultzMB.Why NAD(+) declines during aging: It’s destroyed..Cell Metab2016;23:965-6 PMCID:PMC5088772

[28] MaS,GaleckiA,BurantCF.Cell culture-based profiling across mammals reveals DNA repair and metabolism as determinants of species longevity..Elife2016;5:e19130 PMCID:PMC5148604

[29] MarinTL,MartinM,SmithL.Identification of AMP-activated protein kinase targets by a consensus sequence search of the proteome..BMC Syst Biol2015;9:13 PMCID:PMC4357066

[30] KuanV,GriffithsCJ,WaltonR.DHCR7 mutations linked to higher vitamin D status allowed early human migration to northern latitudes..BMC Evol Biol2013;13:144 PMCID:PMC3708787

[31] MarkKA,BhaumikD,DavisS.Vitamin D promotes protein homeostasis and longevity via the stress response pathway genes skn-1, ire-1, and xbp-1..Cell Rep2016;17:1227-37 PMCID:PMC5689451