Search

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Objective: Genetic factors underpin the narrow intraindividual variability of thyroid function, although precise contributions of environmental vs genetic factors remain uncert ain. We sought to clarify the heritability of thyroid function traits and thyroid peroxidase antibody (TPOAb) positiv ity 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 y ears) from the Brisbane Longitudinal Twin Study (BLTS) was assessed using structural equation modelling. A geno me-wide analysis was conducted on 2832 of these individuals across 7 522 526 SNPs as well as gene-based associa tion analyses. Replication analysis of the association results was performed in the Raine Study (n = 1115) followed by meta-analysis to maximise power for discover y. 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 associat ion 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 es timates of thyroid function traits and TPOAb positivity. GWAS and gene-based association analysis iden tified variants accounting for a component of this heritability. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract: 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract: 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. [Copyright &y& Elsevier]

Abstract: Bone mineral density (BMD) is a complex trait that is the single best predictor of the risk of osteoporotic fractures. Candidate gene and genome-wide association studies have identified genetic variations in approximately 30 genetic loci associated with BMD variation in humans. α-Actinin-3 (ACTN3) is highly expressed in fast skeletal muscle fibres. There is a common null-polymorphism R577X in human ACTN3 that results in complete deficiency of the α-actinin-3 protein in approximately 20% of Eurasians. Absence of α-actinin-3 does not cause any disease phenotypes in muscle because of compensation by α-actinin-2. However, α-actinin-3 deficiency has been shown to be detrimental to athletic sprint/power performance. In this report we reveal additional functions for α-actinin-3 in bone. α-Actinin-3 but not α-actinin-2 is expressed in osteoblasts. The Actn3 −/− mouse displays significantly reduced bone mass, with reduced cortical bone volume (−14%) and trabecular number (−61%) seen by microCT. Dynamic histomorphometry indicated this was due to a reduction in bone formation. In a cohort of postmenopausal Australian women, ACTN3 577XX genotype was associated with lower BMD in an additive genetic model, with the R577X genotype contributing 1.1% of the variance in BMD. Microarray analysis of cultured osteoprogenitors from Actn3 −/− mice showed alterations in expression of several genes regulating bone mass and osteoblast/osteoclast activity, including Enpp1, Opg and Wnt7b. Our studies suggest that ACTN3 likely contributes to the regulation of bone mass through alterations in bone turnover. Given the high frequency of R577X in the general population, the potential role of ACTN3 R577X as a factor influencing variations in BMD in elderly humans warrants further study. [Copyright &y& Elsevier]

Abstract: Osteoporosis is a highly heritable trait that appears to be influenced by multiple genes. Genome-wide linkage studies have highlighted the chromosomal region 3p14-p21 as a quantitative trait locus for BMD. We have previously published evidence suggesting that the ARHGEF3 gene from this region is associated with BMD in women. The product of this gene activates the RHOA GTPase, the gene for which is also located within this region. The aim of this study was to evaluate the influence of genetic polymorphism in RHOA on bone density in women. Sequence variation within the RHOA gene region was determined using 9 single nucleotide polymorphisms (SNPs) in a discovery cohort of 769 female sibs. Of the 9 SNPs, one was found to be monomorphic with the others representing 3 distinct linkage disequilibrium (LD) blocks. Using FBAT software, significant associations were found between two of these LD blocks and BMD Z-score of the spine and hip (P =0.001–0.036). The LD block tagged by the SNP rs17595772 showed maximal association, with the more common G allele at rs17595772 associated with decreased BMD Z-score. Genotyping for rs17595772 in a replication cohort of 780 postmenopausal women confirmed an association with BMD Z-score (P =0.002–0.036). Again, the G allele was found to be associated with a reduced hip and spine BMD Z-score. These results support the implication of the RhoGTPase-RhoGEF pathway in osteoporosis, and suggest that one or more genes in this pathway may be responsible for the linkage observed between 3p14-p21 and BMD. [Copyright &y& Elsevier]

Osteoarthritis (OA), the most prevalent form of arthritis in the elderly, is characterized by the degradation of articular cartilage and has a strong genetic component. Our aim was to identify genetic variants involved in risk of knee OA in women. A pooled genome-wide association scan with the Illumina550 Duo array was performed in 255 controls and 387 cases. Twenty-eight variants with p < 1 × 10-5 were estimated to have probabilities of being false positives ≤0.5 and were genotyped individually in the original samples and in replication cohorts from the UK and the U.S. (599 and 272 cases, 1530 and 258 controls, respectively). The top seven associations were subsequently tested in samples from the Netherlands (306 cases and 584 controls). rs4140564 on chromosome 1 mapping 5′ to both the PTGS2 and PLA2G4A genes was associated with risk of knee OA in all the cohorts studied (overall odds ratio ORMH = 1.55 95% C.I. 1.30-1.85, p < 6.9 × 10-7). Differential allelic expression analysis of PTGS2 with mRNA extracted from the cartilage of joint-replacement surgery OA patients revealed a significant difference in allelic expression (p < 1.0 × 10-6). These results suggest the existence of cis-acting regulatory polymorphisms that are in, or near to, PTGS2 and in modest linkage disequilibrium with rs4140564. Our results and previous studies on the role of the cyclooxygenase 2 enzyme encoded by PTGS2 underscore the importance of this signaling pathway in the pathogenesis of knee OA. [ABSTRACT FROM AUTHOR]

Osteoporosis is a common and debilitating bone disease characterized by low bone mineral density (BMD), a highly heritable and polygenic trait. Genome-wide linkage studies have identified 3p14-p21 as a quantitative trait locus for BMD. The ARHGEF3 gene is situated within this region and was identified as a strong positional candidate. The aim of this study was to evaluate the role of variation in ARHGEF3 on bone density in women. Sequence variation within ARHGEF3 was analyzed with 17 single-nucleotide polymorphisms (SNPs) in a discovery cohort of 769 female sibs. Significant associations were found with family-based association tests between five SNPs and various measures of age-adjusted BMD (p = 0.0007-0.041) with rs7646054 showing maximal association. Analysis of the data with QPDTPHASE suggested that the more common G allele at rs7646054 is associated with decreased age-adjusted BMD. Significant associations were also demonstrated between 3-SNP haplotypes and age-adjusted spine and femoral-neck BMD (p = 0.002 and 0.003, respectively). rs7646054 was then genotyped in a replication cohort, and significant associations with hip and spine BMD were confirmed (p = 0.003-0.038), as well as an association with fracture rate (p = 0.02). Again, the G allele was associated with a decrease in age-adjusted BMD at each site studied. In conclusion, genetic variation in ARHGEF3 plays a role in the determination of bone density in Caucasian women. This data implicates the RhoGTPase-RhoGEF pathway in osteoporosis. [ABSTRACT FROM AUTHOR]

Paget's disease of bone (PDB) is characterised by focal abnormalities of bone remodelling, with increased osteoclastic resorption the primary feature of the disease. Genetic factors have been shown to play an important role in PDB, and genome-wide association studies (GWAS) have identified 7 genetic loci as associated with PDB at the genome-wide level. Expression quantitative trait locus (eQTL) studies using cell types that are directly relevant to the disease of interest are increasingly being used to identify putative effector genes for GWAS loci. We have recently constructed a unique osteoclast-specific eQTL resource using cells differentiated in vitro from 158 subjects for study of the genetics of bone disease. Considering the major role osteoclasts have in PDB, we used this resource to investigate potential genetic regulatory effects for the 7 PDB genome-wide significant loci on genes located within 500 kb of each locus. After correction for multiple testing, we observed statistically significant associations for rs4294134 with expression of the gene STMP1, and rs2458413 with expression of the genes DPYS and DCSTAMP. The eQTL associations observed for rs4294134 with STMP1, and rs2458413 with DCSTAMP were further supported by eQTL data from other tissue types. The product of the STMP1 gene has not been extensively studied, however the DCSTAMP gene has an established role in osteoclast differentiation and the associations seen between rs2458413 and PDB are likely mediated through regulatory effects on this gene. This study highlights the value of eQTL data in determining which genes are relevant to GWAS loci. [ABSTRACT FROM AUTHOR]

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., Competing Interests: Declaration of interests The authors have no conflicts to declare., (Copyright © 2024 The Author(s). Published by Elsevier B.V. All rights reserved.)

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Endocrine Society.)

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/m 2 ). 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., (© 2024. Crown.)

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., (© 2024. The Author(s).)

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., (© 2023. The Author(s).)

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., (© 2023. The Author(s).)

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.

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.

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Endocrine Society. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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.

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.

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Endocrine Society.)

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., (© 2019 by the American Diabetes Association.)

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 (r 2  > 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., (© 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.)

Polycystic ovary syndrome (PCOS) is a disorder characterized by hyperandrogenism, ovulatory dysfunction and polycystic ovarian morphology. Affected women frequently have metabolic disturbances including insulin resistance and dysregulation of glucose homeostasis. PCOS is diagnosed with two different sets of diagnostic criteria, resulting in a phenotypic spectrum of PCOS cases. The genetic similarities between cases diagnosed based on the two criteria have been largely unknown. Previous studies in Chinese and European subjects have identified 16 loci associated with risk of PCOS. We report a fixed-effect, inverse-weighted-variance meta-analysis from 10,074 PCOS cases and 103,164 controls of European ancestry and characterisation of PCOS related traits. We identified 3 novel loci (near PLGRKT, ZBTB16 and MAPRE1), and provide replication of 11 previously reported loci. Only one locus differed significantly in its association by diagnostic criteria; otherwise the genetic architecture was similar between PCOS diagnosed by self-report and PCOS diagnosed by NIH or non-NIH Rotterdam criteria across common variants at 13 loci. Identified variants were associated with hyperandrogenism, gonadotropin regulation and testosterone levels in affected women. Linkage disequilibrium score regression analysis revealed genetic correlations with obesity, fasting insulin, type 2 diabetes, lipid levels and coronary artery disease, indicating shared genetic architecture between metabolic traits and PCOS. Mendelian randomization analyses suggested variants associated with body mass index, fasting insulin, menopause timing, depression and male-pattern balding play a causal role in PCOS. The data thus demonstrate 3 novel loci associated with PCOS and similar genetic architecture for all diagnostic criteria. The data also provide the first genetic evidence for a male phenotype for PCOS and a causal link to depression, a previously hypothesized comorbid disease. Thus, the genetics provide a comprehensive view of PCOS that encompasses multiple diagnostic criteria, gender, reproductive potential and mental health., Competing Interests: Members of the 23andMe Research team are employees of and hold stock or stock options in 23andMe, Inc. GT, UT, KS, and US are employees of deCODE genetics/Amgen Inc. MIM serves on advisory panels for Pfizer and NovoNordisk. MIM has received honoraria from Pfizer, NovoNordisk and EliLilly, and has received research funding from Pfizer, NovoNordisk, EliLilly, AstraZeneca, Sanofi Aventis, Boehringer Ingelheim, Merck, Roche, Janssen, Takeda, and Servier. JL has received consultancy fees from Danone, Metagenics inc., Titus Healthcare, Roche and Euroscreen. CW is a consultant for Novartis and has received UptoDate royalties.

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves' disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets.

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., (Copyright © 2017 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

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.

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.

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.

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.

Dehydroepiandrosterone sulphate (DHEAS) is the most abundant circulating steroid secreted by adrenal glands--yet its function is unknown. Its serum concentration declines significantly with increasing age, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. We conducted a meta-analysis of genome-wide association data with 14,846 individuals and identified eight independent common SNPs associated with serum DHEAS concentrations. Genes at or near the identified loci include ZKSCAN5 (rs11761528; p = 3.15 × 10(-36)), SULT2A1 (rs2637125; p =  2.61 × 10(-19)), ARPC1A (rs740160; p =  1.56 × 10(-16)), TRIM4 (rs17277546; p =  4.50 × 10(-11)), BMF (rs7181230; p = 5.44 × 10(-11)), HHEX (rs2497306; p =  4.64 × 10(-9)), BCL2L11 (rs6738028; p = 1.72 × 10(-8)), and CYP2C9 (rs2185570; p = 2.29 × 10(-8)). These genes are associated with type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins. Several SNPs were associated with changes in gene expression levels, and the related genes are connected to biological pathways linking DHEAS with ageing. This study provides much needed insight into the function of DHEAS., Competing Interests: The authors have declared that no competing interests exist.

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., (2010 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Bone mineral density (BMD) is a heritable complex trait used in the clinical diagnosis of osteoporosis and the assessment of fracture risk. We performed meta-analysis of five genome-wide association studies of femoral neck and lumbar spine BMD in 19,195 subjects of Northern European descent. We identified 20 BMD loci that reached genome-wide significance (GWS; P < 5 x 10(-8)), of which 13 map to regions not previously associated with this trait: 1p31.3 (GPR177), 2p21 (SPTBN1), 3p22 (CTNNB1), 4q21.1 (MEPE), 5q14 (MEF2C), 7p14 (STARD3NL), 7q21.3 (FLJ42280), 11p11.2 (LRP4, ARHGAP1, F2), 11p14.1 (DCDC5), 11p15 (SOX6), 16q24 (FOXL1), 17q21 (HDAC5) and 17q12 (CRHR1). The meta-analysis also confirmed at GWS level seven known BMD loci on 1p36 (ZBTB40), 6q25 (ESR1), 8q24 (TNFRSF11B), 11q13.4 (LRP5), 12q13 (SP7), 13q14 (TNFSF11) and 18q21 (TNFRSF11A). The many SNPs associated with BMD map to genes in signaling pathways with relevance to bone metabolism and highlight the complex genetic architecture that underlies osteoporosis and variation in BMD.

Hip structural analysis (HSA) has been developed over 20 yr, applied extensively in research, and has demonstrated useful outcomes associating bone structural geometry with bone fragility (research-HSA or r-HSA). In 2007, Hologic Inc. (Bedford, MA) incorporated HSA with some modifications as an option for Hologic dual-energy X-ray absorptiometry (DXA) scanners (clinical HSA or c-HSA). This brought HSA from the research environment into the clinical environment. This article reports a comparison of r-HSA and c-HSA implementations using DXA scans from a group of 191 females. Bland-Altman plots at the narrow-neck (NN) HSA region indicated higher r-HSA areal bone mineral density (mean difference: 0.27 g/cm(2); 21.7% [of mean]); cross-sectional area (0.63 cm(2); 18.7%); cross-sectional moment of inertia (0.26 cm(4); 11.1%), and section modulus (0.22 cm(3); 14.5%) compared with c-HSA. The converse was observed for NN subperiosteal width (-0.09 cm; 3.1%). High linear correlations (r(2) > 0.81) were found between r-HSA and c-HSA NN structural geometric outcomes, with the exception of neck shaft angle (r(2) > 0.47). As differences were significant (p < 0.001), slopes and intercepts are provided to enable linear transformations from r-HSA to corresponding c-HSA structural geometric data.