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.

OBJECTIVE

Vitamin D deficiency (25-hydroxyvitamin D [25(OH)D] <50 nmol/L) is commonly reported in both children and adults worldwide, and growing evidence indicates that vitamin D deficiency is associated with many extraskeletal chronic disorders, including the autoimmune diseases type 1 diabetes and multiple sclerosis.

RESEARCH DESIGN AND METHODS

We measured 25(OH)D concentrations in 720 case and 2,610 control plasma samples and genotyped single nucleotide polymorphisms from seven vitamin D metabolism genes in 8,517 case, 10,438 control, and 1,933 family samples. We tested genetic variants influencing 25(OH)D metabolism for an association with both circulating 25(OH)D concentrations and disease status.

RESULTS

Type 1 diabetic patients have lower circulating levels of 25(OH)D than similarly aged subjects from the British population. Only 4.3 and 18.6% of type 1 diabetic patients reached optimal levels (≥75 nmol/L) of 25(OH)D for bone health in the winter and summer, respectively. We replicated the associations of four vitamin D metabolism genes (GC, DHCR7, CYP2R1, and CYP24A1) with 25(OH)D in control subjects. In addition to the previously reported association between type 1 diabetes and CYP27B1 (P = 1.4 × 10−4), we obtained consistent evidence of type 1 diabetes being associated with DHCR7 (P = 1.2 × 10−3) and CYP2R1 (P = 3.0 × 10−3).

CONCLUSIONS

Circulating levels of 25(OH)D in children and adolescents with type 1 diabetes vary seasonally and are under the same genetic control as in the general population but are much lower. Three key 25(OH)D metabolism genes show consistent evidence of association with type 1 diabetes risk, indicating a genetic etiological role for vitamin D deficiency in type 1 diabetes.

Genome-wide association studies have identified many genetic variants associated with complex traits. However, at only a minority of loci have the molecular mechanisms mediating these associations been characterized. In parallel, whilst cis-regulatory patterns of gene expression have been extensively explored, the identification of trans-regulatory effects in humans has attracted less attention. We demonstrate that the Type 2 diabetes and HDL-cholesterol associated cis-acting eQTL of the maternally-expressed transcription factor KLF14 acts as a master trans-regulator of adipose gene expression. Expression levels of genes regulated by this trans-eQTL are highly-correlated with concurrently-measured metabolic traits, and a subset of the trans-genes harbor variants directly-associated with metabolic phenotypes. This trans-eQTL network provides a mechanistic understanding of the effect of the KLF14 locus on metabolic disease risk, providing a potential model for other complex traits.

The corrected QT (QTc) interval is a complex quantitative trait, believed to be influenced by several genetic and environmental factors. It is a strong prognostic indicator of cardiovascular mortality in patients with and without cardiac disease. More than 700 mutations have been described in 12 genes (LQT1-LQT12) involved in congenital long QT syndrome. However, the heritability (genetic contribution) of QTc interval in the general population cannot be adequately explained by these long QT syndrome genes. In order to further investigate the genetic architecture underlying QTc interval in the general population, genome-wide association studies, in which up to one million single nucleotide polymorphisms are assayed in thousands of individuals, are now being employed and have already led to the discovery of variants in seven novel loci and five loci that are known to cause congenital long or short QT syndrome. Here we show that a combined risk score using 11 of these loci explains about 10% of the heritability of QTc. Additional discovery of both common and rare variants will yield further etiological insight and accelerate clinical applications.

Imputation is an extremely valuable tool in conducting and synthesising genome-wide association studies (GWASs). Directly typed SNP quality control (QC) is thought to affect imputation quality. It is, therefore, common practise to use quality-controlled (QCed) data as an input for imputing genotypes. This study aims to determine the effect of commonly applied QC steps on imputation outcomes. We performed several iterations of imputing SNPs across chromosome 22 in a dataset consisting of 3177 samples with Illumina 610k (Illumina, San Diego, CA, USA) GWAS data, applying different QC steps each time. The imputed genotypes were compared with the directly typed genotypes. In addition, we investigated the correlation between alternatively QCed data. We also applied a series of post-imputation QC steps balancing elimination of poorly imputed SNPs and information loss. We found that the difference between the unQCed data and the fully QCed data on imputation outcome was minimal. Our study shows that imputation of common variants is generally very accurate and robust to GWAS QC, which is not a major factor affecting imputation outcome. A minority of common-frequency SNPs with particular properties cannot be accurately imputed regardless of QC stringency. These findings may not generalise to the imputation of low frequency and rare variants.

Background

Telomeres shorten as cells divide. This shortening is compensated by the enzyme telomerase. We evaluated the effect of common variants in the telomerase RNA component (TERC) gene on telomere length (TL) in the population-based Health Aging and Body Composition (Health ABC) Study and in two replication samples (the TwinsUK Study and the Amish Family Osteoporosis Study, AFOS).

Methodology

Five variants were identified in the TERC region by sequence analysis and only one SNP was common (rs2293607, G/A). The frequency of the G allele was 0.26 and 0.07 in white and black, respectively. Testing for association between TL and rs2293607 was performed using linear regression models or variance component analysis conditioning on relatedness among subjects.

Results

The adjusted mean TL was significantly shorter in 665 white carriers of the G allele compared to 887 non-carriers from the Health ABC Study (4.69±0.05 kbp vs. 4.86±0.04 kbp, measured by quantitative PCR, p = 0.005). This association was replicated in another white sample from the TwinsUK Study (6.90±0.03 kbp in 301 carriers compared to 7.06±0.03 kbp in 395 non-carriers, measured by Southern blots, p = 0.009). A similar pattern of association was observed in whites from the family-based AFOS and blacks from the Health ABC cohort, although not statistically significant, possibly due to the lower allele frequency in these populations. Combined analysis using 2,953 white subjects from 3 studies showed a significant association between TL and rs2293607 (β = −0.19±0.04 kbp, p = 0.001).

Conclusion

Our study shows a significant association between a common variant in TERC and TL in humans, suggesting that TERC may play a role in telomere homeostasis.

Objective

Osteoarthritis (OA) and osteoporosis are often considered to lie at opposite ends of a spectrum of bone phenotypes. Lumbar degenerative disc disease (LDD) may be associated with low back pain (LBP) and is similar in many ways to OA. LDD is reported in small studies to be associated with increased spine bone mineral density (BMD). The present work aimed to confirm this association in a large population sample using MRI and explore the relationship further, in particular to determine whether it is mediated genetically.

Methods

A population based sample (N=908, age range 32–74 years) of UK female twins having MRI of the lumbar spine was used in this study. LDD traits and summary measures and their relationship with BMD at the lumbar spine and hip were examined using multivariate multiple regression and maximum likelihood based variance decomposition.

Results

There was a significant positive correlation between LDD and BMD at the lumbar spine and hip, which remained significant after adjustment for confounders. Both traits were highly heritable and the associations between them were mediated genetically.

Conclusions

A clear, significant and independent association of BMD at hip and lumbar spine with LDD was found which is, in part, genetically mediated. The association with the non-axial site, the hip, is of particular interest and suggests a systemic bone effect. This should encourage the search for pleiotropic genes to help in the understanding of the bone–cartilage relationship. Moreover, genetic variants identified could provide novel therapeutic targets in the management of LBP.

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

Author Summary

We measured human eye color to hue and saturation values from high-resolution, digital, full-eye photographs of several thousand Dutch Europeans. This quantitative approach, which is extremely cost-effective, portable, and time efficient, revealed that human eye color varies along more dimensions than the one represented by the blue-green-brown categories studied previously. Our work represents the first genome-wide study of quantitative human eye color. We clearly identified 3 new loci, LYST, 17q25.3, TTC3/DSCR9, in contributing to the natural and subtle eye color variation along multiple dimensions, providing new leads towards a more detailed understanding of the genetic basis of human eye color. Our quantitative prediction model explained over 50% of eye color variance, representing the highest accuracy achieved so far in genomic prediction of human complex and quantitative traits, with relevance for future forensic applications.

Summary

The association of FTO gene variants with body mass index (BMI) and other obesity characteristics is well established. However, uncertainties remain whether the association is present only in young populations and whether it is attributable to body fat mass specifically. We aimed to clarify these two questions in a large sample (N= 4,523 individuals) of middle-aged and older (range 40–80 years) British female twins. The women were assessed for BMI, waist and hip circumference, total lean (LBM) and fat (FBM) body mass. Since the majority of FTO association signals have been reported in a haploblock bordering 52,355–52,408 kb (on chromosome 16q12.2), we examined five genotyped and 43 imputed SNPs mapped to this block. Canonical correlation and other association analyses showed significant and consistent association between the selected SNP and studied body composition phenotypes, with p-values reaching p= 0.000004. Of particular interest, in addition to the expected significant associations between FTO variants and FBM, we also identified significant associations with LBM. These results suggest that the association between FTO variants and body composition phenotypes is present across a wide range of ages, and that FTO appears primarily to affect the amount of body soft tissue, influencing both fat and lean mass.

Background

Strontium ranelate is currently used for osteoporosis. The international, double-blind, randomised, placebo-controlled Strontium ranelate Efficacy in Knee OsteoarthrItis triAl evaluated its effect on radiological progression of knee osteoarthritis.

Methods

Patients with knee osteoarthritis (Kellgren and Lawrence grade 2 or 3, and joint space width (JSW) 2.5–5 mm) were randomly allocated to strontium ranelate 1 g/day (n=558), 2 g/day (n=566) or placebo (n=559). The primary endpoint was radiographical change in JSW (medial tibiofemoral compartment) over 3 years versus placebo. Secondary endpoints included radiological progression, Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) score, and knee pain. The trial is registered (ISRCTN41323372).

Results

The intention-to-treat population included 1371 patients. Treatment with strontium ranelate was associated with smaller degradations in JSW than placebo (1 g/day: −0.23 (SD 0.56) mm; 2 g/day: −0.27 (SD 0.63) mm; placebo: −0.37 (SD 0.59) mm); treatment-placebo differences were 0.14 (SE 0.04), 95% CI 0.05 to 0.23, p<0.001 for 1 g/day and 0.10 (SE 0.04), 95% CI 0.02 to 0.19, p=0.018 for 2 g/day. Fewer radiological progressors were observed with strontium ranelate (p<0.001 and p=0.012 for 1 and 2 g/day). There were greater reductions in total WOMAC score (p=0.045), pain subscore (p=0.028), physical function subscore (p=0.099) and knee pain (p=0.065) with strontium ranelate 2 g/day. Strontium ranelate was well tolerated.

Conclusions

Treatment with strontium ranelate 1 and 2 g/day is associated with a significant effect on structure in patients with knee osteoarthritis, and a beneficial effect on symptoms for strontium ranelate 2 g/day.

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Author Summary

Monozygotic twins look more alike than dizygotic twins or other siblings, and siblings in turn look more alike than unrelated individuals, indicating that human facial morphology has a strong genetic component. We quantitatively assessed human facial shape phenotypes based on statistical shape analyses of facial landmarks obtained from three-dimensional magnetic resonance images of the head. These phenotypes turned out to be highly promising for studying the genetic basis of human facial variation in that they showed high heritability in our twin data. A subsequent genome-wide association study (GWAS) identified five candidate genes affecting facial shape in Europeans: PRDM16, PAX3, TP63, C5orf50, and COL17A1. In addition, our data suggest that genetic variants associated with NSCL/P also influence normal facial shape variation. Overall, this study provides novel and confirmatory links between common DNA variants and normal variation in human facial morphology. Our results also suggest that the high heritability of facial phenotypes seems to be explained by a large number of DNA variants with relatively small individual effect size, a phenomenon well known for other complex human traits, such as adult body height.

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.

Author Summary

Genetic information is transmitted to the cell only through RNA molecules. A special class of RNAs is comprised of the small (up to 30 nucleotide) ones, known to be potent regulators of various cellular processes. At the same time, they have not been as widely studied as messenger RNAs—we do not know how much variation in their sequence and expression level occurs naturally in human populations or how this variability influences other traits. We measured small RNA levels and genetic variability in fat tissue from 131 individuals by high-throughput sequencing. We could associate the expression levels with genetic background of the individuals, as well as changes in metabolic traits. Surprisingly, we found no large scale influence of small RNA variation on mRNA levels, their main regulatory target. Overall, our study is the first to give a quantitative picture of the naturally occurring variation in these important regulatory molecules in human fat tissue.

Background

Type 2 diabetes mellitus (T2DM) has been linked to a state of pre-clinical chronic inflammation resulting from abnormalities in the innate immune pathway. Serum levels of pro-inflammatory cytokines and acute-phase proteins, collectively known as 'inflammatory network', are elevated in the pre-, or early, stages of T2DM and increase with disease progression. Genetic variation can affect the innate immune response to certain environmental factors, and may, therefore, determine an individual's lifetime risk of disease.

Methods

We conducted a cross-sectional study in 6,720 subjects from the TwinsUK Registry to evaluate the association between 18 single nucleotide polymorphisms (SNPs) in five genes (TLR4, IL1A, IL6, TNFA, and CRP) along the innate immunity-related inflammatory pathway and biomarkers of predisposition to T2DM [fasting insulin and glucose, HDL- and LDL- cholesterols, triglycerides (TGs), amyloid-A, sensitive C-reactive protein (sCRP) and vitamin D binding protein (VDBP) and body mass index (BMI)].

Results

Of 18 the SNPs examined for their association with nine metabolic phenotypes of interest, six were significantly associated with five metabolic phenotypes (Bonferroni correction, P ≤ 0.0027). Fasting insulin was associated with SNPs in IL6 and TNFA, serum HDL-C with variants of TNFA and CRP and serum sCRP level with SNPs in CRP. Cross-correlation analysis among the different metabolic factors related to risk of T2DM showed several significant associations. For example, BMI was directly correlated with glucose (r = 0.11), insulin (r = 0.15), sCRP (r = 0.23), LDL-C (r = 0.067) and TGs (r = 0.18) but inversely with HDL-C (r = -0.14). sCRP was also positively correlated (P < 0.0001) with insulin (r = 0.17), amyloid-A (r = 0.39), TGs (r = 0.26), and VDBP (r = 0.36) but inversely with HDL-C (r = -0.12).

Conclusion

Genetic variants in the innate immunity pathway and its related inflammatory cascade is associated with some metabolic risk factors for T2DM; an observation that may provide a rationale for further studying their role as biomarkers for disease early risk prediction.

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Author Summary

Regulation of gene expression is a fundamental cellular process determining a large proportion of the phenotypic variance. Previous studies have identified genetic loci influencing gene expression levels (eQTLs), but the complexity of their tissue-specific properties has not yet been well-characterized. In this study, we perform cis-eQTL analysis in a unique matched co-twin design for three human tissues derived simultaneously from the same set of individuals. The study design allows validation of the substantial discoveries we make in each tissue. We explore in depth the tissue-dependent features of regulatory variants and estimate the proportions of shared and specific effects. We use continuous measures of eQTL sharing to circumvent the statistical power limitations of comparing direct overlap of eQTLs in multiple tissues. In this framework, we demonstrate that 30% of eQTLs are shared among tissues, while 29% are exclusively tissue-specific. Furthermore, we show that the fold change in expression between eQTL genotypic classes differs between tissues. Even among shared eQTLs, we report a substantial proportion (10%–20%) of significant tissue differences in magnitude of these effects. The complexities we highlight here are essential for understanding the impact of regulatory variants on complex traits.

Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.6×10−8), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6×10−13; SOX6, p = 6.4×10−10) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.

Author Summary

BMD and hip geometry are two major predictors of osteoporotic fractures, the most severe consequence of osteoporosis in elderly persons. We performed sex-specific genome-wide association studies (GWAS) for BMD at the lumbar spine and femor neck skeletal sites as well as hip geometric indices (NSA, NL, and NW) in the Framingham Osteoporosis Study and then replicated the top findings in two independent studies. Three novel loci were significant: in women, including chromosome 1p13.2 (RAP1A) for NW; in men, 2q11.2 (TBC1D8) for NSA and 18q11.2 (OSBPL1A) for NW. We confirmed a previously reported region on 8q24.12 (TNFRSF11B/OPG) for lumbar spine BMD in women. In addition, we integrated GWAS signals with eQTL in several tissues and publicly available expression signature profiling in cellular and whole-animal models, and prioritized 16 candidate genes/loci based on their potential involvement in skeletal metabolism. Among three prioritized loci (GPR177, SOX6, and CASR genes) associated with BMD in women, GPR177 and SOX6 have been successfully replicated later in a large-scale meta-analysis, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of using expression profiling to support the candidacy of suggestive GWAS signals that may contain important genes of interest.

Central corneal thickness (CCT), one of the most highly heritable human traits (h2 typically>0.9), is important for the diagnosis of glaucoma and a potential risk factor for glaucoma susceptibility. We conducted genome-wide association studies in five cohorts from Australia and the United Kingdom (total N = 5058). Three cohorts were based on individually genotyped twin collections, with the remaining two cohorts genotyped on pooled samples from singletons with extreme trait values. The pooled sample findings were validated by individual genotyping the pooled samples together with additional samples also within extreme quantiles. We describe methods for efficient combined analysis of the results from these different study designs. We have identified and replicated quantitative trait loci on chromosomes 13 and 16 for association with CCT. The locus on chromosome 13 (nearest gene FOXO1) had an overall meta-analysis p-value for all the individually genotyped samples of 4.6×10−10. The locus on chromosome 16 was associated with CCT with p = 8.95×10−11. The nearest gene to the associated chromosome 16 SNPs was ZNF469, a locus recently implicated in Brittle Cornea Syndrome (BCS), a very rare disorder characterized by abnormal thin corneas. Our findings suggest that in addition to rare variants in ZNF469 underlying CCT variation in BCS patients, more common variants near this gene may contribute to CCT variation in the general population.

Author Summary

Central corneal thickness (CCT) is an important eye measurement. It has been considered as a prognosticator for the development of glaucoma, with a thin cornea potentially increasing the risk of developing a subtype known as open-angle glaucoma. CCT is highly heritable, yet its genetic determinants are poorly characterized. We have revealed two loci near gene FOXO1 and ZNF469 associated with CCT in this multi-stage genome-wide association study examining over 5,000 samples. It is of particular interest that, while rare mutations in ZNF469 cause Brittle Cornea Syndrome, more common variants near this gene also contribute to CCT variation in the general population. Furthermore, given the relation between CCT and glaucoma, results from our CCT studies will implement the search for the disease-susceptibility genes of glaucoma.

Optic nerve assessment is important for many blinding diseases, with cup-to-disc ratio (CDR) assessments commonly used in both diagnosis and progression monitoring of glaucoma patients. Optic disc, cup, rim area and CDR measurements all show substantial variation between human populations and high heritability estimates within populations. To identify loci underlying these quantitative traits, we performed a genome-wide association study in two Australian twin cohorts and identified rs3858145, P = 6.2 × 10−10, near the ATOH7 gene as associated with the mean disc area. ATOH7 is known from studies in model organisms to play a key role in retinal ganglion cell formation. The association with rs3858145 was replicated in a cohort of UK twins, with a meta-analysis of the combined data yielding P = 3.4 × 10−10. Imputation further increased the evidence for association for several SNPs in and around ATOH7 (P = 1.3 × 10−10 to 4.3 × 10−11, top SNP rs1900004). The meta-analysis also provided suggestive evidence for association for the cup area at rs690037, P = 1.5 × 10−7, in the gene RFTN1. Direct sequencing of ATOH7 in 12 patients with optic nerve hypoplasia, one of the leading causes of blindness in children, revealed two novel non-synonymous mutations (Arg65Gly, Ala47Thr) which were not found in 90 unrelated controls (combined Fisher's exact P = 0.0136). Furthermore, the Arg65Gly variant was found to have very low frequency (0.00066) in an additional set of 672 controls.

Recent genome-wide (GW) scans have identified several independent loci affecting human stature, but their contribution through the different skeletal components of height is still poorly understood. We carried out a genome-wide scan in 12,611 participants, followed by replication in an additional 7,187 individuals, and identified 17 genomic regions with GW-significant association with height. Of these, two are entirely novel (rs11809207 in CATSPER4, combined P-value = 6.1×10−8 and rs910316 in TMED10, P-value = 1.4×10−7) and two had previously been described with weak statistical support (rs10472828 in NPR3, P-value = 3×10−7 and rs849141 in JAZF1, P-value = 3.2×10−11). One locus (rs1182188 at GNA12) identifies the first height eQTL. We also assessed the contribution of height loci to the upper- (trunk) and lower-body (hip axis and femur) skeletal components of height. We find evidence for several loci associated with trunk length (including rs6570507 in GPR126, P-value = 4×10−5 and rs6817306 in LCORL, P-value = 4×10−4), hip axis length (including rs6830062 at LCORL, P-value = 4.8×10−4 and rs4911494 at UQCC, P-value = 1.9×10−4), and femur length (including rs710841 at PRKG2, P-value = 2.4×10−5 and rs10946808 at HIST1H1D, P-value = 6.4×10−6). Finally, we used conditional analyses to explore a possible differential contribution of the height loci to these different skeletal size measurements. In addition to validating four novel loci controlling adult stature, our study represents the first effort to assess the contribution of genetic loci to three skeletal components of height. Further statistical tests in larger numbers of individuals will be required to verify if the height loci affect height preferentially through these subcomponents of height.

Author Summary

The first genetic association studies of adult height have confirmed a role of many common variants in influencing human height, but to date, the genetic basis of differences between different skeletal components of height have not been addressed. Here, we take advantage of recent technical and methodological advances to examine the role of common genetic variants on both height and skeletal components of height. By examining nearly 20,000 individuals from the UK and the Netherlands, we provide statistically significant evidence that 17 genomic regions are associated with height, including four novel regions. We also examine, for the first time, the association of these 17 regions with skeletal size measurements of spine, femur, and hip axis length, a measurement of hip geometry known to influence the risk of osteoporotic fractures. We find that some height loci are also associated with these skeletal components, although further statistical tests will be required to verify if these genetic variants act differentially on the individual skeletal measurements. The knowledge generated by this and other studies will not only inform the genetics of human quantitative variation, but will also lead to the potential discovery of many medically important polymorphisms.