Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications

BACKGROUND

Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression.

METHODS

We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites.

FINDINGS

We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2).

INTERPRETATION

The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications.

FUNDING

Funding acknowledgements for each cohort can be found in the Supplementary Note.

Epigenome-wide Association Study Shows Differential DNA Methylation of MDC1, KLF9, and CUTA in Autoimmune Thyroid Disease

CONTEXT

Autoimmune thyroid disease (AITD) includes Graves disease (GD) and Hashimoto disease (HD), which often run in the same family. AITD etiology is incompletely understood: Genetic factors may account for up to 75% of phenotypic variance, whereas epigenetic effects (including DNA methylation [DNAm]) may contribute to the remaining variance (eg, why some individuals develop GD and others HD).

OBJECTIVE

This work aimed to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) comparing GD to HD.

METHODS

Whole-blood DNAm was measured across the genome using the Infinium MethylationEPIC array in 32 Australian patients with GD and 30 with HD (discovery cohort) and 32 Danish patients with GD and 32 with HD (replication cohort). Linear mixed models were used to test for differences in quantile-normalized β values of DNAm between GD and HD and data were later meta-analyzed. Comb-p software was used to identify DMRs.

RESULTS

We identified epigenome-wide significant differences (P < 9E-8) and replicated (P < .05) 2 DMPs between GD and HD (cg06315208 within MDC1 and cg00049440 within KLF9). We identified and replicated a DMR within CUTA (5 CpGs at 6p21.32). We also identified 64 DMPs and 137 DMRs in the meta-analysis.

CONCLUSION

Our study reveals differences in DNAm between GD and HD, which may help explain why some people develop GD and others HD and provide a link to environmental risk factors. Additional research is needed to advance understanding of the role of DNAm in AITD and investigate its prognostic and therapeutic potential.

Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry

BACKGROUND

Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes.

RESULTS

The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata.

CONCLUSIONS

Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes

OBJECTIVES

Osteoarthritis (OA) is a joint disease with a heritable component. Genetic loci identified via genome-wide association studies (GWAS) account for an estimated 26.3% of the disease trait variance in humans. Currently, there is no method for predicting the onset or progression of OA. We describe the first use of the Collaborative Cross (CC), a powerful genetic resource, to investigate knee OA in mice, with follow-up targeted multi-omics analysis of homologous regions of the human genome.

METHODS

We histologically screened 275 mice for knee OA and conducted quantitative trait locus (QTL) mapping in the complete cohort (> 8 months) and the younger onset sub-cohort (8-12 months). Multi-omic analysis of human genetic datasets was conducted to investigate significant loci.

RESULTS

We observed a range of OA phenotypes. QTL mapping identified a genome-wide significant locus on mouse chromosome 19 containing Glis3, the human equivalent of which has been identified as associated with OA in recent GWAS. Mapping the younger onset sub-cohort identified a genome-wide significant locus on chromosome 17. Multi-omic analysis of the homologous region of the human genome (6p21.32) indicated the presence of pleiotropic effects on the expression of the HLA - DPB2 gene and knee OA development risk, potentially mediated through the effects on DNA methylation.

CONCLUSIONS

The significant associations at the 6p21.32 locus in human datasets highlight the value of the CC model of spontaneous OA that we have developed and lend support for an immune role in the disease. Our results in mice also add to the accumulating evidence of a role for Glis3 in OA.

Leveraging osteoclast genetic regulatory data to identify genes with a role in osteoarthritis

There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.

Genetic determinants of thyroid function in children

OBJECTIVE

Genome-wide association studies in adults have identified 42 loci associated with thyroid stimulating hormone (TSH) and 21 loci associated with free thyroxine (FT4) concentrations. While biologically plausible, age-dependent effects have not been assessed. We aimed to study the association of previously identified genetic determinants of TSH and FT4 with TSH and FT4 concentrations in newborns and (pre)school children.

METHODS

We selected participants from three population-based prospective cohorts with data on genetic variants and thyroid function: Generation R (N = 2169 children, mean age 6 years; N = 2388 neonates, the Netherlands), the Avon Longitudinal Study of Parents and Children (ALSPAC; N = 3382, age 7.5 years, United Kingdom), and the Brisbane Longitudinal Twin Study (BLTS; N = 1680, age 12.1 years, Australia). The association of single nucleotide polymorphisms (SNPs) with TSH and FT4 concentrations was studied with multivariable linear regression models. Weighted polygenic risk scores (PRSs) were defined to combine SNP effects.

RESULTS

In childhood, 30/60 SNPs were associated with TSH and 11/31 SNPs with FT4 after multiple testing correction. The effect sizes for AADAT, GLIS3, TM4SF4, and VEGFA were notably larger than in adults. The TSH PRS explained 5.3%-8.4% of the variability in TSH concentrations; the FT4 PRS explained 1.5%-4.2% of the variability in FT4 concentrations. Five TSH SNPs and no FT4 SNPs were associated with thyroid function in neonates.

CONCLUSIONS

The effects of many known thyroid function SNPs are already apparent in childhood and some might be notably larger in children as compared to adults. These findings provide new knowledge about genetic regulation of thyroid function in early life.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Digenic Congenital Hypogonadotropic Hypogonadism Due to Heterozygous GNRH1 p.R31C and AMHR2 p.G445_L453del Variants

A 28-year-old man with congenital hypogonadotropic hypogonadism (CHH) was found to be heterozygous for the GNRH1 p.R31C mutation, reported in the literature as pathogenic and dominant. The same mutation was found in his son at birth, but the testing of the infant at 64 days confirmed the hormonal changes associated with minipuberty. This led to further genetic sequencing of the patient and his son, which found a second variant, AMHR2 p.G445_L453del, in the heterozygous form, reported as pathogenic in the patient but not in his son. This suggests a digenic cause of the patient's CHH. Together, these mutations are postulated to contribute to CHH by the lack of anti-Müllerian hormone (AMH) signalling, leading to the impaired migration of gonadotrophin releasing hormone (GnRH) neurons, the lack of the AMH effect on GnRH secretion, and altered GnRH decapeptide with reduced binding to GnRH receptors. This led us to the conclusion that the observed GNRH1 mutation in the heterozygous state is not certain to be dominant or, at least, exhibits incomplete penetrance and variable expressivity. This report also emphasises the opportunity afforded by the time window of minipuberty in assessing the inherited genetic disorders of hypothalamic function.

Identification of Differentially Expressed Genes and Molecular Pathways Involved in Osteoclastogenesis Using RNA-seq

Osteoporosis is a disease that is characterised by reduced bone mineral density (BMD) and can be exacerbated by the excessive bone resorption of osteoclasts (OCs). Bioinformatic methods, including functional enrichment and network analysis, can provide information about the underlying molecular mechanisms that participate in the progression of osteoporosis. In this study, we harvested human OC-like cells differentiated in culture and their precursor peripheral blood mononuclear cells (PBMCs) and characterised the transcriptome of the two cell types using RNA-sequencing in order to identify differentially expressed genes. Differential gene expression analysis was performed in RStudio using the edgeR package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to identify enriched GO terms and signalling pathways, with inter-connected regions characterised using protein-protein interaction analysis. In this study, we identified 3201 differentially expressed genes using a 5% false discovery rate; 1834 genes were upregulated, whereas 1367 genes were downregulated. We confirmed a significant upregulation of several well-established OC genes including CTSK, DCSTAMP, ACP5, MMP9, ITGB3, and ATP6V0D2. The GO analysis suggested that upregulated genes are involved in cell division, cell migration, and cell adhesion, while the KEGG pathway analysis highlighted oxidative phosphorylation, glycolysis and gluconeogenesis, lysosome, and focal adhesion pathways. This study provides new information about changes in gene expression and highlights key biological pathways involved in osteoclastogenesis.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

DNA Methylation in Autoimmune Thyroid Disease

Graves disease and Hashimoto disease form part of the spectrum of autoimmune thyroid disease (AITD), to which genetic and environmental factors are recognized contributors. Epigenetics provides a potential link between environmental influences, gene expression, and thyroid autoimmunity. DNA methylation (DNAm) is the best studied epigenetic process, and global hypomethylation of leukocyte DNA is reported in several autoimmune disorders. This review summarizes the current understanding of DNAm in AITD. Targeted DNAm studies of blood samples from AITD patients have reported differential DNAm in the promoter regions of several genes implicated in AITD, including TNF, IFNG, IL2RA, IL6, ICAM1, and PTPN22. In many cases, however, the findings await replication and are unsupported by functional studies to support causal roles in AITD pathogenesis. Furthermore, thyroid hormones affect DNAm, and in many studies confounding by reverse causation has not been considered. Recent studies have shown that DNAm patterns in candidate genes including ITGA6, PRKAA2, and DAPK1 differ between AITD patients from regions with different iodine status, providing a potential mechanism for associations between iodine and AITD. Research focus in the field is moving from candidate gene studies to an epigenome-wide approach. Genome-wide methylation studies of AITD patients have demonstrated multiple differentially methylated positions, including some in immunoregulatory genes such as NOTCH1, HLA-DRB1, TNF, and ICAM1. Large, epigenome-wide studies are required to elucidate the pathophysiological role of DNAm in AITD, with the potential to provide novel diagnostic and prognostic biomarkers as well as therapeutic targets.

Shared genetics and causal relationships between migraine and thyroid function traits

BACKGROUND

Epidemiological studies have reported a comorbid relationship between migraine and thyroid dysfunction.

METHODS

We investigated the genetic relationship between migraine and thyroid function traits using genome-wide association study (GWAS) data.

RESULTS

We found a significant genetic correlation (r_g) with migraine for hypothyroidism (r_g = 0.0608), secondary hypothyroidism (r_g = 0.195), free thyroxine (fT4) (r_g = 0.0772), and hyperthyroidism (r_g = -0.1046), but not thyroid stimulating hormone (TSH). Pairwise GWAS analysis revealed two shared loci with TSH and 11 shared loci with fT4. Cross-trait GWAS meta-analysis of migraine identified novel genome-wide significant loci: 17 with hypothyroidism, one with hyperthyroidism, five with secondary hypothyroidism, eight with TSH, and 15 with fT4. Of the genes at these loci, six (RERE, TGFB2, APLF, SLC9B1, SGTB, BTBD16; migraine + hypothyroidism), three (GADD45A, PFDN1, RSPH6A; migraine + TSH), and three (SSBP3, BRD3, TEF; migraine + fT4) were significant in our gene-based analysis (p_{Fisher's combined P-value} < 2.04 × 10^-6). In addition, causal analyses suggested a negative causal relationship between migraine and hyperthyroidism (p = 8.90 × 10^-3) and a positive causal relationship between migraine and secondary hypothyroidism (p = 1.30 × 10^-3).

CONCLUSION

These findings provide strong evidence for genetic correlation and suggest complex causal relationships between migraine and thyroid traits.

Cross-Trait Genetic Analyses Indicate Pleiotropy and Complex Causal Relationships between Headache and Thyroid Function Traits

Epidemiological studies have reported a comorbid relationship between headache and thyroid traits; however, little is known about the shared genetics and causality that contributes to this association. We investigated the genetic overlap and associations between headache and thyroid function traits using genome-wide association study (GWAS) data. We found a significant genetic correlation (rg) with headache and hypothyroidism (rg = 0.09, p = 2.00 × 10−4), free thyroxine (fT4) (rg = 0.08, p = 5.50 × 10−3), and hyperthyroidism (rg = −0.14, p = 1.80 × 10−3), a near significant genetic correlation with secondary hypothyroidism (rg = 0.20, p = 5.24 × 10−2), but not with thyroid stimulating hormone (TSH). Pairwise-GWAS analysis revealed six, 14, four and five shared (pleiotropic) loci with headache and hypothyroidism, hyperthyroidism, secondary hypothyroidism, and fT4, respectively. Cross-trait GWAS meta-analysis identified novel genome-wide significant loci for headache: five with hypothyroidism, three with secondary hypothyroidism, 12 with TSH, and nine with fT4. Of the genes at these loci, six (FAF1, TMX2-CTNND1, AARSD1, PLCD3, ZNF652, and C20orf203; headache-TSH) and six (HMGB1P45, RPL30P1, ZNF462, TMX2-CTNND1, ITPK1, SECISBP2L; headache-fT4) were significant in our gene-based analysis (pFisher’s combined p-value < 2.09 × 10−6). Our causal analysis suggested a positive causal relationship between headache and secondary hypothyroidism (p = 3.64 × 10−4). The results also suggest a positive causal relationship between hypothyroidism and headache (p = 2.45 × 10−3) and a negative causal relationship between hyperthyroidism and headache (p = 1.16 × 10−13). These findings suggest a strong evidence base for a genetic correlation and complex causal relationships between headache and thyroid traits.

Functional Assessment of Calcium-Sensing Receptor Variants Confirms Familial Hypocalciuric Hypercalcaemia

Context

In the clinic it is important to differentiate primary hyperparathyroidism (PHPT) from the more benign, inherited disorder, familial hypocalciuric hypercalcemia (FHH). Since the conditions may sometimes overlap biochemically, identification of calcium-sensing receptor (CASR) gene variants causative of FHH (but not PHPT) is the most decisive diagnostic aid. When novel variants are identified, bioinformatics and functional assessment are required to establish pathogenicity.

Objective

We identified 3 novel CASR transmembrane domain missense variants, Thr699Asn, Arg701Gly, and Thr808Pro, in 3 probands provisionally diagnosed with FHH and examined the variants using bioinformatics and functional analysis.

Methods

Bioinformatics assessment utilized wANNOVAR software. For functional characterization, each variant was cloned into a mammalian expression vector; wild-type and variant receptors were transfected into HEK293 cells, and their expression and cellular localization were assessed by Western blotting and confocal immunofluorescence, respectively. Receptor activation in HEK293 cells was determined using an IP-One ELISA assay following stimulation with Ca⁺⁺ ions.

Results

Bioinformatics analysis of the variants was unable to definitively assign pathogenicity. Compared with wild-type receptor, all variants demonstrated impaired expression of mature receptor reaching the cell surface and diminished activation at physiologically relevant Ca⁺⁺ concentrations.

Conclusion

Three CASR missense variants identified in probands provisionally diagnosed with FHH result in receptor inactivation and are therefore likely causative of FHH. Inactivation may be due to inadequate processing/trafficking of mature receptor and/or conformational changes induced by the variants affecting receptor signaling. This study demonstrates the value of functional studies in assessing genetic variants identified in hypercalcemic patients.

Genome-wide analysis of thyroid function in Australian adolescents highlights SERPINA7 and NCOA3

OBJECTIVE

Genetic factors underpin the narrow intraindividual variability of thyroid function, although precise contributions of environmental vs genetic factors remain uncertain. We sought to clarify the heritability of thyroid function traits and thyroid peroxidase antibody (TPOAb) positivity and identify single nucleotide polymorphisms (SNPs) contributing to the trait variance.

METHODS

Heritability of thyroid-stimulating hormone (TSH), free T4 (fT4), free T3 (fT3) and TPOAb in a cohort of 2854 euthyroid, dizygous and monozygous twins (age range 11.9-16.9 years) from the Brisbane Longitudinal Twin Study (BLTS) was assessed using structural equation modelling. A genome-wide analysis was conducted on 2832 of these individuals across 7 522 526 SNPs as well as gene-based association analyses. Replication analysis of the association results was performed in the Raine Study (n = 1115) followed by meta-analysis to maximise power for discovery.

RESULTS

Heritability of thyroid function parameters in the BLTS was 70.8% (95% CI: 66.7-74.9%) for TSH, 67.5% (59.8-75.3%) for fT4, 59.7% (54.4-65.0%) for fT3 and 48.8% (40.6-56.9%) for TPOAb. The genome-wide association study (GWAS) in the discovery cohort identified a novel association between rs2026401 upstream of NCOA3 and TPOAb. GWAS meta-analysis found associations between TPOAb and rs445219, also near NCOA3, and fT3 and rs12687280 near SERPINA7. Gene-based association analysis highlighted SERPINA7 for fT3 and NPAS3 for fT4.

CONCLUSION

Our findings resolve former contention regarding heritability estimates of thyroid function traits and TPOAb positivity. GWAS and gene-based association analysis identified variants accounting for a component of this heritability.

Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L

CONTEXT

Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited.

OBJECTIVE

To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts.

METHOD

We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed.

RESULTS

We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH.

CONCLUSIONS

This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Familial Dysalbuminemic Hyperthyroxinemia as a Cause for Discordant Thyroid Function Tests

Introduction

Discordant thyroid function tests are routinely encountered in clinical practice. Differential diagnoses include acute thyroxine (T4) ingestion, laboratory interference from heterophilic antibodies, thyroid hormone resistance, thyroid-stimulating hormone (TSH)-secreting pituitary adenomas, and T4 protein binding abnormalities. The impact of abnormal binding proteins may be less recognized since widespread use of free T4 (FT4) assays compared to older total T4 assays.

Case report

A 69-year-old female was referred for assessment of discordant thyroid function tests. Biochemistry since July 2015 showed persistently elevated FT4 levels by immunoassay ranging between 25 to 34 pmol/L with normal or slightly decreased TSH ranging between 0.05 to 2.74 mU/L. The patient was clinically euthyroid on 100 mcg daily of levothyroxine for Hashimoto’s thyroiditis. FT4 measured using liquid chromatography-tandem mass spectrometry (LC-MS/MS) was 19.5 pmol/L. Exome sequencing (confirmed by Sanger sequencing) detected a guanine to adenine substitution at residue 725 of the ALB gene previously associated with dysalbuminemic hyperthyroxinemia. The patient’s daughter had similar thyroid function tests and the same genetic variant. FT4 results from 3 different automated immunoassays showed the Roche Cobas and Siemens Centaur platforms to be most affected by the variant, and Abbott Architect had the best agreement with LC-MS/MS.

Conclusion

Familial dysalbuminemic hyperthyroxinemia is a potential cause of discordant thyroid function tests. Clinicians suspecting protein-binding abnormalities may further investigate using reference methods such as LC-MS/MS and equilibrium dialysis if available. The increasing accessibility of exome sequencing offers a cost-effective method of diagnosing genetic variants that cause discordant thyroid function tests.

Functional Analysis of Calcium-Sensing Receptor Variants Identified in Families Provisionally Diagnosed with Familial Hypocalciuric Hypercalcaemia

Identification of variants in the calcium-sensing receptor (CASR) gene is an important means of distinguishing between familial hypocalciuric hypercalcaemia (FHH) and primary hyperparathyroidism. However, identification and bioinformatics analysis of genetic variants alone is now considered insufficient as definitive proof; additional functional assessment is required to diagnose FHH with certainty. We identified two novel variants, D433Y and C739Y, and one previously reported variant G509R in the CASR of four kindreds provisionally diagnosed with FHH and aimed to functionally characterise these variants to confirm the diagnosis. Variant receptors were cloned as FLAG-tagged constructs into the mammalian expression vector, pcDNA3.1. Wild type and variant receptor constructs were expressed in HEK293 cells and their expression assessed by Western blot analysis and their functionality analysed using an IP-One assay which measures myo-inositol 1-phosphate accumulation following CaSR activation. Western blot analysis showed that the D433Y receptor had diminished mature glycosylated receptor compared with wild type CaSR whereas the G509R receptor had a complete lack of mature receptor. The C739Y receptor was consistently overexpressed. Functional assessment showed the D433Y receptor to be mildly inactivating at physiological calcium concentrations whereas the G509R receptor was inactive at all calcium concentrations. By contrast, the C739Y variant was activating compared to wild type receptor which is inconsistent with it causing FHH. We conclude that functional assessment of CaSR variants using the IP-One assay was useful in the investigation of suspected FHH probands, confirming the D433Y and G509R variants as likely pathogenic/pathogenic, but dismissing the C739Y variant as causing FHH.

Decreased Immunoglobulin G Core Fucosylation, A Player in Antibody-dependent Cell-mediated Cytotoxicity, is Associated with Autoimmune Thyroid Diseases

Autoimmune thyroid diseases (AITD) are the most common group of autoimmune diseases, associated with lymphocyte infiltration and the production of thyroid autoantibodies, like thyroid peroxidase antibodies (TPOAb), in the thyroid gland. Immunoglobulins and cell-surface receptors are glycoproteins with distinctive glycosylation patterns that play a structural role in maintaining and modulating their functions. We investigated associations of total circulating IgG and peripheral blood mononuclear cells glycosylation with AITD and the influence of genetic background in a case-control study with several independent cohorts and over 3,000 individuals in total. The study revealed an inverse association of IgG core fucosylation with TPOAb and AITD, as well as decreased peripheral blood mononuclear cells antennary α1,2 fucosylation in AITD, but no shared genetic variance between AITD and glycosylation. These data suggest that the decreased level of IgG core fucosylation is a risk factor for AITD that promotes antibody-dependent cell-mediated cytotoxicity previously associated with TPOAb levels.

Changes in Thyroid Function Across Adolescence: A Longitudinal Study

OBJECTIVE

There are no large, longitudinal studies of thyroid function across adolescence. The aims of this study were to examine longitudinal trends in thyrotropin (TSH), free triiodothyronine (fT3) and free thyroxine (fT4) and determine age-specific reference ranges.

METHODS

Thyroid function was assessed in 3415 participants in the Brisbane Longitudinal Twin Study at ages 12, 14, and 16, using the Abbott ARCHITECT immunoassay. Longitudinal analyses were adjusted for body mass index and puberty.

RESULTS

In girls, mean fT4 (± SE) increased between age 12 and 14 (by 0.30 ± 0.08 pmol/L; P < 0.001), while remaining unchanged in boys; from age 14 to 16, fT4 increased in both girls (by 0.42 ± 0.07 pmol/L; P < 0.001) and boys (0.64 ± 0.07 pmol/L, P < 0.001). There was a slight increase in fT3 from age 12 to 14 years in girls (by 0.07 ± 0.03 pmol/L; P = 0.042), with a more marked increase in boys (0.29 ± 0.03 pmol/L; P < 0.001), followed by a decrease from age 14 to 16 in both sexes (girls, by 0.53 ± 0.02 pmol/L; P < 0.001; boys, by 0.62 ± 0.03 pmol/L; P < 0.001). From age 12 to 14, TSH showed no significant change in girls or boys, then levels increased from age 14 to 16 in both sexes (in girls, by 4.9%, 95% CI: 2.4%-10.3%, P = 0.020; in boys, by 7.2%, 95% CI: 3.0%-11.6%, P = 0.001). Reference ranges differed substantially from adults, particularly for fT4 and fT3.

CONCLUSIONS

Thyroid function tests in adolescents display complex, sexually dimorphic patterns. Implementation of adolescence-specific reference ranges may be appropriate.

Characterisation of genetic regulatory effects for osteoporosis risk variants in human osteoclasts

Background

Osteoporosis is a complex disease with a strong genetic contribution. A recently published genome-wide association study (GWAS) for estimated bone mineral density (eBMD) identified 1103 independent genome-wide significant association signals. Most of these variants are non-coding, suggesting that regulatory effects may drive many of the associations. To identify genes with a role in osteoporosis, we integrate the eBMD GWAS association results with those from our previous osteoclast expression quantitative trait locus (eQTL) dataset.

Results

We identify sixty-nine significant cis-eQTL effects for eBMD GWAS variants after correction for multiple testing. We detect co-localisation of eBMD GWAS and osteoclast eQTL association signals for 21 of the 69 loci, implicating a number of genes including CCR5, ZBTB38, CPE, GNA12, RIPK3, IQGAP1 and FLCN. Summary-data-based Mendelian Randomisation analysis of the eBMD GWAS and osteoclast eQTL datasets identifies significant associations for 53 genes, with TULP4 presenting as a strong candidate for pleiotropic effects on eBMD and gene expression in osteoclasts. By performing analysis using the GARFIELD software, we demonstrate significant enrichment of osteoporosis risk variants among high-confidence osteoclast eQTL across multiple GWAS P value thresholds. Mice lacking one of the genes of interest, the apoptosis/necroptosis gene RIPK3, show disturbed bone micro-architecture and increased osteoclast number, highlighting a new biological pathway relevant to osteoporosis.

Conclusion

We utilise a unique osteoclast eQTL dataset to identify a number of potential effector genes for osteoporosis risk variants, which will help focus functional studies in this area.

Dysregulated Antibody, Natural Killer Cell and Immune Mediator Profiles in Autoimmune Thyroid Diseases

The pathogenesis of autoimmune thyroid diseases (AITD) is poorly understood and the association between different immune features and the germline variants involved in AITD are yet unclear. We previously observed systemic depletion of IgG core fucosylation and antennary α1,2 fucosylation in peripheral blood mononuclear cells in AITD, correlated with anti-thyroid peroxidase antibody (TPOAb) levels. Fucose depletion is known to potentiate strong antibody-mediated NK cell activation and enhanced target antigen-expressing cell killing. In autoimmunity, this may translate to autoantibody-mediated immune cell recruitment and attack of self-antigen expressing normal tissues. Hence, we investigated the crosstalk between immune cell traits, secreted proteins, genetic variants and the glycosylation patterns of serum IgG, in a multi-omic and cross-sectional study of 622 individuals from the TwinsUK cohort, 172 of whom were diagnosed with AITD. We observed associations between two genetic variants (rs505922 and rs687621), AITD status, the secretion of Desmoglein-2 protein, and the profile of two IgG N-glycan traits in AITD, but further studies need to be performed to better understand their crosstalk in AITD. On the other side, enhanced afucosylated IgG was positively associated with activatory CD335- CD314+ CD158b+ NK cell subsets. Increased levels of the apoptosis and inflammation markers Caspase-2 and Interleukin-1α positively associated with AITD. Two genetic variants associated with AITD, rs1521 and rs3094228, were also associated with altered expression of the thyrocyte-expressed ligands known to recognize the NK cell immunoreceptors CD314 and CD158b. Our analyses reveal a combination of heightened Fc-active IgG antibodies, effector cells, cytokines and apoptotic signals in AITD, and AITD genetic variants associated with altered expression of thyrocyte-expressed ligands to NK cell immunoreceptors. Together, TPOAb responses, dysregulated immune features, germline variants associated with immunoactivity profiles, are consistent with a positive autoreactive antibody-dependent NK cell-mediated immune response likely drawn to the thyroid gland in AITD.

Complement C5a Induces Renal Injury in Diabetic Kidney Disease by Disrupting Mitochondrial Metabolic Agility

The sequelae of diabetes include microvascular complications such as diabetic kidney disease (DKD), which involves glucose-mediated renal injury associated with a disruption in mitochondrial metabolic agility, inflammation, and fibrosis. We explored the role of the innate immune complement component C5a, a potent mediator of inflammation, in the pathogenesis of DKD in clinical and experimental diabetes. Marked systemic elevation in C5a activity was demonstrated in patients with diabetes; conventional renoprotective agents did not therapeutically target this elevation. C5a and its receptor (C5aR1) were upregulated early in the disease process and prior to manifest kidney injury in several diverse rodent models of diabetes. Genetic deletion of C5aR1 in mice conferred protection against diabetes-induced renal injury. Transcriptomic profiling of kidney revealed diabetes-induced downregulation of pathways involved in mitochondrial fatty acid metabolism. Interrogation of the lipidomics signature revealed abnormal cardiolipin remodeling in diabetic kidneys, a cardinal sign of disrupted mitochondrial architecture and bioenergetics. In vivo delivery of an orally active inhibitor of C5aR1 (PMX53) reversed the phenotypic changes and normalized the renal mitochondrial fatty acid profile, cardiolipin remodeling, and citric acid cycle intermediates. In vitro exposure of human renal proximal tubular epithelial cells to C5a led to altered mitochondrial respiratory function and reactive oxygen species generation. These experiments provide evidence for a pivotal role of the C5a/C5aR1 axis in propagating renal injury in the development of DKD by disrupting mitochondrial agility, thereby establishing a new immunometabolic signaling pathway in DKD.

Meta-Analysis of Genomewide Association Studies Reveals Genetic Variants for Hip Bone Geometry

Hip geometry is an important predictor of fracture. We performed a meta-analysis of GWAS studies in adults to identify genetic variants that are associated with proximal femur geometry phenotypes. We analyzed four phenotypes: (i) femoral neck length; (ii) neck-shaft angle; (iii) femoral neck width, and (iv) femoral neck section modulus, estimated from DXA scans using algorithms of hip structure analysis. In the Discovery stage, 10 cohort studies were included in the fixed-effect meta-analysis, with up to 18,719 men and women ages 16 to 93 years. Association analyses were performed with ∼2.5 million polymorphisms under an additive model adjusted for age, body mass index, and height. Replication analyses of meta-GWAS significant loci (at adjusted genomewide significance [GWS], threshold p ≤ 2.6 × 10^-8 ) were performed in seven additional cohorts in silico. We looked up SNPs associated in our analysis, for association with height, bone mineral density (BMD), and fracture. In meta-analysis (combined Discovery and Replication stages), GWS associations were found at 5p15 (IRX1 and ADAMTS16); 5q35 near FGFR4; at 12p11 (in CCDC91); 11q13 (near LRP5 and PPP6R3 (rs7102273)). Several hip geometry signals overlapped with BMD, including LRP5 (chr. 11). Chr. 11 SNP rs7102273 was associated with any-type fracture (p = 7.5 × 10^-5 ). We used bone transcriptome data and discovered several significant eQTLs, including rs7102273 and PPP6R3 expression (p = 0.0007), and rs6556301 (intergenic, chr.5 near FGFR4) and PDLIM7 expression (p = 0.005). In conclusion, we found associations between several genes and hip geometry measures that explained 12% to 22% of heritability at different sites. The results provide a defined set of genes related to biological pathways relevant to BMD and etiology of bone fragility. © 2019 American Society for Bone and Mineral Research.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Expression Quantitative Trait Locus Study of Bone Mineral Density GWAS Variants in Human Osteoclasts

Osteoporosis is a complex disease with a strong genetic component. Genomewide association studies (GWAS) have been very successful at identifying common genetic variants associated with bone parameters. A recently published study documented the results of the largest GWAS for bone mineral density (BMD) performed to date (n = 142,487), identifying 307 conditionally independent single-nucleotide polymorphisms (SNPs) as associated with estimated BMD (eBMD) at the genomewide significance level. The vast majority of these variants are non-coding SNPs. Expression quantitative trait locus (eQTL) studies using disease-specific cell types have increasingly been integrated with the results from GWAS to identify genes through which the observed GWAS associations are likely mediated. We generated a unique human osteoclast-specific eQTL data set using cells differentiated in vitro from 158 participants. We then used this resource to characterize the 307 recently identified BMD GWAS SNPs for association with nearby genes (±500 kb). After correction for multiple testing, 24 variants were found to be significantly associated with the expression of 32 genes in the osteoclast-like cells. Bioinformatics analysis suggested that these variants and those in strong linkage disequilibrium with them are enriched in regulatory regions. Several of the eQTL associations identified are relevant to genes that present strongly as having a role in bone, particularly IQGAP1, CYP19A1, CTNNB1, and COL6A3. Supporting evidence for many of the associations was obtained from publicly available eQTL data sets. We have also generated strong evidence for the presence of a regulatory region on chromosome 15q21.2 relevant to both the GLDN and CYP19A1 genes. In conclusion, we have generated a unique osteoclast-specific eQTL resource and have used this to identify 32 eQTL associations for recently identified BMD GWAS loci, which should inform functional studies of osteoclast biology. © 2018 American Society for Bone and Mineral Research.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

How many cases of disease in a pedigree imply familial disease?

The ability to perform whole-exome and, increasingly, whole-genome sequencing on large numbers of individuals has led to increased efforts to identify rare genetic variants that affect the risk of both common and rare diseases. In such applications, it is important to identify families that are segregating the rare variants of interest. For rare diseases or rare familial forms of common diseases, pedigrees with multiple affected members are clearly harbouring risk variants. For more common diseases, however, it may be unclear whether a family with a few affected members is segregating a familial disease, is the result of multiple sporadic cases, or is a mixture of familial cases and phenocopies. We provide calculations for the probability that a family is harbouring familial disease, presented in general terms that admit working guidelines for selecting families for current sequencing studies. Using examples motivated by our own studies of thyroid cancer and published studies of colorectal cancer, we show that for common diseases, families with exactly two affected first-degree relatives have only a moderate probability of segregating familial disease, but this probability is higher for families with three or more affected relatives, and those families should therefore be prioritised in sequencing studies.

Genome-wide association study meta-analysis for quantitative ultrasound parameters of bone identifies five novel loci for broadband ultrasound attenuation

Osteoporosis is a common and debilitating bone disease that is characterised by low bone mineral density, typically assessed using dual-energy X-ray absorptiometry. Quantitative ultrasound (QUS), commonly utilising the two parameters velocity of sound (VOS) and broadband ultrasound attenuation (BUA), is an alternative technology used to assess bone properties at peripheral skeletal sites. The genetic influence on the bone qualities assessed by QUS remains an under-studied area. We performed a comprehensive genome-wide association study (GWAS) including low-frequency variants (minor allele frequency ≥0.005) for BUA and VOS using a discovery population of individuals with whole-genome sequence (WGS) data from the UK10K project (n = 1268). These results were then meta-analysed with those from two deeply imputed GWAS replication cohorts (n = 1610 and 13 749). In the gender-combined analysis, we identified eight loci associated with BUA and five with VOS at the genome-wide significance level, including three novel loci for BUA at 8p23.1 (PPP1R3B), 11q23.1 (LOC387810) and 22q11.21 (SEPT5) (P = 2.4 × 10-8 to 1.6 × 10-9). Gene-based association testing in the gender-combined dataset revealed eight loci associated with BUA and seven with VOS after correction for multiple testing, with one novel locus for BUA at FAM167A (8p23.1) (P = 1.4 × 10-6). An additional novel locus for BUA was seen in the male-specific analysis at DEFB103B (8p23.1) (P = 1.8 × 10-6). Fracture analysis revealed significant associations between variation at the WNT16 and RSPO3 loci and fracture risk (P = 0.004 and 4.0 × 10-4, respectively). In conclusion, by performing a large GWAS meta-analysis for QUS parameters of bone using a combination of WGS and deeply imputed genotype data, we have identified five novel genetic loci associated with BUA.

Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common- and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.

HABP2 germline variants are uncommon in familial nonmedullary thyroid cancer

Background

The genetic basis of nonsyndromic familial nonmedullary thyroid cancer (FNMTC) is poorly understood. A recent study identified HABP2 as a tumor suppressor gene and identified a germline variant (G534E) in an extended FNMTC kindred. The relevance of this to other FNMTC kindreds is uncertain.

Methods

Sanger sequencing was performed on peripheral blood DNA from probands from 37 Australian FNMTC kindreds to detect the G534E variant. Whole exome data from 59 participants from 20 kindreds were examined for mutations in HABP2 and the thyroid cancer susceptibility genes SRGAP1, NKX2-1, SRRM2 and FOXE1. The population prevalence of the G534E variant in HABP2 was examined in two independent cohorts.

Results

Heterozygosity for the G534E variant in HABP2 was found in 1 of 37 probands (2.7 %), but did not cosegregate with disease in this kindred, being absent in the proband’s affected sister. From whole exome data, pathogenic mutations were not identified in HABP2, SRGAP1, NKX2-1, SRRM2 or FOXE1. Heterozygosity for the G534E variant in HABP2 was present in 7.6 % of Busselton Health Study participants (N = 4634, unknown disease status) and 9.3 % of TwinsUK participants (N = 1195, no history of thyroid cancer).

Conclusions

The G534E variant in HABP2 does not account for the familial nature of NMTC in Australian kindreds, and is common in the general population. Further research is required to elucidate the genetic basis of nonsyndromic FNMTC.

DNA methylation and the social gradient of osteoporotic fracture: A conceptual model

INTRODUCTION

Although there is a documented social gradient for osteoporosis, the underlying mechanism(s) for that gradient remain unknown. We propose a conceptual model based upon the allostatic load theory, to suggest how DNA methylation (DNAm) might underpin the social gradient in osteoporosis and fracture. We hypothesise that social disadvantage is associated with priming of inflammatory pathways mediated by epigenetic modification that leads to an enhanced state of inflammatory reactivity and oxidative stress, and thus places socially disadvantaged individuals at greater risk of osteoporotic fracture.

METHODS/RESULTS

Based on a review of the literature, we present a conceptual model in which social disadvantage increases stress throughout the lifespan, and engenders a proinflammatory epigenetic signature, leading to a heightened inflammatory state that increases risk for osteoporotic fracture in disadvantaged groups that are chronically stressed.

CONCLUSIONS

Our model proposes that, in addition to the direct biological effects exerted on bone by factors such as physical activity and nutrition, the recognised socially patterned risk factors for osteoporosis also act via epigenetic-mediated dysregulation of inflammation. DNAm is a dynamic modulator of gene expression with considerable relevance to the field of osteoporosis. Elucidating the extent to which this epigenetic mechanism transduces the psycho-social environment to increase the risk of osteoporotic fracture may yield novel entry points for intervention that can be used to reduce individual and population-wide risks for osteoporotic fracture. Specifically, an epigenetic evidence-base may strengthen the importance of lifestyle modification and stress reduction programs, and help to reduce health inequities across social groups.

MINI ABSTRACT

Our conceptual model proposes how DNA methylation might underpin the social gradient in osteoporotic fracture. We suggest that social disadvantage is associated with priming of inflammatory signalling pathways, which is mediated by epigenetic modifications, leading to a chronically heightened inflammatory state that places disadvantaged individuals at greater risk of osteoporosis.

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⁻⁰⁹). 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⁻²³). 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⁻⁰⁶, 2.0 × 10⁻⁰⁶ and 2.0 × 10⁻⁰⁶ 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.

Common genetic variants associated with thyroid function may be risk alleles for Hashimoto's disease and Graves' disease

BACKGROUND

Recent studies have identified common genetic variants associated with TSH, free T4 and thyroid peroxidase antibodies, but it is unclear whether these differ between patients with Hashimoto's disease and Graves' disease.

OBJECTIVE

To examine whether 11 common genetic variants differ between Graves' disease and Hashimoto's disease.

PATIENTS AND MEASUREMENTS

We genotyped 11 common variants in a discovery cohort of 203 Australian patients with autoimmune thyroid disease (AITD). Two variants with significant or suggestive associations were analysed in a replication cohort of 384 Danish patients.

RESULTS

For rs753760 (PDE10A), the minor allele frequency in Graves' disease and Hashimoto's disease was 0·38 vs. 0·23, respectively, (P = 6·42 × 10^-4 ) in the discovery cohort, 0·29 vs. 0·24 (P = 0·147) in the replication cohort and 0·32 vs. 0·24 in combined analysis (P = 0·0021; all analyses adjusted for sex). In healthy controls from Busselton, the frequency was 0·29, significantly different from Hashimoto's disease but not Graves' disease. For rs4889009 (MAF gene region), the frequency of the minor G-allele in Graves' disease and Hashimoto's disease was 0·48 vs. 0·36 (P = 0·0156) in the discovery cohort, 0·48 vs. 0·34 (P = 1·83 × 10^-4 ) in the replication cohort and 0·48 vs. 0·35 in the combined analysis (P = 7·53 × 10^-6 ); in controls, the frequency was 0·38, significantly different from Graves' disease but not Hashimoto's disease. After further adjustment for smoking, associations with rs4889009 remained significant, whereas those with rs753760 were not.

CONCLUSION

Common variants in PDE10A and MAF gene regions may influence whether patients with AITD develop Graves' disease or Hashimoto's disease.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

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⁻⁸), with minor allele frequencies of 1.3–23.9%. Novel signals included variants for progesterone (P=7.68 × 10⁻¹²), oestradiol (P=1.63 × 10⁻⁸) and FAI (P=1.50 × 10⁻⁸). A genetic variant near the FSHB gene was identified which influenced both FSH (P=1.74 × 10⁻⁸) and LH (P=3.94 × 10⁻⁹) levels. A separate locus on chromosome 7 was associated with both DHEAS (P=1.82 × 10⁻¹⁴) and progesterone (P=6.09 × 10⁻¹⁴). This study highlights loci that are relevant to reproductive function and suggests overlap in the genetic basis of hormone regulation.

The management of the surgical complications of HELLP syndrome

INTRODUCTION

Complications from HELLP (Haemolysis, Elevated Liver enzymes and Low Platelet) syndrome may present as an emergency to any surgeon. We review the ten-year experience of a tertiary hepatobiliary centre managing HELLP patients. Three selected cases are described to highlight our management strategy and a systematic review of the recent literature is presented.

METHODS

All patients with HELLP syndrome were identified from a prospectively maintained database and their details collated. Subsequently, a detailed search of PubMed was carried out to identify all case series of HELLP syndrome in the literature in the English language since 1999.

RESULTS

On review of 1,002 cases, 10 patients were identified with surgical complications of HELLP syndrome. Seven of these patients had a significant liver injury. Only three of these required surgical intervention for liver injury although four other patients required surgical intervention for other complications. There was no maternal mortality in this series. Review of the literature identified 49 cases in 31 publications. The management approaches of these patients were compared with ours.

CONCLUSIONS

We have presented a large series of patients with surgical complications resulting from HELLP syndrome managed without maternal mortality. This review has confirmed that haemodynamically stable patients with HELLP syndrome associated hepatic rupture can be conservatively treated successfully. However, in unstable patients, perihepatic packing and transfer to a specialist liver unit is recommended.

Influence of ARHGEF3 and RHOA knockdown on ACTA2 and other genes in osteoblasts and osteoclasts

Osteoporosis is a common bone disease that has a strong genetic component. Genome-wide linkage studies have identified the chromosomal region 3p14-p22 as a quantitative trait locus for bone mineral density (BMD). We have previously identified associations between variation in two related genes located in 3p14-p22, ARHGEF3 and RHOA, and BMD in women. In this study we performed knockdown of these genes using small interfering RNA (siRNA) in human osteoblast-like and osteoclast-like cells in culture, with subsequent microarray analysis to identify genes differentially regulated from a list of 264 candidate genes. Validation of selected findings was then carried out in additional human cell lines/cultures using quantitative real-time PCR (qRT-PCR). The qRT-PCR results showed significant down-regulation of the ACTA2 gene, encoding the cytoskeletal protein alpha 2 actin, in response to RHOA knockdown in both osteoblast-like (P<0.001) and osteoclast-like cells (P = 0.002). RHOA knockdown also caused up-regulation of the PTH1R gene, encoding the parathyroid hormone 1 receptor, in Saos-2 osteoblast-like cells (P<0.001). Other findings included down-regulation of the TNFRSF11B gene, encoding osteoprotegerin, in response to ARHGEF3 knockdown in the Saos-2 and hFOB 1.19 osteoblast-like cells (P = 0.003–0.02), and down-regulation of ARHGDIA, encoding the Rho GDP dissociation inhibitor alpha, in response to RHOA knockdown in osteoclast-like cells (P<0.001). These studies identify ARHGEF3 and RHOA as potential regulators of a number of genes in bone cells, including TNFRSF11B, ARHGDIA, PTH1R and ACTA2, with influences on the latter evident in both osteoblast-like and osteoclast-like cells. This adds further evidence to previous studies suggesting a role for the ARHGEF3 and RHOA genes in bone metabolism.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Conditional testing of multiple variants associated with bone mineral density in the FLNB gene region suggests that they represent a single association signal

BACKGROUND

Low bone mineral density (BMD) is a primary risk factor for osteoporosis and is a highly heritable trait, but appears to be influenced by many genes. Genome-wide linkage studies have highlighted the chromosomal region 3p14-p22 as a quantitative trait locus for BMD (LOD 1.1 - 3.5). The FLNB gene, which is thought to have a role in cytoskeletal actin dynamics, is located within this chromosomal region and presents as a strong candidate for BMD regulation. We have previously identified significant associations between four SNPs in the FLNB gene and BMD in women. We have also previously identified associations between five SNPs located 5' of the transcription start site (TSS) and in intron 1 of the FLNB gene and expression of FLNB mRNA in osteoblasts in vitro. The latter five SNPs were genotyped in this study to test for association with BMD parameters in a family-based population of 769 Caucasian women.

RESULTS

Using FBAT, significant associations were seen for femoral neck BMD Z-score with the SNPs rs11720285, rs11130605 and rs9809315 (P = 0.004 - 0.043). These three SNPs were also found to be significantly associated with total hip BMD Z-score (P = 0.014 - 0.026). We then combined the genotype data for these three SNPs with the four SNPs we previously identified as associated with BMD and performed a conditional analysis to determine whether they represent multiple independent associations with BMD. The results from this analysis suggested that these variants represent a single association signal.

CONCLUSIONS

The SNPs identified in our studies as associated with BMD appear to be part of a single association signal between the FLNB gene and BMD in our data. FLNB is one of several genes located in 3p14-p22 that has been identified as significantly associated with BMD in Caucasian women.

Metabolomic markers reveal novel pathways of ageing and early development in human populations

BACKGROUND

Human ageing is a complex, multifactorial process and early developmental factors affect health outcomes in old age.

METHODS

Metabolomic profiling on fasting blood was carried out in 6055 individuals from the UK. Stepwise regression was performed to identify a panel of independent metabolites which could be used as a surrogate for age. We also investigated the association with birthweight overall and within identical discordant twins and with genome-wide methylation levels.

RESULTS

We identified a panel of 22 metabolites which combined are strongly correlated with age (R(2) = 59%) and with age-related clinical traits independently of age. One particular metabolite, C-glycosyl tryptophan (C-glyTrp), correlated strongly with age (beta = 0.03, SE = 0.001, P = 7.0 × 10(-157)) and lung function (FEV1 beta = -0.04, SE = 0.008, P = 1.8 × 10(-8) adjusted for age and confounders) and was replicated in an independent population (n = 887). C-glyTrp was also associated with bone mineral density (beta = -0.01, SE = 0.002, P = 1.9 × 10(-6)) and birthweight (beta = -0.06, SE = 0.01, P = 2.5 × 10(-9)). The difference in C-glyTrp levels explained 9.4% of the variance in the difference in birthweight between monozygotic twins. An epigenome-wide association study in 172 individuals identified three CpG-sites, associated with levels of C-glyTrp (P < 2 × 10(-6)). We replicated one CpG site in the promoter of the WDR85 gene in an independent sample of 350 individuals (beta = -0.20, SE = 0.04, P = 2.9 × 10(-8)). WDR85 is a regulator of translation elongation factor 2, essential for protein synthesis in eukaryotes.

CONCLUSIONS

Our data illustrate how metabolomic profiling linked with epigenetic studies can identify some key molecular mechanisms potentially determined in early development that produce long-term physiological changes influencing human health and ageing.

A meta-analysis of thyroid-related traits reveals novel loci and gender-specific differences in the regulation of thyroid function

Thyroid hormone is essential for normal metabolism and development, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over their life span. In addition, even mild alterations in thyroid function are associated with weight changes, atrial fibrillation, osteoporosis, and psychiatric disorders. To identify novel variants underlying thyroid function, we performed a large meta-analysis of genome-wide association studies for serum levels of the highly heritable thyroid function markers TSH and FT4, in up to 26,420 and 17,520 euthyroid subjects, respectively. Here we report 26 independent associations, including several novel loci for TSH (PDE10A, VEGFA, IGFBP5, NFIA, SOX9, PRDM11, FGF7, INSR, ABO, MIR1179, NRG1, MBIP, ITPK1, SASH1, GLIS3) and FT4 (LHX3, FOXE1, AADAT, NETO1/FBXO15, LPCAT2/CAPNS2). Notably, only limited overlap was detected between TSH and FT4 associated signals, in spite of the feedback regulation of their circulating levels by the hypothalamic-pituitary-thyroid axis. Five of the reported loci (PDE8B, PDE10A, MAF/LOC440389, NETO1/FBXO15, and LPCAT2/CAPNS2) show strong gender-specific differences, which offer clues for the known sexual dimorphism in thyroid function and related pathologies. Importantly, the TSH-associated loci contribute not only to variation within the normal range, but also to TSH values outside the reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings explain, respectively, 5.64% and 2.30% of total TSH and FT4 trait variance, and they improve the current knowledge of the regulation of hypothalamic-pituitary-thyroid axis function and the consequences of genetic variation for hypo- or hyperthyroidism.

Levels of thyroid hormones are tightly regulated by TSH produced in the pituitary, and even mild alterations in their concentrations are strong indicators of thyroid pathologies, which are very common worldwide. To identify common genetic variants associated with the highly heritable markers of thyroid function, TSH and FT4, we conducted a meta-analysis of genome-wide association studies in 26,420 and 17,520 individuals, respectively, of European ancestry with normal thyroid function. Our analysis identified 26 independent genetic variants regulating these traits, several of which are new, and confirmed previously detected polymorphisms affecting TSH (within the PDE8B gene and near CAPZB, MAF/LOC440389, and NR3C2) and FT4 (within DIO1) levels. Gender-specific differences in the genetic effects of several variants for TSH and FT4 levels were identified at several loci, which offer clues to understand the known sexual dimorphism in thyroid function and pathology. Of particular clinical interest, we show that TSH-associated loci contribute not only to normal variation, but also to TSH values outside reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings add to the developing landscape of the regulation of thyroid homeostasis and the consequences of genetic variation for thyroid related diseases.

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.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)<P<5.9 × 10(-4)) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture.