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

Genetics of thyroid function

Recent studies show that subtle variations in thyroid function, including subclinical thyroid dysfunction, and even variation in thyroid function within the normal range, are associated with morbidity and mortality. It is estimated that 40-65% of the inter-individual variation in serum TSH and FT4 levels is determined by genetic factors. To identify these factors, various linkage and candidate gene studies have been performed in the past, which have identified only a few genes. In the last decade, genome-wide association studies identified many new genes, while recent whole-genome sequencing efforts have also been proven to be effective. In the current review, we provide a systematic overview of these studies, including strengths and limitations. We discuss new techniques which will further clarify the genetic basis of thyroid function in the near future, as well as the potential use of these genetic markers in personalizing the management of thyroid disease patients.

Diversity in non-repetitive human sequences not found in the reference genome

Genomes usually contain some non-repetitive sequences that are missing from the reference genome and occur only in a population subset. Such non-repetitive, non-reference (NRNR) sequences have remained largely unexplored in terms of their characterization and downstream analyses. Here we describe 3,791 breakpoint-resolved NRNR sequence variants called using PopIns from whole-genome sequence data of 15,219 Icelanders. We found that over 95% of the 244 NRNR sequences that are 200 bp or longer are present in chimpanzees, indicating that they are ancestral. Furthermore, 149 variant loci are in linkage disequilibrium (r² > 0.8) with a genome-wide association study (GWAS) catalog marker, suggesting disease relevance. Additionally, we report an association (P = 3.8 × 10^-8, odds ratio (OR) = 0.92) with myocardial infarction (23,360 cases, 300,771 controls) for a 766-bp NRNR sequence variant. Our results underline the importance of including variation of all complexity levels when searching for variants that associate with disease.

A molecular census of arcuate hypothalamus and median eminence cell types

The hypothalamic arcuate-median eminence complex (Arc-ME) controls energy balance, fertility, and growth through molecularly distinct cell types, many of which remain unknown. To catalog cell types in an unbiased way, we profiled gene expression in 20,921 individual cells in and around the adult mouse Arc-ME using Drop-seq. We identify 50 transcriptionally distinct Arc-ME cell populations, including a rare tanycyte population at the Arc-ME diffusion barrier, a novel leptin-sensing neuronal population, multiple AgRP and POMC subtypes, and an orexigenic somatostatin neuronal population. We extended Drop-seq to detect dynamic expression changes across relevant physiological perturbations, revealing cell type-specific responses to energy status, including distinctly responsive subtypes of AgRP and POMC neurons. Finally, integrating our data with human GWAS data implicates two previously unknown neuronal subtypes in the genetic control of obesity. This resource will accelerate biological discovery by providing insights into molecular and cell type diversity from which function can be inferred.

Nakalanga Syndrome: Clinical Characteristics, Potential Causes, and Its Relationship with Recently Described Nodding Syndrome

Nakalanga syndrome is a condition that was described in Uganda and various other African countries decades ago. Its features include growth retardation, physical deformities, endocrine dysfunction, mental impairment, and epilepsy, amongst others. Its cause remains obscure. Nodding syndrome is a neurological disorder with some features in common with Nakalanga syndrome, which has been described mainly in Uganda, South Sudan, and Tanzania. It has been considered an encephalopathy affecting children who, besides head nodding attacks, can also present with stunted growth, delayed puberty, and mental impairment, amongst other symptoms. Despite active research over the last years on the pathogenesis of Nodding syndrome, to date, no convincing single cause of Nodding syndrome has been reported. In this review, by means of a thorough literature search, we compare features of both disorders. We conclude that Nakalanga and Nodding syndromes are closely related and may represent the same condition. Our findings may provide new directions in research on the cause underlying this neurological disorder.

Dissecting the genetics of complex traits using summary association statistics

During the past decade, genome-wide association studies (GWAS) have been used to successfully identify tens of thousands of genetic variants associated with complex traits and diseases. These studies have produced extensive repositories of genetic variation and trait measurements across large numbers of individuals, providing tremendous opportunities for further analyses. However, privacy concerns and other logistical considerations often limit access to individual-level genetic data, motivating the development of methods that analyse summary association statistics. Here, we review recent progress on statistical methods that leverage summary association data to gain insights into the genetic basis of complex traits and diseases.

LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis

MOTIVATION

LD score regression is a reliable and efficient method of using genome-wide association study (GWAS) summary-level results data to estimate the SNP heritability of complex traits and diseases, partition this heritability into functional categories, and estimate the genetic correlation between different phenotypes. Because the method relies on summary level results data, LD score regression is computationally tractable even for very large sample sizes. However, publicly available GWAS summary-level data are typically stored in different databases and have different formats, making it difficult to apply LD score regression to estimate genetic correlations across many different traits simultaneously.

RESULTS

In this manuscript, we describe LD Hub - a centralized database of summary-level GWAS results for 173 diseases/traits from different publicly available resources/consortia and a web interface that automates the LD score regression analysis pipeline. To demonstrate functionality and validate our software, we replicated previously reported LD score regression analyses of 49 traits/diseases using LD Hub; and estimated SNP heritability and the genetic correlation across the different phenotypes. We also present new results obtained by uploading a recent atopic dermatitis GWAS meta-analysis to examine the genetic correlation between the condition and other potentially related traits. In response to the growing availability of publicly accessible GWAS summary-level results data, our database and the accompanying web interface will ensure maximal uptake of the LD score regression methodology, provide a useful database for the public dissemination of GWAS results, and provide a method for easily screening hundreds of traits for overlapping genetic aetiologies.

AVAILABILITY AND IMPLEMENTATION

The web interface and instructions for using LD Hub are available at http://ldsc.broadinstitute.org/ CONTACT: jie.zheng@bristol.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Association of timing of menarche with depressive symptoms and depression in adolescence: Mendelian randomisation study

BACKGROUND

Observational studies report associations between early menarche and higher levels of depressive symptoms and depression. However, no studies have investigated whether this association is causal.

AIMS

To determine whether earlier menarche is a causal risk factor for depressive symptoms and depression in adolescence.

METHOD

The associations between a genetic score for age at menarche and depressive symptoms at 14, 17 and 19 years, and depression at 18 years, were examined using Mendelian randomisation analysis techniques.

RESULTS

Using a genetic risk score to indicate earlier timing of menarche, we found that early menarche is associated with higher levels of depressive symptoms at 14 years (odds ratio per risk allele 1.02, 95% CI 1.005-1.04, n = 2404). We did not find an association between the early menarche risk score and depressive symptoms or depression after age 14.

CONCLUSIONS

Our results provide evidence for a causal effect of age at menarche on depressive symptoms at age 14.

Saving the spandrels? Adaptive genomic variation in conservation and fisheries management

As highlighted by many of the papers in this issue, research on the genomic basis of adaptive phenotypic variation in natural populations has made spectacular progress in the past few years, largely due to the advances in sequencing technology and analysis. Without question, the resulting genomic data will improve the understanding of regions of the genome under selection and extend knowledge of the genetic basis of adaptive evolution. What is far less clear, but has been the focus of active discussion, is how such information can or should transfer into conservation practice to complement more typical conservation applications of genetic data. Before such applications can be realized, the evolutionary importance of specific targets of selection relative to the genome-wide diversity of the species as a whole must be evaluated. The key issues for the incorporation of adaptive genomic variation in conservation and management are discussed here, using published examples of adaptive genomic variation associated with specific phenotypes in salmonids and other taxa to highlight practical considerations for incorporating such information into conservation programmes. Scenarios are described in which adaptive genomic data could be used in conservation or restoration, constraints on its utility and the importance of validating inferences drawn from new genomic data before applying them in conservation practice. Finally, it is argued that an excessive focus on preserving the adaptive variation that can be measured, while ignoring the vast unknown majority that cannot, is a modern twist on the adaptationist programme that Gould and Lewontin critiqued almost 40 years ago.

Age at menarche and age at natural menopause in East Asian women: a genome-wide association study

Age at menarche (AM) and age at natural menopause (ANM) are complex traits with a high heritability. Abnormal timing of menarche or menopause is associated with a reduced span of fertility and risk for several age-related diseases including breast, endometrial and ovarian cancer, cardiovascular disease, and osteoporosis. To identify novel genetic loci for AM or ANM in East Asian women and to replicate previously identified loci primarily in women of European ancestry by genome-wide association studies (GWASs), we conducted a two-stage GWAS. Stage I aimed to discover promising novel AM and ANM loci using GWAS data of 8073 women from Shanghai, China. The Stage II replication study used the data from another Chinese GWAS (n = 1230 for AM and n = 1458 for ANM), a Korean GWAS (n = 4215 for AM and n = 1739 for ANM), and de novo genotyping of 2877 additional Chinese women. Previous GWAS-identified loci for AM and ANM were also evaluated. We identified two suggestive menarcheal age loci tagged by rs79195475 at 10q21.3 (beta = -0.118 years, P = 3.4 × 10^-6) and rs1023935 at 4p15.1 (beta = -0.145 years, P = 4.9 × 10^-6) and one menopausal age locus tagged by rs3818134 at 22q12.2 (beta = -0.276 years, P = 8.8 × 10^-6). These suggestive loci warrant a further validation in independent populations. Although limited by low statistical power, we replicated 19 of the 98 menarche loci and 5 of the 20 menopause loci previously identified in women of European ancestry in East Asian women, suggesting a shared genetic architecture for these two traits across populations.

Genome-wide analysis identifies 12 loci influencing human reproductive behavior

The genetic architecture of human reproductive behavior-age at first birth (AFB) and number of children ever born (NEB)-has a strong relationship with fitness, human development, infertility and risk of neuropsychiatric disorders. However, very few genetic loci have been identified, and the underlying mechanisms of AFB and NEB are poorly understood. We report a large genome-wide association study of both sexes including 251,151 individuals for AFB and 343,072 individuals for NEB. We identified 12 independent loci that are significantly associated with AFB and/or NEB in a SNP-based genome-wide association study and 4 additional loci associated in a gene-based effort. These loci harbor genes that are likely to have a role, either directly or by affecting non-local gene expression, in human reproduction and infertility, thereby increasing understanding of these complex traits.

RASopathies Are Associated With Delayed Puberty; Are They Associated With Precocious Puberty Too?

RASopathies, such as Noonan, Costello, and cardio-facio-cutaneous syndromes, are developmental disorders caused by mutations in rat sarcoma-mitogen-activated protein kinase pathway genes. Mutations that cause Noonan syndrome have been associated with delayed puberty. Here we report 4 patients with either Costello or cardio-facio-cutaneous syndrome who developed precocious puberty, suggesting complex regulation of the hypothalamic-pituitary-gonadal axis and the timing of puberty by the rat sarcoma-mitogen-activated protein kinase pathway. Additional study of the timing of puberty among patients with RASopathies is warranted to ascertain the incidence of delayed and precocious puberty in these conditions and to examine genotype-phenotype correlations, which may provide insight into pathways that regulate the timing of puberty.

Epigenetic and genetic components of height regulation

Adult height is a highly heritable trait. Here we identified 31.6 million sequence variants by whole-genome sequencing of 8,453 Icelanders and tested them for association with adult height by imputing them into 88,835 Icelanders. Here we discovered 13 novel height associations by testing four different models including parent-of-origin (|β|=0.4-10.6 cm). The minor alleles of three parent-of-origin signals associate with less height only when inherited from the father and are located within imprinted regions (IGF2-H19 and DLK1-MEG3). We also examined the association of these sequence variants in a set of 12,645 Icelanders with birth length measurements. Two of the novel variants, (IGF2-H19 and TET1), show significant association with both adult height and birth length, indicating a role in early growth regulation. Among the parent-of-origin signals, we observed opposing parental effects raising questions about underlying mechanisms. These findings demonstrate that common variations affect human growth by parental imprinting.

An Ancient Fecundability-Associated Polymorphism Switches a Repressor into an Enhancer of Endometrial TAP2 Expression

Variation in female reproductive traits, such as fertility, fecundity, and fecundability, is heritable in humans, but identifying and functionally characterizing genetic variants associated with these traits has been challenging. Here, we explore the functional significance and evolutionary history of a T/C polymorphism of SNP rs2071473, which we have previously shown is an eQTL for TAP2 and significantly associated with fecundability (time to pregnancy). We replicated the association between the rs2071473 genotype and TAP2 expression by using GTEx data and demonstrated that TAP2 is expressed by decidual stromal cells at the maternal-fetal interface. Next, we showed that rs2071473 is located within a progesterone-responsive cis-regulatory element that functions as a repressor with the T allele and an enhancer with the C allele. Remarkably, we found that this polymorphism arose before the divergence of modern and archaic humans, segregates at intermediate to high frequencies across human populations, and has genetic signatures of long-term balancing selection. This variant has also previously been identified in genome-wide association studies of immune-related disease, suggesting that both alleles are maintained as a result of antagonistic pleiotropy.

microRNAs and the adolescent brain: Filling the knowledge gap

Over two decades ago the discovery of microRNAs (miRNA) broadened our understanding of the diverse molecular pathways mediating post-transcriptional control over gene expression. These small non-coding RNAs dynamically fluctuate, temporally and spatially, throughout the lifespan of all organisms. The fundamental role that miRNAs have in shaping embryonic neurodevelopment provides strong evidence that adolescent brain remodeling could be rooted in the changing miRNA landscape of the cell. Few studies have directly measured miRNA gene expression changes in the brain across pubertal development, and even less is known about the functional impact of those miRNAs on the maturational processes that occur in the developing adolescent brain. This review summarizes miRNA biogenesis and function in the brain in the context of normal (i.e. not diseased) physiology. These landmark studies can guide predictions about the role of miRNAs in facilitating maturation of the adolescent brain. However, there are clear indicators that adolescence/puberty is a unique life stage, suggesting miRNA function during adolescence is distinct from those in any other previously described system.

GWIS: Genome-Wide Inferred Statistics for Functions of Multiple Phenotypes

Here we present a method of genome-wide inferred study (GWIS) that provides an approximation of genome-wide association study (GWAS) summary statistics for a variable that is a function of phenotypes for which GWAS summary statistics, phenotypic means, and covariances are available. A GWIS can be performed regardless of sample overlap between the GWAS of the phenotypes on which the function depends. Because a GWIS provides association estimates and their standard errors for each SNP, a GWIS can form the basis for polygenic risk scoring, LD score regression, Mendelian randomization studies, biological annotation, and other analyses. GWISs can also be used to boost power of a GWAS meta-analysis where cohorts have not measured all constituent phenotypes in the function. We demonstrate the accuracy of a BMI GWIS by performing power simulations and type I error simulations under varying circumstances, and we apply a GWIS by reconstructing a body mass index (BMI) GWAS based on a weight GWAS and a height GWAS. Furthermore, we apply a GWIS to further our understanding of the underlying genetic structure of bipolar disorder and schizophrenia and their relation to educational attainment. Our analyses suggest that the previously reported genetic correlation between schizophrenia and educational attainment is probably induced by the observed genetic correlation between schizophrenia and bipolar disorder and the previously reported genetic correlation between bipolar disorder and educational attainment.

Dramatic acceleration of reproductive aging, contraction of biochemical fecundity and healthspan-lifespan implications of opioid-induced endocrinopathy-FSH/LH ratio and other interrelationships

Whilst disturbances of female reproductive hormones and function are commonplace in opioid dependence, their pathophysiological interrelationships are not well understood. Hormonal levels in females were compared in 77 opioid dependent patients (ODP) and 148 medical controls (MC) including 205 and 364 repeat studies. Significant changes in FSH, LH, oestradiol, testosterone and SBG were noted including power functions with age. The FSH/LH was lower in ODP (P=0.0150) and the ratio inversion point occurred at 28.06±9.36v. 46.26±4.76years, implying a 58% reduction in fertility duration. FSH has been shown to induce ovarian failure and GnRH (controlling LH and FSH) has been shown to regulate longevity systemically. This implies that, far from being benign, these findings explicate the adverse experience of female compared to male ODP, exacerbate opioid-dependent aging amongst females, and informs the care of opioid dependent women, particularly relating to the choice, dose and duration of agonist or antagonist therapy.

Gene set analysis for interpreting genetic studies

Interpretation of genome-wide association study (GWAS) results is lacking behind the discovery of new genetic associations. Consequently, there is an urgent need for data-driven methods for interpreting genetic association studies. Gene set analysis (GSA) can identify aetiologic pathways and functional annotations and may hence point towards novel biological insights. However, despite the growing availability of GSA tools, the sizeable amount of variants identified for a vast number of complex traits, and many irrefutably trait-associated gene sets, the gap between discovery and interpretation remains. More efficient interpretation requires more complete and consistent gene set representations of biological pathways, phenotypes and functional annotations. In this review, I examine different types of gene sets, discuss how inconsistencies in gene set definitions impact GSA, describe how GSA has helped to elucidate biology and outline potential future directions.

Control of the onset of puberty

PURPOSE OF REVIEW

The mechanism of puberty initiation remains an enigma, despite extensive research in the field. Pulsatile pituitary gonadotropin secretion under the guidance of hypothalamic gonadotropin-releasing hormone (GnRH) constitutes a sine qua non for pubertal onset. In turn, the secretion of GnRH in the human hypothalamus is regulated by kisspeptin and its receptor as well as by permissive or opposing signals mediated by neurokinin B and dynorphin acting on their respective receptors. These three supra-GnRH regulators compose the Kisspeptin, Neurokinin B and Dynorhin neurons (KNDy) system, a key player in pubertal onset and progression.

RECENT FINDINGS

The recent discovery that makorin ring finger protein 3 is also involved in puberty initiation provided further insights into the regulation of the KNDy pathway. In fact, the inhibitory (γ-amino butyric acid, neuropeptide Y, and RFamide-related peptide-3) and stimulatory signals (glutamate) acting upstream of KNDy called into question the role of makorin ring finger protein 3 as the gatekeeper of puberty. Meanwhile, the findings that 'neuroestradiol' produced locally and endocrine disruptors from the environment may influence GnRH secretion is intriguing. Finally, epigenetic mechanisms have been implicated in pubertal onset through recently discovered mechanisms.

SUMMARY

The exact molecular machinery underlying puberty initiation in humans is under intensive investigation. In this review, we summarize research evidence in the field, while emphasizing the areas of uncertainty and underlining the impact of current information on the evolving theory regarding this fascinating phenomenon.

Genetic evidence for natural selection in humans in the contemporary United States

Recent findings from molecular genetics now make it possible to test directly for natural selection by analyzing whether genetic variants associated with various phenotypes have been under selection. I leverage these findings to construct polygenic scores that use individuals' genotypes to predict their body mass index, educational attainment (EA), glucose concentration, height, schizophrenia, total cholesterol, and (in females) age at menarche. I then examine associations between these scores and fitness to test whether natural selection has been occurring. My study sample includes individuals of European ancestry born between 1931 and 1953 who participated in the Health and Retirement Study, a representative study of the US population. My results imply that natural selection has been slowly favoring lower EA in both females and males, and are suggestive that natural selection may have favored a higher age at menarche in females. For EA, my estimates imply a rate of selection of about -1.5 mo of education per generation (which pales in comparison with the increases in EA observed in contemporary times). Although they cannot be projected over more than one generation, my results provide additional evidence that humans are still evolving-albeit slowly, especially compared with the rapid changes that have occurred over the past few generations due to cultural and environmental factors.

Contrasting the Genetic Architecture of 30 Complex Traits from Summary Association Data

Variance-component methods that estimate the aggregate contribution of large sets of variants to the heritability of complex traits have yielded important insights into the genetic architecture of common diseases. Here, we introduce methods that estimate the total trait variance explained by the typed variants at a single locus in the genome (local SNP heritability) from genome-wide association study (GWAS) summary data while accounting for linkage disequilibrium among variants. We applied our estimator to ultra-large-scale GWAS summary data of 30 common traits and diseases to gain insights into their local genetic architecture. First, we found that common SNPs have a high contribution to the heritability of all studied traits. Second, we identified traits for which the majority of the SNP heritability can be confined to a small percentage of the genome. Third, we identified GWAS risk loci where the entire locus explains significantly more variance in the trait than the GWAS reported variants. Finally, we identified loci that explain a significant amount of heritability across multiple traits.

Assessing the genetic architecture of epithelial ovarian cancer histological subtypes

Epithelial ovarian cancer (EOC) is one of the deadliest common cancers. The five most common types of disease are high-grade and low-grade serous, endometrioid, mucinous and clear cell carcinoma. Each of these subtypes present distinct molecular pathogeneses and sensitivities to treatments. Recent studies show that certain genetic variants confer susceptibility to all subtypes while other variants are subtype-specific. Here, we perform an extensive analysis of the genetic architecture of EOC subtypes. To this end, we used data of 10,014 invasive EOC patients and 21,233 controls from the Ovarian Cancer Association Consortium genotyped in the iCOGS array (211,155 SNPs). We estimate the array heritability (attributable to variants tagged on arrays) of each subtype and their genetic correlations. We also look for genetic overlaps with factors such as obesity, smoking behaviors, diabetes, age at menarche and height. We estimated the array heritabilities of high-grade serous disease ([Formula: see text] = 8.8 ± 1.1 %), endometrioid ([Formula: see text] = 3.2 ± 1.6 %), clear cell ([Formula: see text] = 6.7 ± 3.3 %) and all EOC ([Formula: see text] = 5.6 ± 0.6 %). Known associated loci contributed approximately 40 % of the total array heritability for each subtype. The contribution of each chromosome to the total heritability was not proportional to chromosome size. Through bivariate and cross-trait LD score regression, we found evidence of shared genetic backgrounds between the three high-grade subtypes: serous, endometrioid and undifferentiated. Finally, we found significant genetic correlations of all EOC with diabetes and obesity using a polygenic prediction approach.

Detection and interpretation of shared genetic influences on 42 human traits

We performed a scan for genetic variants associated with multiple phenotypes by comparing large genome-wide association studies (GWAS) of 42 traits or diseases. We identified 341 loci (at a false discovery rate of 10%) associated with multiple traits. Several loci are associated with multiple phenotypes; for example, a nonsynonymous variant in the zinc transporter SLC39A8 influences seven of the traits, including risk of schizophrenia (rs13107325: log-transformed odds ratio (log OR) = 0.15, P = 2 × 10(-12)) and Parkinson disease (log OR = -0.15, P = 1.6 × 10(-7)), among others. Second, we used these loci to identify traits that have multiple genetic causes in common. For example, variants associated with increased risk of schizophrenia also tended to be associated with increased risk of inflammatory bowel disease. Finally, we developed a method to identify pairs of traits that show evidence of a causal relationship. For example, we show evidence that increased body mass index causally increases triglyceride levels.

Ovarian Physiology and GWAS: Biobanks, Biology, and Beyond

Ovarian function is central to female fertility, and several genome-wide association studies (GWAS) have been carried out to elucidate the genetic background of traits and disorders that reflect and affect ovarian physiology. While GWAS have been successful in reporting numerous genetic associations and highlighting involved pathways relevant to reproductive aging, for ovarian disorders, such as premature ovarian insufficiency and polycystic ovary syndrome, research has lagged behind due to insufficient study sample size. Novel approaches to study design and analysis methods that help to fit GWAS findings into biological context will improve our knowledge about genetics governing ovarian function in fertility and disease, and provide input for clinical tools and better patient management.

Pubertal development in healthy children is mirrored by DNA methylation patterns in peripheral blood

Puberty marks numerous physiological processes which are initiated by central activation of the hypothalamic–pituitary–gonadal axis, followed by development of secondary sexual characteristics. To a large extent, pubertal timing is heritable, but current knowledge of genetic polymorphisms only explains few months in the large inter-individual variation in the timing of puberty. We have analysed longitudinal genome-wide changes in DNA methylation in peripheral blood samples (n = 102) obtained from 51 healthy children before and after pubertal onset. We show that changes in single methylation sites are tightly associated with physiological pubertal transition and altered reproductive hormone levels. These methylation sites cluster in and around genes enriched for biological functions related to pubertal development. Importantly, we identified that methylation of the genomic region containing the promoter of TRIP6 was co-ordinately regulated as a function of pubertal development. In accordance, immunohistochemistry identified TRIP6 in adult, but not pre-pubertal, testicular Leydig cells and circulating TRIP6 levels doubled during puberty. Using elastic net prediction models, methylation patterns predicted pubertal development more accurately than chronological age. We demonstrate for the first time that pubertal attainment of secondary sexual characteristics is mirrored by changes in DNA methylation patterns in peripheral blood. Thus, modulations of the epigenome seem involved in regulation of the individual pubertal timing.

Population variation of LIN28B in the Mediterranean: Novel markers for microgeographic discrimination

OBJECTIVES

The aim of this study is to determine whether the LIN28B gene is differentially distributed in the Mediterranean region through the analysis of the allele distribution of three single nucleotide polymorphisms (SNPs), namely rs7759938, rs314277, and rs221639, in 24 populations. These SNPs have been recently related to the age at menarche, pubertal height growth, peripubertal body mass index, levels of prenatal testosterone exposure, and cancer survival.

METHODS

A total of 1,197 DNA samples were genotyped. The allele frequencies were used to determine the relationship between populations, with data from the 1000 Genomes Project being used for external comparisons. The genotype distributions and the population structure between populations and groups of populations were determined.

RESULTS

The population results indicate a significant degree of variation (F_ST  = 0.043, P < 0.0001). Allele frequencies show significant differences among populations. A hierarchical variance analysis is consistent with a primary differentiation between populations on the North and South coasts of the Mediterranean. This difference is especially evident in the unexpected distribution of the SNP rs221639, which shows one of the highest F_ST (11.5%, P < 0.0001) values described in the Mediterranean region thus far.

CONCLUSION

The population differentiation and the structuring of the genetic variance, in agreement with previous studies, indicate that the SNPs in question are good tools for the study of human populations, even at a microgeographic level. Am. J. Hum. Biol. 28:905-912, 2016. © 2016Wiley Periodicals, Inc.

Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes

BACKGROUND

Specific language impairment (SLI) is a common neurodevelopmental disorder, observed in 5-10 % of children. Family and twin studies suggest a strong genetic component, but relatively few candidate genes have been reported to date. A recent genome-wide association study (GWAS) described the first statistically significant association specifically for a SLI cohort between a missense variant (rs4280164) in the NOP9 gene and language-related phenotypes under a parent-of-origin model. Replications of these findings are particularly challenging because the availability of parental DNA is required.

METHODS

We used two independent family-based cohorts characterised with reading- and language-related traits: a longitudinal cohort (n = 106 informative families) including children with language and reading difficulties and a nuclear family cohort (n = 264 families) selected for dyslexia.

RESULTS

We observed association with language-related measures when modelling for parent-of-origin effects at the NOP9 locus in both cohorts: minimum P = 0.001 for phonological awareness with a paternal effect in the first cohort and minimum P = 0.0004 for irregular word reading with a maternal effect in the second cohort. Allelic and parental trends were not consistent when compared to the original study.

CONCLUSIONS

A parent-of-origin effect at this locus was detected in both cohorts, albeit with different trends. These findings contribute in interpreting the original GWAS report and support further investigations of the NOP9 locus and its role in language-related traits. A systematic evaluation of parent-of-origin effects in genetic association studies has the potential to reveal novel mechanisms underlying complex traits.

UGT2B4 previously implicated in the risk of breast cancer is associated with menarche timing in Ukrainian females

Age at menarche (AAM) is a multifactorial trait that is regulated by dozens environmental and genetic factors. Recent meta-analysis of GWAS showed significant association of 106 loci with AAM. These polymorphisms need replicating in different ethnic populations in order to confirm their association with menarche timing. This study was aimed to replicate 53 polymorphisms that were previously associated with AAM. DNA samples were collected from 416 Ukrainian young females for further genotyping. After data quality control 47 polymorphisms remained for the association analysis using the linear regression model. SNP rs13111134 located in UGT2B4 showed the most significant association with AAM (0.431years per allele A, padj=0.044 after the Bonferroni correction). Polymorphisms rs7589318 in POMC, rs11724758 in FABP2, rs7753051 in IGF2R, rs2288696 in FGFR1 and rs12444979 in GPRC5B may also contribute to menarche timing. However, none of these associations remained significant after the Bonferroni correction for multiple testing. The obtained results provide evidence that UGT2B4, which was previously associated with predisposition to breast cancer, may play a role in the onset of menarche.

Plasma levels of vitamin K and the risk of ischemic heart disease: a Mendelian randomization study

Essentials Vitamin K plays a role in coagulation, and deficiency may promote coronary artery calcification. The role of vitamin K1 in heart disease was assessed using Mendelian randomization in Caucasians. Genetically higher vitamin K1 was associated with a higher risk of ischemic heart disease. Further research elucidating the role of vitamin K1 in ischemic heart disease could be useful.

SUMMARY

Background Vitamin K1 is a nutrient in green leafy vegetables; deficiency may promote coronary artery calcification. Warfarin, an anticoagulant used in secondary prevention of thrombotic events, is a vitamin K antagonist. Thrombotic and coronary events may share risk factors. Objectives To clarify the role of vitamin K1 in ischemic heart disease, the risk of coronary artery disease/myocardial infarction (CAD/MI) was assessed according to genetically determined vitamin K1 levels. Given vitamin K1 is fat soluble, associations with lipids were similarly assessed to assess pleotropic effects via lipids. Methods Separate sample instrumental variable analysis with genetic instruments (Mendelian randomization) was used to obtain an unconfounded estimate of the association of vitamin K1 (based on rs2108622 [CYP4F2], rs4645543 [KCNK9] and rs2192574 [CTNNA2] from a genome-wide association study) with CAD/MI using CARDIoGRAMplusC4D (cases = 64 374; controls = 130 681) and with lipids using Global Lipids Genetics Consortium Results (n = 196 475). Results Vitamin K1 single nucleotide polymorphisms were positively associated with CAD/MI (odds ratio [OR], 1.17 per unit [nmol L(-1) ] of natural log-transformed genetically predicted vitamin K1 ; 95% confidence interval [CI], 1.08-1.26), but not with inverse normal transformed low-density lipoprotein cholesterol (-0.0003; 95% CI, -0.03 to 0.03), high-density lipoprotein cholesterol (0.02; 95% CI, -0.01 to 0.05) or triglycerides (-0.01; 95% CI, -0.04 to 0.02). Considering only rs2108622, which is functionally relevant to vitamin K1 , the association for CAD/MI was stronger (OR, 1.21; 95% CI, 1.08-1.36). Conclusions Vitamin K may cause CAD/MI; whether vitamin K or other determinants of coagulation could be relevant to primary prevention might be worth considering.

DNA methylation of LINE-1 and Alu repetitive elements in relation to sex hormones and pubertal timing in Mexican-American children

BACKGROUND

The molecular mechanisms linking environmental exposures to earlier pubertal development are not well characterized. Epigenetics may play an important role, but data on the relationship between epigenetic marks and puberty, particularly in humans, is limited.

METHODS

We used pyrosequencing to measure Alu and long interspersed nucleotide elements (LINE-1) methylation in DNA isolated from whole blood samples collected from newborns and 9-y-old children (n = 266). Tanner staging was completed six times between ages 9 and 12 y to determine pubertal status, and hormone levels were measured in 12-y-old boys.

RESULTS

Among girls, we observed a suggestive trend of increased odds of breast and pubic hair development with higher Alu and LINE-1 methylation in 9-y-old blood, respectively. The strongest association identified was an inverse association of LINE-1 methylation in 9-y-old girls with odds of experiencing menarche by age 12 (OR (95% CI): 0.63 (0.46, 0.87); P = 0.005). We observed a consistent inverse relationship for Alu and LINE-1 methylation at 9 y with luteinizing hormone (LH), testosterone and follicle-stimulating hormone levels in boys but it was only significant between LINE-1 and LH.

CONCLUSION

DNA methylation of Alu and LINE-1 may be involved in puberty initiation and development. This relationship should be confirmed in future studies.

MANAGEMENT OF ENDOCRINE DISEASE: Reversible hypogonadotropic hypogonadism

Congenital hypogonadotropic hypogonadism (CHH) is characterized by lack of puberty and infertility. Traditionally, it has been considered a life-long condition yet cases of reversibility have been described wherein patients spontaneously recover function of the reproductive axis following treatment. Reversibility occurs in both male and female CHH cases and appears to be more common (~10-15%) than previously thought. These reversal patients span a range of GnRH deficiency from mild to severe and many reversal patients harbor mutations in genes underlying CHH. However, to date there are no clear factors for predicting reversible CHH. Importantly, recovery of reproductive axis function may not be permanent. Thus, CHH is not always life-long and the incidence of reversal warrants periodic treatment withdrawal with close monitoring and follow-up. Reversible CHH highlights the importance of environmental (epigenetic) factors such as sex steroid treatment on the reproductive axis in modifying the phenotype. This review provides an overview and an update on what is known about this phenomenon.

SLCO1B1*5 polymorphism (rs4149056) is associated with chemotherapy-induced amenorrhea in premenopausal women with breast cancer: a prospective cohort study

Background

Because inheritance is recognized as playing a role in age at menarche and natural menopause, the development of chemotherapy-induced amenorrhea (CIA) might depend on inherited genetic factors; however, studies that explore such a correlation are few and have received scant attention. Given the importance of this topic we conducted a comprehensive genotype study in young women (≤45 years) with early-stage breast cancer.

Methods

Our approach tested the effect of variant polymorphisms in drug metabolism enzymes (DMEs) using a predesigned pharmacogenomics panel (TaqMan® OpenArray®, Life Technologies GmbH, Darmstadt, Germany) in premenopausal women (n = 50). Patients received contemporary chemotherapy; in all cases a cyclophosphamide-based regimen with a dose of at least 500 mg/m² for six cycles. CIA was considered to be present in women with no resumption of menstrual bleeding within 12 months after completion of chemotherapy or goserelin.

Results

Twenty-six patients (52 %) showed CIA during follow-up whereas 24 women (48 %) remained premenopausal. Of all the DMEs studied, only the SLCO1B1*5 (rs4149056) genotype was associated with the development of CIA (P = 0.017). Of the 26 patients who were homozygous for the T/T allele SLCO1B1*5, 18 (69.2 %) developed CIA compared with 8 (30.8 %) of the 22 patients who were heterozygous (C/T allele). The association of heterozygous SLCO1B1*5 allele (OR 0.038; 95%CI: 0.05–0.92) with a lower risk of developing CIA remained significant in a binary logistic regression analysis that include age, SLCO1B1*5 allele variants, and goserelin therapy.

Conclusions

Patient age and SLCO1B1*5 allele variants predict the likelihood of young women with breast cancer developing CIA.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2373-3) contains supplementary material, which is available to authorized users.

Identification of Common Genetic Variants Influencing Spontaneous Dizygotic Twinning and Female Fertility

Spontaneous dizygotic (DZ) twinning occurs in 1%-4% of women, with familial clustering and unknown physiological pathways and genetic origin. DZ twinning might index increased fertility and has distinct health implications for mother and child. We performed a GWAS in 1,980 mothers of spontaneous DZ twins and 12,953 control subjects. Findings were replicated in a large Icelandic cohort and tested for association across a broad range of fertility traits in women. Two SNPs were identified (rs11031006 near FSHB, p = 1.54 × 10(-9), and rs17293443 in SMAD3, p = 1.57 × 10(-8)) and replicated (p = 3 × 10(-3) and p = 1.44 × 10(-4), respectively). Based on ∼90,000 births in Iceland, the risk of a mother delivering twins increased by 18% for each copy of allele rs11031006-G and 9% for rs17293443-C. A higher polygenic risk score (PRS) for DZ twinning, calculated based on the results of the DZ twinning GWAS, was significantly associated with DZ twinning in Iceland (p = 0.001). A higher PRS was also associated with having children (p = 0.01), greater lifetime parity (p = 0.03), and earlier age at first child (p = 0.02). Allele rs11031006-G was associated with higher serum FSH levels, earlier age at menarche, earlier age at first child, higher lifetime parity, lower PCOS risk, and earlier age at menopause. Conversely, rs17293443-C was associated with later age at last child. We identified robust genetic risk variants for DZ twinning: one near FSHB and a second within SMAD3, the product of which plays an important role in gonadal responsiveness to FSH. These loci contribute to crucial aspects of reproductive capacity and health.

Recent advances in the understanding and management of delayed puberty

Delayed puberty, especially in boys, is a common presentation in paediatrics. Recent advances have improved our understanding of the neuroendocrine, genetic and environmental factors controlling pubertal development, and hence inform the pathophysiology of delayed puberty. The discovery of kisspeptin signalling through its receptor identified neuroendocrine mechanisms controlling the gonadotrophin-releasing hormone (GnRH) pulse generator at the onset of puberty. Genetic mechanisms from single gene mutations to single nucleotide polymorphism associated with delayed puberty are being identified. Environmental factors, including nutritional factors and endocrine disruptors, have also been implicated in changes in secular trends and abnormal timing of puberty. Despite these advances, the key clinical question is to distinguish delayed puberty associated with an underlying pathology or hypogonadism from constitutional delay in growth and puberty, which remains challenging as biochemical tests are not always discriminatory. The diagnostic accuracies of newer investigations, including 36-hour luteinising hormone releasing hormone (LHRH) tests, GnRH-agonist tests, antimullerian hormone and inhibin-B, require further evaluation. Sex hormone replacement remains the main available treatment for delayed puberty, the choice of which is largely dictated by clinical practice and availability of the various sex steroid preparations. Spontaneous reversal of hypogonadism has been reported in boys with idiopathic hypogonadotrophic hypogonadism after a period of sex steroid treatment, highlighting the importance of reassessment at the end of pubertal induction. Novel therapies with a more physiological basis such as gonadotrophins or kisspeptin-agonist are being investigated for the management of hypogonadotrophic hypogonadism. Careful clinical assessment and appreciation of the normal physiology remain the key approach to patients with delayed puberty.

New Perspectives on Genomic Imprinting, an Essential and Multifaceted Mode of Epigenetic Control in the Developing and Adult Brain

Mammalian evolution entailed multiple innovations in gene regulation, including the emergence of genomic imprinting, an epigenetic regulation leading to the preferential expression of a gene from its maternal or paternal allele. Genomic imprinting is highly prevalent in the brain, yet, until recently, its central roles in neural processes have not been fully appreciated. Here, we provide a comprehensive survey of adult and developmental brain functions influenced by imprinted genes, from neural development and wiring to synaptic function and plasticity, energy balance, social behaviors, emotions, and cognition. We further review the widespread identification of parental biases alongside monoallelic expression in brain tissues, discuss their potential roles in dosage regulation of key neural pathways, and suggest possible mechanisms underlying the dynamic regulation of imprinting in the brain. This review should help provide a better understanding of the significance of genomic imprinting in the normal and pathological brain of mammals including humans.

Polymorphisms and genes associated with puberty in heifers

Puberty onset is a multifactorial process influenced by genetic determinants and environmental conditions, especially nutritional status. Genes, genetic variations, and regulatory networks compose the molecular basis of achieving puberty. In this article, we reviewed the discovery of multiple polymorphisms and genes associated with heifer puberty phenotypes and discuss the opportunities to use this evolving knowledge of genetic determinants for breeding early pubertal Bos indicus-influenced cattle. The discovery of polymorphisms and genes was mainly achieved through candidate gene studies, quantitative trait loci analyses, genome-wide association studies, and recently, global gene expression studies (transcriptome). These studies are recapitulated and summarized in the current review.

Physical and neurobehavioral determinants of reproductive onset and success

The ages of puberty, first sexual intercourse and first birth signify the onset of reproductive ability, behavior and success, respectively. In a genome-wide association study of 125,667 UK Biobank participants, we identify 38 loci associated (P < 5 × 10(-8)) with age at first sexual intercourse. These findings were taken forward in 241,910 men and women from Iceland and 20,187 women from the Women's Genome Health Study. Several of the identified loci also exhibit associations (P < 5 × 10(-8)) with other reproductive and behavioral traits, including age at first birth (variants in or near ESR1 and RBM6-SEMA3F), number of children (CADM2 and ESR1), irritable temperament (MSRA) and risk-taking propensity (CADM2). Mendelian randomization analyses infer causal influences of earlier puberty timing on earlier first sexual intercourse, earlier first birth and lower educational attainment. In turn, likely causal consequences of earlier first sexual intercourse include reproductive, educational, psychiatric and cardiometabolic outcomes.

Pubertal development and prostate cancer risk: Mendelian randomization study in a population-based cohort

BACKGROUND

Epidemiological studies have observed a positive association between an earlier age at sexual development and prostate cancer, but markers of sexual maturation in boys are imprecise and observational estimates are likely to suffer from a degree of uncontrolled confounding. To obtain causal estimates, we examined the role of pubertal development in prostate cancer using genetic polymorphisms associated with Tanner stage in adolescent boys in a Mendelian randomization (MR) approach.

METHODS

We derived a weighted genetic risk score for pubertal development, combining 13 SNPs associated with male Tanner stage. A higher score indicated a later puberty onset. We examined the association of this score with prostate cancer risk, stage and grade in the UK-based ProtecT case-control study (n = 2,927), and used the PRACTICAL consortium (n = 43,737) as a replication sample.

RESULTS

In ProtecT, the puberty genetic score was inversely associated with prostate cancer grade (odds ratio (OR) of high- vs. low-grade cancer, per tertile of the score: 0.76; 95 % CI, 0.64-0.89). In an instrumental variable estimation of the causal OR, later physical development in adolescence (equivalent to a difference of one Tanner stage between pubertal boys of the same age) was associated with a 77 % (95 % CI, 43-91 %) reduced odds of high Gleason prostate cancer. In PRACTICAL, the puberty genetic score was associated with prostate cancer stage (OR of advanced vs. localized cancer, per tertile: 0.95; 95 % CI, 0.91-1.00) and prostate cancer-specific mortality (hazard ratio amongst cases, per tertile: 0.94; 95 % CI, 0.90-0.98), but not with disease grade.

CONCLUSIONS

Older age at sexual maturation is causally linked to a reduced risk of later prostate cancer, especially aggressive disease.

Discordant Haplotype Sequencing Identifies Functional Variants at the 2q33 Breast Cancer Risk Locus

The findings from genome-wide association studies hold enormous potential for novel insight into disease mechanisms. A major challenge in the field is to map these low-risk association signals to their underlying functional sequence variants (FSV). Simple sequence study designs are insufficient, as the vast numbers of statistically comparable variants and a limited knowledge of noncoding regulatory elements complicate prioritization. Furthermore, large sample sizes are typically required for adequate power to identify the initial association signals. One important question is whether similar sample sizes need to be sequenced to identify the FSVs. Here, we present a proof-of-principle example of an extreme discordant design to map FSVs within the 2q33 low-risk breast cancer locus. Our approach employed DNA sequencing of a small number of discordant haplotypes to efficiently identify candidate FSVs. Our results were consistent with those from a 2,000-fold larger, traditional imputation-based fine-mapping study. To prioritize further, we used expression-quantitative trait locus analysis of RNA sequencing from breast tissues, gene regulation annotations from the ENCODE consortium, and functional assays for differential enhancer activities. Notably, we implicate three regulatory variants at 2q33 that target CASP8 (rs3769823, rs3769821 in CASP8, and rs10197246 in ALS2CR12) as functionally relevant. We conclude that nested discordant haplotype sequencing is a promising approach to aid mapping of low-risk association loci. The ability to include more efficient sequencing designs into mapping efforts presents an opportunity for the field to capitalize on the potential of association loci and accelerate translation of association signals to their underlying FSVs. Cancer Res; 76(7); 1916-25. ©2016 AACR.

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.

Analyzing the Role of MicroRNAs in Schizophrenia in the Context of Common Genetic Risk Variants

IMPORTANCE

The recent implication of 108 genomic loci in schizophrenia marked a great advancement in our understanding of the disease. Against the background of its polygenic nature there is a necessity to identify how schizophrenia risk genes interplay. As regulators of gene expression, microRNAs (miRNAs) have repeatedly been implicated in schizophrenia etiology. It is therefore of interest to establish their role in the regulation of schizophrenia risk genes in disease-relevant biological processes.

OBJECTIVE

To examine the role of miRNAs in schizophrenia in the context of disease-associated genetic variation.

DESIGN, SETTING, AND PARTICIPANTS

The basis of this study was summary statistics from the largest schizophrenia genome-wide association study meta-analysis to date (83 550 individuals in a meta-analysis of 52 genome-wide association studies) completed in 2014 along with publicly available data for predicted miRNA targets. We examined whether schizophrenia risk genes were more likely to be regulated by miRNA. Further, we used gene set analyses to identify miRNAs that are regulators of schizophrenia risk genes.

MAIN OUTCOMES AND MEASURES

Results from association tests for miRNA targetomes and related analyses.

RESULTS

In line with previous studies, we found that similar to other complex traits, schizophrenia risk genes were more likely to be regulated by miRNAs (P < 2 × 10-16). Further, the gene set analyses revealed several miRNAs regulating schizophrenia risk genes, with the strongest enrichment for targets of miR-9-5p (P = .0056 for enrichment among the top 1% most-associated single-nucleotide polymorphisms, corrected for multiple testing). It is further of note that MIR9-2 is located in a genomic region showing strong evidence for association with schizophrenia (P = 7.1 × 10-8). The second and third strongest gene set signals were seen for the targets of miR-485-5p and miR-137, respectively.

CONCLUSIONS AND RELEVANCE

This study provides evidence for a role of miR-9-5p in the etiology of schizophrenia. Its implication is of particular interest as the functions of this neurodevelopmental miRNA tie in with established disease biology: it has a regulatory loop with the fragile X mental retardation homologue FXR1 and regulates dopamine D2 receptor density.

Constitutional delay of puberty: presentation and inheritance pattern in 48 familial cases

Background

The mechanism that initiates the onset of puberty is largely unknown but the age of onset is mainly under genetic control and influenced by environmental factors including nutrition. Familial forms of constitutional delay of puberty (CDP) suggest the involvement of genetic factors. The purpose of this study is to describe the presentation and the mode of inheritance of CDP in a series of familial cases.

Methods

A retrospective, single center study was carried out over 10 years on 48 probands (14 girls and 34 boys) from 48 families seen for CDP with a familial component.

Results

Of the 48 probands, 46 (96 %) had at least one affected 1^st degree relatives and 2 (4 %, 2 boys) had only 2^nd degree relatives affected.

In girls, 11 families (79 %) exhibited exclusive maternal inheritance, 1 (7 %) paternal inheritance and 2 (14 %) both maternal and paternal inheritance. In boys, 14 families (41 %) exhibited exclusive maternal inheritance, 12 (35 %) paternal inheritance and 8 (24 %) both maternal and paternal inheritance.

In the boys with bilineal inheritance, the ages at onset of puberty (16 ± 1.41 years) and at evaluation (16.05 ± 2.47 years) were higher than in those with unilineal inheritance (15.25 ± 0.35 and 15.1 ± 0.42 years respectively), but the difference was not significant.

Conclusions

In girls exclusive maternal inheritance seems to be the major mode of inheritance whereas for boys the mode of inheritance was almost equally maternal, paternal or bilineal. Clinical phenotype of boys with bilineal inheritance seems to be more severe, but the difference did not reach statistical significance, perhaps because of the small sample size. This greater severity of the phenotype in boys with bilineal inheritance is likely due to inheriting different puberty timing genes from each parent. Future research should be directed at identifying such genes.

Age at menarche and the risk of diabetic microvascular complications in patients with type 1 diabetes

AIMS/HYPOTHESIS

The aim of this study was to evaluate the relationship among age at onset of diabetes, age at onset of menarche and risk of diabetic nephropathy and laser-treated retinopathy in type 1 diabetes.

METHODS

Data related to age at menarche were collected through questionnaires and were available for 1,304 women who participated in the Finnish Diabetic Nephropathy Study (FinnDiane). A possible association between age at menarche and diabetic nephropathy and retinopathy was investigated.

RESULTS

There was an inverse relationship between the age at onset of diabetes and age at menarche: the younger the age at onset of diabetes, the higher the age at menarche (p < 0.0001). A non-linear relationship between the age of menarche and risk of diabetic microvascular complications was found in patients with diabetes onset before menarche, but there was no such association in patients with diabetes onset after menarche. Women with delayed menarche (> mean age + 2 years) had a 2.30 (95% CI 1.27, 4.17; p < 0.006) times higher risk of nephropathy compared with the women who underwent menarche at the mean age ± 2 years. Delayed menarche also increased the risk of retinopathy (OR 2.34 [95% CI 1.36, 4.01]). After excluding patients with nephropathy, the OR for retinopathy was 2.11 (95% CI 1.15, 3.90). Earlier menarche (< mean age - 2 years) did not have any effect on this risk.

CONCLUSIONS/INTERPRETATION

Delayed menarche was associated with an increased risk of diabetic nephropathy and retinopathy, whereas early menarche was not. Delayed menarche may be used as a new tool to identify women at risk of diabetic microvascular complications.

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.