Genetic meta-analysis of twin birth weight shows high genetic correlation with singleton birth weight.

Birth weight (BW) is an important predictor of newborn survival and health and has associations with many adult health outcomes, including cardiometabolic disorders, autoimmune diseases and mental health. On average, twins have a lower BW than singletons as a result of a different pattern of fetal growth and shorter gestational duration. Therefore, investigations into the genetics of BW often exclude data from twins, leading to a reduction in sample size and remaining ambiguities concerning the genetic contribution to BW in twins. In this study, we carried out a genome-wide association meta-analysis of BW in 42 212 twin individuals and found a positive correlation of beta values (Pearson's r = 0.66, 95% confidence interval [CI]: 0.47-0.77) with 150 previously reported genome-wide significant variants for singleton BW. We identified strong positive genetic correlations between BW in twins and numerous anthropometric traits, most notably with BW in singletons (genetic correlation [rg] = 0.92, 95% CI: 0.66-1.18). Genetic correlations of BW in twins with a series of health-related traits closely resembled those previously observed for BW in singletons. Polygenic scores constructed from a genome-wide association study on BW in the UK Biobank demonstrated strong predictive power in a target sample of Dutch twins and singletons. Together, our results indicate that a similar genetic architecture underlies BW in twins and singletons and that future genome-wide studies might benefit from including data from large twin registers.

Quickly home again: patients' experiences of early discharge after minor stroke.

BACKGROUND: Patients with minor stroke experience some of the same issues as patients experiencing stroke of increased severity such as fatigue, anxiety and cognitive symptoms. It is current practice that patients with minor stroke receive accelerated treatment and care, yet studies indicate that patients find it difficult to return to their everyday lives after being discharged. We aimed to explore how patients with minor stroke experience the transitional period from the hospital through the first 2-4 weeks after an accelerated care pathway with discharge within 72 hours after stroke onset. METHODS: A qualitative study consisting of semi-structured interviews with 11 patients experiencing first-time stroke 2-4 weeks after discharge. RESULTS: The patients struggled to identify themselves as having had a stroke. They strived to find a new everyday life, but were challenged by existential concerns, mental fatigue and the fear of having a stroke again. Unresolved questions and misunderstandings arose, and the patients expressed a need for health professionals to support them and discuss unclear issues after discharge. Patients searched for others with similar issues in order to find a new sense of self. CONCLUSION: Patients with minor stroke struggle with everyday life after discharge. There is a need for support after discharge from healthcare professionals with specialised knowledge of stroke. Patients also requested an opportunity to meet other patients with minor stroke.

Cohort Differences in the Associations of Selected Candidate Genes With Risk of All-Cause Mortality at Advanced Ages.

Considerable efforts have been made to identify the genetic basis of human longevity, with only limited progress. One important drawback of current genetic studies is the limited knowledge of gene-environment interaction. Using 2 cohorts of long-lived individuals born in 1905 and 1915 in Denmark, we performed survival analysis to estimate risk of mortality for major candidate genes of aging and longevity and their cohort effects. Through statistical modeling that combines individual genetic and survival information with cohort-specific survival data, we estimated the relative risks of mortality from ages 95 to 103 years associated with genetic variants in apolipoprotein E (APOE), forkhead box class O3a, clusterin, and phosphatidylinositol binding clathrin assembly protein. Our analysis estimated a decreased risk of carrying the APOE$\varepsilon $4 allele (change in risk = -0.403, 95% confidence interval (CI): -0.831, 0.021; P = 0.040) in men of the later cohort, although the allele itself was harmful to survival across sexes (relative risk = 1.161, 95% CI: 1.027, 1.345; P = 0.026). We also estimated a cohort effect of increased risk for the minor allele of rs3851179 in phosphatidylinositol binding clathrin assembly protein with borderline significance (change in risk = 0.165, 95% CI: -0.010, 0.331; P = 0.052) in women. Our estimated significant cohort effect on APOE$\varepsilon $4 is indicative of the interplay between the gene and the changing environment that modulates survival at extreme ages.

DNA methylome profiling in identical twin pairs discordant for body mass index.

OBJECTIVE: Body mass index (BMI) serves as an important measurement of obesity and adiposity, which are highly correlated with cardiometabolic diseases. Although high heritability has been estimated, the identified genetic variants by genetic association studies only explain a small proportion of BMI variation. As an active effort for further exploring the molecular basis of BMI variation, large-scale epigenome-wide association studies have been conducted but with limited number of loci reported, perhaps due to poorly controlled confounding factors, including genetic factors. Being genetically identical, monozygotic twins discordant for BMI are ideal subjects for analyzing the epigenetic association between DNA methylation and BMI, providing perfect control on their genetic makeups largely responsible for BMI variation. SUBJECTS: We performed an epigenome-wide association study on BMI using 30 identical twin pairs (15 male and 15 female pairs) with age ranging from 39 to 72 years and degree of BMI discordance ranging from 3-7.5 kg/m2. Methylation data from whole blood samples were collected using the reduced representation bisulfite sequencing technique. RESULTS: After adjusting for blood cell composition and clinical variables, we identified 136 CpGs with p-value < 1e-4, 30 CpGs with p < 1e-05 but no CpGs reached genome-wide significance. Genomic region-based analysis found 11 differentially methylated regions harboring coding and non-coding genes some of which were validated by gene expression analysis on independent samples. CONCLUSIONS: Our DNA methylation sequencing analysis on identical twins provides new references for the epigenetic regulation on BMI and obesity.

On the power of epigenome-wide association studies using a disease-discordant twin design.

Motivation: Many studies have investigated the association between DNA methylation alterations and disease occurrences using two design paradigms, traditional case-control and disease-discordant twins. In the disease-discordant twin design, the affected twin serves as the case and the unaffected twin serves as the control. Theoretically the twin design takes advantage of controlling for the shared genetic make-up, but it is still highly debatable if and how much researchers may benefit from such a design over the traditional case-control design. Results: In this study, we investigate and compare the power of both designs with simulations. A liability threshold model was used assuming that identical twins share the same genetic contribution with respect to the liability of complex human diseases. Varying ranges of parameters have been used to ensure that the simulation is close to real-world scenarios. Our results reveal that the disease-discordant twin design implies greater statistical power over the traditional case-control design. For diseases with moderate and high heritability (>0.3), the disease-discordant twin design allows for large sample size reductions compared to the ordinary case-control design. Our simulation results indicate that the discordant twin design is indeed a powerful tool for epigenetic association studies. Availability and implementation: Computer scripts are available at https://github.com/zickyls/EWAS-Twin-Simulation. Supplementary information: Supplementary data are available at Bioinformatics online.

White blood cell mitochondrial DNA copy number is decreased in rheumatoid arthritis and linked with risk factors. A twin study.

Low mitochondrial DNA copy number (mtDNA CN) has been associated with e.g. cancer, cardiovascular and autoimmune diseases. We aimed to study a potential association between mtDNA CN and rheumatoid arthritis (RA). The relative quantity of mitochondrial DNA compared to nuclear DNA was measured in peripheral white blood cells from 149 RA affected twin pairs and 1321 non-affected twin pairs. Multiple regression analysis including RA discordant twin pairs was performed in order to separate specific effects of RA and familial RA predisposition using non-RA affected twin pairs as reference group. In addition, we performed a twin pair level analysis including only RA discordant twin pairs evaluating the effect of cell type, auto antibodies and RA genetic risk factors. Both the RA twins and their non-affected co-twins had significantly lower mtDNA CN than non-affected twins (-28.7 and -23.1 mtDNA CN, respectively). Adjusting for cell count attenuated these differences (-23.1 mtDNA CN and -20.1 mtDNA CN respectively). Within RA discordant twin pairs PTPN22(T) positive RA twins had a significantly lower amount than their co-twins (-16.3 mtDNA CN). PTPN22(T) had no effect among twins from non-affected twin pairs. MtDNA CN is significantly lower in persons with established RA and in predisposed non-affected RA co-twins suggesting that mitochondrial variation may be involved in the RA disease pathways. Our results also suggest that the RA associated genetic risk factor, PTPN22(T), further decreases the mtDNA CN, but only in carriers with established RA.

Epigenome-wide exploratory study of monozygotic twins suggests differentially methylated regions to associate with hand grip strength.

Hand grip strength is a measure of muscular strength and is used to study age-related loss of physical capacity. In order to explore the biological mechanisms that influence hand grip strength variation, an epigenome-wide association study (EWAS) of hand grip strength in 672 middle-aged and elderly monozygotic twins (age 55-90 years) was performed, using both individual and twin pair level analyses, the latter controlling the influence of genetic variation. Moreover, as measurements of hand grip strength performed over 8 years were available in the elderly twins (age 73-90 at intake), a longitudinal EWAS was conducted for this subsample. No genome-wide significant CpG sites or pathways were found, however two of the suggestive top CpG sites were mapped to the COL6A1 and CACNA1B genes, known to be related to muscular dysfunction. By investigating genomic regions using the comb-p algorithm, several differentially methylated regions in regulatory domains were identified as significantly associated to hand grip strength, and pathway analyses of these regions revealed significant pathways related to the immune system, autoimmune disorders, including diabetes type 1 and viral myocarditis, as well as negative regulation of cell differentiation. The genes contributing to the immunological pathways were HLA-B, HLA-C, HLA-DMA, HLA-DPB1, MYH10, ERAP1 and IRF8, while the genes implicated in the negative regulation of cell differentiation were IRF8, CEBPD, ID2 and BRCA1. In conclusion, this exploratory study suggests hand grip strength to associate with differentially methylated regions enriched in immunological and cell differentiation pathways, and hence merits further investigations.

Surfactant protein-D, a potential mediator of inflammation in axial spondyloarthritis.

Objectives: Surfactant protein-D (SP-D), an innate immune defence molecule of the collectin family, is expressed in lungs and additional extrapulmonary epithelia. SP-D has immune modulatory and anti-microbial effects depending on its oligomerization. The ratio of high molecular weight (HMW): low molecular weight (LMW) SP-D in serum is mainly determined by the Met11Thr polymorphism (SNP rs721917). We aimed to study the SP-D serum level and the molecular size distribution in patients with untreated axial spondyloarthritis (axSpA) as compared with control subjects. Methods: Thirty-four patients with disease modifier untreated axSpA according to the ASAS criteria, age 19-63 years, disease duration 3.9 (2.2-5.6) years were included. Demographics, smoking habits, HLA-B27 status, ASDAS, BASDAI, BASFI, BASMI and visual analogue scale scores were recorded. SP-D in serum was measured by ELISA. DNA was isolated from whole blood and single nucleotide polymorphism rs721917 was genotyped. SP-D molecular size distribution was determined using gel filtration chromatography. Results: SP-D in serum did not differ between patients with axSpA and healthy controls, 1177 (869, 1536) vs 910 (494, 1682) (P = 0.35) and SP-D did not correlate with disease activity. However, the HMW/LMW ratio of SP-D in serum was significantly lower in axSpA, 0.38 (0.18, 0.53) compared with controls 1.49 (0.37, 3.24) when adjusting for the Met11Thr polymorphism, gender, age, BMI and smoking (P = 0.0004). There was no correlation between HMW/LMW ratio and CRP or composite diseases outcome measures. Conclusion: We suggest that predominance of LMW oligomeric variants of SP-D may enhance local or systemic inflammatory responses in axSpA.

Stratification by smoking status reveals an association of CHRNA5-A3-B4 genotype with body mass index in never smokers.

We previously used a single nucleotide polymorphism (SNP) in the CHRNA5-A3-B4 gene cluster associated with heaviness of smoking within smokers to confirm the causal effect of smoking in reducing body mass index (BMI) in a Mendelian randomisation analysis. While seeking to extend these findings in a larger sample we found that this SNP is associated with 0.74% lower body mass index (BMI) per minor allele in current smokers (95% CI -0.97 to -0.51, P = 2.00 × 10(-10)), but also unexpectedly found that it was associated with 0.35% higher BMI in never smokers (95% CI +0.18 to +0.52, P = 6.38 × 10(-5)). An interaction test confirmed that these estimates differed from each other (P = 4.95 × 10(-13)). This difference in effects suggests the variant influences BMI both via pathways unrelated to smoking, and via the weight-reducing effects of smoking. It would therefore be essentially undetectable in an unstratified genome-wide association study of BMI, given the opposite association with BMI in never and current smokers. This demonstrates that novel associations may be obscured by hidden population sub-structure. Stratification on well-characterized environmental factors known to impact on health outcomes may therefore reveal novel genetic associations.

Differentially Methylated Genomic Regions in Birth-Weight Discordant Twin Pairs.

Poor nutrition during critical growth phases may alter the structural and physiologic development of vital organs thus "programming" the susceptibility to adult-onset diseases and disease-related health conditions. Epigenome-wide association studies have been performed in birth-weight discordant twin pairs to find evidence for such "programming" effects, but no significant results emerged. We further investigated this issue using a new computational approach: Instead of probing single genomic sites for significant alterations in epigenetic marks, we scan for differentially methylated genomic regions. Whole genome DNA methylation levels were measured in whole blood from 150 pairs of adult identical twins discordant for birth-weight. Intrapair differential DNA methylation was associated with qualitative (large or small) and quantitative (percentage) birth-weight discordance at each genomic site using regression models adjusting for age and sex. Based on the regression results, genomic regions with consistent alteration patterns of DNA methylation were located and tested for significant robustness using computational permutation tests. This yielded an interesting genomic region on chromosome 1, which is significantly differentially methylated for quantitative birth-weight discordance. The region covers two genes (TYW3 and CRYZ) both reportedly associated with metabolism. We conclude that prenatal conditions for birth-weight discordance may result in persistent epigenetic modifications potentially affecting even adult health.

G×E Interaction Influences Trajectories of Hand Grip Strength.

Age-related decline in grip strength predicts later life disability, frailty, lower well-being and cognitive change. While grip strength is heritable, genetic influence on change in grip strength has been relatively ignored, with non-shared environmental influence identified as the primary contributor in a single longitudinal study. The extent to which gene-environment interplay, particularly gene-environment interactions, contributes to grip trajectories has yet to be examined. We considered longitudinal grip strength measurements in seven twin studies of aging in the Interplay of Genes and Environment across Multiple Studies consortium. Growth curve parameters were estimated for same-sex pairs, aged 34-99 (N = 10,681). Fisher's test for mixture distribution of within-monozygotic twin-pair differences (N = 1724) was performed on growth curve parameters. We observed significant gene-environment interaction on grip strength trajectories. Finally, we compared the variability of within-pair differences of growth curve parameters by APOE haplotypes. Though not statistically significant, the results suggested that APOE ɛ2ɛ2/ɛ2ɛ3 haplotypes might buffer environmental influences on grip strength trajectories.

Gene-Environment Interplay in Physical, Psychological, and Cognitive Domains in Mid to Late Adulthood: Is APOE a Variability Gene?

Despite emerging interest in gene-environment interaction (GxE) effects, there is a dearth of studies evaluating its potential relevance apart from specific hypothesized environments and biometrical variance trends. Using a monozygotic within-pair approach, we evaluated evidence of G×E for body mass index (BMI), depressive symptoms, and cognition (verbal, spatial, attention, working memory, perceptual speed) in twin studies from four countries. We also evaluated whether APOE is a 'variability gene' across these measures and whether it partly represents the 'G' in G×E effects. In all three domains, G×E effects were pervasive across country and gender, with small-to-moderate effects. Age-cohort trends were generally stable for BMI and depressive symptoms; however, they were variable-with both increasing and decreasing age-cohort trends-for different cognitive measures. Results also suggested that APOE may represent a 'variability gene' for depressive symptoms and spatial reasoning, but not for BMI or other cognitive measures. Hence, additional genes are salient beyond APOE.

DCAF4, a novel gene associated with leucocyte telomere length.

BACKGROUND: Leucocyte telomere length (LTL), which is fashioned by multiple genes, has been linked to a host of human diseases, including sporadic melanoma. A number of genes associated with LTL have already been identified through genome-wide association studies. The main aim of this study was to establish whether DCAF4 (DDB1 and CUL4-associated factor 4) is associated with LTL. In addition, using ingenuity pathway analysis (IPA), we examined whether LTL-associated genes in the general population might partially explain the inherently longer LTL in patients with sporadic melanoma, the risk for which is increased with ultraviolet radiation (UVR). RESULTS: Genome-wide association (GWA) meta-analysis and de novo genotyping of 20 022 individuals revealed a novel association (p=6.4×10(-10)) between LTL and rs2535913, which lies within DCAF4. Notably, eQTL analysis showed that rs2535913 is associated with decline in DCAF4 expressions in both lymphoblastoid cells and sun-exposed skin (p=4.1×10(-3) and 2×10(-3), respectively). Moreover, IPA revealed that LTL-associated genes, derived from GWA meta-analysis (N=9190), are over-represented among genes engaged in melanoma pathways. Meeting increasingly stringent p value thresholds (p<0.05, <0.01, <0.005, <0.001) in the LTL-GWA meta-analysis, these genes were jointly over-represented for melanoma at p values ranging from 1.97×10(-169) to 3.42×10(-24). CONCLUSIONS: We uncovered a new locus associated with LTL in the general population. We also provided preliminary findings that suggest a link of LTL through genetic mechanisms with UVR and melanoma in the general population.

Change in Depression Symptomatology and Cognitive Function in Twins: A 10-Year Follow-Up Study.

A complex interrelation exists between change in depression symptomatology and cognitive decline. Studies indicate either that depression is a direct risk factor for cognitive change over time, or vice versa. Longitudinal twin studies provide the possibility to unravel cause and effect of correlated traits. Here, we have applied twin modeling approaches to shed light on the genetic correlation between both level and change of depression symptomatology and cognitive functioning, and to further explore the bidirectionality of any such correlation using assessments of both phenotypes at two occasions 10 years apart. The study included 2,866 Danish twins with a mean age of 56.8 years at intake (range: 45-68 years). Of these, 1,267 were intact pairs. A total number of 1,582 twins (55%), of whom 557 were intact pairs, participated in the follow-up survey. We found stable cross-sectional heritability estimates of approximately 60% for general cognitive abilities and 30% for affective depressive symptoms. There was a considerable decline in the mean cognitive performance over 10 years, whereas the mean affective depression symptoms score was stable and with no genetic contribution to any individual change. Additionally, we saw a small but significant cross-trait correlation at both occasions (-0.11 and -0.09, respectively), but cross-trait cross-occasion analysis revealed no evidence that either of the two traits predicts the other over a 10-year interval. Thus, our study was not able to detect any causal association between change in depressive symptomatology and cognitive decline in middle-aged and elderly people over a 10-year interval.

A twin study of the trough plasma steady-state concentration of metformin.

OBJECTIVE: The aim of this study was to determine the intrapair similarity in trough steady-state plasma concentrations of metformin in monozygotic and dizygotic twin pairs. METHODS: We included 16 twin pairs (eight monozygotic and eight dizygotic twin pairs) for this study after contacting 524 twin pairs. They were dosed with metformin to steady state (1 g twice daily) for 6 days and on day 7, the trough concentration of metformin was determined 12 h after the last dose. RESULTS: There was no strong intrapair similarity in trough steady-state plasma concentrations of metformin in either dizygotic or monozygotic twin pairs. CONCLUSION: The trough steady-state plasma concentration of metformin does not appear to be tightly genetically regulated. The interpretation of this finding is limited by the small sample size.

Twin methodology in epigenetic studies.

Since the final decades of the last century, twin studies have made a remarkable contribution to the genetics of human complex traits and diseases. With the recent rapid development in modern biotechnology of high-throughput genetic and genomic analyses, twin modelling is expanding from analysis of diseases to molecular phenotypes in functional genomics especially in epigenetics, a thriving field of research that concerns the environmental regulation of gene expression through DNA methylation, histone modification, microRNA and long non-coding RNA expression, etc. The application of the twin method to molecular phenotypes offers new opportunities to study the genetic (nature) and environmental (nurture) contributions to epigenetic regulation of gene activity during developmental, ageing and disease processes. Besides the classical twin model, the case co-twin design using identical twins discordant for a trait or disease is becoming a popular and powerful design for epigenome-wide association study in linking environmental exposure to differential epigenetic regulation and to disease status while controlling for individual genetic make-up. It can be expected that novel uses of twin methods in epigenetic studies are going to help with efficiently unravelling the genetic and environmental basis of epigenomics in human complex diseases.

DNA Methylation Changes in the IGF1R Gene in Birth Weight Discordant Adult Monozygotic Twins.

Low birth weight (LBW) can have an impact on health outcomes in later life, especially in relation to pre-disposition to metabolic disease. Several studies suggest that LBW resulting from restricted intrauterine growth leaves a footprint on DNA methylation in utero, and this influence likely persists into adulthood. To investigate this further, we performed epigenome-wide association analyses of blood DNA methylation using Infinium HumanMethylation450 BeadChip profiles in 71 adult monozygotic (MZ) twin pairs who were extremely discordant for birth weight. A signal mapping to the IGF1R gene (cg12562232, p = 2.62 × 10(-8)), was significantly associated with birth weight discordance at a genome-wide false-discovery rate (FDR) of 0.05. We pursued replication in three additional independent datasets of birth weight discordant MZ pairs and observed the same direction of association, but the results were not significant. However, a meta-analysis across the four independent samples, in total 216 birth-weight discordant MZ twin pairs, showed a significant positive association between birth weight and DNA methylation differences at IGF1R (random-effects meta-analysis p = .04), and the effect was particularly pronounced in older twins (random-effects meta-analysis p = .008, 98 older birth-weight discordant MZ twin pairs). The results suggest that severe intra-uterine growth differences (birth weight discordance >20%) are associated with methylation changes in the IGF1R gene in adulthood, independent of genetic effects.

Epigenetic signature of birth weight discordance in adult twins.

BACKGROUND: A low birth weight has been extensively related to poor adult health outcomes. Birth weight can be seen as a proxy for environmental conditions during prenatal development. Identical twin pairs discordant for birth weight provide an extraordinary model for investigating the association between birth weight and adult life health while controlling for not only genetics but also postnatal rearing environment. We performed an epigenome-wide profiling on blood samples from 150 pairs of adult monozygotic twins discordant for birth weight to look for molecular evidence of epigenetic signatures in association with birth weight discordance. RESULTS: Our association analysis revealed no CpG site with genome-wide statistical significance (FDR < 0.05) for either qualitative (larger or smaller) or quantitative discordance in birth weight. Even with selected samples of extremely birth weight discordant twin pairs, no significant site was found except for 3 CpGs that displayed age-dependent intra-pair differential methylation with FDRs 0.014 (cg26856578, p = 3.42e-08), 0.0256 (cg15122603, p = 1.25e-07) and 0.0258 (cg16636641, p = 2.05e-07). Among the three sites, intra-pair differential methylation increased with age for cg26856578 but decreased with age for cg15122603 and cg16636641. There was no genome-wide statistical significance for sex-dependent effects on intra-pair differential methylation in either the whole samples or the extremely discordant twins. CONCLUSIONS: Genome-wide DNA methylation profiling did not reveal epigenetic signatures of birth weight discordance although some sites displayed age-dependent intra-pair differential methylation in the extremely discordant twin pairs.

Surfactant protein D is a candidate biomarker for subclinical tobacco smoke-induced lung damage.

Variation in surfactant protein D (SP-D) is associated with lung function in tobacco smoke-induced chronic respiratory disease. We hypothesized that the same association exists in the general population and could be used to identify individuals sensitive to smoke-induced lung damage. The association between serum SP-D (sSP-D) and expiratory lung function was assessed in a cross-sectional design in a Danish twin population (n = 1,512, 18-72 yr old). The adjusted heritability estimates for expiratory lung function, associations between SP-D gene (SFTPD) single-nucleotide polymorphisms or haplotypes, and expiratory lung function were assessed using twin study methodology and mixed-effects models. Significant inverse associations were evident between sSP-D and the forced expiratory volume in 1 s and forced vital capacity in the presence of current tobacco smoking but not in nonsmokers. The two SFTPD single-nucleotide polymorphisms, rs1923536 and rs721917, and haplotypes, including these single-nucleotide polymorphisms or rs2243539, were inversely associated with expiratory lung function in interaction with smoking. In conclusion, SP-D is phenotypically and genetically associated with lung function measures in interaction with tobacco smoking. The obtained data suggest sSP-D as a candidate biomarker in risk assessments for subclinical tobacco smoke-induced lung damage. The data and derived conclusion warrant confirmation in a longitudinal population following chronic obstructive pulmonary disease initiation and development.

Genome-wide meta-analysis points to CTC1 and ZNF676 as genes regulating telomere homeostasis in humans.

Leukocyte telomere length (LTL) is associated with a number of common age-related diseases and is a heritable trait. Previous genome-wide association studies (GWASs) identified two loci on chromosomes 3q26.2 (TERC) and 10q24.33 (OBFC1) that are associated with the inter-individual LTL variation. We performed a meta-analysis of 9190 individuals from six independent GWAS and validated our findings in 2226 individuals from four additional studies. We confirmed previously reported associations with OBFC1 (rs9419958 P = 9.1 × 10(-11)) and with the telomerase RNA component TERC (rs1317082, P = 1.1 × 10(-8)). We also identified two novel genomic regions associated with LTL variation that map near a conserved telomere maintenance complex component 1 (CTC1; rs3027234, P = 3.6 × 10(-8)) on chromosome17p13.1 and zinc finger protein 676 (ZNF676; rs412658, P = 3.3 × 10(-8)) on 19p12. The minor allele of rs3027234 was associated with both shorter LTL and lower expression of CTC1. Our findings are consistent with the recent observations that point mutations in CTC1 cause short telomeres in both Arabidopsis and humans affected by a rare Mendelian syndrome. Overall, our results provide novel insights into the genetic architecture of inter-individual LTL variation in the general population.