Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche

Pubertal development and regulation

Puberty marks the end of childhood and is a period when individuals undergo physiological and psychological changes to achieve sexual maturation and fertility. The hypothalamic-pituitary-gonadal axis controls puberty and reproduction and is tightly regulated by a complex network of excitatory and inhibitory factors. This axis is active in the embryonic and early postnatal stages of life and is subsequently restrained during childhood, and its reactivation culminates in puberty initiation. The mechanisms underlying this reactivation are not completely known. The age of puberty onset varies between individuals and the timing of puberty initiation is associated with several health outcomes in adult life. In this Series paper, we discuss pubertal markers, epidemiological trends of puberty initiation over time, and the mechanisms whereby genetic, metabolic, and other factors control secretion of gonadotropin-releasing hormone to determine initiation of puberty.

Significant impact of miRNA-target gene networks on genetics of human complex traits

The impact of microRNA (miRNA) on the genetics of human complex traits, especially in the context of miRNA-target gene networks, has not been fully assessed. Here, we developed a novel analytical method, MIGWAS, to comprehensively evaluate enrichment of genome-wide association study (GWAS) signals in miRNA–target gene networks. We applied the method to the GWAS results of the 18 human complex traits from >1.75 million subjects, and identified significant enrichment in rheumatoid arthritis (RA), kidney function, and adult height (P < 0.05/18 = 0.0028, most significant enrichment in RA with P = 1.7 × 10⁻⁴). Interestingly, these results were consistent with current literature-based knowledge of the traits on miRNA obtained through the NCBI PubMed database search (adjusted P = 0.024). Our method provided a list of miRNA and target gene pairs with excess genetic association signals, part of which included drug target genes. We identified a miRNA (miR-4728-5p) that downregulates PADI2, a novel RA risk gene considered as a promising therapeutic target (rs761426, adjusted P = 2.3 × 10⁻⁹). Our study indicated the significant impact of miRNA–target gene networks on the genetics of human complex traits, and provided resources which should contribute to drug discovery and nucleic acid medicine.

Sex-specific regulation of weight and puberty by the Lin28/let-7 axis

Growth and pubertal timing differ in boys and girls. Variants in/near LIN28B associate with age at menarche (AAM) in genome-wide association studies and some AAM-related variants associate with growth in a sex-specific manner. Sex-specific growth patterns in response to Lin28b perturbation have been detected in mice, and overexpression of Lin28a has been shown to alter pubertal timing in female mice. To investigate further how Lin28a and Lin28b affect growth and puberty in both males and females, we evaluated Lin28b loss-of-function (LOF) mice and Lin28a gain-of-function (GOF) mice. Because both Lin28a and Lin28b can act via the conserved microRNA let-7, we also examined let-7 GOF mice. As reported previously, Lin28b LOF led to lighter body weights only in male mice while Lin28a GOF yielded heavier mice of both sexes. Let-7 GOF mice weighed less than controls, and males were more affected than females. Timing of puberty was assessed by vaginal opening (VO) and preputial separation (PS). Male Lin28b LOF and male let-7 GOF, but not female, mice displayed alteration of pubertal timing, with later PS than controls. In contrast, both male and female Lin28a GOF mice displayed late onset of puberty. Together, these data point toward a complex system of regulation by Lin28a, Lin28b, and let-7, in which Lin28b and let-7 can impact both puberty and growth in a sex-specific manner, raising the possibility that this pathway may contribute to differential regulation of male and female growth and puberty in humans.

New loci for body fat percentage reveal link between adiposity and cardiometabolic disease risk

To increase our understanding of the genetic basis of adiposity and its links to cardiometabolic disease risk, we conducted a genome-wide association meta-analysis of body fat percentage (BF%) in up to 100,716 individuals. Twelve loci reached genome-wide significance (P<5 × 10(-8)), of which eight were previously associated with increased overall adiposity (BMI, BF%) and four (in or near COBLL1/GRB14, IGF2BP1, PLA2G6, CRTC1) were novel associations with BF%. Seven loci showed a larger effect on BF% than on BMI, suggestive of a primary association with adiposity, while five loci showed larger effects on BMI than on BF%, suggesting association with both fat and lean mass. In particular, the loci more strongly associated with BF% showed distinct cross-phenotype association signatures with a range of cardiometabolic traits revealing new insights in the link between adiposity and disease risk.

The genetics of pubertal timing in the general population: recent advances and evidence for sex-specificity

PURPOSE OF REVIEW

This article overviews advances in the genetics of puberty based on studies in the general population, describes evidence for sex-specific genetic effects on pubertal timing, and briefly reviews possible mechanisms mediating sexually dimorphic genetic effects.

RECENT FINDINGS

Pubertal timing is highly polygenic, and many loci are conserved among ethnicities. A number of identified loci underlie both pubertal timing and related traits such as height and BMI. It is increasingly apparent that understanding the factors modulating the onset of puberty is important because the timing of this developmental stage is associated with a wider range of adult health outcomes than previously appreciated. Although most of the genetic effects underlying the timing of puberty are common between boys and girls, some effects show sex-specificity and many are epigenetically modulated. Several potential mechanisms, including hormone-independent ones, may be responsible for observed sex differences.

SUMMARY

Studies of pubertal timing in the general population have provided new knowledge about the genetic architecture of this complex trait. Increasing attention paid to sex-specific effects may provide key insights into the sexual dimorphism in pubertal timing and even into the associations between puberty and adult health risks by identifying common underlying biological pathways.

Withdrawn: Discovering Genes Essential to the Hypothalamic Regulation of Human Reproduction Using a Human Disease Model: Adjusting to Life in the "-Omics" Era

The neuroendocrine regulation of reproduction is an intricate process requiring the exquisite coordination of an assortment of cellular networks, all converging on the GnRH neurons. These neurons have a complex life history, migrating mainly from the olfactory placode into the hypothalamus, where GnRH is secreted and acts as the master regulator of the hypothalamic-pituitary-gonadal axis. Much of what we know about the biology of the GnRH neurons has been aided by discoveries made using the human disease model of isolated GnRH deficiency (IGD), a family of rare Mendelian disorders that share a common failure of secretion and/or action of GnRH causing hypogonadotropic hypogonadism. Over the last 30 years, research groups around the world have been investigating the genetic basis of IGD using different strategies based on complex cases that harbor structural abnormalities or single pleiotropic genes, endogamous pedigrees, candidate gene approaches as well as pathway gene analyses. Although such traditional approaches, based on well-validated tools, have been critical to establish the field, new strategies, such as next-generation sequencing, are now providing speed and robustness, but also revealing a surprising number of variants in known IGD genes in both patients and healthy controls. Thus, before the field moves forward with new genetic tools and continues discovery efforts, we must reassess what we know about IGD genetics and prepare to hold our work to a different standard. The purpose of this review is to: 1) look back at the strategies used to discover the "known" genes implicated in the rare forms of IGD; 2) examine the strengths and weaknesses of the methodologies used to validate genetic variation; 3)substantiate the role of known genes in the pathophysiology of the disease; and 4) project forward as we embark upon a widening use of these new and powerful technologies for gene discovery. (Endocrine Reviews 36: 603-621, 2015).

A two-phase procedure for non-normal quantitative trait genetic association study

BACKGROUND

The nonparametric trend test (NPT) is well suitable for identifying the genetic variants associated with quantitative traits when the trait values do not satisfy the normal distribution assumption. If the genetic model, defined according to the mode of inheritance, is known, the NPT derived under the given genetic model is optimal. However, in practice, the genetic model is often unknown beforehand. The NPT derived from an uncorrected model might result in loss of power. When the underlying genetic model is unknown, a robust test is preferred to maintain satisfactory power.

RESULTS

We propose a two-phase procedure to handle the uncertainty of the genetic model for non-normal quantitative trait genetic association study. First, a model selection procedure is employed to help choose the genetic model. Then the optimal test derived under the selected model is constructed to test for possible association. To control the type I error rate, we derive the joint distribution of the test statistics developed in the two phases and obtain the proper size.

CONCLUSIONS

The proposed method is more robust than existing methods through the simulation results and application to gene DNAH9 from the Genetic Analysis Workshop 16 for associated with Anti-cyclic citrullinated peptide antibody further demonstrate its performance.

Genome-wide association study identifies common and low-frequency variants at the AMH gene locus that strongly predict serum AMH levels in males

Anti-Müllerian hormone (AMH) is an essential messenger of sexual differentiation in the foetus and is an emerging biomarker of postnatal reproductive function in females. Due to a paucity of adequately sized studies, the genetic determinants of circulating AMH levels are poorly characterized. In samples from 2815 adolescents aged 15 from the ALSPAC study, we performed the first genome-wide association study of serum AMH levels across a set of ∼9 m '1000 Genomes Reference Panel' imputed genetic variants. Genetic variants at the AMH protein-coding gene showed considerable allelic heterogeneity, with both common variants [rs4807216 (P(Male) = 2 × 10(-49), Beta: ∼0.9 SDs per allele), rs8112524 (P(Male) = 3 × 10(-8), Beta: ∼0.25)] and low-frequency variants [rs2385821 (P(Male) = 6 × 10(-31), Beta: ∼1.2, frequency 3.6%)] independently associated with apparently large effect sizes in males, but not females. For all three SNPs, we highlight mechanistic links to AMH gene function and demonstrate highly significant sex interactions (P(Het) 0.0003-6.3 × 10(-12)), culminating in contrasting estimates of trait variance explained (24.5% in males versus 0.8% in females). Using these SNPs as a genetic proxy for AMH levels, we found no evidence in additional datasets to support a biological role for AMH in complex traits and diseases in men.

Genetic evidence that lower circulating FSH levels lengthen menstrual cycle, increase age at menopause and impact female reproductive health

STUDY QUESTION

How does a genetic variant in the FSHB promoter, known to alter FSH levels, impact female reproductive health?

SUMMARY ANSWER

The T allele of the FSHB promoter polymorphism (rs10835638; c.-211G>T) results in longer menstrual cycles and later menopause and, while having detrimental effects on fertility, is protective against endometriosis.

WHAT IS KNOWN ALREADY

The FSHB promoter polymorphism (rs10835638; c.-211G>T) affects levels of FSHB transcription and, as a result, circulating levels of FSH. FSH is required for normal fertility and genetic variants at the FSHB locus are associated with age at menopause and polycystic ovary syndrome (PCOS).

STUDY DESIGN, SIZE, DURATION

We used cross-sectional data from the UK Biobank to look at associations between the FSHB promoter polymorphism and reproductive traits, and performed a genome-wide association study (GWAS) for length of menstrual cycle.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We included white British individuals aged 40-69 years in 2006-2010, in the May 2015 release of genetic data from UK Biobank. We tested the FSH-lowering T allele of the FSHB promoter polymorphism (rs10835638; c.-211G>T) for associations with 29, mainly female, reproductive phenotypes in up to 63 350 women and 56 608 men. We conducted a GWAS in 9534 individuals to identify genetic variants associated with length of menstrual cycle.

MAIN RESULTS AND THE ROLE OF CHANCE

The FSH-lowering T allele of the FSHB promoter polymorphism (rs10835638; MAF 0.16) was associated with longer menstrual cycles [0.16 SD (c. 1 day) per minor allele; 95% confidence interval (CI) 0.12-0.20; P = 6 × 10(-16)], later age at menopause (0.13 years per minor allele; 95% CI 0.04-0.22; P = 5.7 × 10(-3)), greater female nulliparity [odds ratio (OR) = 1.06; 95% CI 1.02-1.11; P = 4.8 × 10(-3)] and lower risk of endometriosis (OR = 0.79; 95% CI 0.69-0.90; P = 4.1 × 10(-4)). The FSH-lowering T allele was not associated with other female reproductive illnesses or conditions in our study and we did not replicate associations with male infertility or PCOS. In the GWAS for menstrual cycle length, only variants near the FSHB gene reached genome-wide significance (P < 5 × 10(-9)).

LIMITATIONS, REASONS FOR CAUTION

The data included might be affected by recall bias. Cycle length was not available for 25% of women still cycling (1% did not answer, 6% did not know and for 18% cycle length was recorded as 'irregular'). Women with a cycle length recorded were aged over 40 and were approaching menopause; however, we did not find evidence that this affected the results. Many of the groups with illnesses had relatively small sample sizes and so the study may have been under-powered to detect an effect.

WIDER IMPLICATIONS OF THE FINDINGS

We found a strong novel association between a genetic variant that lowers FSH levels and longer menstrual cycles, at a locus previously robustly associated with age at menopause. The variant was also associated with nulliparity and endometriosis risk. These findings should now be verified in a second independent group of patients. We conclude that lifetime differences in circulating levels of FSH between individuals can influence menstrual cycle length and a range of reproductive outcomes, including menopause timing, infertility, endometriosis and PCOS.

STUDY FUNDING/COMPETING INTERESTS

None.

TRIAL REGISTRATION NUMBER

Not applicable.

Puberty, Socioeconomic Status, and Depression in Girls: Evidence for Gene × Environment Interactions

In the current study, we tested for Gene × Environment interactions in the association between pubertal timing and adolescent depression by examining how socioeconomic factors might moderate age at menarche's relation with depressive symptoms. Participants comprised 630 female twin and sibling pairs from the National Longitudinal Study of Adolescent Health. Consistent with previous studies, results showed that genetic predispositions toward later menarche were associated with fewer depressive symptoms and that genetic predispositions toward earlier menarche were associated with more depressive symptoms. However, this pattern was subtle and evident only in girls from higher socioeconomic backgrounds. Although girls from lower socioeconomic families showed the highest overall levels of depression, their symptoms appeared unrelated to timing of physical development through either a genetic or an environmental path.

Epigenetic regulation of puberty via Zinc finger protein-mediated transcriptional repression

In primates, puberty is unleashed by increased GnRH release from the hypothalamus following an interval of juvenile quiescence. GWAS implicates Zinc finger (ZNF) genes in timing human puberty. Here we show that hypothalamic expression of several ZNFs decreased in agonadal male monkeys in association with the pubertal reactivation of gonadotropin secretion. Expression of two of these ZNFs, GATAD1 and ZNF573, also decreases in peripubertal female monkeys. However, only GATAD1 abundance increases when gonadotropin secretion is suppressed during late infancy. Targeted delivery of GATAD1 or ZNF573 to the rat hypothalamus delays puberty by impairing the transition of a transcriptional network from an immature repressive epigenetic configuration to one of activation. GATAD1 represses transcription of two key puberty-related genes, KISS1 and TAC3, directly, and reduces the activating histone mark H3K4me2 at each promoter via recruitment of histone demethylase KDM1A. We conclude that GATAD1 epitomizes a subset of ZNFs involved in epigenetic repression of primate puberty.

Molecular insights into the aetiology of female reproductive ageing

As age at pubertal onset declines and age at first pregnancy increases, the mechanisms that regulate female reproductive lifespan become increasingly relevant to population health. The timing of menarche and menopause can have profound effects not only on fertility but also on the risk of diseases such as type 2 diabetes mellitus, cardiovascular disease and breast cancer. Genetic studies have identified dozens of highly penetrant rare mutations associated with reproductive disorders, and also ∼175 common genetic variants associated with the timing of puberty or menopause. These findings, alongside other functional studies, have highlighted a diverse range of mechanisms involved in reproductive ageing, implicating core biological processes such as cell cycle regulation and energy homeostasis. The aim of this article is to review the contribution of such genetic findings to our understanding of the molecular regulation of reproductive timing, as well as the biological basis of the epidemiological links between reproductive ageing and disease risk.

Discovering Genes Essential to the Hypothalamic Regulation of Human Reproduction Using a Human Disease Model: Adjusting to Life in the "-Omics" Era

The neuroendocrine regulation of reproduction is an intricate process requiring the exquisite coordination of an assortment of cellular networks, all converging on the GnRH neurons. These neurons have a complex life history, migrating mainly from the olfactory placode into the hypothalamus, where GnRH is secreted and acts as the master regulator of the hypothalamic-pituitary-gonadal axis. Much of what we know about the biology of the GnRH neurons has been aided by discoveries made using the human disease model of isolated GnRH deficiency (IGD), a family of rare Mendelian disorders that share a common failure of secretion and/or action of GnRH causing hypogonadotropic hypogonadism. Over the last 30 years, research groups around the world have been investigating the genetic basis of IGD using different strategies based on complex cases that harbor structural abnormalities or single pleiotropic genes, endogamous pedigrees, candidate gene approaches as well as pathway gene analyses. Although such traditional approaches, based on well-validated tools, have been critical to establish the field, new strategies, such as next-generation sequencing, are now providing speed and robustness, but also revealing a surprising number of variants in known IGD genes in both patients and healthy controls. Thus, before the field moves forward with new genetic tools and continues discovery efforts, we must reassess what we know about IGD genetics and prepare to hold our work to a different standard. The purpose of this review is to: 1) look back at the strategies used to discover the "known" genes implicated in the rare forms of IGD; 2) examine the strengths and weaknesses of the methodologies used to validate genetic variation; 3) substantiate the role of known genes in the pathophysiology of the disease; and 4) project forward as we embark upon a widening use of these new and powerful technologies for gene discovery.

The epigenomic landscape of African rainforest hunter-gatherers and farmers

The genetic history of African populations is increasingly well documented, yet their patterns of epigenomic variation remain uncharacterized. Moreover, the relative impacts of DNA sequence variation and temporal changes in lifestyle and habitat on the human epigenome remain unknown. Here we generate genome-wide genotype and DNA methylation profiles for 362 rainforest hunter-gatherers and sedentary farmers. We find that the current habitat and historical lifestyle of a population have similarly critical impacts on the methylome, but the biological functions affected strongly differ. Specifically, methylation variation associated with recent changes in habitat mostly concerns immune and cellular functions, whereas that associated with historical lifestyle affects developmental processes. Furthermore, methylation variation—particularly that correlated with historical lifestyle—shows strong associations with nearby genetic variants that, moreover, are enriched in signals of natural selection. Our work provides new insight into the genetic and environmental factors affecting the epigenomic landscape of human populations over time.

Genetic and environmental factors affect genome-wide patterns of epigenetic variation. Here, the authors show that while current habitat and historical lifestyle impact the methylome of rainforest hunter-gatherers and sedentary farmers, the biological functions affected and the degree of genetic control differ.

Uncovering the Genetic Architectures of Quantitative Traits

The aim of a genome-wide association study (GWAS) is to identify loci in the human genome affecting a phenotype of interest. This review summarizes some recent work on conceptual and methodological aspects of GWAS. The average effect of gene substitution at a given causal site in the genome is the key estimand in GWAS, and we argue for its fundamental importance. Implicit in the definition of average effect is a linear model relating genotype to phenotype. The fraction of the phenotypic variance ascribable to polymorphic sites with nonzero average effects in this linear model is called the heritability, and we describe methods for estimating this quantity from GWAS data. Finally, we show that the theory of compressed sensing can be used to provide a sharp estimate of the sample size required to identify essentially all sites contributing to the heritability of a given phenotype.

Shared genetic aetiology of puberty timing between sexes and with health-related outcomes

Understanding of the genetic regulation of puberty timing has come largely from studies of rare disorders and population-based studies in women. Here, we report the largest genomic analysis for puberty timing in 55,871 men, based on recalled age at voice breaking. Analysis across all genomic variants reveals strong genetic correlation (0.74, P=2.7 × 10⁻⁷⁰) between male and female puberty timing. However, some loci show sex-divergent effects, including directionally opposite effects between sexes at the SIM1/MCHR2 locus (P_{heterogeneity}=1.6 × 10⁻¹²). We find five novel loci for puberty timing (P<5 × 10⁻⁸), in addition to nine signals in men that were previously reported in women. Newly implicated genes include two retinoic acid-related receptors, RORB and RXRA, and two genes reportedly disrupted in rare disorders of puberty, LEPR and KAL1. Finally, we identify genetic correlations that indicate shared aetiologies in both sexes between puberty timing and body mass index, fasting insulin levels, lipid levels, type 2 diabetes and cardiovascular disease.

Past studies on genetics of puberty relied on rare disorders or age of menarche in women. Here, Day et al. examine puberty timing in men by the age of voice breaking, and find some loci with sexually dimorphic effects and genetic architectures shared with other health conditions.

An atlas of genetic correlations across human diseases and traits

Identifying genetic correlations between complex traits and diseases can provide useful etiological insights and help prioritize likely causal relationships. The major challenges preventing estimation of genetic correlation from genome-wide association study (GWAS) data with current methods are the lack of availability of individual-level genotype data and widespread sample overlap among meta-analyses. We circumvent these difficulties by introducing a technique-cross-trait LD Score regression-for estimating genetic correlation that requires only GWAS summary statistics and is not biased by sample overlap. We use this method to estimate 276 genetic correlations among 24 traits. The results include genetic correlations between anorexia nervosa and schizophrenia, anorexia and obesity, and educational attainment and several diseases. These results highlight the power of genome-wide analyses, as there currently are no significantly associated SNPs for anorexia nervosa and only three for educational attainment.

Partitioning heritability by functional annotation using genome-wide association summary statistics

Recent work has demonstrated that some functional categories of the genome contribute disproportionately to the heritability of complex diseases. Here we analyze a broad set of functional elements, including cell type-specific elements, to estimate their polygenic contributions to heritability in genome-wide association studies (GWAS) of 17 complex diseases and traits with an average sample size of 73,599. To enable this analysis, we introduce a new method, stratified LD score regression, for partitioning heritability from GWAS summary statistics while accounting for linked markers. This new method is computationally tractable at very large sample sizes and leverages genome-wide information. Our findings include a large enrichment of heritability in conserved regions across many traits, a very large immunological disease-specific enrichment of heritability in FANTOM5 enhancers and many cell type-specific enrichments, including significant enrichment of central nervous system cell types in the heritability of body mass index, age at menarche, educational attainment and smoking behavior.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Bivariate Genome-Wide Association Study Implicates ATP6V1G1 as a Novel Pleiotropic Locus Underlying Osteoporosis and Age at Menarche

OBJECTIVE

Age at menarche (AAM) is determined by the overall duration of endocrine-tissue sex hormone exposure levels. Osteoporosis, the most common metabolic bone disease, is characterized primarily by reduced bone mineral density (BMD) and an increased risk of low trauma fractures. Bone was an endocrine organ regulating the synthesis and secretion of sex steroid hormones. The mutual dependence between bone and gonads underscore the importance of genetic approaches to identify novel pleiotropic genetic factors coregulating BMD and AAM. In this study, we performed a bivariate genome-wide association study (GWAS) to explore novel ethnic common loci and/or genes that may influence both AAM and BMD.

METHODS

We analyzed genotyping data available for 826 unrelated Chinese subjects using genome-wide human genotyping arrays. After quality control, a total of 702 413 single-nucleotide polymorphisms (SNPs) were tested for association using a bivariate linear regression model. The interesting SNPs were replicated in three independent cohorts including 1728 unrelated Caucasians, 709 African-Americans, and 408 Hispanic-Americans.

RESULTS

We found four SNPs (rs10817638, rs7851259, rs10982287, and rs4979427), located upstream of the ATP6V1G1 gene, were bivariately associated with hip BMD-AAM (P = 4.90 × 10(-7), P = 1.07 × 10(-6), P = 1.28 × 10(-5), and P = 5.42 × 10(-5), respectively). These four SNPs were replicated in African-Americans, with corresponding values of P = 1.95 × 10(-2), P = 3.18 × 10(-2), P = 2.57 × 10(-2), and P = 3.64 × 10(-2), respectively. rs10817638 and rs10982287 were further replicated in Caucasians (P = 1.76 × 10(-2) and P = 9.42 × 10(-3), respectively) and Hispanic-Americans (P = 8.37 × 10(-3) and P = 1.52 × 10(-3), respectively). Meta-analyses yielded stronger association signals for rs10817638 and rs10982287 with combined values of P = 3.02 × 10(-9) and P = 3.49 × 10(-9), respectively.

CONCLUSIONS

Our study implicated ATP6V1G1 as a novel pleiotropic gene underlying variation of both BMD and AAM. The findings enhance our knowledge of genetic associations between BMD and AAM and provide a rationale for subsequent functional studies of these implicated genes in the pathophysiology of diseases/traits, such as osteoporosis and AAM.

Season of birth is associated with birth weight, pubertal timing, adult body size and educational attainment: a UK Biobank study

Season of birth, a marker of in utero vitamin D exposure, has been associated with a wide range of health outcomes. Using a dataset of ∼450,000 participants from the UK Biobank study, we aimed to assess the impact of this seasonality on birth weight, age at menarche, adult height and body mass index (BMI). Birth weight, age at menarche and height, but not BMI, were highly significantly associated with season of birth. Individuals born in summer (June-July-August) had higher mean birth weight (P = 8 × 10^-10), later pubertal development (P = 1.1 × 10^-45) and taller adult height (P = 6.5 × 10^-9) compared to those born in all other seasons. Concordantly, those born in winter (December-January-February) showed directionally opposite differences in these outcomes. A secondary comparison of the extreme differences between months revealed higher odds ratios [95% confidence intervals (CI)] for low birth weight in February vs. September (1.23 [1.15-1.32], P = 4.4 × 10^-10), for early puberty in September vs. July (1.22 [1.16-1.28], P = 7.3 × 10^-15) and for short stature in December vs. June (1.09 [1.03-1.17], P = 0.006). The above associations were also seen with total hours of sunshine during the second trimester, but not during the first three months after birth. Additional associations were observed with educational attainment; individuals born in autumn vs. summer were more likely to continue in education post age 16 years (P = 1.1 × 10^-91) or attain a degree-level qualification (P = 4 × 10^-7). However, unlike other outcomes, an abrupt difference was seen between those born in August vs. September, which flank the start of the school year. Our findings provide support for the 'fetal programming' hypothesis, refining and extending the impact that season of birth has on childhood growth and development. Whilst other mechanisms may contribute to these associations, these findings are consistent with a possible role of in utero vitamin D exposure.

Parent-of-Origin Effects of the APOB Gene on Adiposity in Young Adults

Loci identified in genome-wide association studies (GWAS) of cardio-metabolic traits account for a small proportion of the traits' heritability. To date, most association studies have not considered parent-of-origin effects (POEs). Here we report investigation of POEs on adiposity and glycemic traits in young adults. The Jerusalem Perinatal Family Follow-Up Study (JPS), comprising 1250 young adults and their mothers was used for discovery. Focusing on 18 genes identified by previous GWAS as associated with cardio-metabolic traits, we used linear regression to examine the associations of maternally- and paternally-derived offspring minor alleles with body mass index (BMI), waist circumference (WC), fasting glucose and insulin. We replicated and meta-analyzed JPS findings in individuals of European ancestry aged ≤50 belonging to pedigrees from the Framingham Heart Study, Family Heart Study and Erasmus Rucphen Family study (total N≅4800). We considered p<2.7x10^-4 statistically significant to account for multiple testing. We identified a common coding variant in the 4^th exon of APOB (rs1367117) with a significant maternally-derived effect on BMI (β = 0.8; 95%CI:0.4,1.1; p = 3.1x10^-5) and WC (β = 2.7; 95%CI:1.7,3.7; p = 2.1x10^-7). The corresponding paternally-derived effects were non-significant (p>0.6). Suggestive maternally-derived associations of rs1367117 were observed with fasting glucose (β = 0.9; 95%CI:0.3,1.5; p = 4.0x10^-3) and insulin (ln-transformed, β = 0.06; 95%CI:0.03,0.1; p = 7.4x10^-4). Bioinformatic annotation for rs1367117 revealed a variety of regulatory functions in this region in liver and adipose tissues and a 50% methylation pattern in liver only, consistent with allelic-specific methylation, which may indicate tissue-specific POE. Our findings demonstrate a maternal-specific association between a common APOB variant and adiposity, an association that was not previously detected in GWAS. These results provide evidence for the role of regulatory mechanisms, POEs specifically, in adiposity. In addition this study highlights the benefit of utilizing family studies for deciphering the genetic architecture of complex traits.

To date, genetic variants identified in large-scale genetic studies using recent technical and methodological advances explain only a small proportion of the genetic basis of obesity, diabetes and other cardiovascular risk factors. These studies were typically conducted in samples of unrelated individuals. Here we utilize a family-based approach to identify genetic variants associated with obesity-related traits. Specifically, we examined the separate contribution of maternally- vs. paternally-inherited common genetic variants to these traits. By examining 1250 young adults and their mothers from Jerusalem, we show that a specific genetic variant, rs1367117, located in the APOB gene on chromosome 2 is related to body mass index and waist circumference when inherited from mother and not from father. This maternal effect is not restricted to Jerusalemites, but is also seen in a large sample of individuals of European descent from independent family studies worldwide. Our findings provide support of the role of complex genetic mechanisms in obesity, and highlight the benefit of utilizing family studies for uncovering genetic pathways underlying common risk factors and diseases.

The Cannabinoid Receptor 2 Q63R Variant Modulates the Relationship between Childhood Obesity and Age at Menarche

Background

The ovary is an important site where gene variants modulate pubertal timing. The cannabinoid receptor 2 (CB2) is expressed in the ovary, plays a role in folliculogenesis and ovulation, and can be modulated by estrogens. Obesity is strictly associated with early menarche and is characterized by sex hormone and endocannabinoid derangement.

Aim

In this study, we investigated the role of the CB2 receptor in determining the age at menarche in obese girls.

Methods

We studied a cohort of 240 obese girls (age 11.9±3 years; BMI z-score 2.8±0.8). The age at menarche (if it had already occurred) was recorded at the time of the visit or via phonecall. The CNR2 rs35761398 polymorphism, which leads to the CB2 Q63R variant, was detected by the TaqMan assay.

Results

In total, 105 patients were homozygous for the R63-coding allele (RR), 113 were QR and 22 were QQ. Variance analysis revealed a significantly earlier age of menarche in subjects carrying the Q63 allele, which was also found after adjusting for BMI z-score (11±1.2 vs. 11.6±1.2 years, p = 0.0003). Logistic regression analysis demonstrated that patients homozygous for the Q allele had a 2.2-fold higher risk (odds ratio = 2.2; CI1.1–3.4; p = 0.02) of presenting with an early menarche (age at menarche <12 years).

Conclusion

We demonstrated for the first time the association between the CB2 Q63R functional variant and the age at menarche in a cohort of Italian obese girls.

Comprehensive genetic assessment of the ESR1 locus identifies a risk region for endometrial cancer

Excessive exposure to estrogen is a well-established risk factor for endometrial cancer (EC), particularly for cancers of endometrioid histology. The physiological function of estrogen is primarily mediated by estrogen receptor alpha, encoded by ESR1. Consequently, several studies have investigated whether variation at the ESR1 locus is associated with risk of EC, with conflicting results. We performed comprehensive fine-mapping analyses of 3633 genotyped and imputed single nucleotide polymorphisms (SNPs) in 6607 EC cases and 37 925 controls. There was evidence of an EC risk signal located at a potential alternative promoter of the ESR1 gene (lead SNP rs79575945, P=1.86×10(-5)), which was stronger for cancers of endometrioid subtype (P=3.76×10(-6)). Bioinformatic analysis suggests that this risk signal is in a functionally important region targeting ESR1, and eQTL analysis found that rs79575945 was associated with expression of SYNE1, a neighbouring gene. In summary, we have identified a single EC risk signal located at ESR1, at study-wide significance. Given SNPs located at this locus have been associated with risk for breast cancer, also a hormonally driven cancer, this study adds weight to the rationale for performing informed candidate fine-scale genetic studies across cancer types.

Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index

We propose a method (GREML-LDMS) to estimate heritability for human complex traits in unrelated individuals using whole-genome sequencing data. We demonstrate using simulations based on whole-genome sequencing data that ∼97% and ∼68% of variation at common and rare variants, respectively, can be captured by imputation. Using the GREML-LDMS method, we estimate from 44,126 unrelated individuals that all ∼17 million imputed variants explain 56% (standard error (s.e.) = 2.3%) of variance for height and 27% (s.e. = 2.5%) of variance for body mass index (BMI), and we find evidence that height- and BMI-associated variants have been under natural selection. Considering the imperfect tagging of imputation and potential overestimation of heritability from previous family-based studies, heritability is likely to be 60-70% for height and 30-40% for BMI. Therefore, the missing heritability is small for both traits. For further discovery of genes associated with complex traits, a study design with SNP arrays followed by imputation is more cost-effective than whole-genome sequencing at current prices.

Impact of gene polymorphisms of gonadotropins and their receptors on human reproductive success

Gonadotropins and their receptors' genes carry several single-nucleotide polymorphisms resulting in endocrine genotypes modulating reproductive parameters, diseases, and lifespan leading to important implications for reproductive success and potential relevance during human evolution. Here we illustrate common genotypes of the gonadotropins and gonadotropin receptors' genes and their clinical implications in phenotypes relevant for reproduction such as ovarian cycle length, age of menopause, testosterone levels, polycystic ovary syndrome, and cancer. We then discuss their possible role in human reproduction and adaptation to the environment. Gonadotropins and their receptors' variants are differently distributed among human populations. Some hints suggest that they may be the result of natural selection that occurred in ancient times, increasing the individual chance of successful mating, pregnancy, and effective post-natal parental cares. The gender-related differences in the regulation of the reproductive endocrine systems imply that many of these genotypes may lead to sex-dependent effects, increasing the chance of mating and reproductive success in one sex at the expenses of the other sex. Also, we suggest that sexual conflicts within the FSH and LH-choriogonadotropin receptor genes contributed to maintain genotypes linked to subfertility among humans. Because the distribution of polymorphic markers results in a defined geographical pattern due to human migrations rather than natural selection, these polymorphisms may have had only a weak impact on reproductive success. On the contrary, such genotypes could acquire relevant consequences in the modern, developed societies in which parenthood attempts often occur at a later age, during a short, suboptimal reproductive window, making clinical fertility treatments necessary.

Evidence of a Causal Association Between Insulinemia and Endometrial Cancer: A Mendelian Randomization Analysis

BACKGROUND

Insulinemia and type 2 diabetes (T2D) have been associated with endometrial cancer risk in numerous observational studies. However, the causality of these associations is uncertain. Here we use a Mendelian randomization (MR) approach to assess whether insulinemia and T2D are causally associated with endometrial cancer.

METHODS

We used single nucleotide polymorphisms (SNPs) associated with T2D (49 variants), fasting glucose (36 variants), fasting insulin (18 variants), early insulin secretion (17 variants), and body mass index (BMI) (32 variants) as instrumental variables in MR analyses. We calculated MR estimates for each risk factor with endometrial cancer using an inverse-variance weighted method with SNP-endometrial cancer associations from 1287 case patients and 8273 control participants.

RESULTS

Genetically predicted higher fasting insulin levels were associated with greater risk of endometrial cancer (odds ratio [OR] per standard deviation = 2.34, 95% confidence internal [CI] = 1.06 to 5.14, P = .03). Consistently, genetically predicted higher 30-minute postchallenge insulin levels were also associated with endometrial cancer risk (OR = 1.40, 95% CI = 1.12 to 1.76, P = .003). We observed no associations between genetic risk of type 2 diabetes (OR = 0.91, 95% CI = 0.79 to 1.04, P = .16) or higher fasting glucose (OR = 1.00, 95% CI = 0.67 to 1.50, P = .99) and endometrial cancer. In contrast, endometrial cancer risk was higher in individuals with genetically predicted higher BMI (OR = 3.86, 95% CI = 2.24 to 6.64, P = 1.2x10(-6)).

CONCLUSION

This study provides evidence to support a causal association of higher insulin levels, independently of BMI, with endometrial cancer risk.

Expert consensus document: European Consensus Statement on congenital hypogonadotropic hypogonadism--pathogenesis, diagnosis and treatment

Congenital hypogonadotropic hypogonadism (CHH) is a rare disorder caused by the deficient production, secretion or action of gonadotropin-releasing hormone (GnRH), which is the master hormone regulating the reproductive axis. CHH is clinically and genetically heterogeneous, with >25 different causal genes identified to date. Clinically, the disorder is characterized by an absence of puberty and infertility. The association of CHH with a defective sense of smell (anosmia or hyposmia), which is found in ∼50% of patients with CHH is termed Kallmann syndrome and results from incomplete embryonic migration of GnRH-synthesizing neurons. CHH can be challenging to diagnose, particularly when attempting to differentiate it from constitutional delay of puberty. A timely diagnosis and treatment to induce puberty can be beneficial for sexual, bone and metabolic health, and might help minimize some of the psychological effects of CHH. In most cases, fertility can be induced using specialized treatment regimens and several predictors of outcome have been identified. Patients typically require lifelong treatment, yet ∼10-20% of patients exhibit a spontaneous recovery of reproductive function. This Consensus Statement summarizes approaches for the diagnosis and treatment of CHH and discusses important unanswered questions in the field.

Study of Five Pubertal Transition-Related Gene Polymorphisms as Risk Factors for Premature Coronary Artery Disease in a Chinese Han Population

Background

Recently, single nucleotide polymorphisms (SNPs) (DLK-rs10144321, SIX6-rs1254337, MKRN3-rs12148769, LIN28B-rs7759938, and KCNK9-rs1469039) were found to be strongly associated with age at menarche. Recent studies also suggested that age at menarche is a heritable trait and is associated with risks for obesity, type 2 diabetes mellitus (T2DM), cardiovascular disease, and all-cause mortality. Since an association between these five SNPs and premature coronary artery disease (CAD) has never been reported, we investigated whether these SNPs are associated with premature CAD and its severity in a Chinese Han population.

Methods

We enrolled 432 consecutive patients including 198 with premature CAD (<55 years in men and <65 years in women) and 234 controls. All subjects were genotyped for the five SNPs by the PCR-ligase detection reaction method. The associations between these SNPs and premature CAD and its severity were analyzed.

Results

The following genotypes were identified: GG, AG, and AA at rs10144321 and rs12148769; TT, AT, and AA at rs1254337; CC, CT, and TT at rs1469039; and TT and CT at rs7759938. Significant differences in genotype distribution frequencies at rs1254337 were found between controls and patients with premature CAD (P<0.05). No associations were found between the five SNPs and the severity of coronary lesions (all P>0.05). Compared with controls, patients with premature CAD had a higher prevalence of T2DM and dyslipidemia, and the proportion of patients with T2DM rose significantly with an increase in the number of stenosed coronary vessels (all P<0.05). After adjustment for the clinical parameters in multivariable analysis, three factors were identified that significantly increased the risk of premature CAD: the AA genotype at rs1254337 (OR: 2.388, 95% CI: 1.190–4.792, P = 0.014), male gender (OR: 1.565, 95% CI: 1.012–2.420, P = 0.044), and T2DM (OR 2.252, 95% CI: 1.233–4.348, P = 0.015).

Conclusions

Among the five pubertal transition-related gene polymorphisms, we identified an association between rs1254337 and premature CAD in a Chinese Han population.

PLAG1 and NCAPG-LCORL in livestock

A recent progress on stature genetics has revealed simple genetic architecture in livestock animals in contrast to that in humans. PLAG1 and/or NCAPG‐LCORL, both of which are known as a locus for adult human height, have been detected for association with body weight/height in cattle and horses, and for selective sweep in dogs and pigs. The findings indicate a significant impact of these loci on mammalian growth or body size and usefulness of the natural variants for selective breeding. However, association with an unfavorable trait, such as late puberty or risk for a neuropathic disease, was also reported for the respective loci, indicating an importance to discriminate between causality and association. Here I review the recent findings on quantitative trait loci (QTL) for stature in livestock animals, mainly focusing on the PLAG1 and NCAPG‐LCORL loci. I also describe our recent efforts to identify the causative variation for the third major locus for carcass weight in Japanese Black cattle.

Rare coding variants and X-linked loci associated with age at menarche

More than 100 loci have been identified for age at menarche by genome-wide association studies; however, collectively these explain only ∼3% of the trait variance. Here we test two overlooked sources of variation in 192,974 European ancestry women: low-frequency protein-coding variants and X-chromosome variants. Five missense/nonsense variants (in ALMS1/LAMB2/TNRC6A/TACR3/PRKAG1) are associated with age at menarche (minor allele frequencies 0.08-4.6%; effect sizes 0.08-1.25 years per allele; P<5 × 10(-8)). In addition, we identify common X-chromosome loci at IGSF1 (rs762080, P=9.4 × 10(-13)) and FAAH2 (rs5914101, P=4.9 × 10(-10)). Highlighted genes implicate cellular energy homeostasis, post-transcriptional gene silencing and fatty-acid amide signalling. A frequently reported mutation in TACR3 for idiopathic hypogonatrophic hypogonadism (p.W275X) is associated with 1.25-year-later menarche (P=2.8 × 10(-11)), illustrating the utility of population studies to estimate the penetrance of reportedly pathogenic mutations. Collectively, these novel variants explain ∼0.5% variance, indicating that these overlooked sources of variation do not substantially explain the 'missing heritability' of this complex trait.

Does salt have a permissive role in the induction of puberty?

Puberty is starting earlier than ever before and there are serious physiological and sociological implications as a result of this development. Current research has focused on the potential role of high caloric, and commensurate high adiposity, contributions to early puberty. However, girls with normal BMI also appear to be initiating puberty earlier. Westernized diets, in addition to being high in fat and sugar, are also high in salt. To date, no research has investigated a link between elevated salt and the reproductive axis. We hypothesize that a high salt diet can result in an earlier onset of puberty through three mechanisms that are not mutually exclusive. (1) High salt activates neurokinin B, a hormone that is involved in both the reproductive axis and salt regulation, and this induces kisspeptin release and ultimate activation of the reproductive axis. (2) Vasopressin released in response to high salt acts on vasopressin receptors expressed on kisspeptin neurons in the anteroventral periventricular nucleus, thereby stimulating gonadotropin releasing hormone and subsequently luteinizing hormone secretion. (3) Salt induces metabolic changes that affect the reproductive axis. Specifically, salt acts indirectly to modulate adiposity, ties in with the obesity epidemic, and further compounds the pathologic effects of obesity. Our overall hypothesis offers an additional cause behind the induction of puberty and provides testable postulates to determine the mechanism of potential salt-mediated affects on puberty.

Neuroendocrine control of the onset of puberty

This chapter is based on the Geoffrey Harris Memorial Lecture presented at the 8th International Congress of Neuroendocrinology, which was held in Sydney, August 2014. It provides the development of our understanding of the neuroendocrine control of puberty since Harris proposed in his 1955 monograph (Harris, 1955) that "a major factor responsible for puberty is an increased rate of release of pituitary gonadotrophin" and posited "that a neural (hypothalamic) stimulus, via the hypophysial portal vessels, may be involved." Emphasis is placed on the neurobiological mechanisms governing puberty in highly evolved primates, although an attempt is made to reverse translate a model for the timing of puberty in man and monkey to non-primate species.

Targeted resequencing of the pericentromere of chromosome 2 linked to constitutional delay of growth and puberty

Constitutional delay of growth and puberty (CDGP) is the most common cause of pubertal delay. CDGP is defined as the proportion of the normal population who experience pubertal onset at least 2 SD later than the population mean, representing 2.3% of all adolescents. While adolescents with CDGP spontaneously enter puberty, they are at risk for short stature, decreased bone mineral density, and psychosocial problems. Genetic factors contribute heavily to the timing of puberty, but the vast majority of CDGP cases remain biologically unexplained, and there is no definitive test to distinguish CDGP from pathological absence of puberty during adolescence. Recently, we published a study identifying significant linkage between a locus at the pericentromeric region of chromosome 2 (chr 2) and CDGP in Finnish families. To investigate this region for causal variation, we sequenced chr 2 between the genomic coordinates of 79-124 Mb (genome build GRCh37) in the proband and affected parent of the 13 families contributing most to this linkage signal. One gene, DNAH6, harbored 6 protein-altering low-frequency variants (< 6% in the Finnish population) in 10 of the CDGP probands. We sequenced an additional 135 unrelated Finnish CDGP subjects and utilized the unique Sequencing Initiative Suomi (SISu) population reference exome set to show that while 5 of these variants were present in the CDGP set, they were also present in the Finnish population at similar frequencies. Additional variants in the targeted region could not be prioritized for follow-up, possibly due to gaps in sequencing coverage or lack of functional knowledge of non-genic genomic regions. Thus, despite having a well-characterized sample collection from a genetically homogeneous population with a large population-based reference sequence dataset, we were unable to pinpoint variation in the linked region predisposing delayed puberty. This study highlights the difficulties of detecting genetic variants under linkage regions for complex traits and suggests that advancements in annotation of gene function and regulatory regions of the genome will be critical for solving the genetic background of complex phenotypes like CDGP.

A new pathway in the control of the initiation of puberty: the MKRN3 gene

Pubertal timing is influenced by complex interactions among genetic, nutritional, environmental, and socioeconomic factors. The role of MKRN3, an imprinted gene located in the Prader-Willi syndrome critical region (chromosome 15q11-13), in pubertal initiation was first described in 2013 after the identification of deleterious MKRN3 mutations in five families with central precocious puberty (CPP) using whole-exome sequencing analysis. Since then, additional loss-of-function mutations of MKRN3 have been associated with the inherited premature sexual development phenotype in girls and boys from different ethnic groups. In all of these families, segregation analysis clearly demonstrated autosomal dominant inheritance with complete penetrance, but with exclusive paternal transmission, consistent with the monoallelic expression of MKRN3 (a maternally imprinted gene). Interestingly, the hypothalamic Mkrn3 mRNA expression pattern in mice correlated with a putative inhibitory input on puberty initiation. Indeed, the initiation of puberty depends on a decrease in factors that inhibit the release of GnRH combined with an increase in stimulatory factors. These recent human and animal findings suggest that MKRN3 plays an inhibitory role in the reproductive axis to represent a new pathway in pubertal regulation.

Copy number variations and cognitive phenotypes in unselected populations

IMPORTANCE

The association of copy number variations (CNVs), differing numbers of copies of genetic sequence at locations in the genome, with phenotypes such as intellectual disability has been almost exclusively evaluated using clinically ascertained cohorts. The contribution of these genetic variants to cognitive phenotypes in the general population remains unclear.

OBJECTIVE

To investigate the clinical features conferred by CNVs associated with known syndromes in adult carriers without clinical preselection and to assess the genome-wide consequences of rare CNVs (frequency ≤0.05%; size ≥250 kilobase pairs [kb]) on carriers' educational attainment and intellectual disability prevalence in the general population.

DESIGN, SETTING, AND PARTICIPANTS

The population biobank of Estonia contains 52,000 participants enrolled from 2002 through 2010. General practitioners examined participants and filled out a questionnaire of health- and lifestyle-related questions, as well as reported diagnoses. Copy number variant analysis was conducted on a random sample of 7877 individuals and genotype-phenotype associations with education and disease traits were evaluated. Our results were replicated on a high-functioning group of 993 Estonians and 3 geographically distinct populations in the United Kingdom, the United States, and Italy.

MAIN OUTCOMES AND MEASURES

Phenotypes of genomic disorders in the general population, prevalence of autosomal CNVs, and association of these variants with educational attainment (from less than primary school through scientific degree) and prevalence of intellectual disability.

RESULTS

Of the 7877 in the Estonian cohort, we identified 56 carriers of CNVs associated with known syndromes. Their phenotypes, including cognitive and psychiatric problems, epilepsy, neuropathies, obesity, and congenital malformations are similar to those described for carriers of identical rearrangements ascertained in clinical cohorts. A genome-wide evaluation of rare autosomal CNVs (frequency, ≤0.05%; ≥250 kb) identified 831 carriers (10.5%) of the screened general population. Eleven of 216 (5.1%) carriers of a deletion of at least 250 kb (odds ratio [OR], 3.16; 95% CI, 1.51-5.98; P = 1.5e-03) and 6 of 102 (5.9%) carriers of a duplication of at least 1 Mb (OR, 3.67; 95% CI, 1.29-8.54; P = .008) had an intellectual disability compared with 114 of 6819 (1.7%) in the Estonian cohort. The mean education attainment was 3.81 (P = 1.06e-04) among 248 (≥250 kb) deletion carriers and 3.69 (P = 5.024e-05) among 115 duplication carriers (≥1 Mb). Of the deletion carriers, 33.5% did not graduate from high school (OR, 1.48; 95% CI, 1.12-1.95; P = .005) and 39.1% of duplication carriers did not graduate high school (OR, 1.89; 95% CI, 1.27-2.8; P = 1.6e-03). Evidence for an association between rare CNVs and lower educational attainment was supported by analyses of cohorts of adults from Italy and the United States and adolescents from the United Kingdom.

CONCLUSIONS AND RELEVANCE

Known pathogenic CNVs in unselected, but assumed to be healthy, adult populations may be associated with unrecognized clinical sequelae. Additionally, individually rare but collectively common intermediate-size CNVs may be negatively associated with educational attainment. Replication of these findings in additional population groups is warranted given the potential implications of this observation for genomics research, clinical care, and public health.

Genomic imprinting: A missing piece of the Multiple Sclerosis puzzle?

Evidence for parent-of-origin effects in complex diseases such as Multiple Sclerosis (MS) strongly suggests a role for epigenetic mechanisms in their pathogenesis. In this review, we describe the importance of accounting for parent-of-origin when identifying new risk variants for complex diseases and discuss how genomic imprinting, one of the best-characterized epigenetic mechanisms causing parent-of-origin effects, may impact etiology of complex diseases. While the role of imprinted genes in growth and development is well established, the contribution and molecular mechanisms underlying the impact of genomic imprinting in immune functions and inflammatory diseases are still largely unknown. Here we discuss emerging roles of imprinted genes in the regulation of inflammatory responses with a particular focus on the Dlk1 cluster that has been implicated in etiology of experimental MS-like disease and Type 1 Diabetes. Moreover, we speculate on the potential wider impact of imprinting via the action of imprinted microRNAs, which are abundantly present in the Dlk1 locus and predicted to fine-tune important immune functions. Finally, we reflect on how unrelated imprinted genes or imprinted genes together with non-imprinted genes can interact in so-called imprinted gene networks (IGN) and suggest that IGNs could partly explain observed parent-of-origin effects in complex diseases. Unveiling the mechanisms of parent-of-origin effects is therefore likely to teach us not only about the etiology of complex diseases but also about the unknown roles of this fascinating phenomenon underlying uneven genetic contribution from our parents. This article is part of a Directed Issue entitled: Epigenetics dynamics in development and disease.

Widespread seasonal gene expression reveals annual differences in human immunity and physiology

Seasonal variations are rarely considered a contributing component to human tissue function or health, although many diseases and physiological process display annual periodicities. Here we find more than 4,000 protein-coding mRNAs in white blood cells and adipose tissue to have seasonal expression profiles, with inverted patterns observed between Europe and Oceania. We also find the cellular composition of blood to vary by season, and these changes, which differ between the United Kingdom and The Gambia, could explain the gene expression periodicity. With regards to tissue function, the immune system has a profound pro-inflammatory transcriptomic profile during European winter, with increased levels of soluble IL-6 receptor and C-reactive protein, risk biomarkers for cardiovascular, psychiatric and autoimmune diseases that have peak incidences in winter. Circannual rhythms thus require further exploration as contributors to various aspects of human physiology and disease.

Molecular basis of imprinting disorders affecting chromosome 14: lessons from murine models

Uniparental inheritance of chromosome 14q32 causes developmental failure during gestation and early postnatal development due to mis-expression of a cluster of imprinted genes under common epigenetic control. Two syndromes associated with chromosome 14q32 abnormalities have been described, Kagami–Ogata and Temple syndromes. Both of these syndromes are characterised by specific impairments of intrauterine development, placentation and early postnatal survival. Such abnormalities arise because the processes of intrauterine growth and postnatal adaptation are critically modulated by the dosage of imprinted genes in the chromosome 14q32 cluster. Much of our understanding of how the imprinted genes in this cluster are regulated, as well as their individual functions in the molecular pathways controlling growth and postnatal adaptation, has come from murine models. Mouse chromosome 12qF1 contains an imprinted region syntenic to human chromosome 14q32, collectively referred to as the Dlk1–Dio3 cluster. In this review, we will summarise the wealth of information derived from animal models of chromosome 12 imprinted gene mis-regulation, and explore the relationship between the functions of individual genes and the phenotypic result of their mis-expression. As there is often a considerable overlap between the functions of genes in the Dlk1–Dio3 cluster, we propose that the expression dosage of these genes is controlled by common regulatory mechanisms to co-ordinate the timing of growth and postnatal adaptation. While the diseases associated with mis-regulated chromosome 14 imprinting are rare, studies carried out in mice on the functions of the affected genes as well as their normal regulatory mechanisms have revealed new mechanistic pathways for the control of growth and survival in early life.

A shared genetic basis for self-limited delayed puberty and idiopathic hypogonadotropic hypogonadism

CONTEXT

Delayed puberty (DP) is a common issue and, in the absence of an underlying condition, is typically self limited. Alhough DP seems to be heritable, no specific genetic cause for DP has yet been reported. In contrast, many genetic causes have been found for idiopathic hypogonadotropic hypogonadism (IHH), a rare disorder characterized by absent or stalled pubertal development.

OBJECTIVE

The objective of this retrospective study, conducted at academic medical centers, was to determine whether variants in IHH genes contribute to the pathogenesis of DP.

SUBJECTS AND OUTCOME MEASURES

Potentially pathogenic variants in IHH genes were identified in two cohorts: 1) DP family members of an IHH proband previously found to have a variant in an IHH gene, with unaffected family members serving as controls, and 2) DP individuals with no family history of IHH, with ethnically matched control subjects drawn from the Exome Aggregation Consortium.

RESULTS

In pedigrees with an IHH proband, the proband's variant was shared by 53% (10/19) of DP family members vs 12% (4/33) of unaffected family members (P = .003). In DP subjects with no family history of IHH, 14% (8/56) had potentially pathogenic variants in IHH genes vs 5.6% (1 907/33 855) of controls (P = .01). Potentially pathogenic variants were found in multiple DP subjects for the genes IL17RD and TAC3.

CONCLUSIONS

These findings suggest that variants in IHH genes can contribute to the pathogenesis of self-limited DP. Thus, at least in some cases, self-limited DP shares an underlying pathophysiology with IHH.

Genetic Regulation of Puberty Timing in Humans

Understanding the regulation of puberty timing has relevance to developmental and human biology and to the pathogenesis of various diseases. Recent large-scale genome-wide association studies on puberty timing and adult height, body mass index (BMI) and central body shape provide evidence for shared biological mechanisms that regulate these traits. There is a substantial genetic overlap between age at menarche in women and BMI, with almost invariable directional consistency with the epidemiological associations between earlier menarche and higher BMI. By contrast, the genetic loci identified for age at menarche are largely distinct from those identified for central body shape, while alleles that confer earlier menarche can be associated with taller or shorter adult height. The findings of population-based studies on age at menarche show increasing relevance for other studies of rare monogenic disorders and enrich our understanding of the mechanisms that regulate the timing of puberty and reproductive function.

Cohort Profile Update: the TRacking Adolescents' Individual Lives Survey (TRAILS)

TRAILS consists of a population cohort (N = 2230) and a clinical cohort (N = 543), both of which were followed from about age 11 years onwards. To date, the population cohort has been assessed five times over a period of 11 years, with retention rates ranging between 80% and 96%. The clinical cohort has been assessed four times over a period of 8 years, with retention rates ranging between 77% and 85%. Since the IJE published a cohort profile on the TRAILS in 2008, the participants have matured from adolescents into young adults. The focus shifted from parents and school to entry into the labour market and family formation, including offspring. Furthermore, psychiatric diagnostic interviews were administered, the database was linked to a Psychiatric Case Registry, and the availability of genome-wide SNP variations opened the door to genome-wide association studies regarding a wide range of (endo)phenotypes. With some delay, TRAILS data are available to researchers outside the TRAILS consortium without costs; access can be obtained by submitting a publication proposal (see www.trails.nl).

Pubertal onset in girls is strongly influenced by genetic variation affecting FSH action

Age at pubertal onset varies substantially in healthy girls. Although genetic factors are responsible for more than half of the phenotypic variation, only a small part has been attributed to specific genetic polymorphisms identified so far. Follicle-stimulating hormone (FSH) stimulates ovarian follicle maturation and estradiol synthesis which is responsible for breast development. We assessed the effect of three polymorphisms influencing FSH action on age at breast deveopment in a population-based cohort of 964 healthy girls. Girls homozygous for FSHR -29AA (reduced FSH receptor expression) entered puberty 7.4 (2.5–12.4) months later than carriers of the common variants FSHR -29GG+GA, p = 0.003. To our knowledge, this is the strongest genetic effect on age at pubertal onset in girls published to date.