Genome-wide association study of sexual maturation in males and females highlights a role for body mass and menarche loci in male puberty

Sex-specific developmental gene expression atlas unveils dimorphic gene networks in C. elegans

Sex-specific traits and behaviors emerge during development by the acquisition of unique properties in the nervous system of each sex. However, the genetic events responsible for introducing these sex-specific features remain poorly understood. In this study, we create a comprehensive gene expression atlas of pure populations of hermaphrodites and males of the nematode Caenorhabditis elegans across development. We discover numerous differentially expressed genes, including neuronal gene families like transcription factors, neuropeptides, and G protein-coupled receptors. We identify INS-39, an insulin-like peptide, as a prominent male-biased gene expressed specifically in ciliated sensory neurons. We show that INS-39 serves as an early-stage male marker, facilitating the effective isolation of males in high-throughput experiments. Through complex and sex-specific regulation, ins-39 plays pleiotropic sexually dimorphic roles in various behaviors, while also playing a shared, dimorphic role in early life stress. This study offers a comparative sexual and developmental gene expression database for C. elegans. Furthermore, it highlights conserved genes that may underlie the sexually dimorphic manifestation of different human diseases.

Association of obesity and menarche SNPs and interaction with environmental factors on precocious puberty

BACKGROUND

Obesity is an important cause for the precocious or early puberty. However, the association between obesity-related loci and the risk of precocious puberty as well as the effect of gene-environment interaction are unclear, especially in the Chinese children population.

METHODS

This was a case-control study using baseline data from two cohorts and hospital cases in China. 15 SNPs loci and several environmental factors were included in the analysis of 1201 participants. Chi-square test and logistic regression were used to analyze the association between SNPs and precocious puberty. Additionally, exploratory factor analysis was conducted on 13 environmental variables, and then to explore their interaction with genes on precocious puberty.

RESULTS

The effect allele C of rs571312, and G of rs12970134 MC4R were associated with precocious puberty in girls with obesity. Regarding the gene-environment interaction, we found that when girls were in the high socioeconomic status, the rs571312 (OR: 3.996; 95% CI: 1.694-9.423) and rs12970134 (OR: 3.529; 95% CI: 1.452-8.573) risk genotypes had a greater effect on precocious puberty.

CONCLUSIONS

The obesity risk gene polymorphisms MC4R rs571312 and rs12970134 were associated with precocious puberty in Chinese girls with obesity, and girls with risk genotypes and high socioeconomic status should be given extra attention.

IMPACT

This is the first study that identified the association between rs571312 and rs12970134 of MC4R gene and precocious puberty in Chinese children. We found that when girls were in the high socioeconomic status, the risk genotypes of rs571312 and rs12970134 had a greater effect on precocious puberty. The results of this study have great public health implications. It is recommended that girls who are in high socioeconomic status and have a high genetic risk for early sexual maturity should closely monitor their pubertal development and consider early intervention strategies.

Maternal Age at Menarche Genes Determines Fetal Growth Restriction Risk

We aimed to explore the potential link of maternal age at menarche (mAAM) gene polymorphisms with risk of the fetal growth restriction (FGR). This case (FGR)-control (FGR free) study included 904 women (273 FGR and 631 control) in the third trimester of gestation examined/treated in the Departments of Obstetrics. For single nucleotide polymorphism (SNP) multiplex genotyping, 50 candidate loci of mAAM were chosen. The relationship of mAAM SNPs and FGR was appreciated by regression procedures (logistic/model-based multifactor dimensionality reduction [MB-MDR]) with subsequent in silico assessment of the assumed functionality pithy of FGR-related loci. Three mAAM-appertain loci were FGR-linked to genes such as KISS1 (rs7538038) (effect allele G-odds ratio (OR)allelic = 0.63/pperm = 0.0003; ORadditive = 0.61/pperm = 0.001; ORdominant = 0.56/pperm = 0.001), NKX2-1 (rs999460) (effect allele A-ORallelic = 1.37/pperm = 0.003; ORadditive = 1.45/pperm = 0.002; ORrecessive = 2.41/pperm = 0.0002), GPRC5B (rs12444979) (effect allele T-ORallelic = 1.67/pperm = 0.0003; ORdominant = 1.59/pperm = 0.011; ORadditive = 1.56/pperm = 0.009). The haplotype ACA FSHB gene (rs555621*rs11031010*rs1782507) was FRG-correlated (OR = 0.71/pperm = 0.05). Ten FGR-implicated interworking models were founded for 13 SNPs (pperm ≤ 0.001). The rs999460 NKX2-1 and rs12444979 GPRC5B interplays significantly influenced the FGR risk (these SNPs were present in 50% of models). FGR-related mAAM-appertain 15 polymorphic variants and 350 linked SNPs were functionally momentous in relation to 39 genes participating in the regulation of hormone levels, the ovulation cycle process, male gonad development and vitamin D metabolism. Thus, this study showed, for the first time, that the mAAM-appertain genes determine FGR risk.

Metabolic control of puberty: 60 years in the footsteps of Kennedy and Mitra's seminal work

An individual's nutritional status has a powerful effect on sexual maturation. Puberty onset is delayed in response to chronic energy insufficiency and is advanced under energy abundance. The consequences of altered pubertal timing for human health are profound. Late puberty increases the chances of cardiometabolic, musculoskeletal and neurocognitive disorders, whereas early puberty is associated with increased risks of adult obesity, type 2 diabetes mellitus, cardiovascular diseases and various cancers, such as breast, endometrial and prostate cancer. Kennedy and Mitra's trailblazing studies, published in 1963 and using experimental models, were the first to demonstrate that nutrition is a key factor in puberty onset. Building on this work, the field has advanced substantially in the past decade, which is largely due to the impressive development of molecular tools for experimentation and population genetics. In this Review, we discuss the latest advances in basic and translational sciences underlying the nutritional and metabolic control of pubertal development, with a focus on perspectives and future directions.

Association Between Markers of Adiposity During Childhood and Puberty Onset in Latino Girls

CONTEXT

Prepubertal adiposity is associated with earlier puberty. It is unclear when this association starts, if all adiposity markers are similarly associated, and whether all pubertal milestones are similarly affected.

OBJECTIVE

To evaluate the association between different adiposity markers during childhood and the timing of different pubertal milestones in Latino girls.

DESIGN, SETTING, AND PARTICIPANTS

Longitudinal follow-up of 539 female participants of the Chilean Growth and Obesity Cohort recruited from childcare centers (mean age 3.5 years) from the southeast area of Santiago, Chile. Participants were singletons born between 2002 and 2003 within the normal birthweight range. Since 2006, a trained dietitian measured weight, height, waist circumference (WC) and skinfolds to estimate body mass index (BMI) Centers for Disease Control and Prevention percentiles, central obesity, percentage of fat mass (%FM), and fat mass index (FMI, fat mass/height2).

MAIN OUTCOME

Since 2009, sexual maturation was assessed every 6 months to assess age at (1) thelarche, (2) pubarche, (3) menarche, and (4) peak height velocity (PHV).

RESULTS

At thelarche, 12.5% were obese and 2% had central obesity. The median age of pubarche, menarche, and PHV were all associated with markers of adiposity at different time points during childhood whereas thelarche only with %FM and FMI. Adiposity clusters models showed that children with trajectories of high WC, %FM, and FMI during childhood were related with earlier thelarche, pubarche, menarche, and PHV but BMI trajectories only with menarche and PHV.

CONCLUSIONS

Higher WC, %FM, and FMI were associated with earlier age at thelarche, pubarche, menarche, and PHV. The effect of BMI was less consistent.

Maternal Age at Menarche Gene Polymorphisms Are Associated with Offspring Birth Weight

In this study, the association between maternal age at menarche (AAM)-related polymorphisms and offspring birth weight (BW) was studied. The work was performed on a sample of 716 pregnant women and their newborns. All pregnant women underwent genotyping of 50 SNPs of AAM candidate genes. Regression methods (linear and Model-Based Multifactor Dimensionality Reduction (MB-MDR)) with permutation procedures (the indicator p_perm was calculated) were used to identify the correlation between SNPs and newborn weight (transformed BW values were analyzed) and in silico bioinformatic examination was applied to assess the intended functionality of BW-associated loci. Four AAM-related genetic variants were BW-associated including genes such as POMC (rs7589318) (β_additive = 0.202/p_perm = 0.015), KDM3B (rs757647) (β_recessive = 0.323/p_perm = 0.005), INHBA (rs1079866) (β_additive = 0.110/p_perm = 0.014) and NKX2-1 (rs999460) (β_recessive = -0.176/p_perm = 0.015). Ten BW-significant models of interSNPs interactions (p_perm ≤ 0.001) were identified for 20 polymorphisms. SNPs rs7538038 KISS1, rs713586 RBJ, rs12324955 FTO and rs713586 RBJ-rs12324955 FTO two-locus interaction were included in the largest number of BW-associated models (30% models each). BW-associated AAM-linked 22 SNPs and 350 proxy loci were functionally related to 49 genes relevant to pathways such as the hormone biosynthesis/process and female/male gonad development. In conclusion, maternal AMM-related genes polymorphism is associated with the offspring BW.

Functional genomics of reproduction in pigs: Are we there yet?

Reproductive failure is the main reason for culling females in swine herds and is both a financial and sustainability issue. Because reproductive traits are complex and lowly to moderately heritable, genomic selection within populations can achieve substantial genetic gain in reproductive efficiency. A better understanding of the physiological components affecting the expression of these traits will facilitate greater understanding of the genes affecting reproductive traits and is necessary to improve and optimize management strategies to maximize reproductive success of gilts and sows. Large-scale genotyping with single-nucleotide polymorphism (SNP) arrays are used for genome-wide association studies (GWAS) and have facilitated identification of positional candidate genes. Transcriptomic data can be used to weight SNP for GWAS and could lead to previously unidentified candidate genes. Resequencing and fine mapping of candidate genes are necessary to identify putative functional variants and some of these have been incorporated into new genotyping arrays. Sequence imputation and genotype by sequence are newer strategies that could reveal novel functional mutations. In this study, these approaches are discussed. Advantages and limitations are highlighted where additional research is needed.

Endocrine Health and Health Care Disparities in the Pediatric and Sexual and Gender Minority Populations: An Endocrine Society Scientific Statement

Endocrine care of pediatric and adult patients continues to be plagued by health and health care disparities that are perpetuated by the basic structures of our health systems and research modalities, as well as policies that impact access to care and social determinants of health. This scientific statement expands the Society's 2012 statement by focusing on endocrine disease disparities in the pediatric population and sexual and gender minority populations. These include pediatric and adult lesbian, gay, bisexual, transgender, queer, intersex, and asexual (LGBTQIA) persons. The writing group focused on highly prevalent conditions-growth disorders, puberty, metabolic bone disease, type 1 (T1D) and type 2 (T2D) diabetes mellitus, prediabetes, and obesity. Several important findings emerged. Compared with females and non-White children, non-Hispanic White males are more likely to come to medical attention for short stature. Racially and ethnically diverse populations and males are underrepresented in studies of pubertal development and attainment of peak bone mass, with current norms based on European populations. Like adults, racial and ethnic minority youth suffer a higher burden of disease from obesity, T1D and T2D, and have less access to diabetes treatment technologies and bariatric surgery. LGBTQIA youth and adults also face discrimination and multiple barriers to endocrine care due to pathologizing sexual orientation and gender identity, lack of culturally competent care providers, and policies. Multilevel interventions to address these disparities are required. Inclusion of racial, ethnic, and LGBTQIA populations in longitudinal life course studies is needed to assess growth, puberty, and attainment of peak bone mass. Growth and development charts may need to be adapted to non-European populations. In addition, extension of these studies will be required to understand the clinical and physiologic consequences of interventions to address abnormal development in these populations. Health policies should be recrafted to remove barriers in care for children with obesity and/or diabetes and for LGBTQIA children and adults to facilitate comprehensive access to care, therapeutics, and technological advances. Public health interventions encompassing collection of accurate demographic and social needs data, including the intersection of social determinants of health with health outcomes, and enactment of population health level interventions will be essential tools.

The Congenital and Acquired Mechanisms Implicated in the Etiology of Central Precocious Puberty

The etiology of central precocious puberty (CPP) is multiple and heterogeneous, including congenital and acquired causes that can be associated with structural or functional brain alterations. All causes of CPP culminate in the premature pulsatile secretion of hypothalamic GnRH and, consequently, in the premature reactivation of hypothalamic-pituitary-gonadal axis. The activation of excitatory factors or suppression of inhibitory factors during childhood represent the 2 major mechanisms of CPP, revealing a delicate balance of these opposing neuronal pathways. Hypothalamic hamartoma (HH) is the most well-known congenital cause of CPP with central nervous system abnormalities. Several mechanisms by which hamartoma causes CPP have been proposed, including an anatomical connection to the anterior hypothalamus, autonomous neuroendocrine activity in GnRH neurons, trophic factors secreted by HH, and mechanical pressure applied to the hypothalamus. The importance of genetic and/or epigenetic factors in the underlying mechanisms of CPP has grown significantly in the last decade, as demonstrated by the evidence of genetic abnormalities in hypothalamic structural lesions (eg, hamartomas, gliomas), syndromic disorders associated with CPP (Temple, Prader-Willi, Silver-Russell, and Rett syndromes), and isolated CPP from monogenic defects (MKRN3 and DLK1 loss-of-function mutations). Genetic and epigenetic discoveries involving the etiology of CPP have had influence on the diagnosis and familial counseling providing bases for potential prevention of premature sexual development and new treatment targets in the future. Global preventive actions inducing healthy lifestyle habits and less exposure to endocrine-disrupting chemicals during the lifespan are desirable because they are potentially associated with CPP.

Shared genetic basis between reproductive behaviors and anxiety-related disorders

Reproductive behaviors are associated with risks for psychiatric disorders. Reproductive phenotypes are moderately heritable and have genetic overlaps with risks for psychiatric disorders. However, the genetic and causal relationships between anxiety-related disorders or specific anxiety disorders and reproductive phenotypes remain unknown. We utilized large-scale genome-wide association study (GWAS) results (n = 9537-542,901) for five reproductive phenotypes [age at menarche, age at first sexual intercourse (AFS), age at first birth (AFB), number of children ever born (NEB), and age at menopause] and five anxiety-related disorders [panic disorder, anxiety disorders from the ANGST and the UK biobank (UKBB), posttraumatic stress disorder (PTSD) and obsessive-compulsive disorder (OCD)]. To assess genetic correlations and causal associations, linkage disequilibrium score regression and Mendelian randomization analyses, respectively, were performed. We found that AFS and AFB were negatively correlated with anxiety disorders ANGST (AFS: r_g ± SE = -0.28 ± 0.08, p = 6.00 × 10^-4; AFB: -0.45 ± 0.11, p = 3.26 × 10^-5), anxiety disorders UKBB (AFS: -0.18 ± 0.03, p = 9.64 × 10^-9; AFB; -0.25 ± 0.03, p = 2.90 × 10^-13) and PTSD (AFS: -0.42 ± 0.12, p = 4.00 × 10^-4; AFB: -0.44 ± 0.12, p = 2.00 × 10^-4) and positively correlated with OCD (AFS: 0.25 ± 0.05, p = 2.46 × 10^-6; AFB: 0.25 ± 0.05, p = 3.92 × 10^-7). Conversely, NEB was negatively correlated with OCD (-0.28 ± 0.08, p = 6.00 × 10^-4). We revealed bidirectional effects between earlier AFS and AFB and anxiety disorders (odds ratios: OR_{earlier AFS→Anxiety} = 1.64, p = 2.27 × 10^-8; OR_{earlier AFB→Anxiety} = 1.15, p = 2.28 × 10^-3; OR_{Anxiety→earlier AFS} = 1.02, p = 6.62 × 10^-8; OR_{Anxiety→earlier AFB} = 1.08, p = 1.60 × 10^-4). In contrast, we observed unidirectional effects of later AFS and AFB on OCD (OR_{later AFS→OCD} = 2.18, p = 2.16 × 10^-6; OR_{later AFB→OCD} = 1.22, p = 0.016). We suggest that those who have earlier sexual debut and childbirth are prone to risk for anxiety disorders and vice versa, while those who have later sexual debut and childbirth are genetically prone to risk for OCD. Our findings further support revising the diagnostic criteria (DSM-5) such that OCD is independent from anxiety disorders.

The Association Between Puberty Timing and Body Mass Index in a Longitudinal Setting: The Contribution of Genetic Factors

We analyzed the contribution of genetic factors on the association between puberty timing and body mass index (BMI) using longitudinal data and two approaches: (i) genetic twin design and (ii) polygenic scores (PGS) of obesity indices. Our data were derived from Finnish cohorts: 9080 twins had information on puberty timing and BMI and 2468 twins also had genetic data. Early puberty timing was moderately associated with higher BMI in childhood in both boys and girls; in adulthood these correlations were weaker and largely disappeared after adjusting for childhood BMI. The largest proportion of these correlations was attributable to genetic factors. The higher PGSs of BMI and waist circumference were associated with earlier timing of puberty in girls, whereas weaker associations were found in boys. Early puberty is not an independent risk factor for adult obesity but rather reflects the association between puberty timing and childhood BMI contributed by genetic predisposition.

Fitness, Food, and Biomarkers: Characterizing Body Composition in 19,634 Early Adolescents

Adolescent obesity persists as a major concern, especially in Central and Eastern Europe, yet evidence gaps exist regarding the pivotal early adolescent years. Our objective was to provide a comprehensive picture using a holistic approach of measured anthropometry in early adolescence, including body composition, cardiorespiratory fitness (CRF), and reported lifestyle characteristics. We aimed to elucidate potential sex/gender differences throughout and associations to biomarkers of disease risk for obese adolescents. Methods: Trained nurses measured 19,634 early adolescents (12−14-year-olds), we collected parental reports, and, for obese adolescents, fasting blood samples in four major Polish cities using a cross-sectional developmental design. Results: 24.7% boys and 18.6% girls were overweight/obese, and 2886 had BMI ≥ 90th percentile. With increasing age, there was greater risk of obesity among boys (p for trend = 0.001) and a decreasing risk of thinness for girls (p for trend = 0.01). Contrary to debate, we found BMI (continuous) was a useful indicator of measured fat mass (FM). There were 38.6% with CRF in the range of poor/very poor and was accounted for primarily by FM in boys, rather than BMI, and systolic blood pressure in girls. Boys, in comparison to girls, engaged more in sports (t = 127.26, p < 0.0001) and consumed more fast food (t = 188.57, p < 0.0001) and sugar-sweetened beverages (167.46, p < 0.0001). Uric acid, a potential marker for prediabetes, was strongly related to BMI in the obese subsample for both boys and girls. Obese girls showed signs of undernutrition. Conclusion: these findings show that overweight/obesity is by far a larger public health problem than thinness in early adolescence and is characterized differentially by sex/gender. Moreover, poor CRF in this age, which may contribute to life course obesity and disease, highlights the need for integrated and personalized intervention strategies taking sex/gender into account.

Genome-Wide Association Study and Polygenic Risk Scores of Serum DHEAS Levels in a Chilean Children Cohort

CONTEXT

Adrenarche reflects the developmental growth of the adrenal zona reticularis, which produces increasing adrenal androgen secretion (eg, dehydroepiandrosterone [DHEA]/dehydroepiandrosterone sulfate [DHEAS]) from approximately age 5 to 15 years.

OBJECTIVE

We hypothesized that the study of the genetic determinants associated with variations in serum DHEAS during adrenarche might detect genetic variants influencing the rate or timing of this process.

METHODS

Genome-wide genotyping was performed in participants of the Chilean pediatric Growth and Obesity Chilean Cohort Study (GOCS) cohort (n = 788). We evaluated the genetic determinants of DHEAS levels at the genome-wide level and in targeted genes associated with steroidogenesis. To corroborate our findings, we evaluated a polygenic risk score (PRS) for age at pubarche, based on the discovered variants, in children from the same cohort.

RESULTS

We identified one significant variant at the genome-wide level in the full cohort, close to the GALR1 gene (P = 3.81 × 10-8). In addition, variants suggestive of association (P < 1 × 10-5) were observed in PRLR, PITX1, PTPRD, NR1H4, and BCL11B. Stratifying by sex, we found variants suggestive of association in SERBP1 and CAMTA1/VAMP3 for boys and near ZNF98, TRPC6, and SULT2A1 for girls. We also found significant reductions in age at pubarche in those children with higher PRS for greater DHEAS based on these newly identified variants.

CONCLUSION

Our results disclose one variant associated with DHEAS concentrations at the level of genome-wide association study significance, and several variants with a suggestive association that may be involved in the genetic regulation of adrenarche.

Regional variation in growth status. The Peruvian health and optimist growth study

OBJECTIVE

This study aims to (1) investigate differences in weight, body mass index (BMI), and waist circumference (WC) among Peruvian children and adolescents living in three areas located at different altitudes; (2) compare age- and sex-specific height, weight, and BMI within each site with US reference data.

METHODS

We sampled 8753 subjects (4130 boys), aged 4 to 17 years from sea level, rainforest, and high-altitude. Height, weight, and WC were measured and BMI was calculated. Analysis of variance was used to compare variables across geographic regions, and the Hoff and Blackburn procedure was used to compare the Peruvian results with US reference data.

RESULTS

Participants living at sea level were taller, heavier, had greater BMI and WC relative to those living at high-altitude and in the rainforest. Peruvian schoolchildren of both sexes from the three geographical areas were shorter and lighter than their American peers. Boys and girls living in the rainforest and at high-altitude had lower BMI, whereas WC values of American schoolchildren are higher than those of the Peruvian children by age and gender.

CONCLUSIONS

Peruvians living at different altitudes differ in their growth indicators (height, weight, BMI, and WC), with significant differences between those living at sea level relative to their peers from other regions. Further, Peruvian schoolchildren of both sexes from the three geographical areas significantly differ from their US counterparts.

Transcriptome Analyses Identify Potential Key microRNAs and Their Target Genes Contributing to Ovarian Reserve

Female endocrinological symptoms, such as premature ovarian inefficiency (POI) are caused by diminished ovarian reserve and chemotherapy. The etiology of POI remains unknown, but this can lead to infertility. This has accelerated the search for master regulator genes or other molecules that contribute as enhancers or silencers. The impact of regulatory microRNAs (miRNAs) on POI has gained attention; however, their regulatory function in this condition is not well known. RNA sequencing was performed at four stages, 2-(2 W), 6-(6 W), 15-(15 W), and 20-(20 W) weeks, on ovarian tissue samples and 5058 differentially expressed genes (DEGs) were identified. Gene expression and enrichment were analyzed based on the gene ontology and KEGG databases, and their association with other proteins was assessed using the STRING database. Gene set enrichment analysis was performed to identify the key target genes. The DEGs were most highly enriched in 6 W and 15 W groups. Figla, GDF9, Nobox, and Pou51 were significantly in-creased at 2 W compared with levels at 6 W and 20 W, whereas the expression of Foxo1, Inha, and Taf4b was significantly de-creased at 20 W. Ccnd2 and Igf1 expression was maintained at similar levels in each stage. In total, 27 genes were upregulated and 26 genes interacted with miRNAs; moreover, stage-specific upregulated and downregulated interactions were demonstrated. Increased and decreased miRNAs were identified at each stage in the ovaries. The constitutively expressed genes, Ccnd2 and Igf1, were identified as the major targets of many miRNAs (p < 0.05), and Fshr and Foxo3 interacted with miRNAs, namely mmu-miR-670-3p and mmu-miR-153-3p. miR-26a-5p interacted with Piwil2, and its target genes were downregulated in the 20 W mouse ovary. In this study, we aimed to identify key miRNAs and their target genes encompassing the reproductive span of mouse ovaries using mRNA and miRNA sequencing. These results indicated that gene sets are regulated in the reproductive stage-specific manner via interaction with miRNAs. Furthermore, consistent expression of Ccnd2 and Igf1 is considered crucial for the ovarian reserve and is regulated by many interactive miRNAs.

A genome-wide association study of interhemispheric theta EEG coherence: implications for neural connectivity and alcohol use behavior

Aberrant connectivity of large-scale brain networks has been observed among individuals with alcohol use disorders (AUDs) as well as in those at risk, suggesting deficits in neural communication between brain regions in the liability to develop AUD. Electroencephalographical (EEG) coherence, which measures the degree of synchrony between brain regions, may be a useful measure of connectivity patterns in neural networks for studying the genetics of AUD. In 8810 individuals (6644 of European and 2166 of African ancestry) from the Collaborative Study on the Genetics of Alcoholism (COGA), we performed a Multi-Trait Analyses of genome-wide association studies (MTAG) on parietal resting-state theta (3-7 Hz) EEG coherence, which previously have been associated with AUD. We also examined developmental effects of GWAS findings on trajectories of neural connectivity in a longitudinal subsample of 2316 adolescent/young adult offspring from COGA families (ages 12-30) and examined the functional and clinical significance of GWAS variants. Six correlated single nucleotide polymorphisms located in a brain-expressed lincRNA (ENSG00000266213) on chromosome 18q23 were associated with posterior interhemispheric low theta EEG coherence (3-5 Hz). These same variants were also associated with alcohol use behavior and posterior corpus callosum volume, both in a subset of COGA and in the UK Biobank. Analyses in the subsample of COGA offspring indicated that the association of rs12954372 with low theta EEG coherence occurred only in females, most prominently between ages 25 and 30 (p < 2 × 10-9). Converging data provide support for the role of genetic variants on chromosome 18q23 in regulating neural connectivity and alcohol use behavior, potentially via dysregulated myelination. While findings were less robust, genome-wide associations were also observed with rs151174000 and parieto-frontal low theta coherence, rs14429078 and parieto-occipital interhemispheric high theta coherence, and rs116445911 with centro-parietal low theta coherence. These novel genetic findings highlight the utility of the endophenotype approach in enhancing our understanding of mechanisms underlying addiction susceptibility.

Delayed Puberty in Girls with Primary Amenorrhea: A Report of Cases

Background:Female puberty starts when the pituitary hormone producing follicle-stimulating hormone (FSH) and luteinizing hormone (LH), which will stimulate the ovaries to produce estrogen. Delayed puberty with primary amenorrhea in female is the lack of breast development followed by the absence of menses 3 years after the initiation of breast development. Sex chromosomes have an important role in determining the sex, germ cell differentiation of foetus, and reproductive functions of an offspring, thus, sex chromosomal aberrations frequently cause primary amenorrheaCase presentation: We report two delayed puberty cases with the chief complain of primary amenorrhea. Both cases showed hypoplasia of uterus and ovaries on pelvic imaging and hormonal assay showed low of FSH. The first case was gonadal dysgenesis with 46,XX karyotype and low level of estrogen and the second case was a turner syndrome with 45,X karyotype and normal level of estrogen. Conclusion:This study reported delayed puberty with primary amenorrhea cases due to different chromosomal aberration pattern which have similar clinical features. Therefore, cytogenetic examination is needed for any primary amenorrhea cases in order to accomplish the confirmatory diagnosis and for the clinicians to make a correct intervention and treatment. 

Large-scale collaboration in ENIGMA-EEG: A perspective on the meta-analytic approach to link neurological and psychiatric liability genes to electrophysiological brain activity

BACKGROUND AND PURPOSE

The ENIGMA-EEG working group was established to enable large-scale international collaborations among cohorts that investigate the genetics of brain function measured with electroencephalography (EEG). In this perspective, we will discuss why analyzing the genetics of functional brain activity may be crucial for understanding how neurological and psychiatric liability genes affect the brain.

METHODS

We summarize how we have performed our currently largest genome-wide association study of oscillatory brain activity in EEG recordings by meta-analyzing the results across five participating cohorts, resulting in the first genome-wide significant hits for oscillatory brain function located in/near genes that were previously associated with psychiatric disorders. We describe how we have tackled methodological issues surrounding genetic meta-analysis of EEG features. We discuss the importance of harmonizing EEG signal processing, cleaning, and feature extraction. Finally, we explain our selection of EEG features currently being investigated, including the temporal dynamics of oscillations and the connectivity network based on synchronization of oscillations.

RESULTS

We present data that show how to perform systematic quality control and evaluate how choices in reference electrode and montage affect individual differences in EEG parameters.

CONCLUSION

The long list of potential challenges to our large-scale meta-analytic approach requires extensive effort and organization between participating cohorts; however, our perspective shows that these challenges are surmountable. Our perspective argues that elucidating the genetic of EEG oscillatory activity is a worthwhile effort in order to elucidate the pathway from gene to disease liability.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Maternal Age at Menarche and Pubertal Timing in Boys and Girls: A Cohort Study From Chongqing, China

PURPOSE

This study explored the association of maternal age at menarche (AAM) with pubertal timing among girls and boys in Chongqing, China.

METHODS

Pubertal development of 1,237 children (542 girls and 695 boys) were examined half-yearly through inspection and palpation from April 2014 to June 2019. Characteristics of parents and maternal AAM were collected by a parental questionnaire at baseline. Maternal AAM was used both as a continuous and a categorical variable in Cox regression models.

RESULTS

A total of 1,198 children (528 girls and 670 boys) were included in the study. In the simple Cox model, earlier maternal AAM was associated with girls' earlier menarche, breast and pubic hair development, and boys' first ejaculation, testicular development, and genital development. When adjusting for children's body mass index z-scores (BMIz) and socioeconomic covariates, we found that girls whose mothers had early AAM had a higher risk of earlier onset of menarche (hazard ratio [HR]: .922, 95% confidence interval [CI]: .852-.998 for continuous maternal AAM, HR: 1.297, 95% CI: 1.041-1.616 for maternal AAM ≤13 years), and boys whose mother achieved menarche earlier experienced a higher risk of earlier onset of first ejaculation (HR: .896, 95% CI: .830-.968). Children's BMIz were related to all nine pubertal milestones. Parental education and relationship, birth weight, parity, and family type were also associated with pubertal timing.

CONCLUSIONS

Earlier maternal AAM was related to earlier pubertal timing in both girls and boys in Chongqing, especially girls' age at menarche and boys' first ejaculation. Children's BMIz was the most consistent factor for pubertal timing. Children's BMIz and socioeconomic conditions had greater influence on most pubertal milestones than maternal AAM.

Exploring the midline soft tissue surface changes from 12 to 15 years of age in three distinct country population cohorts

INTRODUCTION

Several studies have highlighted differences in the facial features in a White European population. Genetics appear to have a major influence on normal facial variation, and environmental factors are likely to have minor influences on face shape directly or through epigenetic mechanisms.

AIM

The aim of this longitudinal cohort study is to determine the rate of change in midline facial landmarks in three distinct homogenous population groups (Finnish, Latvian, and Welsh) from 12.8 to 15.3 years of age. This age range covers the pubertal growth period for the majority of boys and girls.

METHODS

A cohort of children aged 12 were monitored for facial growth in three countries [Finland (n = 60), Latvia (n = 107), and Wales (n = 96)]. Three-dimensional facial surface images were acquired (using either laser or photogrammetric methods) at regular intervals (6-12 months) for 4 years. Ethical approval was granted in each country. Nine midline landmarks were identified and the relative spatial positions of these surface landmarks were measured relative to the mid-endocanthion (men) over a 4-year period.

RESULTS

This study reports the children who attended 95 per cent of all scanning sessions (Finland 48 out of 60; Latvia 104 out of 107; Wales 50 out of 96). Considerable facial variation is seen for all countries and sexes. There are clear patterns of growth that show different magnitudes at different age groups for the different country groups, sexes, and facial parameters. The greatest single yearly growth rate (5.4 mm) was seen for Welsh males for men-pogonion distance at 13.6 years of age. Males exhibit greater rates of growth compared to females. These variations in magnitude and timings are likely to be influenced by genetic ancestry as a result of population migration.

CONCLUSION

The midline points are a simple and valid method to assess the relative spatial positions of facial surface landmarks. This study confirms previous reports on the subtle differences in facial shapes and sizes of male and female children in different populations and also highlights the magnitudes and timings of growth for various midline landmark distances to the men point.

Candidate Genes for Age at Menarche Are Associated With Uterine Leiomyoma

Age at menarche (AAM) is an important marker of the pubertal development and function of the hypothalamic-pituitary-ovarian system. It was reported as a possible factor for a risk of uterine leiomyoma (UL). However, while more than 350 loci for AAM have been determined by genome-wide association studies (GWASs) to date, no studies of these loci for their association with UL have been conducted so far. In this study, we analyzed 52 candidate loci for AAM for possible association with UL in a sample of 569 patients and 981 controls. The results of the study suggested that 23 out of the 52 studied polymorphisms had association with UL. Locus rs7759938 LIN28B was individually associated with the disease according to the dominant model. Twenty loci were associated with UL within 11 most significant models of intergenic interactions. Nine loci involved in 16 most significant models of interactions between single-nucleotide polymorphism (SNP), induced abortions, and chronic endometritis were associated with UL. Among the 23 loci associated with UL, 16 manifested association also with either AAM (7 SNPs) or height and/or body mass index (BMI) (13 SNPs). The above 23 SNPs and 514 SNPs linked to them have non-synonymous, regulatory, and expression quantitative trait locus (eQTL) significance for 35 genes, which play roles in the pathways related to development of the female reproductive organs and hormone-mediated signaling [false discovery rate (FDR) ≤ 0.05]. This is the first study reporting associations of candidate genes for AAM with UL.

Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS

We hereby provide the initial portrait of lincNORS, a spliced lincRNA generated by the MIR193BHG locus, entirely distinct from the previously described miR-193b-365a tandem. While inducible by low O₂ in a variety of cells and associated with hypoxia in vivo, our studies show that lincNORS is subject to multiple regulatory inputs, including estrogen signals. Biochemically, this lincRNA fine-tunes cellular sterol/steroid biosynthesis by repressing the expression of multiple pathway components. Mechanistically, the function of lincNORS requires the presence of RALY, an RNA-binding protein recently found to be implicated in cholesterol homeostasis. We also noticed the proximity between this locus and naturally occurring genetic variations highly significant for sterol/steroid-related phenotypes, in particular the age of sexual maturation. An integrative analysis of these variants provided a more formal link between these phenotypes and lincNORS, further strengthening the case for its biological relevance.

Common Genetic Influences on Age at Pubertal Voice Change and BMI in Male Twins

The present study aimed to explore secular trends in age at voice change (AVC), estimate heritability of AVC and investigate to what extent common genes influence the association between AVC and body mass index (BMI) in South Korean males. The sample of 955 male twins consisted of 241 pairs and 118 co-twin missing monozygotic (MZ) twins, 82 pairs and 50 co-twin missing dizygotic (DZ) twins and 141 male members of opposite-sex DZ twins who participated in telephone surveys in the South Korean Twin Registry. AVC was asked of twins during the surveys. The mean (SD) age of the sample was 18.92 (2.42) years (range: 16.00-29.25 years). The birth years of the twins were divided into two groups (1988-1993, 1994-2001). Kaplan-Meyer survival analyses were conducted to compute the mean age of AVC in the total sample as well as to test mean differences between the two birth cohorts. Maximum likelihood twin correlations and univariate and bivariate model-fitting analyses were performed. The mean AVC in the total sample was 14.19 (95% CI [14.09, 14.29]) years. The mean AVC significantly declined from 14.38 to 14.02 years from 1988 to 2001, confirming downward trends in AVC in recent years. Heritability for AVC was .59 (95% CI [.50, .67]), which was within the range reported in most Western twin studies. Although the phenotypic correlation between AVC and BMI was modest (r = -.14; 95% CI [-.07, -.21]), it was entirely mediated by common genes, similar to what has been found in females in prior twin studies. In conclusion, the present twin study underscores the importance of genetic influences on pubertal timing and its association with BMI in South Korean males.

Adipokines: A gear shift in puberty

In this review, we discuss the role of adipokines in the onset of puberty in children with obesity during adrenarche and gonadarche and provide a clear and detailed overview of the biological processes of two major players, leptin and adiponectin. Adipokines, especially leptin and adiponectin, seem to induce an early onset of puberty in girls and boys with obesity by affecting the hypothalamic-pituitary-gonadal (HPG) axis. Moreover, adipokines and their receptors are expressed in the gonads, suggesting a role in sexual maturation and reproduction. All in all, adipokines may be a clue in understanding mechanisms underlying the onset of puberty in childhood obesity and puberty onset variability.

Racial differences in prostate cancer: does timing of puberty play a role?

The burden of prostate cancer has a remarkably disproportionate distribution across racial groups. For example, in the USA, African Americans are twice as likely as individuals of European ancestry to develop or die from prostate cancer, and have a more aggressive disease nature at diagnosis. In contrast, Asian American men have the lowest incidence and mortality rates of prostate cancer. That considerable racial disparities exist even in the subclinical stage of prostate cancer among young men in their 20–30s suggests that patterns of prostate carcinogenesis start to diverge even earlier, perhaps during puberty, when the prostate matures at its most rapid rate. Mendelian randomisation studies have provided strong population-based evidence supporting the hypothesis that earlier onset of puberty increases the risk of prostate cancer—particularly of high grade—and prostate cancer-specific mortality later in life, observations which correspond to the epidemiology of the disease in African Americans. Notably, African American boys initiate genital development ~1 year earlier and thus go through longer periods of pubertal maturation compared with European American boys. In this perspective, bringing together existing evidence, we point to puberty as a potential critical window of increased susceptibility to prostate carcinogenesis that could account for the marked prevailing racial differences in the burden of prostate cancer.

Lin28 gene and mammalian puberty

Lin28a and Lin28b, homologs of the Caenorhabditis elegans Lin28 gene, play important roles in cell pluripotency, reprogramming, and tumorigenicity. Recently, genome-wide association and transgenic studies showed that Lin28a and/or Lin28b gene were involved in the onset of mammalian puberty, the stage representing the attainment of reproduction capacity; however, the detailed mechanism of these genes in mammalian puberty remains largely unknown. The present paper reviews the research progress on the roles of Lin28a/b genes in the onset of mammalian puberty by analyzing the results coming from gene expression patterns, mutations, and transgenic studies, and put forward possible pathways for further studies on their roles in animal reproduction.

Age at onset of different pubertal signs in boys and girls and differential DNA methylation at age 10 and 18 years: an epigenome-wide follow-up study

STUDY QUESTION

Is the age of onset of pubertal markers related to subsequent changes in DNA methylation (DNAm)?

SUMMARY ANSWER

We identified 273 cytosine-phosphate-guanine (CpG) dinucleotides in girls and 67 CpGs in boys that were related to puberty and that were replicable in two other investigations.

WHAT IS KNOWN ALREADY

Previously, 457 CpGs (not gender-specific) and 347 (in girls) and 50 (in boys), respectively, were found to be associated with puberty, according to investigations of studies from Denmark (20 girls and 31 boys) and North America (30 girls and 25 boys).

STUDY DESIGN SIZE DURATION

The study was based on a birth cohort of 1456 participants born in 1989/90, with follow-up at age 10 and 18 years.

PARTICIPANTS/MATERIALS SETTING METHODS

The follow-up included 470 participants with information on DNAm and age of pubertal onset (244 girls and 226 boys). Age of pubertal onset was ascertained retrospectively at age 18 years. Using the Pubertal Development Scale, both genders were asked about ages of onset of growth spurt, body hair growth and skin changes. Ages at voice deepening and growth of facial hair were inquired from boys; ages at breast development and menarche from girls. Blood samples were collected at 10 and 18 years of age. DNA was extracted using a standard salting out procedure. The methylation level for each CpG site was assessed using one of two different platforms. DNAm was measured by a ratio of intensities denoted as β values for each CpG site. After quality control, 349 455 CpG sites were available for analysis. M values were calculated (log2(β/(1-β)) to approximate a normal distribution, and their levels were adjusted for blood cell proportions. Linear mixed models were applied to test the association between age of pubertal markers and repeated measurement of DNAm at 10 and 18 years.

MAIN RESULTS AND THE ROLE OF CHANCE

In girls, a total of 63 019 CpGs statistically significantly changed after occurrence of any of the five pubertal events and 13 487 were changed subsequent to all five events: the respective number is boys were 3072 and 301. To further exclude false-positive findings, we investigated which CpGs were replicable in prior studies from Denmark or North America, resulting in 273 replicable CpG in girls and 67 CpGs in boys (236 and 68 genes, respectively). Most identified genes are known to be related to biological processes of puberty; however, genetic polymorphisms of only four of these genes were previously linked to pubertal markers in humans.

LIMITATIONS REASONS FOR CAUTION

The relative age of pubertal onset to the age of DNAm measurements does not allow causal inference, since DNAm at an earlier age may have affected the pubertal age or pubertal age may have altered later DNAm. This investigation concentrates on autosomes. CpGs on X and Y chromosomes are not included in the current study.

WIDER IMPLICATIONS OF THE FINDINGS

Assessment of biological processes involved in pubertal transitions should include epigenetic information. Differential DNAm related to puberty needs to be investigated to determine whether it can act as an early marker for adult diseases known to be associated with puberty.

STUDY FUNDING/COMPETING INTERESTS

This work was supported by NIH grants R03HD092776 (Epigenetic characterization of pubertal transitions) and R01AI121226. The 10-year follow-up of this study was funded by National Asthma Campaign, UK (Grant No 364), and the 18-year follow-up by a grant from the National Heart and Blood Institute (R01 HL082925). The authors have no conflicts to report.

A Polygenic Risk Score Suggests Shared Genetic Architecture of Voice Break With Early Markers of Pubertal Onset in Boys

CONTEXT

Voice break, as a landmark of advanced male puberty in genome-wide association studies (GWAS), has revealed that pubertal timing is a highly polygenic trait. Although voice break is easily recorded in large cohorts, it holds quite low precision as a marker of puberty. In contrast, gonadarche and pubarche are early and clinically well-defined measures of puberty onset.

OBJECTIVE

To determine whether a polygenic risk score (PRS) of alleles that confer risk for voice break associates with age at gonadarche (AAG) and age at pubarche (AAP) in Chilean boys.

EXPERIMENTAL DESIGN

Longitudinal study.

SUBJECTS AND METHODS

401 boys from the Growth and Obesity Chilean Cohort Study (n = 1194; 49.2% boys).

MAIN OUTCOME MEASURES

Biannual clinical pubertal staging including orchidometry. AAG and AAP were estimated by censoring methods. Genotyping was performed using the Multi-Ethnic Global Array (Illumina). Using GWAS summary statistics from the UK-Biobank, 29 significant and independent single nucleotide polymorphisms associated with age at voice break were extracted. Individual PRS were computed as the sum of risk alleles weighted by the effect size.

RESULTS

The PRS was associated with AAG (β=0.01, P = 0.04) and AAP (β=0.185, P = 0.0004). In addition, boys within the 20% highest PRS experienced gonadarche and pubarche 0.55 and 0.67 years later than those in the lowest 20%, respectively (P = 0.013 and P = 0.007).

CONCLUSIONS

Genetic variants identified in large GWAS on age at VB significantly associate with age at testicular growth and pubic hair development, suggesting that these events share a genetic architecture across ethnically distinct populations.

Age of puberty and Sleep duration: Observational and Mendelian randomization study

Earlier age of puberty has detrimental consequences for many aspects of health. Here, for the first time, we assessed the association of earlier puberty with sleep duration observationally and with validation using Mendelian Randomization. In the “Children of 1997” birth cohort (n = 8,327), we used adjusted multivariable logistic regression to assess the associations of each clinically assessed marker of earlier puberty with self-report sleep duration in adolescence. Using two-sample MR, we assessed the effect of earlier puberty timing based on 203 single nucleotide polymorphisms applied to genome wide association studies of sleep duration in adults (n = 335,410). In “Children of 1997”, cross-sectionally, older age of menarche was associated with longer (9+ hours) sleep duration [odds ratio (OR) 1.11, 95% confidence interval (CI) 1.01 to 1.21] at 13.5 years. The other earlier puberty markers were unrelated to sleep duration. Using inverse variance weighting, later of age at menarche increased adult sleep duration [0.020 per category, 95% CI 0.006 to 0.034]. This study demonstrated a causal effect of age at menarche on adult sleep duration, since age of menarche also affects obesity, our novel finding may be relevant to the observed relation of sleep duration with obesity and poor health.

HOPS: a quantitative score reveals pervasive horizontal pleiotropy in human genetic variation is driven by extreme polygenicity of human traits and diseases

Horizontal pleiotropy, where one variant has independent effects on multiple traits, is important for our understanding of the genetic architecture of human phenotypes. We develop a method to quantify horizontal pleiotropy using genome-wide association summary statistics and apply it to 372 heritable phenotypes measured in 361,194 UK Biobank individuals. Horizontal pleiotropy is pervasive throughout the human genome, prominent among highly polygenic phenotypes, and enriched in active regulatory regions. Our results highlight the central role horizontal pleiotropy plays in the genetic architecture of human phenotypes. The HOrizontal Pleiotropy Score (HOPS) method is available on Github at https://github.com/rondolab/HOPS.

Effects of the Timing of Sex-Steroid Exposure in Adolescence on Adult Health Outcomes

CONTEXT

Variation in pubertal timing is associated with a wide range of adult risks and outcomes, but it is unclear whether these associations are causal, and it is largely unknown whether these associations can be modified by treatment.

EVIDENCE ACQUISITION

We conducted PubMed searches to identify Mendelian randomization (MR) studies on the influence of pubertal timing on adult health and studies on sex-steroid treatment of the following conditions associated with reduced reproductive endocrine function in adolescence: constitutional delay, Turner syndrome, and Klinefelter syndrome.

EVIDENCE SYNTHESIS

Results of MR studies suggest that earlier pubertal timing increases body mass index; increases risk for breast, ovarian, endometrial, and prostate cancers; elevates fasting glucose levels and blood pressure; impairs lung capacity and increases risk for asthma; leads to earlier sexual intercourse and first birth; decreases time spent in education; and increases depressive symptoms in adolescence. Later pubertal timing appears to lower bone mineral density (BMD). Although studies of constitutional delay have not shown that sex-steroid treatment alters adult height or BMD, studies of girls with Turner syndrome and boys with Klinefelter syndrome suggest that earlier initiation of sex-steroid treatment improves physical and neurocognitive outcomes.

CONCLUSIONS

Despite having some limitations, MR studies suggest that pubertal timing causally influences many adult conditions and disease risks. Studies of Turner syndrome and Klinefelter syndrome suggest that earlier sex-steroid exposure may have short- and long-term benefits. The mechanisms underlying these findings and the effects of trends and treatments affecting pubertal timing remain to be determined.

Interdisciplinary Work Is Essential for Research on Puberty: Complexity and Dynamism in Action

Puberty is associated with changes in behavior and psychosocial well-being, and is important in lifelong health. We present five different facets regarding interdisciplinary research that are important to puberty. A short history of philosophical issues instrumental in promoting early interdisciplinary research is first presented. We discuss then what is hard and what is easy about interdisciplinary research, the purpose of which is to alert scientists to challenges and opportunities for interdisciplinary research on puberty. Readers then are introduced to advances and obstacles in interdisciplinary research on development. Recommendations for tailoring graduate education toward interdisciplinarity are introduced. Finally, issues related to publication, education of scientists, and policy makers are described. The report concludes with a discussion of funding and policy issues.

The Early Growth Genetics (EGG) and EArly Genetics and Lifecourse Epidemiology (EAGLE) consortia: design, results and future prospects

The impact of many unfavorable childhood traits or diseases, such as low birth weight and mental disorders, is not limited to childhood and adolescence, as they are also associated with poor outcomes in adulthood, such as cardiovascular disease. Insight into the genetic etiology of childhood and adolescent traits and disorders may therefore provide new perspectives, not only on how to improve wellbeing during childhood, but also how to prevent later adverse outcomes. To achieve the sample sizes required for genetic research, the Early Growth Genetics (EGG) and EArly Genetics and Lifecourse Epidemiology (EAGLE) consortia were established. The majority of the participating cohorts are longitudinal population-based samples, but other cohorts with data on early childhood phenotypes are also involved. Cohorts often have a broad focus and collect(ed) data on various somatic and psychiatric traits as well as environmental factors. Genetic variants have been successfully identified for multiple traits, for example, birth weight, atopic dermatitis, childhood BMI, allergic sensitization, and pubertal growth. Furthermore, the results have shown that genetic factors also partly underlie the association with adult traits. As sample sizes are still increasing, it is expected that future analyses will identify additional variants. This, in combination with the development of innovative statistical methods, will provide detailed insight on the mechanisms underlying the transition from childhood to adult disorders. Both consortia welcome new collaborations. Policies and contact details are available from the corresponding authors of this manuscript and/or the consortium websites.

Links between age at menarche, antral follicle count, and body mass index in African American and European American women

OBJECTIVE

To examine the relationships between age at menarche, antral follicle count (AFC), and body mass index (BMI) in a multi-ethnic population of women.

DESIGN

Community-based, cross-sectional study.

SETTING

Academic setting.

PATIENT(S)

A total of 245 African American women and 273 European American women, aged 25-45 years, with regular menstrual cycles and no reproductive disorders. The ethnicity of these women was self-reported and genetically validated.

INTERVENTION(S)

The AFCs were measured by transvaginal ultrasound during the early follicular phase. Anthropometric measurements were taken, and age at menarche was gathered by questionnaire.

MAIN OUTCOME MEASURE(S)

Determination of the associations between age of menarche and adult AFC and BMI.

RESULT(S)

Earlier age of menarche was associated with both higher BMIs and higher AFCs in adulthood, with control for female age. The antral follicle difference between early (<12 years) vs. late (≥15 years) initiation of menarche in both white and black women was +3.81 and +3.34 follicles, respectively, which is equivalent to an approximately 20% difference in AFC.

CONCLUSION(S)

This study provides the first evidence that timing of menarche may influence AFC. Because of limited studies on African American women, this work provides additional needed data and may enhance our ability to prospectively screen and better treat various diseases associated with the female reproductive lifespan.

Association of genetic polymorphisms with age at menarche in Russian women

OBJECTIVES

Examine the association of genetic polymorphisms with age at menarche (AAM) in Russian women.

STUDY DESIGN

A total of 1613 Russian females were recruited for the study. Fifty two polymorphisms were analyzed for their association with AAM, height, and BMI. The associations were analyzed assuming the additive, dominant, and recessive models and using the log-linear regression as implemented in PLINK v. 2.050. The 2-, 3-, and 4-loci models of gene-gene interactions were analyzed using the MB-MDR method and validated by the permutation test.

MAIN OUTCOME MEASURES

Genetic polymorphism rs6438424 3q13.32 was independently associated with AAM in Russian women. In addition, 14 SNPs were determined as possible contributors to this trait through gene-gene interactions.

RESULTS

The obtained results suggest that 14 out of 52 studied polymorphisms may contribute to AAM in Russian women. The rs6438424 3q13.32 polymorphism was associated with AAM according to both additive and dominant models (р_perm = 0.005). In total 12 two-, three-, and four-locus models of gene-gene interactions were determined as contributing to AAM (p_perm ≤ 0.006). Nine of the 14 AAM-associated SNPs are also associated with height and BMI (p_perm ≤ 0.003). Among 14 AAM-associated SNPs (a priori all having regulatory significance), the highest regulatory potential was determined for rs4633 COMT, rs2164808 POMC, rs2252673INSR, rs6438424 3q13.32, and rs10769908 STK33. Eleven loci are cis-eQTL and affect expression of 14 genes in various tissues and organs (FDR < 0.05). The neuropeptide-encoding genes were overrepresented among the AAM-associated genes (p_bonf = 0.039).

CONCLUSIONS

The rs6438424 polymorphism is independently associated with AAM in Russian females in this study. The other 14 SNPs manifest this association through gene-gene interactions.

Assessing the genetic correlations between early growth parameters and bone mineral density: A polygenic risk score analysis

OBJECTIVE

The relationships between early growth parameters and bone mineral density (BMD) remain elusive now. In this study, we performed a large scale polygenic risk score (PRS) analysis to evaluate the potential impact of early growth parameters on the variations of BMD.

METHODS

We used 2286 Caucasian subjects as cohort 1 and 3404 Framingham Heart Study (FHS) subjects as cohort 2 in this study. BMD at ulna & radius, hip and spine were measured using dual energy X-ray absorptiometry. BMD values were adjusted for age, sex, height and weight as covariates. Genome-wide single-nucleotide polymorphism (SNP) genotyping of the 2286 Caucasian subjects was performed using Affymetrix Human SNP Array 6.0. The GWAS datasets of early growth parameters were driven from the Early Growth Genetics Consortium, including birth weight (BW), birth head circumference (BHC), childhood body mass index (CBMI), pubertal height growth related indexes and tanner stage. Polygenic Risk Score (PRSice) and linkage disequilibrium (LD) score regression analysis were conducted to assess the genetic correlation between early growth parameters and BMD.

RESULTS

We detected significant genetic correlations in cohort 1, such as total spine BMD vs. CBMI (p value = 1.51 × 10^-4, rg = 0.4525), right ulna and radius BMD vs. CBMI (p value = 1.51 × 10^-4, rg = 0.4399) and total body BMD vs. tanner stage (p value = 7.00 × 10^-4, rg = -0.0721). For cohort 2, significant correlations were observed for total spine BMD vs. height change standard deviation score (SDS) between 8 years and adult (denoted as PGF + PGM) (p value = 3.97 × 10^-4, rg = -0.1425), femoral neck BMD vs. the timing of peak height velocity by looking at the height change SDS between age 14 years and adult (denoted as PTF + PTM) (p value = 7.04 × 10^-4, rg = -0.2185), and total spine BMD vs. PTF + PTM (p value = 6.86 × 10^-4, rg = -0.2180).

CONCLUSION

Our study results suggest that some early growth parameters could affect the variations of BMD.

Facial Genetics: A Brief Overview

Historically, craniofacial genetic research has understandably focused on identifying the causes of craniofacial anomalies and it has only been within the last 10 years, that there has been a drive to detail the biological basis of normal-range facial variation. This initiative has been facilitated by the availability of low-cost hi-resolution three-dimensional systems which have the ability to capture the facial details of thousands of individuals quickly and accurately. Simultaneous advances in genotyping technology have enabled the exploration of genetic influences on facial phenotypes, both in the present day and across human history. There are several important reasons for exploring the genetics of normal-range variation in facial morphology.     - Disentangling the environmental factors and relative parental biological contributions to heritable traits can help to answer the age-old question "why we look the way that we do?"     - Understanding the etiology of craniofacial anomalies; e.g., unaffected family members of individuals with non-syndromic cleft lip/palate (nsCL/P) have been shown to differ in terms of normal-range facial variation to the general population suggesting an etiological link between facial morphology and nsCL/P.     - Many factors such as ancestry, sex, eye/hair color as well as distinctive facial features (such as, shape of the chin, cheeks, eyes, forehead, lips, and nose) can be identified or estimated using an individual's genetic data, with potential applications in healthcare and forensics.     - Improved understanding of historical selection and adaptation relating to facial phenotypes, for example, skin pigmentation and geographical latitude.     - Highlighting what is known about shared facial traits, medical conditions and genes.

Genetics of pubertal timing

PURPOSE OF REVIEW

To summarize advances in the genetics underlying variation in normal pubertal timing, precocious puberty, and delayed puberty, and to discuss mechanisms by which genes may regulate pubertal timing.

RECENT FINDINGS

Genome-wide association studies have identified hundreds of loci that affect pubertal timing in the general population in both sexes and across ethnic groups. Single genes have been implicated in both precocious and delayed puberty. Potential mechanisms for how these genetic loci influence pubertal timing may include effects on the development and function of the GnRH neuronal network and the responsiveness of end-organs.

SUMMARY

There has been significant progress in identifying genetic loci that affect normal pubertal timing, and the first single-gene causes of precocious and delayed puberty are being described. How these genes influence pubertal timing remains to be determined.

The genetic underpinnings of variation in ages at menarche and natural menopause among women from the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) Study: A trans-ethnic meta-analysis

Current knowledge of the genetic architecture of key reproductive events across the female life course is largely based on association studies of European descent women. The relevance of known loci for age at menarche (AAM) and age at natural menopause (ANM) in diverse populations remains unclear. We investigated 32 AAM and 14 ANM previously-identified loci and sought to identify novel loci in a trans-ethnic array-wide study of 196,483 SNPs on the MetaboChip (Illumina, Inc.). A total of 45,364 women of diverse ancestries (African, Hispanic/Latina, Asian American and American Indian/Alaskan Native) in the Population Architecture using Genomics and Epidemiology (PAGE) Study were included in cross-sectional analyses of AAM and ANM. Within each study we conducted a linear regression of SNP associations with self-reported or medical record-derived AAM or ANM (in years), adjusting for birth year, population stratification, and center/region, as appropriate, and meta-analyzed results across studies using multiple meta-analytic techniques. For both AAM and ANM, we observed more directionally consistent associations with the previously reported risk alleles than expected by chance (p-values_binomial≤0.01). Eight densely genotyped reproductive loci generalized significantly to at least one non-European population. We identified one trans-ethnic array-wide SNP association with AAM and two significant associations with ANM, which have not been described previously. Additionally, we observed evidence of independent secondary signals at three of six AAM trans-ethnic loci. Our findings support the transferability of reproductive trait loci discovered in European women to women of other race/ethnicities and indicate the presence of additional trans-ethnic associations both at both novel and established loci. These findings suggest the benefit of including diverse populations in future studies of the genetic architecture of female growth and development.

Early pubertal timing is associated with lower sperm concentration in college students

To study the associations between pubertal timing and semen quality and reproductive hormones, 680 volunteers were recruited from universities in Chongqing, China. Pubertal timing was obtained using a questionnaire. The main measurements were five routine semen parameters and six reproductive hormones. After adjusting for potential confounders, we found that early pubertal timing was associated with lower sperm concentration. An one-year increase in age of peak height velocity was associated with a 4.7% (95% confidence interval [CI] = 1.0 to 8.6) increase in sperm concentration. An one-year increase in age of first spermatorrhea was associated with a 6.4% increase in sperm concentration and a 2.9% decrease in semen volume (95% CI = 1.7 to 11.3, -5.5 to -0.3; respectively). Regarding reproductive hormones, an one-year increase in age of height spurt and peak height velocity was associated with a 6.5% and a 6.7% decrease in estrogen (95% CI = -9.8 to -3.0, -10.4 to -2.8; respectively). While an one-year increase in age of height spurt was associated with higher follicle-stimulating hormone (% change = 2.6, 95% CI = 0.2 to 4.7). This was the first report that has suggested that early pubertal timing is associated with lower sperm concentration. However, further study is still needed to validate this association and fully elucidate the mechanism behind it.

Putative Effects of Obesity on Linear Growth and Puberty

Childhood obesity is a major public health problem that has grown to epidemic proportions throughout the world. Obesity is influenced by genetic and environmental factors. The nutritional status plays an important role in growth and body weight regulation. Excess adiposity during childhood can affect the process of growth and puberty. Obese children are frequently tall for their age, with accelerated epiphyseal growth plate maturation despite low growth hormone levels. Several regulatory hormones may affect the process of linear growth in the constellation of obesity, as high levels of insulin and leptin are observed in obese children. Leptin can act as a skeletal growth factor, with a direct effect on skeletal growth centers. The finding that overweight children, especially girls, tend to mature earlier than lean children has led to the hypothesis that the degree of body fatness may trigger the neuroendocrine events that lead to the onset of puberty. Leptin receptors have been identified in the hypothalamus, as well as in gonadotrope cells, ovarian follicular cells, and Leydig cells. The increased leptin and androgen levels seen in obese children may be implicated in their earlier onset of puberty and accelerated pubertal growth. This review is focused on the interaction between childhood obesity and growth and pubertal processes.
.

Muscle strength characteristics of the hamstrings and quadriceps in players from a high-level youth football (soccer) Academy

The objective of this study was to investigate knee muscle strength characteristics in players from a high-level youth football Academy. In total, 110 players (aged 8-15 years) underwent muscle strength assessments carried out by a research physiotherapist using a computer-linked hand-held dynamometer. Results indicated that isometric hamstrings and quadriceps strength increased with age, whereas the isometric hamstring to quadriceps (H/Q) ratio decreased with age. A number of youth football players (n = 20; 18%; 95% CI: 11-27%) demonstrated isometric H/Q ratios of less than 0.60, as well as muscle strength asymmetries between limbs for the hamstrings (n = 40, 36%; 95% CI: 27-46%) and quadriceps (n = 51, 46%; 95% CI 37-56%), potentially increasing injury risk. This study provides new evidence that the isometric H/Q ratio reduces with advancing age during adolescence which may have important implications for junior athlete development and long-term injury prevention in football.

Contributions of Function-Altering Variants in Genes Implicated in Pubertal Timing and Body Mass for Self-Limited Delayed Puberty

CONTEXT

Self-limited delayed puberty (DP) is often associated with a delay in physical maturation, but although highly heritable the causal genetic factors remain elusive. Genome-wide association studies of the timing of puberty have identified multiple loci for age at menarche in females and voice break in males, particularly in pathways controlling energy balance.

OBJECTIVE/MAIN OUTCOME MEASURES

We sought to assess the contribution of rare variants in such genes to the phenotype of familial DP.

DESIGN/PATIENTS

We performed whole-exome sequencing in 67 pedigrees (125 individuals with DP and 35 unaffected controls) from our unique cohort of familial self-limited DP. Using a whole-exome sequencing filtering pipeline one candidate gene [fat mass and obesity-associated gene (FTO)] was identified. In silico, in vitro, and mouse model studies were performed to investigate the pathogenicity of FTO variants and timing of puberty in FTO+/- mice.

RESULTS

We identified potentially pathogenic, rare variants in genes in linkage disequilibrium with genome-wide association studies of age at menarche loci in 283 genes. Of these, five genes were implicated in the control of body mass. After filtering for segregation with trait, one candidate, FTO, was retained. Two FTO variants, found in 14 affected individuals from three families, were also associated with leanness in these patients with DP. One variant (p.Leu44Val) demonstrated altered demethylation activity of the mutant protein in vitro. Fto+/- mice displayed a significantly delayed timing of pubertal onset (P < 0.05).

CONCLUSIONS

Mutations in genes implicated in body mass and timing of puberty in the general population may contribute to the pathogenesis of self-limited DP.

On the origin of obesity: identifying the biological, environmental and cultural drivers of genetic risk among human populations

Genetic predisposition to obesity presents a paradox: how do genetic variants with a detrimental impact on human health persist through evolutionary time? Numerous hypotheses, such as the thrifty genotype hypothesis, attempt to explain this phenomenon yet fail to provide a justification for the modern obesity epidemic. In this critical review, we appraise existing theories explaining the evolutionary origins of obesity and explore novel biological and sociocultural agents of evolutionary change to help explain the modern-day distribution of obesity-predisposing variants. Genetic drift, acting as a form of 'blind justice,' may randomly affect allele frequencies across generations while gene pleiotropy and adaptations to diverse environments may explain the rise and subsequent selection of obesity risk alleles. As an adaptive response, epigenetic regulation of gene expression may impact the manifestation of genetic predisposition to obesity. Finally, exposure to malnutrition and disease epidemics in the wake of oppressive social systems, culturally mediated notions of attractiveness and desirability, and diverse mating systems may play a role in shaping the human genome. As an important first step towards the identification of important drivers of obesity gene evolution, this review may inform empirical research focused on testing evolutionary theories by way of population genetics and mathematical modelling.

Gene expression profiling of puberty-associated genes reveals abundant tissue and sex-specific changes across postnatal development

The timing of human puberty is highly variable, sexually dimorphic, and associated with adverse health outcomes. Over 20 genes carrying rare mutations have been identified in known pubertal disorders, many of which encode critical components of the hypothalamic-pituitary-gonadal (HPG) axis. Recent genome-wide association studies (GWAS) have identified more than 100 candidate genes at loci associated with age at menarche or voice breaking in males. We know little about the spatial, temporal or postnatal expression patterns of the majority of these puberty-associated genes. Using a high-throughput and sensitive microfluidic quantitative PCR strategy, we profiled the gene expression patterns of the mouse orthologs of 178 puberty-associated genes in male and female mouse HPG axis tissues, the pineal gland, and the liver at five postnatal ages spanning the pubertal transition. The most dynamic gene expression changes were observed prior to puberty in all tissues. We detected known and novel tissue-enhanced gene expression patterns, with the hypothalamus expressing the largest number of the puberty-associated genes. Notably, over 40 puberty-associated genes in the pituitary gland showed sex-biased gene expression, most of which occurred peri-puberty. These sex-biased genes included the orthologs of candidate genes at GWAS loci that show sex-discordant effects on pubertal timing. Our findings provide new insight into the expression of puberty-associated genes and support the possibility that the pituitary plays a role in determining sex differences in the timing of puberty.

Neonatal Overnutrition Increases Testicular Size and Expression of Luteinizing Hormone β-Subunit in Peripubertal Male Rats

Proper nutrition is important for growth and development. Maturation of the reproductive axis and the timing of pubertal onset can be delayed when insufficient nutrition is available, or possibly advanced with nutritional abundance. The childhood obesity epidemic has been linked to a secular trend in advanced puberty in some populations. The increase in circulating leptin that occurs in association with obesity has been suggested to act as a signal that an adequate nutritional status exists for puberty to occur, allowing activation of central mechanisms. However, obesity-associated hyperleptinemia is linked to decreased leptin sensitivity, at least in adults. Here, we analyzed whether neonatal overnutrition modifies the response to an increase in leptin in peripubertal male rats, as previously demonstrated in females. Wistar rats were raised in litters of 4 (neonatal overnutrition) or 12 pups (controls) per dam. Leptin was administered sc (3 µg/g body weight) at postnatal day 35 and the rats killed 45 min or 2 h later. Postnatal overfeeding resulted in increased body weight and circulating leptin levels; however, we found no overweight-related changes in the mRNA levels of neuropeptides involved in metabolism or reproduction. In contrast, pituitary expression of luteinizing hormone (LH) beta-subunit was increased in overweight rats, as was testicular weight. There were no basal differences between L4 and L12 males or in their response to leptin administration in pSTAT3 levels in the hypothalamus at either 45 min or 2 h. In contrast, pJAK2 was found to be higher at 45 min in L4 compared to L12 males regardless of leptin treatment, while at 2 h it was higher in L4 leptin-treated males compared to L12 leptin-treated males, as well as L4 vehicle-treated rats. There were no changes in response to leptin administration in the expression of the neuropeptides analyzed. However, serum LH levels rose only in L4 males in response to leptin, but with no change in testosterone levels. In conclusion, the advancement in pubertal onset in males with neonatal overnutrition does not appear to be related to overt modifications in the central response to exogenous leptin during the peripubertal period.

The POZ/BTB and AT-Hook Containing Zinc Finger 1 (PATZ1) Transcription Regulator: Physiological Functions and Disease Involvement

PATZ1 is a zinc finger protein, belonging to the POZ domain Krüppel-like zinc finger (POK) family of architectural transcription factors, first discovered in 2000 by three independent groups. Since that time accumulating evidences have shown its involvement in a variety of biological processes (i.e., embryogenesis, stemness, apoptosis, senescence, proliferation, T-lymphocyte differentiation) and human diseases. Here we summarize these studies with a focus on the PATZ1 emerging and controversial role in cancer, where it acts as either a tumor suppressor or an oncogene. Finally, we give some insight on clinical perspectives using PATZ1 as a prognostic marker and therapeutic target.

Association between age at menarche and cardiovascular disease: A systematic review on risk and potential mechanisms

Age at menarche (AAM) has been reported to be associated with the risk of cardiovascular disease (CVD), but the shape of and the mechanisms behind this association remain unclear. We reviewed the data on the association between AAM and different subtypes of CVD, and used shared genetic loci to identify possible mechanisms underlying this association using shared genetic association. We searched the databases of PubMed, Web of Science and Embase through to April 2017. We included articles with any clinically manifest CVD endpoint and for any ethnicity. We identified single nucleotide polymorphisms (SNPs) for AAM in genome-wide association studies (GWAS) in Caucasians through PubMed and HuGE Navigator, and searched whether these SNPs or any of their proxies were associated with any CVD-related trait. Eight studies in Caucasian populations reported an inverse linear relation between AAM and CVD risk, whereas one large study reported a significant U-shaped relation between them. Data from Asian populations were contradictory and inconclusive. In total, 122 AAM SNPs were identified at a genome-wide significance level (p<5×10^-8). Of those, 18 were also associated with various CVD-related traits, primarily body mass index (BMI), obesity, and height. In conclusion, early AAM and possibly also late AAM increase the risk of CVD in Caucasian populations. Weight and height may be part of the mechanism underlying the relation between AAM and CVD risk in Caucasians. Data on other ethnicities are too limited for meaningful analysis and conclusions.

Pleiotropic Effects of Trait-Associated Genetic Variation on DNA Methylation: Utility for Refining GWAS Loci

Most genetic variants identified in genome-wide association studies (GWASs) of complex traits are thought to act by affecting gene regulation rather than directly altering the protein product. As a consequence, the actual genes involved in disease are not necessarily the most proximal to the associated variants. By integrating data from GWAS analyses with those from genetic studies of regulatory variation, it is possible to identify variants pleiotropically associated with both a complex trait and measures of gene regulation. In this study, we used summary-data-based Mendelian randomization (SMR), a method developed to identify variants pleiotropically associated with both complex traits and gene expression, to identify variants associated with complex traits and DNA methylation. We used large DNA methylation quantitative trait locus (mQTL) datasets generated from two different tissues (blood and fetal brain) to prioritize genes for >40 complex traits with robust GWAS data and found considerable overlap with the results of SMR analyses performed with expression QTL (eQTL) data. We identified multiple examples of variable DNA methylation associated with GWAS variants for a range of complex traits, demonstrating the utility of this approach for refining genetic association signals.