Genome-wide association study using family-based cohorts identifies the WLS and CCDC170/ESR1 loci as associated with bone mineral density.

BACKGROUND: Osteoporosis is a common and debilitating bone disease that is characterised by a low bone mineral density (BMD), a highly heritable trait. Genome-wide association studies (GWAS) have proven to be very successful in identifying common genetic variants associated with BMD adjusted for age, gender and weight, however a large portion of the genetic variance for this trait remains unexplained. There is evidence to suggest significant genetic correlation between body size traits and BMD. It has also recently been suggested that unintended bias can be introduced as a result of adjusting a phenotype for a correlated trait. We performed a GWAS meta-analysis in two populations (total n = 6,696) using BMD data adjusted for only age and gender, in an attempt to identify genetic variants associated with BMD including those that may have potential pleiotropic effects on BMD and body size traits. RESULTS: We observed a single variant, rs2566752, associated with spine BMD at the genome-wide significance level in the meta-analysis (P = 3.36 × 10(-09)). Logistic regression analysis also revealed an association between rs2566752 and fracture rate in one of our study cohorts (P = 0.017, n = 5,654). This is an intronic variant located in the wntless Wnt ligand secretion mediator (WLS) gene (1p31.3), a known BMD locus which encodes an integral component of the Wnt ligand secretion pathway. Bioinformatics analyses of variants in moderate LD with rs2566752 produced strong evidence for a regulatory role for the variants rs72670452, rs17130567 and rs1430738. Expression quantitative trait locus (eQTL) analysis suggested that the variants rs12568456 and rs17130567 are associated with expression of the WLS gene in whole blood, cerebellum and temporal cortex brain tissue (P = 0.034-1.19 × 10(-23)). Gene-wide association testing using the VErsatile Gene-based Association Study 2 (VEGAS2) software revealed associations between the coiled-coil domain containing 170 (CCDC170) gene, located adjacent to the oestrogen receptor 1 (ESR1) gene, and BMD at the spine, femoral neck and total hip sites (P = 1.0 × 10(-06), 2.0 × 10(-06) and 2.0 × 10(-06) respectively). CONCLUSIONS: Genetic variation at the WLS and CCDC170/ESR1 loci were found to be significantly associated with BMD adjusted for only age and gender at the genome-wide level in this meta-analysis.

Genome-wide association study with 1000 genomes imputation identifies signals for nine sex hormone-related phenotypes.

Genetic factors contribute strongly to sex hormone levels, yet knowledge of the regulatory mechanisms remains incomplete. Genome-wide association studies (GWAS) have identified only a small number of loci associated with sex hormone levels, with several reproductive hormones yet to be assessed. The aim of the study was to identify novel genetic variants contributing to the regulation of sex hormones. We performed GWAS using genotypes imputed from the 1000 Genomes reference panel. The study used genotype and phenotype data from a UK twin register. We included 2913 individuals (up to 294 males) from the Twins UK study, excluding individuals receiving hormone treatment. Phenotypes were standardised for age, sex, BMI, stage of menstrual cycle and menopausal status. We tested 7,879,351 autosomal SNPs for association with levels of dehydroepiandrosterone sulphate (DHEAS), oestradiol, free androgen index (FAI), follicle-stimulating hormone (FSH), luteinizing hormone (LH), prolactin, progesterone, sex hormone-binding globulin and testosterone. Eight independent genetic variants reached genome-wide significance (P<5 × 10(-8)), with minor allele frequencies of 1.3-23.9%. Novel signals included variants for progesterone (P=7.68 × 10(-12)), oestradiol (P=1.63 × 10(-8)) and FAI (P=1.50 × 10(-8)). A genetic variant near the FSHB gene was identified which influenced both FSH (P=1.74 × 10(-8)) and LH (P=3.94 × 10(-9)) levels. A separate locus on chromosome 7 was associated with both DHEAS (P=1.82 × 10(-14)) and progesterone (P=6.09 × 10(-14)). This study highlights loci that are relevant to reproductive function and suggests overlap in the genetic basis of hormone regulation.

Whole-genome sequence-based analysis of thyroid function.

Normal thyroid function is essential for health, but its genetic architecture remains poorly understood. Here, for the heritable thyroid traits thyrotropin (TSH) and free thyroxine (FT4), we analyse whole-genome sequence data from the UK10K project (N=2,287). Using additional whole-genome sequence and deeply imputed data sets, we report meta-analysis results for common variants (MAF≥1%) associated with TSH and FT4 (N=16,335). For TSH, we identify a novel variant in SYN2 (MAF=23.5%, P=6.15 × 10(-9)) and a new independent variant in PDE8B (MAF=10.4%, P=5.94 × 10(-14)). For FT4, we report a low-frequency variant near B4GALT6/SLC25A52 (MAF=3.2%, P=1.27 × 10(-9)) tagging a rare TTR variant (MAF=0.4%, P=2.14 × 10(-11)). All common variants explain ≥20% of the variance in TSH and FT4. Analysis of rare variants (MAF<1%) using sequence kernel association testing reveals a novel association with FT4 in NRG1. Our results demonstrate that increased coverage in whole-genome sequence association studies identifies novel variants associated with thyroid function.

Copy number variation of the APC gene is associated with regulation of bone mineral density.

INTRODUCTION: Genetic studies of osteoporosis have commonly examined SNPs in candidate genes or whole genome analyses, but insertions and deletions of DNA, collectively called copy number variations (CNVs), also comprise a large amount of the genetic variability between individuals. Previously, SNPs in the APC gene have been strongly associated with femoral neck and lumbar spine volumetric bone mineral density in older men. In addition, familial adenomatous polyposis patients carrying heterozygous mutations in the APC gene have been shown to have significantly higher mean bone mineral density than age- and sex-matched controls suggesting the importance of this gene in regulating bone mineral density. We examined CNV within the APC gene region to test for association with bone mineral density. METHODS: DNA was extracted from venous blood, genotyped using the Human Hap610 arrays and CNV determined from the fluorescence intensity data in 2070 Caucasian men and women aged 47.0 ± 13.0 (mean ± SD) years, to assess the effects of the CNV on bone mineral density at the forearm, spine and total hip sites. RESULTS: Data for covariate adjusted bone mineral density from subjects grouped by APC CNV genotype showed significant difference (P=0.02-0.002). Subjects with a single copy loss of APC had a 7.95%, 13.10% and 13.36% increase in bone mineral density at the forearm, spine and total hip sites respectively, compared to subjects with two copies of the APC gene. CONCLUSIONS: These data support previous findings of APC regulating bone mineral density and demonstrate that a novel CNV of the APC gene is significantly associated with bone mineral density in Caucasian men and women.

The use of genome-wide eQTL associations in lymphoblastoid cell lines to identify novel genetic pathways involved in complex traits.

The integrated analysis of genotypic and expression data for association with complex traits could identify novel genetic pathways involved in complex traits. We profiled 19,573 expression probes in Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs) from 299 twins and correlated these with 44 quantitative traits (QTs). For 939 expressed probes correlating with more than one QT, we investigated the presence of eQTL associations in three datasets of 57 CEU HapMap founders and 86 unrelated twins. Genome-wide association analysis of these probes with 2.2 m SNPs revealed 131 potential eQTLs (1,989 eQTL SNPs) overlapping between the HapMap datasets, five of which were in cis (58 eQTL SNPs). We then tested 535 SNPs tagging the eQTL SNPs, for association with the relevant QT in 2,905 twins. We identified nine potential SNP-QT associations (P<0.01) but none significantly replicated in five large consortia of 1,097-16,129 subjects. We also failed to replicate previous reported eQTL associations with body mass index, plasma low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides levels derived from lymphocytes, adipose and liver tissue. Our results and additional power calculations suggest that proponents may have been overoptimistic in the power of LCLs in eQTL approaches to elucidate regulatory genetic effects on complex traits using the small datasets generated to date. Nevertheless, larger tissue-specific expression data sets relevant to specific traits are becoming available, and should enable the adoption of similar integrated analyses in the near future.

A locus on chromosome 1p36 is associated with thyrotropin and thyroid function as identified by genome-wide association study.

Thyroid hormones are key regulators of cellular growth, development, and metabolism, and thyroid disorders are a common cause of ill health in the community. Circulating concentrations of thyrotropin (TSH), thyroxine (T4) and triiodothyronine (T3) have a strong heritable component and are thought to be under polygenic control, but the genes responsible are mostly unknown. In order to identify genetic loci associated with these metabolic phenotypes, we performed a genome-wide association study of 2,120,505 SNPs in 2014 female twins from the TwinsUK study and found a significant association between rs10917469 on chromosome 1p36.13 and serum TSH (p = 3.2 × 10(-8)). The association of rs10917469 with serum TSH was replicated (p = 2.0 × 10(-4)) in an independent community-based sample of 1154 participants in the Busselton Health Study. This SNP is located near CAPZB, which might be a regulator of TSH secretion and thus of pituitary-thyroid axis function. Twenty-nine percent of white individuals carry the variant, and the difference in mean TSH concentrations between wild-type individuals and those homozygous for the minor G allele was 0.5 mU/l, which is likely to be clinically relevant. We also provide evidence of suggestive association (p < 5.0 × 10(-6)) of other SNPs with serum TSH, free T4, and free T3 concentrations, and these SNPs might be good targets for further studies. These results advance understanding of the genetic basis of pituitary-thyroid axis function and metabolic regulation.

Common sequence variation in FLNB regulates bone structure in women in the general population and FLNB mRNA expression in osteoblasts in vitro.

Previous data from our group indicate that BMD is linked to chromosome 3p14-p21. Because the filamin B (FLNB gene resides in this region, is the cause of skeletal dysplasias, and was identified among the top genes in our bioinformatics analysis, we hypothesized a role for FLNB in the regulation of bone structure in the general population. Using a tag single nucleotide polymorphism (SNP) approach, a family study of 767 female sibs in which the 3p14-p21 linkage with BMD was previously shown was examined. FLNB variants showing a BMD association were tested in two additional data sets, a study of 1085 UK female twins and a population study (CAIFOS) of 1315 Australian women. Genotype-expression studies were performed in 96 human osteoblast lines to examine the variants in vitro. rs7637505, rs9822918, rs2177153, and rs2001972 showed association with femoral neck (p = 0.0002-0.02) in the family-based study. The twin study provided further support for an association between rs7637505 and femoral neck and spine BMD (p = 0.02-0.03). The CAIFOS study further suggested an association between rs2177153 and rs9822918 and femoral neck BMD (p = 0.004-0.03). Prevalent fractures were increased in carriers of the A allele of rs2177153 (p = 0.009). In vitro studies showed association between rs11130605, itself in strong LD with rs7637505, and FLNB mRNA expression. These findings suggest common variants in FLNB have effects on bone structure in women. Although the location of variants having effects is not entirely consistent, variation at the 5' end of the gene may reflect effects on levels of FLNB transcription efficiency.

Genetic loci linked to pituitary-thyroid axis set points: a genome-wide scan of a large twin cohort.

OBJECTIVE: Previous studies have shown that circulating concentrations of TSH, free T4, and free T3 are genetically regulated, but the genes responsible remain largely unknown. The aim of this study was to identify genetic loci associated with these parameters. DESIGN: We performed a multipoint, nonparametric genome-wide linkage scan of 613 female dizygotic twin pairs. All subjects were euthyroid (TSH 0.4-4.0 mU/liter) with negative thyroid peroxidase antibodies and no history of thyroid disease. The genome scan comprised 737 microsatellite markers supplemented with dinucleotide markers. Data were analyzed using residualized thyroid hormone data after adjustment for age, smoking, and body mass index. RESULTS: Multipoint linkage analysis gave linkage peaks for free T4 on chromosome 14q13 and 18q21 [logarithm of odds (LOD) 2.4-3.2]; TSH on chromosomes 2q36, 4q32, and 9q34 (LOD 2.1-3.2); and free T3 on chromosomes 7q36, 8q22, and 18q21 (LOD 2.0-2.3). CONCLUSIONS: This study has identified eight genomic locations with linkage of LOD of 2.0 or greater. These results should enable targeted positional candidate and positional cloning studies to advance our understanding of genetic control of the pituitary-thyroid axis.

The association between physical activity in leisure time and leukocyte telomere length.

BACKGROUND: Physical inactivity is an important risk factor for many aging-related diseases. Leukocyte telomere dynamics (telomere length and age-dependent attrition rate) are ostensibly a biological indicator of human aging. We therefore tested the hypothesis that physical activity level in leisure time (over the past 12 months) is associated with leukocyte telomere length (LTL) in normal healthy volunteers. METHODS: We studied 2401 white twin volunteers, comprising 2152 women and 249 men, with questionnaires on physical activity level, smoking status, and socioeconomic status. Leukocyte telomere length was derived from the mean terminal restriction fragment length and adjusted for age and other potential confounders. RESULTS: Leukocyte telomere length was positively associated with increasing physical activity level in leisure time (P< .001); this association remained significant after adjustment for age, sex, body mass index, smoking, socioeconomic status, and physical activity at work. The LTLs of the most active subjects were 200 nucleotides longer than those of the least active subjects (7.1 and 6.9 kilobases, respectively; P= .006). This finding was confirmed in a small group of twin pairs discordant for physical activity level (on average, the LTL of more active twins was 88 nucleotides longer than that of less active twins; P= .03). CONCLUSIONS: A sedentary lifestyle (in addition to smoking, high body mass index, and low socioeconomic status) has an effect on LTL and may accelerate the aging process. This provides a powerful message that could be used by clinicians to promote the potentially antiaging effect of regular exercise.

Higher serum vitamin D concentrations are associated with longer leukocyte telomere length in women.

BACKGROUND: Vitamin D is a potent inhibitor of the proinflammatory response and thereby diminishes turnover of leukocytes. Leukocyte telomere length (LTL) is a predictor of aging-related disease and decreases with each cell cycle and increased inflammation. OBJECTIVE: The objective of the study was to examine whether vitamin D concentrations would attenuate the rate of telomere attrition in leukocytes, such that higher vitamin D concentrations would be associated with longer LTL. DESIGN: Serum vitamin D concentrations were measured in 2160 women aged 18-79 y (mean age: 49.4) from a large population-based cohort of twins. LTL was measured by using the Southern blot method. RESULTS: Age was negatively correlated with LTL (r = -0.40, P < 0.0001). Serum vitamin D concentrations were positively associated with LTL (r = 0.07, P = 0.0010), and this relation persisted after adjustment for age (r = 0.09, P < 0.0001) and other covariates (age, season of vitamin D measurement, menopausal status, use of hormone replacement therapy, and physical activity; P for trend across tertiles = 0.003). The difference in LTL between the highest and lowest tertiles of vitamin D was 107 base pairs (P = 0.0009), which is equivalent to 5.0 y of telomeric aging. This difference was further accentuated by increased concentrations of C-reactive protein, which is a measure of systemic inflammation. CONCLUSION: Our findings suggest that higher vitamin D concentrations, which are easily modifiable through nutritional supplementation, are associated with longer LTL, which underscores the potentially beneficial effects of this hormone on aging and age-related diseases.

Sex and ethnic differences in the association of ASPN, CALM1, COL2A1, COMP, and FRZB with genetic susceptibility to osteoarthritis of the knee.

OBJECTIVE: To assess whether the association of genetic polymorphisms with osteoarthritis (OA) in other populations could be replicated in a large, multicenter, mixed-sex, case-control study of clinical knee OA. METHODS: Genetic polymorphisms in OA candidate genes were genotyped in 298 men and 305 women, ages 50-86 years, all of whom had a diagnosis of knee OA as assessed clinically and radiographically, and in 300 male and 299 female control subjects matched for age and ethnicity. Allele and haplotype frequencies for 5 genes (ASPN, CALM1, COL2A1, COMP, and FRZB) previously tested for association with hip and/or knee OA in other populations were compared between patients and control subjects, analyzing men and women separately. RESULTS: The same FRZB 2-marker single-nucleotide polymorphism (SNP) haplotype associated with hip OA in other populations of Caucasian women was shown to increase the risk of knee OA among the women (but not the men) in the current study (odds ratio [OR] 2.87, P < 0.04). The CALM1 SNP, which affects the risk of hip OA in Japanese individuals, was not shown to be associated with susceptibility to OA in men or women. COL2A1 haplotypes were demonstrated to be associated with a decreased risk of knee OA in men (OR 0.68, P < 0.005) but not in women. COMP haplotypes that were associated with susceptibility to knee OA were different in men and women (P < 0.014 and P < 0.032, respectively). A meta-analysis of these data and those from previously published reports indicated a strong association between the FRZB G324 allele (P < 0.0003) and suggested that an ASPN allele is protective against the risk of knee OA in Caucasians (P < 0.02). CONCLUSION: Our results indicate that genetic polymorphisms affecting knee OA vary between populations (Japanese versus Caucasian) and sexes and indicate a role for ASPN, COMP, FRZB, and COL2A1 in Caucasians.

cMYB is involved in the regulation of fetal hemoglobin production in adults.

A quantitative trait locus (QTL) controlling HbF levels has previously been mapped to chromosome 6q23 in an Asian-Indian kindred with beta thalassemia and heterocellular hereditary persistence of fetal hemoglobin (HPFH). Five protein-coding genes, ALDH8A1, HBS1L, cMYB, AHI1, and PDE7B reside in this 1.5-megabase (Mb) candidate interval of 6q23. To direct sequencing efforts we compared the expression profiles of these 5 genes between 12 individuals with elevated and 14 individuals with normal HbF levels during adult erythropoiesis by real-time quantitative reverse transcription-polymerase chain reaction (RT-PCR). Two genes, cMYB and HBS1L, demonstrated simultaneous transcriptional down-regulation in individuals with elevated HbF levels. Transfection of K562 cells encoding human cDNA of cMYB and HBS1L genes showed that, although overexpression of ectopic cMYB inhibited gamma-globin gene expression, overexpression of HBS1L had no effect. Low levels of cMYB were associated with low cell expansions, accelerated erythroid maturation, and higher number of macrophages in erythroid cell culture. These observations suggest that differences in the intrinsic levels of cMYB may account for some of the variation in adult HbF levels. The possible mechanism of cMYB influencing gamma- to beta-globin switching is discussed.

Reproducible genetic associations between candidate genes and clinical knee osteoarthritis in men and women.

OBJECTIVE: Osteoarthritis (OA) is recognized to have a genetic component, and in this study, we aimed to replicate in a case-control study of men and women with clinical knee OA genetic associations in 12 candidate genes previously reported to be associated with OA. METHODS: Twenty-five single-nucleotide polymorphisms were genotyped in 298 men and 305 women ages 50-86 who were diagnosed as having knee OA, as assessed both clinically and radiographically, and in 297 men and 299 women matched for age and ethnicity (controls). Standardized anteroposterior radiographs of the knee in extension were performed on each of the cases, and all cases met the American College of Rheumatology criteria for OA of the knee. Genotype and haplotype frequencies in cases and controls were compared separately in men and women. The 12 genes tested were AACT, ADAM12, BMP2, CD36, CILP, COX2, ESR1, NCOR2, OPG, TNA, TNFAIP6, and VDR. RESULTS: Eight of the candidate genes were associated in women and 5 in men, and only 3 genes (TNFAIP6, NCOR2, and CD36) were not significantly associated in either sex. The strongest associations in terms of odds ratios (ORs) were a haplotype in ADAM12 (OR 7.1 [95% confidence interval (95% CI) 3.3-33.8]) and a haplotype in ESR1 (OR 3.6 [95% CI 1.18-10.98]) in women. The same ADAM12 haplotype (OR 2.54 [95% CI 1.2-5.4]) and a haplotype in the CILP gene (OR 0.38 [95% CI 0.23-0.62]) were the strongest associations in men. CONCLUSION: We found that genes previously identified by their association with subclinical features of knee OA or progression were also associated with clinical knee OA. These genetic associations may identify individuals at a particularly high risk of developing knee OA.

Linkage and potential association of obesity-related phenotypes with two genes on chromosome 12q24 in a female dizygous twin cohort.

Obesity is a multifactorial disorder with a complex phenotype. It is a significant risk factor for diabetes and hypertension. We assessed obesity-related traits in a large cohort of twins and performed a genome-wide linkage scan and positional candidate analysis to identify genes that play a role in regulating fat mass and distribution in women. Dizygous female twin pairs from 1,094 pedigrees were studied (mean age 47.0+/-11.5 years (range 18-79 years)). Nonparametric multipoint linkage analyses showed linkage for central fat mass to 12q24 (141 cM) with LOD 2.2 and body mass index to 8q11 (67 cM) with LOD 1.3, supporting previously established linkage data. Novel areas of suggestive linkage were for total fat percentage at 6q12 (LOD 2.4) and for total lean mass at 2q37 (LOD 2.4). Data from follow-up fine mapping in an expanded cohort of 1243 twin pairs reinforced the linkage for central fat mass to 12q24 (LOD 2.6; 143 cM) and narrowed the -1 LOD support interval to 22 cM. In all, 45 single-nucleotide polymorphisms (SNPs) from 26 positional candidate genes within the 12q24 interval were then tested for association in a cohort of 1102 twins. Single-point Monks-Kaplan analysis provided evidence of association between central fat mass and SNPs in two genes - PLA2G1B (P = 0.0067) and P2RX4 (P = 0.017). These data provide replication and refinement of the 12q24 obesity locus and suggest that genes involved in phospholipase and purinoreceptor pathways may regulate fat accumulation and distribution.

Mapping genetic loci that determine leukocyte telomere length in a large sample of unselected female sibling pairs.

Telomeres play a central role in cellular senescence and cancer pathobiology and are associated with age-related diseases such as atherosclerosis and dementia. Telomere length varies between individuals of the same age, is influenced by DNA-damaging factors such as oxidative stress, and is heritable. We performed a quantitative-trait linkage analysis using an approximate 10-cM genomewide map for mean leukocyte terminal-restriction fragment (TRF) lengths measured by Southern blotting, in 2,050 unselected women aged 18-80 years, comprising 1,025 complete dizygotic twin pairs. Heritability of mean batch-adjusted TRF was 36% (95% confidence interval [CI] 18%-48%), with a large common environmental effect of 49% (95% CI 40%-58%). Significant linkage was observed on chromosome 14 (LOD 3.9) at 14q23.2, and suggestive linkage at 10q26.13 (LOD 2.4) and 3p26.1 (LOD 2.7). This is the first report of loci, mapped in a sample of healthy individuals, that influence mean telomere variation in humans.

Identification of QTLs for serum lipid levels in a female sib-pair cohort: a novel application to improve the power of two-locus linkage analysis.

Using a novel approach for a two-locus model that provides a greatly increased power to detect multiple quantitative trait loci (QTLs) in simulated data, we identified in a sample of 961 female sib-pairs, three genome-wide significant QTLs for apolipoprotein A1 on chromosomes 8p21.1-q13.1 (LOD score 3.71), 9q21.32-33.1 (LOD score 3.28) and 10p15.1-p13 (LOD score 5.51), two for lipoprotein (a) on chromosomes 6q25.2-q27 (LOD score 10.18) and 21q21.1-q21.3 (LOD score 4.57) and two for triglycerides on chromosomes 4q28.3-32.1 (LOD score 3.71) and 5q23.1-q32 (LOD score 3.60). The two-locus ordered-subset analysis has led to the confirmation of known and likely identification of novel regions linked to serum lipid levels that would have otherwise been missed and deserves wider application in linkage analyses of quantitative traits. Given the relative lack of power for the sample sizes commonly used in human genetics linkage studies, minor QTL effects often go undetected and those that are detected will be upwardly biased. We show through simulation that the discrepancy between the real and estimated QTL-effects is often likely to generate an unpredictable source of false-negative errors, using multi-locus models, reducing the power to detect multiple QTLs through oligogenic linkage analysis. The successful simultaneous modelling of the identified QTLs in a multi-locus context helps to eliminate false positives and increases the power to detect linkages, adding compelling evidence that they are likely to be reliable QTLs for these lipid traits.