Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Jan 12;2(1):100084.
doi: 10.1016/j.xgen.2021.100084. Epub 2022 Jan 13.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Margaret A Taub  1 Matthew P Conomos  2 Rebecca Keener  3 Kruthika R Iyer  4 Joshua S Weinstock  5   6 Lisa R Yanek  7 John Lane  8 Tyne W Miller-Fleming  9 Jennifer A Brody  10 Laura M Raffield  11 Caitlin P McHugh  2 Deepti Jain  2 Stephanie M Gogarten  2 Cecelia A Laurie  2 Ali Keramati  12 Marios Arvanitis  13 Albert V Smith  5   6 Benjamin Heavner  2 Lucas Barwick  14 Lewis C Becker  7 Joshua C Bis  10 John Blangero  15 Eugene R Bleecker  16   17 Esteban G Burchard  18   19 Juan C Celedón  20 Yen Pei C Chang  21 Brian Custer  22   23 Dawood Darbar  24 Lisa de las Fuentes  25 Dawn L DeMeo  26   27 Barry I Freedman  28 Melanie E Garrett  29   30 Mark T Gladwin  31 Susan R Heckbert  32   33 Bertha A Hidalgo  34 Marguerite R Irvin  34 Talat Islam  35 W Craig Johnson  36 Stefan Kaab  37   38 Lenore Launer  39 Jiwon Lee  40 Simin Liu  41 Arden Moscati  42 Kari E North  43 Patricia A Peyser  44 Nicholas Rafaels  45 Christine Seidman  46 Daniel E Weeks  47   48 Fayun Wen  49 Marsha M Wheeler  50 L Keoki Williams  51 Ivana V Yang  45 Wei Zhao  44 Stella Aslibekyan  34 Paul L Auer  52 Donald W Bowden  53 Brian E Cade  54   27 Zhanghua Chen  35 Michael H Cho  26 L Adrienne Cupples  55   56 Joanne E Curran  15 Michelle Daya  45 Ranjan Deka  57 Celeste Eng  18 Tasha E Fingerlin  58   59 Xiuqing Guo  60 Lifang Hou  61 Shih-Jen Hwang  62 Jill M Johnsen  63   64 Eimear E Kenny  65   42 Albert M Levin  66 Chunyu Liu  56   67 Ryan L Minster  47 Take Naseri  68   69 Mehdi Nouraie  31 Muagututi'a Sefuiva Reupena  70 Ester C Sabino  71 Jennifer A Smith  44 Nicholas L Smith  32   33 Jessica Lasky Su  26   27 James G Taylor  49 Marilyn J Telen  29   72 Hemant K Tiwari  73 Russell P Tracy  74 Marquitta J White  18 Yingze Zhang  31 Kerri L Wiggins  10 Scott T Weiss  26   27 Ramachandran S Vasan  75   56 Kent D Taylor  60 Moritz F Sinner  37   38 Edwin K Silverman  26   27 M Benjamin Shoemaker  76 Wayne H-H Sheu  77 Frank Sciurba  78 David A Schwartz  45 Jerome I Rotter  79 Daniel Roden  80 Susan Redline  54   81 Benjamin A Raby  82   83 Bruce M Psaty  84 Juan M Peralta  15 Nicholette D Palmer  53 Sergei Nekhai  49 Courtney G Montgomery  85 Braxton D Mitchell  21   86 Deborah A Meyers  16   17 Stephen T McGarvey  69 NHLBI CARE NetworkAngel Cy Mak  18 Ruth Jf Loos  42   87 Rajesh Kumar  88 Charles Kooperberg  89 Barbara A Konkle  63   64 Shannon Kelly  22   90 Sharon Lr Kardia  44 Robert Kaplan  91 Jiang He  92 Hongsheng Gui  51 Frank D Gilliland  35 Bruce D Gelb  93 Myriam Fornage  94   95 Patrick T Ellinor  96 Mariza de Andrade  97 Adolfo Correa  98 Yii-Der Ida Chen  60 Eric Boerwinkle  99 Kathleen C Barnes  45 Allison E Ashley-Koch  29   30 Donna K Arnett  100 NHLBI Trans-Omics for Precision Medicine (TOPMed) ConsortiumTOPMed Hematology and Hemostasis Working GroupTOPMed Structural Variation Working GroupCathy C Laurie  2 Goncalo Abecasis  5   101 Deborah A Nickerson  50 James G Wilson  102 Stephen S Rich  103 Daniel Levy  56   67 Ingo Ruczinski  1 Abraham Aviv  104 Thomas W Blackwell  5   6 Timothy Thornton  105 Jeff O'Connell  106   107 Nancy J Cox  108 James A Perry  21 Mary Armanios  109 Alexis Battle  3   110 Nathan Pankratz  8 Alexander P Reiner  111   89 Rasika A Mathias  7
Affiliations

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Margaret A Taub et al. Cell Genom. .

Abstract

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

PubMed Disclaimer

Conflict of interest statement

Declaration of Interests: The authors declare the following competing interests: J.C.C. has received research materials from GlaxoSmithKline and Merck (inhaled steroids) and Pharmavite (vitamin D and placebo capsules) to provide medications free of cost to participants in NIH-funded studies, unrelated to the current work. B.I.F. is a consultant for AstraZeneca Pharmaceuticals and RenalytixAI L.W. is on the advisory board for GlaxoSmithKline and receives grant funding from NIAID, NHLBI, and NIDDK, NIH I.V.Y. is a consultant for ElevenP15 S.A. receives equity and salary from 23andMe, Inc. M.H.C. receives grant support from GlaxoSmithKline S.T.W. receives royalties from UpToDate E.K.S. received grant support from GlaxoSmithKline and Bayer in the past three years. B.M.P. serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. F.D.M. is supported by grants from NIH/NHLBI (HL139054,HL091889,HL132523,HL130045,HL098112,HL056177), the NIH/NIEHS (ES006614), the NIH/NIAID (AI126614), and the NIH/ Office of Director (OD023282). Vifor Pharmaceuticals provided medicine and additional funding to support recruitment for HL130045. Dr. Martinez is a council member for the Council for the Developing Child. P.T.E. is supported by a grant from Bayer AG to the Broad Institute focused on the genetics and therapeutics of cardiovascular diseases. Dr. Ellinor has also served on advisory boards or consulted for Bayer AG, Quest Diagnostics, and Novartis. K.C.B. receives royalties from UpToDate G.A. is an employee of Regeneron Pharmaceuticals and owns stock and stock options for Regeneron Pharmaceuticals. A.M. is an employee of Regeneron Pharmaceuticals and owns stock and stock options for Regeneron Pharmaceuticals. A.B. is a consultant for Third Rock Ventures, LLC and holds stock in Google, Inc. D.A.S. is the founder and chief scientific officer of Eleven P15, a company focused on the early diagnosis and treatment of pulmonary fibrosis

Figures

None
Graphical abstract
Figure 1
Figure 1
Genome-wide Manhattan plot (A) Pie chart showing population groups based on HARE for samples included in analysis: European (green, n = 51,654), African (orange, n = 29,260), Hispanic/Latino (purple, n = 18,019), Asian (red, n = 5,683), and Other/Mixed/Unknown (gray, n = 4,506). (B) Trans-population genome-wide tests for association using 163 million sequence-identified variants on n = 109,122 samples with sequence-generated telomere length from TOPMed. All loci had a peak p < 5 × 10−9 in the pooled trans-population analysis. Previously reported loci for TL are indicated in red, and loci newly associated in the present study are indicated in blue. Note the shift in scale above the y axis break; no peak variants had a p value within the y axis break.
Figure 2
Figure 2
LocusZoom plots for multi-hit loci and TINF2 (A) LocusZoom plots for all loci with >1 sentinel variant. Linkage disequilibrium (LD) was calculated from the set of samples used in the analysis with respect to the peak variant in the pooled trans-population primary analysis, thereby reflecting LD patterns specific to the TOPMed samples. For each figure, the peak sentinel variant from the pooled trans-population analysis is indexed and labeled in purple, and all of the independent variants identified through the iterative conditional approach are labeled in green and highlighted with green dotted lines. (B) LocusZoom plots for 4 population groups for the TINF2 locus. (C) Forest plots displaying effect sizes and standard errors, as well as minor allele frequencies, by population group for the 3 sentinel variants in TINF2. See also Table S2.
Figure 3
Figure 3
Fine-mapping of multiple OBFC1 signals (A) LocusZoom plot of the OBFC1 locus, where green dotted lines indicate each independent signal, as in Figure 2. (B) Roadmap Epigenomics Consortium data in hg19 coordinates for skeletal muscle tissue, Primary T CD4+ memory cells from peripheral blood, and primary T CD8+ naive cells from peripheral blood (Roadmap samples E108, E037, and E047, respectively; data were not available for sun-exposed skin). The ChromHMM state model is shown for the 18-state auxiliary model. The state model suggests the primary (rs9420907), secondary (rs111447985), and tertiary (rs112163720) signals are in the promoter region, while the quaternary signal (rs10883948) is in an enhancer region in all Roadmap blood cell types but is transcriptional for peripheral blood monocytes and CD19+ B cells. (C–E) GWAS and eQTL results for the primary (C), tertiary (D), and quaternary (E) signals. The top panels are the GWAS summary statistics from the primary, and iterative conditional analyses that were used to perform colocalization analysis (secondary signal was rare and not available for colocalization). Bottom panels are eQTLs for OBFC1 in the indicated tissue from GTEx. The GTEx eQTLs for these tissues do not colocalize with one another (PPH4 < 4.4 × 10−7), and each signal did not significantly colocalize in the other tissues. LD was calculated from the pooled trans-population samples with respect to the sentinel (black diamond). See also Figures S4 and S5 and Table S5.

References

    1. Aviv A., Shay J.W. Reflections on telomere dynamics and ageing-related diseases in humans. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2018;373:20160436. - PMC - PubMed
    1. McNally E.J., Luncsford P.J., Armanios M. Long telomeres and cancer risk: the price of cellular immortality. J. Clin. Invest. 2019;129:3474–3481. - PMC - PubMed
    1. Codd V., Nelson C.P., Albrecht E., Mangino M., Deelen J., Buxton J.L., Hottenga J.J., Fischer K., Esko T., Surakka I., et al. Identification of seven loci affecting mean telomere length and their association with disease. Nat. Genet. 2013;45:422–427. - PMC - PubMed
    1. Codd V., Mangino M., van der Harst P., Braund P.S., Kaiser M., Beveridge A.J., Rafelt S., Moore J., Nelson C., Soranzo N., et al. Wellcome Trust Case Control Consortium Common variants near TERC are associated with mean telomere length. Nat. Genet. 2010;42:197–199. - PMC - PubMed
    1. Delgado D.A., Zhang C., Chen L.S., Gao J., Roy S., Shinkle J., Sabarinathan M., Argos M., Tong L., Ahmed A., et al. Genome-wide association study of telomere length among South Asians identifies a second RTEL1 association signal. J. Med. Genet. 2018;55:64–71. - PMC - PubMed

Grants and funding