Polygenic basis and biomedical consequences of telomere length variation.

Telomeres, the end fragments of chromosomes, play key roles in cellular proliferation and senescence. Here we characterize the genetic architecture of naturally occurring variation in leukocyte telomere length (LTL) and identify causal links between LTL and biomedical phenotypes in 472,174 well-characterized UK Biobank participants. We identified 197 independent sentinel variants associated with LTL at 138 genomic loci (108 new). Genetically determined differences in LTL were associated with multiple biological traits, ranging from height to bone marrow function, as well as several diseases spanning neoplastic, vascular and inflammatory pathologies. Finally, we estimated that, at the age of 40 years, people with an LTL >1 s.d. shorter than the population mean had a 2.5-year-lower life expectancy compared with the group with ≥1 s.d. longer LDL. Overall, we furnish new insights into the genetic regulation of LTL, reveal wide-ranging influences of LTL on physiological traits, diseases and longevity, and provide a powerful resource available to the global research community.

Rare loss-of-function mutations of PTGIR are enriched in fibromuscular dysplasia.

AIMS: Fibromuscular dysplasia (FMD) and spontaneous coronary artery dissection (SCAD) are related, non-atherosclerotic arterial diseases mainly affecting middle-aged women. Little is known about their physiopathological mechanisms. We aimed to identify rare genetic causes to elucidate molecular mechanisms implicated in FMD and SCAD. METHODS AND RESULTS: We analysed 29 exomes that included familial and sporadic FMD. We identified one rare loss-of-function variant (LoF) (frequencygnomAD = 0.000075) shared by two FMD sisters in the prostaglandin I2 receptor gene (PTGIR), a key player in vascular remodelling. Follow-up was conducted by targeted or Sanger sequencing (1071 FMD and 363 SCAD patients) or lookups in exome (264 FMD) or genome sequences (480 SCAD), all independent and unrelated. It revealed four additional LoF allele carriers, in addition to several rare missense variants, among FMD patients, and two LoF allele carriers among SCAD patients, including one carrying a rare splicing mutation (c.768 + 1C>G). We used burden test to test for enrichment in patients compared to gnomAD controls, which detected a putative enrichment in FMD (PTRAPD = 8 × 10-4), but not a significant enrichment (PTRAPD = 0.12) in SCAD. The biological effects of variants on human prostaclycin receptor (hIP) signalling and protein expression were characterized using transient overexpression in human cells. We confirmed the LoFs (Q163X and P17RfsX6) and one missense (L67P), identified in one FMD and one SCAD patient, to severely impair hIP function in vitro. CONCLUSIONS: Our study shows that rare genetic mutations in PTGIR are enriched among FMD patients and found in SCAD patients, suggesting a role for prostacyclin signalling in non-atherosclerotic stenosis and dissection.

Spontaneous Coronary Artery Dissection: Insights on Rare Genetic Variation From Genome Sequencing.

BACKGROUND: Spontaneous coronary artery dissection (SCAD) occurs when an epicardial coronary artery is narrowed or occluded by an intramural hematoma. SCAD mainly affects women and is associated with pregnancy and systemic arteriopathies, particularly fibromuscular dysplasia. Variants in several genes, such as those causing connective tissue disorders, have been implicated; however, the genetic architecture is poorly understood. Here, we aim to better understand the diagnostic yield of rare variant genetic testing among a cohort of SCAD survivors and to identify genes or gene sets that have a significant enrichment of rare variants. METHODS: We sequenced a cohort of 384 SCAD survivors from the United Kingdom, alongside 13 722 UK Biobank controls and a validation cohort of 92 SCAD survivors. We performed a research diagnostic screen for pathogenic variants and exome-wide and gene-set rare variant collapsing analyses. RESULTS: The majority of patients within both cohorts are female, 29% of the study cohort and 14% validation cohort have a remote arteriopathy. Four cases across the 2 cohorts had a diagnosed connective tissue disorder. We identified pathogenic or likely pathogenic variants in 7 genes (PKD1, COL3A1, SMAD3, TGFB2, LOX, MYLK, and YY1AP1) in 14/384 cases in the study cohort and in 1/92 cases in the validation cohort. In our rare variant collapsing analysis, PKD1 was the highest-ranked gene, and several functionally plausible genes were enriched for rare variants, although no gene achieved study-wide statistical significance. Gene-set enrichment analysis suggested a role for additional genes involved in renal function. CONCLUSIONS: By studying the largest sequenced cohort of SCAD survivors, we demonstrate that, based on current knowledge, only a small proportion have a pathogenic variant that could explain their disease. Our findings strengthen the overlap between SCAD and renal and connective tissue disorders, and we highlight several new genes for future validation.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Adult height and risk of 50 diseases: a combined epidemiological and genetic analysis.

BACKGROUND: Adult height is associated with risk of several diseases, but the breadth of such associations and whether these associations are primary or due to confounding are unclear. We examined the association of adult height with 50 diseases spanning multiple body systems using both epidemiological and genetic approaches, the latter to identify un-confounded associations and possible underlying mechanisms. METHODS: We examined the associations for adult height (using logistic regression adjusted for potential confounders) and genetically determined height (using a two-sample Mendelian randomisation approach with height-associated genetic variants as instrumental variables) in 417,434 individuals of white ethnic background participating in the UK Biobank. We undertook pathway analysis of height-associated genes to identify biological processes that could link height and specific diseases. RESULTS: Height was associated with 32 diseases and genetically determined height associated with 12 diseases. Of these, 11 diseases showed a concordant association in both analyses, with taller height associated with reduced risks of coronary artery disease (odds ratio per standard deviation (SD) increase in height ORepi = 0.80, 95% CI 0.78-0.81; OR per SD increase in genetically determined height ORgen = 0.86, 95% CI 0.82-0.90), hypertension (ORepi = 0.83, 95% CI 0.82-0.84; ORgen = 0.88, 95% CI 0.85-0.91), gastro-oesophageal reflux disease (ORepi = 0.85, 95% CI 0.84-0.86; ORgen = 0.94, 95% CI 0.92-0.97), diaphragmatic hernia (ORepi = 0.81, 95% CI 0.79-0.82; ORgen = 0.91, 95% CI 0.88-0.94), but increased risks of atrial fibrillation (ORepi = 1.42, 95% CI 1.38-1.45; ORgen = 1.33, 95% CI 1.26-1.40), venous thromboembolism (ORepi = 1.18, 95% CI 1.16-1.21; ORgen = 1.15, 95% CI 1.11-1.19), intervertebral disc disorder (ORepi = 1.15, 95% CI 1.13-1.18; ORgen = 1.14, 95% CI 1.09-1.20), hip fracture (ORepi = 1.19, 95% CI 1.12-1.26; ORgen = 1.27, 95% CI 1.17-1.39), vasculitis (ORepi = 1.15, 95% CI 1.11-1.19; ORgen = 1.20, 95% CI 1.14-1.28), cancer overall (ORepi = 1.09, 95% CI 1.08-1.11; ORgen = 1.06, 95% CI 1.04-1.08) and breast cancer (ORepi = 1.08, 95% CI 1.06-1.10; ORgen = 1.07, 95% CI 1.03-1.11). Pathway analysis showed multiple height-associated pathways associating with individual diseases. CONCLUSIONS: Adult height is associated with risk of a range of diseases. We confirmed previously reported height associations for coronary artery disease, atrial fibrillation, venous thromboembolism, intervertebral disc disorder, hip fracture and cancer and identified potential novel associations for gastro-oesophageal reflux disease, diaphragmatic hernia and vasculitis. Multiple biological mechanisms affecting height may affect the risks of these diseases.

Systematic Evaluation of Pleiotropy Identifies 6 Further Loci Associated With Coronary Artery Disease.

BACKGROUND: Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseases or traits. OBJECTIVES: This study sought to systematically test if genetic variants identified for non-CAD diseases/traits also associate with CAD and to undertake a comprehensive analysis of the extent of pleiotropy of all CAD loci. METHODS: In discovery analyses involving 42,335 CAD cases and 78,240 control subjects we tested the association of 29,383 common (minor allele frequency >5%) single nucleotide polymorphisms available on the exome array, which included a substantial proportion of known or suspected single nucleotide polymorphisms associated with common diseases or traits as of 2011. Suggestive association signals were replicated in an additional 30,533 cases and 42,530 control subjects. To evaluate pleiotropy, we tested CAD loci for association with cardiovascular risk factors (lipid traits, blood pressure phenotypes, body mass index, diabetes, and smoking behavior), as well as with other diseases/traits through interrogation of currently available genome-wide association study catalogs. RESULTS: We identified 6 new loci associated with CAD at genome-wide significance: on 2q37 (KCNJ13-GIGYF2), 6p21 (C2), 11p15 (MRVI1-CTR9), 12q13 (LRP1), 12q24 (SCARB1), and 16q13 (CETP). Risk allele frequencies ranged from 0.15 to 0.86, and odds ratio per copy of the risk allele ranged from 1.04 to 1.09. Of 62 new and known CAD loci, 24 (38.7%) showed statistical association with a traditional cardiovascular risk factor, with some showing multiple associations, and 29 (47%) showed associations at p < 1 × 10-4 with a range of other diseases/traits. CONCLUSIONS: We identified 6 loci associated with CAD at genome-wide significance. Several CAD loci show substantial pleiotropy, which may help us understand the mechanisms by which these loci affect CAD risk.

Screening in silico predicted remotely acting NF1 gene regulatory elements for mutations in patients with neurofibromatosis type 1.

Neurofibromatosis type 1 (NF1), a neuroectodermal disorder, is caused by germline mutations in the NF1 gene. NF1 affects approximately 1/3,000 individuals worldwide, with about 50% of cases representing de novo mutations. Although the NF1 gene was identified in 1990, the underlying gene mutations still remain undetected in a small but obdurate minority of NF1 patients. We postulated that in these patients, hitherto undetected pathogenic mutations might occur in regulatory elements far upstream of the NF1 gene. In an attempt to identify such remotely acting regulatory elements, we reasoned that some of them might reside within DNA sequences that (1) have the potential to interact at distance with the NF1 gene and (2) lie within a histone H3K27ac-enriched region, a characteristic of active enhancers. Combining Hi-C data, obtained by means of the chromosome conformation capture technique, with data on the location and level of histone H3K27ac enrichment upstream of the NF1 gene, we predicted in silico the presence of two remotely acting regulatory regions, located, respectively, approximately 600 kb and approximately 42 kb upstream of the NF1 gene. These regions were then sequenced in 47 NF1 patients in whom no mutations had been found in either the NF1 or SPRED1 gene regions. Five patients were found to harbour DNA sequence variants in the distal H3K27ac-enriched region. Although these variants are of uncertain pathological significance and still remain to be functionally characterized, this approach promises to be of general utility for the detection of mutations underlying other inherited disorders that may be caused by mutations in remotely acting regulatory elements.

Exploring the somatic NF1 mutational spectrum associated with NF1 cutaneous neurofibromas.

Neurofibromatosis type-1 (NF1), caused by heterozygous inactivation of the NF1 tumour suppressor gene, is associated with the development of benign and malignant peripheral nerve sheath tumours (MPNSTs). Although numerous germline NF1 mutations have been identified, relatively few somatic NF1 mutations have been described in neurofibromas. Here we have screened 109 cutaneous neurofibromas, excised from 46 unrelated NF1 patients, for somatic NF1 mutations. NF1 mutation screening (involving loss-of-heterozygosity (LOH) analysis, multiplex ligation-dependent probe amplification and DNA sequencing) identified 77 somatic NF1 point mutations, of which 53 were novel. LOH spanning the NF1 gene region was evident in 25 neurofibromas, but in contrast to previous data from MPNSTs, it was absent at the TP53, CDKN2A and RB1 gene loci. Analysis of DNA/RNA from neurofibroma-derived Schwann cell cultures revealed NF1 mutations in four tumours whose presence had been overlooked in the tumour DNA. Bioinformatics analysis suggested that four of seven novel somatic NF1 missense mutations (p.A330T, p.Q519P, p.A776T, p.S1463F) could be of functional/clinical significance. Functional analysis confirmed this prediction for p.S1463F, located within the GTPase-activating protein-related domain, as this mutation resulted in a 150-fold increase in activated GTP-bound Ras. Comparison of the relative frequencies of the different types of somatic NF1 mutation observed with those of their previously reported germline counterparts revealed significant (P=0.001) differences. Although non-identical somatic mutations involving either the same or adjacent nucleotides were identified in three pairs of tumours from the same patients (P<0.0002), no association was noted between the type of germline and somatic NF1 lesion within the same individual.

Comparative analysis of germline and somatic microlesion mutational spectra in 17 human tumor suppressor genes.

Mutations associated with tumorigenesis may either arise somatically or can be inherited through the germline. We performed a comparison of somatic, germline, shared (found in both soma and germline) and somatic recurrent mutational spectra for 17 human tumor suppressor genes, which focused upon missense single base-pair substitutions and microdeletions/microinsertions. Somatic and germline mutational spectra were similar in relation to C.G>T.A transitions but differed with respect to the frequency of A.T>G.C, A.T>T.A, and C.G>A.T substitutions. Shared missense mutations were characterized by higher mutability rates, greater physicochemical differences between wild-type and mutant residues, and a tendency to occur in evolutionarily conserved residues and within CpG/CpHpG oligonucleotides. Mononucleotide runs (≥4 bp) were identified as hotspots for shared microdeletions/microinsertions. Both germline and somatic microdeletions/microinsertions were found to be significantly overrepresented within the "indel-hotspot" motif, GTAAGT. Using a naïve Bayes' classifier trained to discriminate between five missense mutation groups, 63% of mutations in our dataset were on average correctly recognized. Applying this classifier to an independent dataset of probable driver mutations, we concluded that ∼50% of these somatic missense mutations possess features consistent with their being either shared or recurrent, suggesting that a disproportionate number of such lesions are likely to be drivers of tumorigenesis.