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

Association of birth weight, childhood obesity, and age at menarche with the risk of ovarian dysfunction: A mendelian randomization study

BACKGROUND

Observational studies have revealed associations between birth weight, childhood obesity, age at menarche, and ovarian dysfunction. However, these studies are susceptible to unavoidable confounding factors, leading to ongoing debates regarding their conclusions and making causal relationships challenging to infer. In light of these challenges, Mendelian randomization was employed in this study to investigate the causal relationships between birth weight, childhood obesity, age at menarche, and ovarian dysfunction.

METHODS

This study employed a two-sample Mendelian randomization approach using genetic variation as instrumental variables to investigate causal relationships. Genetic variation data were sourced from summary data of genome-wide association studies in European populations. Instrumental variables were selected based on the principles of Mendel's three assumptions. The study utilized the inverse variance weighted method to assess the relationships between birth weight, childhood obesity, age at menarche, and ovarian dysfunction. Supplementary analyses were conducted using MR-Egger regression, the weighted median method, and the weighted median mode to complement the IVW results. Furthermore, the study conducted heterogeneity, horizontal pleiotropy, and sensitivity analyses to evaluate the robustness of the results.

RESULTS

Based on the inverse variance weighted method, it was found that there exists a causal relationship between childhood obesity (OR = 1.378, 95% CI: 1.113∼1.705, p = 0.003), age at menarche (OR = 0.639, 95% CI: 0.468∼0.871, p = 0.005), and ovarian dysfunction, while no causal relationship was observed between birth weight and ovarian dysfunction. Heterogeneity tests, multiplicity tests, and leave-one-out sensitivity analyses did not detect any heterogeneity or multiplicity effects in the estimated impact of these three exposure factors on the risk of ovarian dysfunction.

CONCLUSIONS

This study represents the first evidence suggesting a potential causal relationship between childhood obesity, age at menarche, and ovarian dysfunction. Childhood obesity was found to increase the risk of ovarian dysfunction, while a later age at menarche was associated with a reduced risk of ovarian dysfunction.

Genetic Prediction of Osteoporosis by Anti-Müllerian Hormone Levels and Reproductive Factors in Women: A Mendelian Randomization Study

Previous observational studies have suggested that anti-Müllerian hormone (AMH) and reproductive factors are linked to reduced bone mineral density (BMD) and an increased risk of osteoporosis (OP) in women. However, related studies are limited, and these traditional observational studies may be subject to residual confounders and reverse causation, while also lacking a more comprehensive observation of various reproductive factors. Univariate and multivariate two-sample Mendelian randomization analyses were conducted to determine the causal associations of AMH levels and six reproductive factors with BMD and OP, using the random-effects inverse-variance weighted method. Heterogeneity was assessed using Cochran's Q-statistic, and sensitivity analyses were performed to identify causal correlations. Age at menarche (AAM) was negatively associated with total body BMD (TB-BMD) in females aged 45-60 and over 60 years, as well as with heel bone mineral density (eBMD). Conversely, age at natural menopause (ANM) was positively associated with TB-BMD in the same age ranges and with eBMD. ANM was only causally associated with self-reported OP and showed no significant correlation with definitively diagnosed OP. Neither AMH level nor other reproductive factors were significantly associated with a genetic predisposition to BMD at any age and OP. Later AAM and earlier ANM are significantly genetically causally associated with decreased BMD but not with OP. AMH levels, length of menstrual cycle, age at first birth, age at last birth, and number of live births, in terms of genetic backgrounds, are not causally related to BMD or OP.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Methylome analysis in girls with idiopathic central precocious puberty

BACKGROUND

Genetic and environmental factors are implicated in many developmental processes. Recent evidence, however, has suggested that epigenetic changes may also influence the onset of puberty or the susceptibility to a wide range of diseases later in life. The present study aims to investigate changes in genomic DNA methylation profiles associated with pubertal onset analyzing human peripheral blood leukocytes from three different groups of subjects: 19 girls with central precocious puberty (CPP), 14 healthy prepubertal girls matched by age and 13 healthy pubertal girls matched by pubertal stage. For this purpose, the comparisons were performed between pre- and pubertal controls to identify changes in normal pubertal transition and CPP versus pre- and pubertal controls.

RESULTS

Analysis of methylation changes associated with normal pubertal transition identified 1006 differentially methylated CpG sites, 86% of them were found to be hypermethylated in prepubertal controls. Some of these CpG sites reside in genes associated with the age of menarche or transcription factors involved in the process of pubertal development. Analysis of methylome profiles in CPP patients showed 65% and 55% hypomethylated CpG sites compared with prepubertal and pubertal controls, respectively. In addition, interestingly, our results revealed the presence of 43 differentially methylated genes coding for zinc finger (ZNF) proteins. Gene ontology and IPA analysis performed in the three groups studied revealed significant enrichment of them in some pathways related to neuronal communication (semaphorin and gustation pathways), estrogens action, some cancers (particularly breast and ovarian) or metabolism (particularly sirtuin).

CONCLUSIONS

The different methylation profiles of girls with normal and precocious puberty indicate that regulation of the pubertal process in humans is associated with specific epigenetic changes. Differentially methylated genes include ZNF genes that may play a role in developmental control. In addition, our data highlight changes in the methylation status of genes involved in signaling pathways that determine the migration and function of GnRH neurons and the onset of metabolic and neoplastic diseases that may be associated with CPP in later life.

Puberty classifications in beef heifers are moderately to highly heritable and associated with candidate genes related to cyclicity and timing of puberty

Introduction: Pubertal attainment is critical to reproductive longevity in heifers. Previously, four heifer pubertal classifications were identified according to attainment of blood plasma progesterone concentrations > 1 ng/ml: 1) Early; 2) Typical; 3) Start-Stop; and 4) Non-Cycling. Early and Typical heifers initiated and maintained cyclicity, Start-Stop started and then stopped cyclicity and Non-Cycling never initiated cyclicity. Start-Stop heifers segregated into Start-Stop-Discontinuous (SSD) or Start-Stop-Start (SSS), with SSD having similar phenotypes to Non-Cycling and SSS to Typical heifers. We hypothesized that these pubertal classifications are heritable, and loci associated with pubertal classifications could be identified by genome wide association studies (GWAS). Methods: Heifers (n = 532; 2017 - 2022) genotyped on the Illumina Bovine SNP50 v2 or GGP Bovine 100K SNP panels were used for variant component estimation and GWAS. Heritability was estimated using a univariate Bayesian animal model. Results: When considering pubertal classifications: Early, Typical, SSS, SSD, and Non-Cycling, pubertal class was moderately heritable (0.38 ± 0.08). However, when heifers who initiated and maintained cyclicity were compared to those that did not cycle (Early+Typical vs. SSD+Non-Cycling) heritability was greater (0.59 ± 0.19). A GWAS did not identify single nucleotide polymorphisms (SNPs) significantly associated with pubertal classifications, indicating puberty is a polygenic trait. A candidate gene approach was used, which fitted SNPs within or nearby a set of 71 candidate genes previously associated with puberty, PCOS, cyclicity, regulation of hormone secretion, signal transduction, and methylation. Eight genes/regions were associated with pubertal classifications, and twenty-two genes/regions were associated with whether puberty was attained during the trial. Additionally, whole genome sequencing (WGS) data on 33 heifers were aligned to the reference genome (ARS-UCD1.2) to identify variants in FSHR, a gene critical to pubertal attainment. Fisher's exact test determined if FSHR SNPs segregated by pubertal classification. Two FSHR SNPs that were not on the bovine SNP panel were selected for additional genotyping and analysis, and one was associated with pubertal classifications and whether they cycled during the trial. Discussion: In summary, these pubertal classifications are moderately to highly heritable and polygenic. Consequently, genomic tools to inform selection/management of replacement heifers would be useful if informed by SNPs associated with cyclicity and early pubertal attainment.

Integrating Genetic and Transcriptomic Data to Identify Genes Underlying Obesity Risk Loci

Genome-wide association studies (GWAS) have identified numerous body mass index (BMI) loci. However, most underlying mechanisms from risk locus to BMI remain unknown. Leveraging omics data through integrative analyses could provide more comprehensive views of biological pathways on BMI. We analyzed genotype and blood gene expression data in up to 5,619 samples from the Framingham Heart Study (FHS). Using 3,992 single nucleotide polymorphisms (SNPs) at 97 BMI loci and 20,692 transcripts within 1 Mb, we performed separate association analyses of transcript with BMI and SNP with transcript (PBMI and PSNP, respectively) and then a correlated meta-analysis between the full summary data sets (PMETA). We identified transcripts that met Bonferroni-corrected significance for each omic, were more significant in the correlated meta-analysis than each omic, and were at least nominally associated with BMI in FHS data. Among 308 significant SNP-transcript-BMI associations, we identified seven genes (NT5C2, GSTM3, SNAPC3, SPNS1, TMEM245, YPEL3, and ZNF646) in five association regions. Using an independent sample of blood gene expression data, we validated results for SNAPC3 and YPEL3. We tested for generalization of these associations in hypothalamus, nucleus accumbens, and liver and observed significant (PMETA<0.05 & PMETA Collapse

Key Words Collapse

MESH Headings Collapse

Grants

• R01 HL143885 NHLBI NIH HHS
• R01 HL163262 NHLBI NIH HHS
• R01 HL142302 NHLBI NIH HHS
• R01 HL151152 NHLBI NIH HHS
• N01HC25195 NHLBI NIH HHS
• R01 HG010297 NHGRI NIH HHS
• R01 HL092577 NHLBI NIH HHS
• P20 MD000170 NIMHD NIH HHS
• 75N92019D00031 NHLBI NIH HHS
• R01 HD057194 NICHD NIH HHS
• UL1 TR000371 NCATS NIH HHS
• R01 DK122503 NIDDK NIH HHS

Collapse

Affiliation(s)

Hanfei Xu Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA
Shreyash Gupta Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Ian Dinsmore Department of Genomic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Abbey Kollu Department of Psychology and Neuroscience, University of North Carolina, 235 E. Cameron Avenue, Chapel Hill, NC, 27599, USA
Anne Marie Cawley Marsico Lung Institute, University of North Carolina, 125 Mason Farm Rd, Chapel Hill, NC, 27599, USA
Mohammad Y. Anwar Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Hung-Hsin Chen Institute of Biomedical Sciences, Academia Sinica, No. 128, Section 2, Academia Rd., Taipei, Nangang District, 115201, Taiwan Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Lauren E. Petty Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Sudha Seshadri Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA Department of Neurology, School of Medicine, Boston University, 85 East Concord Street, Boston, MA, 02118, USA Glenn Biggs Institute for Alzheimer’s & Neurodegenerative Diseases, UT Health San Antonio, 8300 Floyd Curl Drive, San Antonio, TX, 78229, USA
Misa Graff Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Piper Below Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Jennifer A. Brody Department of Medicine, Cardiovascular Health Research Unit, University of Washington, 1730 Minor Ave, Seattle, WA, 98101, USA
Geetha Chittoor Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Susan P. Fisher-Hoch Department of Epidemiology, School of Public Health, UT Health Houston, Regional Academic Health Center, One West University Blvd, Brownsville, TX, 78520, USA
Nancy L. Heard-Costa Framingham Heart Study, 73 Mt Wayte Ave, Framingham, MA, 01702, USA Department of Neurology, Chobanian & Avedisian School of Medicine, Boston University, 72 E Concord St, Boston, MA, 02118, USA
Daniel Levy Population Sciences Branch, National Heart, Lung, and Blood Institute of the National Institutes of Health, 6701 Rockledge Drive, Bethesda, MD, 20892, USA
Honghuang Lin Department of Medicine, University of Massachusetts Chan Medical School, 55 N Lake Ave, Worcester, MA, 01655, USA
Ruth JF. Loos Charles Bronfman Institute for Personalized Medicine at Mount Sinai, Icahn School of Medicine at Mount Sinai, 1 Gustave L. Levy Pl, New York, NY, 10029, USA Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3A, 2200, Copenhagen, Denmark
Joseph B. Mccormick Department of Epidemiology, School of Public Health, UT Health Houston, Regional Academic Health Center, One West University Blvd, Brownsville, TX, 78520, USA
Jerome I. Rotter Department of Pediatrics, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, 1124 West Carson Street, Torrance, CA, 90502, USA
Tooraj Mirshahi Department of Genomic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Christopher D. Still Center for Obesity and Metabolic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Anita Destefano Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA Department of Neurology, School of Medicine, Boston University, 85 East Concord Street, Boston, MA, 02118, USA
L. Adrienne Cupples Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA
Karen L Mohlke Department of Genetics, School of Medicine, University of North Carolina, Chapel Hill, NC, 27599, USA
Kari E. North Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Anne E. Justice Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Ching-Ti Liu Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA

Collapse

Genetic parameters, GWAS and selection perspective on gestation length in 16 French cattle breeds

In this paper, we present a comprehensive study of gestation length (GL) in 16 cattle breeds by using large genotype and animal record databases. Data included over 20 million gestations since 2000 and genotypes from one million calves. The study addressed the GL variability within and between breeds, estimation of its direct and maternal heritability coefficients, association with fitness and several economic traits, and QTL detection. The breed average GL varied from 279.7 to 294.4 d, in Holstein and Blonde d'Aquitaine breeds, respectively. Standard deviations per breed were similar and ranged from 5.2 to 5.8 d. Direct heritability (i.e., for GL defined as a trait of the calf) was moderate to high (h2 = 0.40 to 0.67), whereas the maternal heritability was low (0.04 to 0.06). Extreme breeding values for GL were strongly associated with a higher mortality during the first 2 d of life and were associated with milk production of dams for dairy breeds and precocity of females. Finally, several QTL were detected affecting GL with cumulated effects up to a few days, and at least 2 QTL were found to be shared between different breeds. Our study highlights the risks that would be associated with selection toward a reduced gestation length. Further genomic studies are needed to identify the causal variants, and their association with juvenile mortality and other economic traits.

Menarche-a journey into womanhood: age at menarche and health-related outcomes in East Asians

STUDY QUESTION

Are there associations of age at menarche (AAM) with health-related outcomes in East Asians?

SUMMARY ANSWER

AAM is associated with osteoporosis, Type 2 diabetes (T2D), glaucoma, and uterine fibroids, as demonstrated through observational studies, polygenic risk scores, genetic correlations, and Mendelian randomization (MR), with additional findings indicating a causal effect of BMI and T2D on earlier AAM.

WHAT IS KNOWN ALREADY

Puberty timing is linked to adult disease risk, but research predominantly focuses on European populations, with limited studies in other groups.

STUDY DESIGN, SIZE, DURATION

We performed an AAM genome-wide association study (GWAS) with 57 890 Han Taiwanese females and examined the association between AAM and 154 disease outcomes using the Taiwanese database. Additionally, we examined genetic correlations between AAM and 113 diseases and 67 phenotypes using Japanese GWAS summary statistics.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We performed AAM GWAS and gene-based GWAS studies to obtain summary statistics and identify potential AAM-related genes. We applied phenotype, polygenic risk scores, and genetic correlation analyses of AAM to explore health-related outcomes, using multivariate regression and linkage disequilibrium score regression analyses. We also explored potential bidirectional causal relationships between AAM and related outcomes through univariable and multivariable MR analyses.

MAIN RESULTS AND THE ROLE OF CHANCE

Fifteen lead single-nucleotide polymorphisms and 24 distinct genes were associated with AAM in Taiwan. AAM was genetically associated with later menarche and menopause, greater height, increased osteoporosis risk, but lower BMI, and reduced risks of T2D, glaucoma, and uterine fibroids in East Asians. Bidirectional MR analyses indicated that higher BMI/T2D causally leads to earlier AAM.

LIMITATIONS, REASONS FOR CAUTION

Our findings were specific to Han Taiwanese individuals, with genetic correlation analyses conducted in East Asians. Further research in other ethnic groups is necessary.

WIDER IMPLICATIONS OF THE FINDINGS

Our study provides insights into the genetic architecture of AAM and its health-related outcomes in East Asians, highlighting causal links between BMI/T2D and earlier AAM, which may suggest potential prevention strategies for early puberty.

STUDY FUNDING/COMPETING INTEREST(S)

The work was supported by China Medical University, Taiwan (CMU110-S-17, CMU110-S-24, CMU110-MF-49, CMU111-SR-158, CMU111-MF-105, CMU111-MF-21, CMU111-S-35, CMU112-SR-30, and CMU112-MF-101), the China Medical University Hospital, Taiwan (DMR-111-062, DMR-111-153, DMR-112-042, DMR-113-038, and DMR-113-103), and the Ministry of Science and Technology, Taiwan (MOST 111-2314-B-039-063-MY3, MOST 111-2314-B-039-064-MY3, MOST 111-2410-H-039-002-MY3, and NSTC 112-2813-C-039-036-B). The funders had no influence on the data collection, analyses, or conclusions of the study. No conflict of interests to declare.

TRIAL REGISTRATION NUMBER

N/A.

Impact of insomnia on ovarian cancer risk and survival: a Mendelian randomization study

BACKGROUND

Insomnia is the most common sleep disorder in patients with epithelial ovarian cancer (EOC). We investigated the causal association between genetically predicted insomnia and EOC risk and survival through a two-sample Mendelian randomization (MR) study.

METHODS

Insomnia was proxied using genetic variants identified in a genome-wide association study (GWAS) meta-analysis of UK Biobank and 23andMe. Using genetic associations with EOC risk and overall survival from the Ovarian Cancer Association Consortium (OCAC) GWAS in 66,450 women (over 11,000 cases with clinical follow-up), we performed Iterative Mendelian Randomization and Pleiotropy (IMRP) analysis followed by a set of sensitivity analyses. Genetic associations with survival and response to treatment in ovarian cancer study of The Cancer Genome Atlas (TCGA) were estimated controlling for chemotherapy and clinical factors.

FINDINGS

Insomnia was associated with higher risk of endometrioid EOC (OR = 1.60, 95% CI 1.05-2.45) and lower risk of high-grade serous EOC (HGSOC) and clear cell EOC (OR = 0.79 and 0.48, 95% CI 0.63-1.00 and 0.27-0.86, respectively). In survival analysis, insomnia was associated with shorter survival of invasive EOC (OR = 1.45, 95% CI 1.13-1.87) and HGSOC (OR = 1.4, 95% CI 1.04-1.89), which was attenuated after adjustment for body mass index and reproductive age. Insomnia was associated with reduced survival in TCGA HGSOC cases who received standard chemotherapy (OR = 2.48, 95% CI 1.13-5.42), but was attenuated after adjustment for clinical factors.

INTERPRETATION

This study supports the impact of insomnia on EOC risk and survival, suggesting treatments targeting insomnia could be pivotal for prevention and improving patient survival.

FUNDING

National Institutes of Health, National Cancer Institute. Full funding details are provided in acknowledgments.

A comprehensive genome-wide cross-trait analysis of sexual factors and uterine leiomyoma

Age at first sexual intercourse (AFS) and lifetime number of sexual partners (NSP) may influence the pathogenesis of uterine leiomyoma (UL) through their associations with hormonal concentrations and uterine infections. Leveraging summary statistics from large-scale genome-wide association studies conducted in European ancestry for each trait (NAFS = 214,547; NNSP = 370,711; NUL = 302,979), we observed a significant negative genomic correlation for UL with AFS (rg = -0.11, P = 7.83×10-4), but not with NSP (rg = 0.01, P = 0.62). Four specific genomic regions were identified as contributing significant local genetic correlations to AFS and UL, including one genomic region further identified for NSP and UL. Partitioning SNP-heritability with cell-type-specific annotations, a close clustering of UL with both AFS and NSP was identified in immune and blood-related components. Cross-trait meta-analysis revealed 15 loci shared between AFS/NSP and UL, including 7 novel SNPs. Univariable two-sample Mendelian randomization (MR) analysis suggested no evidence for a causal association between genetically predicted AFS/NSP and risk of UL, nor vice versa. Multivariable MR adjusting for age at menarche or/and age at natural menopause revealed a significant causal effect of genetically predicted higher AFS on a lower risk of UL. Such effect attenuated to null when age at first birth was further included. Utilizing participant-level data from the UK Biobank, one-sample MR based on genetic risk scores yielded consistent null findings among both pre-menopausal and post-menopausal females. From a genetic perspective, our study demonstrates an intrinsic link underlying sexual factors (AFS and NSP) and UL, highlighting shared biological mechanisms rather than direct causal effects. Future studies are needed to elucidate the specific mechanisms involved in the shared genetic influences and their potential impact on UL development.

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.

Genetic correlation and Mendelian randomization analyses support causal relationships between dietary habits and age at menarche

Dietary habits are essential in the mean age at menarche (AAM). However, the causal relationship between these factors remains unclear. Therefore, this study aimed to elucidate the genetic relationship between dietary habits and AAM. Genetic summary statistics for dietary habits were obtained from the UK Biobank. GWAS summary data for AAM was obtained from the ReproGen Consortium. Linkage disequilibrium score regression was used to test genetic correlations between dietary habits and AAM. The Mendelian randomization (MR) analyses used the inverse-variance weighted method. Genetic correlations with AAM were identified for 29 candi-date dietary habits, such as milk type (skimmed, semi-skimmed, full cream; coefficient = 0.2704, Pldsc = 1.13 × 10-14). MR evaluations revealed that 19 dietary habits were associated with AAM, including bread type (white vs. any other; OR 1.71, 95% CI 1.28-2.29, Pmr = 3.20 × 10-4), tablespoons of cooked vegetables (OR 0.437, 95% CI 0.29-0.67; Pmr = 1.30 × 10-4), and cups of coffee per day (OR 0.72, 95% CI 0.57-0.92, Pmr = 8.31 × 10-3). These results were observed to be stable under the sensitivity analysis. Our study provides potential insights into the genetic mechanisms underlying AAM and evidence that dietary habits are associated with AAM.

Hormone and reproductive factors and risk of systemic lupus erythematosus: a Mendelian randomized study

Systemic lupus erythematosus (SLE) is an autoimmune and inflammatory disease with a risk associated with hormonal and reproductive factors. However, the potential causal effects between these factors and SLE remain unclear. A two-sample Mendelian randomization study was conducted using the published summary data from the genome-wide association study database. Five independent genetic variants associated with hormonal and reproductive factors were selected as instrumental variables: age at menarche, age at natural menopause, estradiol, testosterone, and follistatin. To estimate the causal relationship between these exposure factors and disease outcome, we employed the inverse-variance weighted, weighted median, and MR-Egger methods. In addition, we carried out multiple sensitivity analyses to validate model assumptions. Inverse variance weighted showed that there was a causal association between circulating follistatin and SLE risk (OR = 1.38, 95% CI 1.03 to 1.86, P = 0.033). However, no evidence was found that correlation between AAM (OR = 1.04, 95% CI 0.77 to 1.40, P = 0.798), ANM (OR = 0.99, 95% CI 0.92 to 1.06, P = 0.721), E2 (OR = 1.40, 95% CI 0.14 to 13.56, P = 0.772), T (OR = 1.25, 95% CI 0.70 to 2.28, P = 0.459), and SLE risk. Our study revealed that elevated circulating follistatin associates with an increased risk of SLE. This finding suggests that the regulatory signals mediated by circulating follistatin may provide a potential mechanism relevant to the treatment of SLE.

Variable parallelism in the genomic basis of age at maturity across spatial scales in Atlantic Salmon

Complex traits often exhibit complex underlying genetic architectures resulting from a combination of evolution from standing variation, hard and soft sweeps, and alleles of varying effect size. Increasingly, studies implicate both large-effect loci and polygenic patterns underpinning adaptation, but the extent that common genetic architectures are utilized during repeated adaptation is not well understood. Sea age or age at maturation represents a significant life history trait in Atlantic Salmon (Salmo salar), the genetic basis of which has been studied extensively in European Atlantic populations, with repeated identification of large-effect loci. However, the genetic basis of sea age within North American Atlantic Salmon populations remains unclear, as does the potential for a parallel trans-Atlantic genomic basis to sea age. Here, we used a large single-nucleotide polymorphism (SNP) array and low-coverage whole-genome resequencing to explore the genomic basis of sea age variation in North American Atlantic Salmon. We found significant associations at the gene and SNP level with a large-effect locus (vgll3) previously identified in European populations, indicating genetic parallelism, but found that this pattern varied based on both sex and geographic region. We also identified nonrepeated sets of highly predictive loci associated with sea age among populations and sexes within North America, indicating polygenicity and low rates of genomic parallelism. Despite low genome-wide parallelism, we uncovered a set of conserved molecular pathways associated with sea age that were consistently enriched among comparisons, including calcium signaling, MapK signaling, focal adhesion, and phosphatidylinositol signaling. Together, our results indicate parallelism of the molecular basis of sea age in North American Atlantic Salmon across large-effect genes and molecular pathways despite population-specific patterns of polygenicity. These findings reveal roles for both contingency and repeated adaptation at the molecular level in the evolution of life history variation.

Sex- and time-specific associations of obesity with glycaemic traits: A two-step multivariate Mendelian randomization study

AIM

To assess the sex- and time-specific causal effects of obesity-related anthropometric traits on glycaemic traits.

MATERIALS AND METHODS

We used univariate and multivariate Mendelian randomization to assess the causal associations of anthropometric traits (gestational variables, birth weight, childhood body mass index [BMI], BMI, waist-to-hip ratio [WHR], BMI-adjusted WHR [WHRadj BMI]) with fasting glucose and insulin in Europeans from the Early Growth Genetics Consortium (n ≤ 298 142), the UK Biobank, the Genetic Investigation of Anthropometric Traits Consortium (n ≤ 697 734; females: n ≤ 434 794; males: n ≤ 374 754) and the Meta-Analyses of Glucose and Insulin-related traits Consortium (n ≤ 151 188; females: n ≤ 73 089; males: n ≤ 67 506), adjusting for maternal genetic effects, smoking, alcohol consumption, and age at menarche.

RESULTS

We observed a null association for gestational variables, a negative association for birth weight, and positive associations for childhood BMI and adult traits (BMI, WHR, and WHRadj BMI). In female participants, increased birth weight causally decreased fasting insulin (betaIVW , -0.07, 95% confidence interval [CI] -0.11 to -0.03; p = 1.92 × 10-3 ), but not glucose levels, which was annulled by adjusting for age at menarche. In male participants, increased birth weight causally decreased fasting glucose (betainverse-variance-weighted (IVW) , -0.07, 95% CI -0.11 to -0.03; p = 3.22 × 10-4 ), but not insulin levels. In time-specific analyses, independent effects of birth weight were absent in female participants, and were more pronounced in male participants. Independent effects of childhood BMI were attenuated in both sexes; independent effects of adult traits differed by sex.

CONCLUSIONS

Our findings provide evidence for causal and independent effects of sex- and time-specific anthropometric traits on glycaemic variables, and highlight the importance of considering multiple obesity exposures at different time points in the life course.

Altered pubertal timing in 7q11.23 copy number variations and associated genetic mechanisms

Pubertal timing, including age at menarche (AAM), is a heritable trait linked to lifetime health outcomes. Here, we investigate genetic mechanisms underlying AAM by combining genome-wide association study (GWAS) data with investigations of two rare genetic conditions clinically associated with altered AAM: Williams syndrome (WS), a 7q11.23 hemideletion characterized by early puberty; and duplication of the same genes (7q11.23 Duplication syndrome [Dup7]) characterized by delayed puberty. First, we confirm that AAM-derived polygenic scores in typically developing children (TD) explain a modest amount of variance in AAM (R2 = 0.09; p = 0.04). Next, we demonstrate that 7q11.23 copy number impacts AAM (WS < TD < Dup7; p = 1.2x10-8, η2 = 0.45) and pituitary volume (WS < TD < Dup7; p = 3x10-5, ηp2 = 0.2) with greater effect sizes. Finally, we relate an AAM-GWAS signal in 7q11.23 to altered expression in postmortem brains of STAG3L2 (p = 1.7x10-17), a gene we also find differentially expressed with 7q11.23 copy number (p = 0.03). Collectively, these data explicate the role of 7q11.23 in pubertal onset, with STAG3L2 and pituitary development as potential mediators.

Gene co-expression patterns in Atlantic salmon adipose tissue provide a molecular link among seasonal changes, energy balance and age at maturity

Sexual maturation in many fishes requires a major physiological change that involves a rapid transition between energy storage and usage. In Atlantic salmon, this transition for the initiation of maturation is tightly controlled by seasonality and requires a high-energy status. Lipid metabolism is at the heart of this transition since lipids are the main energy storing molecules. The balance between lipogenesis (lipid accumulation) and lipolysis (lipid use) determines energy status transitions. A genomic region containing a transcription co-factor of the Hippo pathway, vgll3, is the main determinant of maturation timing in Atlantic salmon. Interestingly, vgll3 acts as an inhibitor of adipogenesis in mice and its genotypes are potentially associated with seasonal heterochrony in lipid storage and usage in juvenile Atlantic salmon. Here, we explored changes in expression of more than 300 genes directly involved in the processes of adipogenesis, lipogenesis and lipolysis, as well as the Hippo pathway in the adipose tissue of immature and mature Atlantic salmon with distinct vgll3 genotypes. We found molecular evidence consistent with a scenario in which immature males with different vgll3 genotypes exhibit contrasting seasonal dynamics in their lipid profiles. We also identified components of the Hippo signalling pathway as potential major drivers of vgll3 genotype-specific differences in adipose tissue gene expression. This study demonstrates the importance of adipose gene expression patterns for directly linking environmental changes with energy balance and age at maturity through genetic factors bridging lipid metabolism, seasonality and sexual maturation.

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.

Identifying causal associations between women's reproductive traits and risk of schizophrenia: a multivariate validated two-sample Mendelian randomization analysis

BACKGROUND

A significant association between women's reproductive traits and the risk of schizophrenia (SCZ) has been discovered, but the causalities remain unclear. We designed a two-sample univariate Mendelian randomization (MR) study using female-specific SNPs collected from a large-scale genome-wide association study as a genetic tool to explore the causal effect of female reproductive traits on the risk of SCZ, and conducted a multivariate MR study to re-validate the above findings.

METHODS

From extensive genome-wide association studies (GWAS) of people with European ancestry (n = 176,881 to 418,758 individuals), summary-level data on five female reproductive variables were extracted. Summary-level information on SCZ was taken from a GWAS meta-analysis involving 320,404 people with European ancestry. The inverse variance weighting estimations for both univariable MR (UVMR) and multivariable MR (MVMR) were presented as the primary results. MR-Egger, weighted median, simple mode, and weighted mode regression methods for UVMR, and MVMR-Egger, MVMR-Lasso, and MVMR-median methods for MVMR were used for sensitivity analyses.

RESULTS

The UVMR produced compelling proof for a connection between genetically predicted later age at first sexual intercourse (AFS) (OR, 0.632; 95% CI, 0.512-0.777; P < 0.01) and decreased SCZ risk. Pleiotropy analysis of the AFS-SCZ association confirmed the robustness of the MR results (P > 0.05). Consistent, substantial causal effects of AFS (OR, 0.592; 95%CI, 0.407-0.862; P < 0.01) on the risk of SCZ were demonstrated after adjusting for body mass index, years of schooling, and smoking initiation using MVMR.

CONCLUSIONS

Our findings provide convincing evidence that early AFS is a risk factor for SCZ. SCZ risk may be decreased by raising awareness of reproductive healthcare for women.

The causal effects of age at menarche and age at menopause on sepsis: A two-sample Mendelian randomization analysis

OBJECTIVES

To determine whether the age at menarche (AAM) and the age at menopause (ANM) are causally related to the development of sepsis.

METHODS

We performed a two-sample Mendelian randomization (MR) analysis by utilizing summary statistics from genome-wide association study (GWAS) datasets for both the exposure and outcome variables. Single nucleotide polymorphisms (SNPs) that exhibited significant associations with AAM and ANM were chosen as instrumental variables to estimate the causal effects on sepsis. Our study employed a variety of methods, including MR-Egger regression, weighted median estimation, inverse variance weighting, a simple model, and a weighted model. Odds ratios (ORs) along with their corresponding 95% confidence intervals (CIs) were used as the primary indicators for assessing causality. Furthermore, we conducted sensitivity analyses to explore the presence of genetic heterogeneity and validate the robustness of the tools employed.

RESULT

Our analysis revealed a significant negative causal relationship between AAM and the risk of sepsis (IVW: OR = 0.870, 95% CI = 0.793-0.955, P = 0.003). However, our Mendelian randomization (MR) analysis did not yield sufficient evidence to support a causal link between ANM and sepsis (IVW: OR = 0.987, 95% CI = 0.971-1.004, P = 0.129).

CONCLUSIONS

Our findings suggest that an earlier AAM may be associated with an increased risk of sepsis. However, we did not find sufficient evidence to support a causal relationship between ANM and sepsis.

The impact of reproductive factors on the metabolic profile of females from menarche to menopause

We explore the relation between age at menarche, parity and age at natural menopause with 249 metabolic traits in over 65,000 UK Biobank women using multivariable regression, Mendelian randomization and negative control (parity only). Older age of menarche is related to a less atherogenic metabolic profile in multivariable regression and Mendelian randomization, which is largely attenuated when accounting for adult body mass index. In multivariable regression, higher parity relates to more particles and lipids in VLDL, which are not observed in male negative controls. In multivariable regression and Mendelian randomization, older age at natural menopause is related to lower concentrations of inflammation markers, but we observe inconsistent results for LDL-related traits due to chronological age-specific effects. For example, older age at menopause is related to lower LDL-cholesterol in younger women but slightly higher in older women. Our findings support a role of reproductive traits on later life metabolic profile and provide insights into identifying novel markers for the prevention of adverse cardiometabolic outcomes in women.

Causal relationship between female reproductive factors, sex hormones and uterine leiomyoma: a Mendelian randomization study

RESEARCH QUESTION

Are the observed associations between female reproductive factors and sex hormones with the risk of uterine leiomyoma truly causal associations?

DESIGN

The putative causal relationships between female reproductive factors and sex hormones with uterine leiomyoma were investigated using two-sample Mendelian randomization. Statistics on exposure-associated genetic variants were obtained from genome-wide association studies (GWAS). The uterine leiomyoma GWAS from the FinnGen and FibroGENE consortia were used as outcome data for discovery and replication analyses, respectively. Results were pooled by meta-analysis. Sensitivity analyses ensured robustness of the Mendelian randomization analysis.

RESULTS

When FinnGen GWAS were used as outcome data, a causal relationship was found between age at menarche (OR 0.84, P < 0.0001), age at menopause (OR 1.08, P < 0.0001), number of live births (OR 0.25, P < 0.001) and total testosterone levels (OR 0.90, P < 0.001) with the risk of uterine leiomyoma. When FibroGENE GWAS were used as outcome data, Mendelian randomization results for age at menopause, the number of live births and total testosterone levels were replicated. In the meta-analysis, a later age at menopause (OR 1.08, P < 0.0001) was associated with an increased risk of uterine leiomyoma. A higher number of live births (OR 0.25, P < 0.0001) and higher total testosterone levels (OR 0.90, P < 0.0001) were associated with a decreased risk of uterine leiomyoma.

CONCLUSIONS

A causal relationship between later age at menopause, lower number of live births and lower total testosterone levels with increased risk of uterine leiomyoma was found.

Cloning, tissue expression and imprinting status analysis of the NDN gene in Dolang sheep

BACKGROUND

Genomic imprinting refers to expressing parent-specific genes in mammalian diploid cells. The NDN gene is maternally imprinted in humans and mice and correlates with the timing of puberty. This study aimed to investigate its imprinting status and its relationship with the onset of puberty in Dolang sheep.

METHODS AND RESULTS

In this study, cloning and sequencing obtained the NDN gene cDNA sequence of 1082 bp of Dolang sheep, coding for 325 amino acids. Similarity analysis and phylogenetic tree showed that the NDN gene conformed to the law of speciation and was highly conserved among mammals. RT-qPCR results showed the highest expression of NDN mRNA was found in the hypothalamus at puberty, and the expression was significantly increased and then significantly decreased from prepuberty to postpuberty in the hypothalamus, pituitary, and ovary and oviduct. Based on expressed single nucleotide polymorphism (SNP), the NDN gene was expressed monoallelically in the tissues of adult and neonatal umbilical cords, and the expressed allele was paternally inherited. The NDN promoter region of 3400 bp was obtained by cloning and identified in monoallelic-expressing tissues (hypothalamus, ovary, spleen) as a differentially methylated region (DMR).

CONCLUSION

These findings will enrich the number of imprinted genes in sheep and suggest that the NDN gene could be a candidate gene for studying puberty initiation in Dolang sheep.

The impact of age-specific childhood body-mass index on adult cardiometabolic traits: a Mendelian randomization study

Objective

To evaluate the causal relationship between childhood body-mass index (BMI) at different ages and adult cardiometabolic traits.

Methods

We retrieved genetic instrument variables (IVs) for exposures (standardized BMI at newborn, infant, toddler and late childhood), cardiometabolic traits and potential confounders or mediators (adult BMI, SHBG, testosterone and age at menarche) from the corresponding genome-wide association analysis. We performed univariate and multivariable Mendelian randomization (MR) to dissect associations between age-specific childhood BMI and adult cardiometabolic outcomes. Odds ratio was used to present the direction of the causal association.

Results

In univariate MR, higher newborn BMI was causally associated with reduced risk for type 2 diabetes in women. Late childhood BMI was associated with increased risk for female diabetes and coronary artery disease (CAD), myocardial infarction (MI), and chronic kidney disease (CKD) in general population. Among these associations, only association between late childhood BMI with MI remained significant after adjusting for adult male BMI and sex hormones, (OR = 1.120, 95% CI 1.023-1.226, p = 0.014). Besides, in multivariable MR, we found evidence for causal association between newborn BMI with reduced risk for CAD (OR = 0.862, 95% CI 0.751-0.989, p = 0.034) and MI (OR = 0.864, 95% CI 0.752-0.991, p = 0.037) in men. No obvious impact of infant or toddler BMI was identified on the above-mentioned diseases. For continuous cardiometabolic traits, in all age epochs except infant, higher BMI was associated with increased level of fasting glucose in women.

Conclusion

BMI at birth and late childhood exerts different impact on adult cardiometabolic diseases, while BMI at infant and toddler ages is not causally associated with these outcomes. The effect of childhood BMI may be influenced by sex disparities.

Mendelian randomization identifies circulating proteins as biomarkers for age at menarche and age at natural menopause

Age at menarche (AAM) and age at natural menopause (ANM) are highly heritable traits and have been linked to various health outcomes. We aimed to identify circulating proteins associated with altered ANM and AAM using an unbiased two-sample Mendelian randomization (MR) and colocalization approach. By testing causal effects of 1,271 proteins on AAM, we identified 22 proteins causally associated with AAM in MR, among which 13 proteins (GCKR, FOXO3, SEMA3G, PATE4, AZGP1, NEGR1, LHB, DLK1, ANXA2, YWHAB, DNAJB12, RMDN1 and HPGDS) colocalized. Among 1,349 proteins tested for causal association with ANM using MR, we identified 19 causal proteins among which 7 proteins (CPNE1, TYMP, DNER, ADAMTS13, LCT, ARL and PLXNA1) colocalized. Follow-up pathway and gene enrichment analyses demonstrated links between AAM-related proteins and obesity and diabetes, and between AAM and ANM-related proteins and various types of cancer. In conclusion, we identified proteomic signatures of reproductive ageing in women, highlighting biological processes at both ends of the reproductive lifespan.

Long-term health consequences of central precocious/early puberty (CPP) and treatment with Gn-RH analogue: a short update

BACKGROUND

The relationship between precocious or early puberty and its treatment has received significant research attention, yielding diverse outcomes. This short review aims to comprehensively analyze and summarize research articles to elucidate the potential link between precocious or early pubertal onset (CPP) and crucial health factors.

METHODS

We conducted a systematic review of studies published from -January 2000 to March 2023, sourced from databases of Medline, PubMed, Google Scholar and Web of Science. We assessed the relationship between CPP and final adult height (FHt), bone health, reproductive function, body mass index, metabolic and cardiovascular abnormalities, and increased cancer risk.

RESULTS

Upon reviewing and analyzing selected studies, the following key findings emerged: (a) treating CPP in girls before age 6-7 and in boys before age 9 improves FHt; (b) bone mineral density (BMD) decreases during GnRHa treatment but normalizes afterward, with no lasting effects on peak bone mass during puberty; (c) GnRH treatment does not negatively affect menstrual cycles; however, untreated CPP increases the risk of premature or early-onset menopause; (d) the incidence of PCOS/hyperandrogenemia may be slightly elevated in women with a history of CPP, but overall reproductive function remains largely unaffected; (e) earlier thelarche and menarche may enhance susceptibility to breast carcinogenesis; (f) CPP contributes to an increased risk of obesity and type 2 diabetes in both genders; (g) early menarche may slightly increase the risk of coronary heart disease and ischemic strokes and (h) early pubertal timing increases the risk of depression and anxiety disorders.

CONCLUSION

Monitoring and early diagnosis of these conditions are of paramount importance for successful management.

Differential Expression of RNAseq Imprinted Genes from Bovine Females Before and After Puberty

The productivity of beef cows depends on early reproduction traits such as puberty and has an economic impact on the efficiency of production system. Imprinted genes modulate many important endocrine processes such as growth, the onset of puberty and maternal reproductive and behavior. The role of imprinted genes in puberty is a challenging subject since they show the reciprocal role of maternal and paternal genomes in progeny. Although, there are evidences of the involvement of imprint genes in puberty in human, the role of this type of genes in the onset of puberty in cattle has not been studied yet. Here we examined the expression of 27 imprinted genes in pre and post puberty in a bovine model to find differentially expressed imprinted genes in maternal-paternal purebreds and reciprocal crosses across eight tissues and discussed the task of these genes in this crucial process of development and in onset of puberty. DLK1 and MKRN3 that previously described as cause of the central precocious puberty (CPP) in human were differentially expressed in this study. Functional annotation analysis of differentially imprinted genes in different tissues showed significant biological processes of cellular response to growth factor stimulus, response to growth factor, response to parathyroid hormone, developmental growth and the importance of alternative splicing. The results of this study have implications in understanding the role of imprinted genes in the onset of puberty in cattle.

DNA methylation sites in early adulthood characterised by pubertal timing and development: a twin study

BACKGROUND

Puberty is a highly heritable and variable trait, with environmental factors having a role in its eventual timing and development. Early and late pubertal onset are both associated with various diseases developing later in life, and epigenetic characterisation of pubertal timing and development could lead to important insights. Blood DNA methylation, reacting to both genotype and environment, has been associated with puberty; however, such studies are relatively scarce. We investigated peripheral blood DNA methylation profiles (using Illumina 450 K and EPIC platforms) of 1539 young adult Finnish twins associated with pubertal development scale (PDS) at ages 12 and 14 as well as pubertal age (PA).

RESULTS

Fixed effect meta-analysis of the two platforms on 347,521 CpGs in common identified 58 CpG sites associated (p < 1 × 10-5) with either PDS or PA. All four CpGs associated with PA and 45 CpGs associated with PDS were sex-specific. Thirteen CpGs had a high heritability (h2: 0.51-0.98), while one CpG site (mapped to GET4) had a high shared environmental component accounting for 68% of the overall variance in methylation at the site. Utilising twin discordance analysis, we found 6 CpG sites (5 associated with PDS and 1 with PA) that had an environmentally driven association with puberty. Furthermore, genes with PDS- or PA-associated CpGs were consistently linked to various developmental processes and diseases such as breast, prostate and ovarian cancer, while methylation quantitative trait loci of associated CpG sites were enriched in immune pathways developing during puberty.

CONCLUSIONS

By identifying puberty-associated DNA methylation sites and examining the effects of sex, environment and genetics, we shed light on the intricate interplay between environment and genetics in the context of puberty. Through our comprehensive analysis, we not only deepen the understanding of the significance of both genetic and environmental factors in the complex processes of puberty and its timing, but also gain insights into potential links with disease risks.

An atlas of genetic determinants of forearm fracture

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Analysis of the Imprinting Status and Expression of the MAGEL2 Gene During Initiation at Puberty in the Dolang Sheep

Genomic imprinting refers to the expression of parent-specific genes in diploid mammalian cells. MAGEL2 gene is a maternally imprinted gene that has been identified in mice and humans and is associated with the onset of puberty. The purpose of this study was to investigate its imprinting status and its relationship with the onset of puberty in Dolang sheep. The sequence of 3734 bp cDNA of MAGEL2 in Dolang sheep was obtained by cloning and sequencing, encoding 1173 amino acids. The results of the nucleotide and amino acid similarity analysis showed that it was highly conserved among different mammalian species. The MAGEL2 gene was expressed monoallelically in the tissues of adult and neonatal umbilical cords, and the expressed allele was paternally inherited. Real Time quantitative PCR (RT-qPCR) results showed that the MAGEL2 gene was highly expressed in the hypothalamus and pituitary gland, increased significantly from prepuberty to puberty, and decreased significantly after puberty. This study suggests that MAGEL2 is a paternally expressed and maternally imprinted gene in Dolang sheep, which may be involved in the initiation of puberty in Dolang sheep. This study provides a theoretical basis for further research on the mechanism of the imprinted gene MAGEL2 regulating the onset of puberty in sheep, and provides a new idea for the future research on the mechanism of onset of puberty in sheep.

Variations in gender identity and sexual orientation of university students

Background

Previous studies have shown that a small percentage of people in the general population have atypical gender identity and/or sexual orientation.

Aim

This study aimed to explore variations in gender identity and sexual orientation in university students and determine genetic factors associated with these variations.

Methods

Deviations from complete gender congruence and exclusive heterosexual orientation in 736 Japanese university students were quantitatively assessed with self-assessment questionnaires. Next, we conducted genetic tests for 80 participants who showed relatively low gender identity scores and/or atypical sexual orientation. These genetic tests consisted of repeat number analysis of the androgen receptor gene (AR) and a SKAT-O: an optimal unified sequence kernel association test, which is an exome-based rare variant association study. The results of the genetic tests were compared with the Japanese reference data and the results of our 637 control samples.

Outcomes

We calculated the gender identity and sexual orientation scores of all participants and analyzed the molecular data of 80 selected participants.

Results

The gender identity scores of 736 participants were broadly distributed: only ~15% of natal males and ~5% of natal females had the maximum score that corresponds to complete gender congruence. The sexual orientation scores also varied: ~80% of natal males and ~60% of natal females showed exclusive heterosexual orientation. We found no association between gender characteristics and AR repeat numbers. The SKAT-O showed that rare damaging variants of TDRP and 3 other genes were more common in the 80 participants than in the control group.

Clinical Implications

Our data support the view that gender is a phenotypic continuum rather than a binary trait.

Strength and Limitations

This study quantitatively assessed the gender characteristics of a large cohort of university students. Moreover, we conducted systematic screening for genetic factors associated with gender variations. The weaknesses of the study were the limited analytic power of the questionnaires, the relatively small sample for molecular analyses, and incomplete clinical information and relatively advanced ages of the control group.

Conclusion

This study revealed significant variations in gender identity and sexual orientation in university students, which may be partly associated with variants in TDRP or other genes.

Effect of women's fertility and sexual development on epigenetic clock: Mendelian randomization study

BACKGROUND AND OBJECTIVES

In observational studies, women's fertility and sexual development traits may have implications for DNA methylation patterns, and pregnancy-related risk factors can also affect maternal DNA methylation patterns. The aim of our study is to disentangle any potential causal associations between women's fertility and sexual development traits and epigenetic clocks, as well as to search for probable mediators by using the Mendelian randomization (MR) method.

METHODS

Instrumental variables for exposures, mediators, and outcomes were adopted from genome-wide association studies data of European ancestry individuals. The potential causal relationship between women's fertility and sexual development traits and four epigenetic clocks were evaluated by inverse variance weighted method and verified by other two methods. Furthermore, we employed multivariable MR (MVMR) adjusting for hypertension, hyperglycemia, BMI changes, and insomnia. Then, combining the MVMR results and previous research, we performed two-step MR to explore the mediating effects of BMI, AFS, and AFB. Multiple sensitivity analyses were further performed to verify the robustness of our findings.

RESULTS

Leveraging two-sample MR analysis, we observed statistically significant associations between earlier age at first birth (AFB) with a higher HannumAge, PhenoAge and GrimAge acceleration(β = - 0.429, 95% CI [- 0.781 to - 0.077], p = 0.017 for HannumAge; β = - 0.571, 95% CI [- 1.006 to - 0.136], p = 0.010 for PhenoAge, and β = - 1.136, 95% CI [- 1.508 to - 0.765], p = 2.03E-09 for GrimAge respectively) and age at first sexual intercourse (AFS) with a higher HannumAge and GrimAge acceleration(β = - 0.175, 95% CI [- 0.336 to - 0.014], p = 0.033 for HannumAge; β = - 0.210, 95% CI [- 0.350 to - 0.070], p = 0.003 for GrimAge, respectively). Further analyses indicated that BMI, AFB and AFS played mediator roles in the path from women's fertility and sexual development traits to epigenetic aging.

CONCLUSIONS

Our study suggested that AFS and AFB are associated with epigenetic aging. These findings may prove valuable in informing the development of prevention strategies and interventions targeted towards women's fertility and sexual development experiences and their relationship with epigenetic aging-related diseases.

Establishing the relationships between adiposity and reproductive factors: a multivariable Mendelian randomization analysis

BACKGROUND

Few studies have investigated associations between adiposity and reproductive factors using causal methods, both of which have a number of consequences on women's health. Here we assess whether adiposity at different points in the lifecourse affects reproductive factors differently and independently, and the plausibility of the impact of reproductive factors on adiposity.

METHODS

We used genetic data from UK Biobank (273,238 women) and other consortia (EGG, GIANT, ReproGen and SSGAC) for eight reproductive factors: age at menarche, age at menopause, age at first birth, age at last birth, number of births, being parous, age first had sexual intercourse and lifetime number of sexual partners, and two adiposity traits: childhood and adulthood body size. We applied multivariable Mendelian randomization to account for genetic correlation and to estimate the causal effects of childhood and adulthood adiposity, independently of each other, on reproductive factors. Additionally, we estimated the effects of reproductive factors, independently of other relevant reproductive factors, on adulthood adiposity.

RESULTS

We found a higher childhood body size leads to an earlier age at menarche, and an earlier age at menarche leads to a higher adulthood body size. Furthermore, we find contrasting and independent effects of childhood and adulthood body size on age at first birth (beta 0.22 SD (95% confidence interval: 0.14, 0.31) vs - 2.49 (- 2.93, - 2.06) per 1 SD increase), age at last birth (0.13 (0.06,0.21) vs - 1.86 (- 2.23, - 1.48) per 1 SD increase), age at menopause (0.17 (0.09, 0.25) vs - 0.99 (- 1.39, - 0.59) per 1 SD increase), and likelihood of having children (Odds ratio 0.97 (0.95, 1.00) vs 1.20 (1.06, 1.37) per 1 SD increase).

CONCLUSIONS

Our findings demonstrate the importance of considering a lifecourse approach when investigating the inter-relationships between adiposity measures and reproductive events, as well as the use of 'age specific' genetic instruments when evaluating lifecourse hypotheses in a Mendelian randomization framework.

Health outcomes of age at menarche in European women: a two-sample Mendelian randomization study

BACKGROUND

Observational studies have shown an association between age at menarche (AAM) and the risk of gynecological diseases. However, the causality cannot be determined due to residual confounding.

METHODS

We conducted a Mendelian randomization (MR) study to evaluate the causal effect of AAM on several gynecological diseases, including endometriosis, female infertility, pre-eclampsia or eclampsia, uterine fibroids, breast cancer, ovarian cancer, and endometrial cancer. Single nucleotide polymorphisms were used as genetic instruments. The inverse variance weighted method was used as the primary approach and several other MR models were conducted for comparison. Cochran's Q test, Egger's intercept test, and leave-one-out analysis were conducted for sensitivity analysis. Radial MR analysis was conducted when detecting the existence of heterogeneity.

RESULTS

After Bonferroni correction and thorough sensitivity analysis, we observed a robust causal effect of AAM on endometrial cancer (odds ratio: 0.80; 95% confidence interval: 0.72-0.89; P = 4.61 × 10-5) and breast cancer (odds ratio: 0.94; 95% confidence interval: 0.90-0.98; P = .003). Sensitivity analysis found little evidence of horizontal pleiotropy. The inverse variance weighted method also detected weak evidence of associations of AAM with endometriosis and pre-eclampsia or eclampsia.

CONCLUSIONS

This MR study demonstrated a causal effect of AAM on gynecological diseases, especially for breast cancer and endometrial cancer, which indicates AAM might be a promising index to use for disease screening and prevention in clinical practice. Key messages What is already known on this topic - Observational studies have reported associations between age at menarche (AAM) and a variety of gynecological diseases but the causality has not been determined. What this study adds - This Mendelian randomization study demonstrated that AAM causally affects the risk of breast cancer and endometrial cancer. How this study might affect research, practice, or policy - The findings of our study imply that AAM could be a candidate marker for early screening of populations at higher risk of breast cancer and endometrial cancer.

Global endometrial DNA methylation analysis reveals insights into mQTL regulation and associated endometriosis disease risk and endometrial function

Endometriosis is a leading cause of pain and infertility affecting millions of women globally. Herein, we characterize variation in DNA methylation (DNAm) and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genotype data and methylation in endometrial samples from 984 deeply-phenotyped participants. We estimate that 15.4% of the variation in endometriosis is captured by DNAm and identify significant differences in DNAm profiles associated with stage III/IV endometriosis, endometriosis sub-phenotypes and menstrual cycle phase, including opening of the window for embryo implantation. Menstrual cycle phase was a major source of DNAm variation suggesting cellular and hormonally-driven changes across the cycle can regulate genes and pathways responsible for endometrial physiology and function. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independent cis-mQTLs including 51 associated with risk of endometriosis, highlighting candidate genes contributing to disease risk. Our work provides functional evidence for epigenetic targets contributing to endometriosis risk and pathogenesis. Data generated serve as a valuable resource for understanding tissue-specific effects of methylation on endometrial biology in health and disease.

Large-scale exome sequence analysis identifies sex- and age-specific determinants of obesity

Obesity contributes substantially to the global burden of disease and has a significant heritable component. Recent large-scale exome sequencing studies identified several genes in which rare, protein-coding variants have large effects on adult body mass index (BMI). Here we extended such work by performing sex-stratified associations in the UK Biobank study (N∼420,000). We identified genes in which rare heterozygous loss-of-function increases adult BMI in women (DIDO1, PTPRG, and SLC12A5) and in men (SLTM), with effect sizes up to ∼8 kg/m2. This is complemented by analyses implicating rare variants in OBSCN and MADD for recalled childhood adiposity. The known functions of these genes, as well as findings of common variant genome-wide pathway enrichment analyses, suggest a role for neuron death, apoptosis, and DNA damage response mechanisms in the susceptibility to obesity across the life-course. These findings highlight the importance of considering sex-specific and life-course effects in the genetic regulation of obesity.

Bidirectional two-sample mendelian randomization analysis identifies causal associations between age at menarche and sarcopenia-related traits

BACKGROUND

Observational studies have shown that the age of menarche is associated with sarcopenia, but confounding factors make the causal relationship difficult to infer.

OBJECTIVE

Therefore, we conducted a bidirectional two-sample Mendelian randomized (MR) analysis to evaluate the potential causal relationship between age at menarche and sarcopenia-related traits (hand grip strength, lean mass, walking pace).

METHODS

We obtained the latest aggregate statistics from the Genome-wide association studies (GWAS) database on the age of menarche of 182,416 participants from ReproGen, the appendicular lean mass of 244,730 participants from EMBL's European Bioinformatics Institute, the left-hand grip strength of 401,026 participants, the right-hand grip strength of 461,089 participants and the usual walking pace of 459,915 participants from the UK Biobank. The inverse variance weighting (IVW) method and other MR methods were used to evaluate the bidirectional causal relationship between the age of menarche and sarcopenia.

RESULTS

The forward MR results showed that the age of menarche predicted by the gene was positively correlated with left-hand grip strength (IVWβ=0.041, P = 2.00 × 10-10), right-hand grip strength (IVWβ=0.053, P = 1.97 × 10-18), appendicular lean mass (IVWβ=0.012, P = 4.38 × 10-13) and usual walking pace (IVWβ=0.033, P = 1.62 × 10-8).In the reverse MR analysis, we also found that the usual walking pace was positively correlated with the age of menarche predicted by genes (IVWβ=0.532, P = 1.65 × 10-4). Still, there was no causal relationship between grip strength and appendicular lean mass and the age at menarche.

CONCLUSION

Our results show that earlier menarche will increase the risk of sarcopenia. In addition, people with higher muscle function tend to have menarche later. These findings may provide a reference for prevention strategies and interventions for menarche in advance and sarcopenia.

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.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Novel MKRN3 Missense Mutations Associated With Central Precocious Puberty Reveal Distinct Effects on Ubiquitination

CONTEXT

Loss-of-function mutations in the maternally imprinted genes, MKRN3 and DLK1, are associated with central precocious puberty (CPP). Mutations in MKRN3 are the most common known genetic etiology of CPP.

OBJECTIVE

This work aimed to screen patients with CPP for MKRN3 and DLK1 mutations and analyze the effects of identified mutations on protein function in vitro.

METHODS

Participants included 84 unrelated children with CPP (79 girls, 5 boys) and, when available, their first-degree relatives. Five academic medical institutions participated. Sanger sequencing of MKRN3 and DLK1 5' upstream flanking and coding regions was performed on DNA extracted from peripheral blood leukocytes. Western blot analysis was performed to assess protein ubiquitination profiles.

RESULTS

Eight heterozygous MKRN3 mutations were identified in 9 unrelated girls with CPP. Five are novel missense mutations, 2 were previously identified in patients with CPP, and 1 is a frameshift variant not previously associated with CPP. No pathogenic variants were identified in DLK1. Girls with MKRN3 mutations had an earlier age of initial pubertal signs and higher basal serum luteinizing hormone and follicle-stimulating hormone compared to girls with CPP without MRKN3 mutations. Western blot analysis revealed that compared to wild-type MKRN3, mutations within the RING finger domain reduced ubiquitination whereas the mutations outside this domain increased ubiquitination.

CONCLUSION

MKRN3 mutations were present in 10.7% of our CPP cohort, consistent with previous studies. The novel identified mutations in different domains of MKRN3 revealed different patterns of ubiquitination, suggesting distinct molecular mechanisms by which the loss of MRKN3 results in early pubertal onset.

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.

Clinical and Genetic Characterization of Familial Central Precocious Puberty

CONTEXT

Central precocious puberty (CPP) can have a familial form in approximately one-quarter of the children. The recognition of this inherited condition increased after the identification of autosomal dominant CPP with paternal transmission caused by mutations in the MKRN3 and DLK1 genes.

OBJECTIVE

We aimed to characterize the inheritance and estimate the prevalence of familial CPP in a large multiethnic cohort; to compare clinical and hormonal features, as well as treatment response to GnRH analogs (GnRHa), in children with distinct modes of transmission; and to investigate the genetic basis of familial CPP.

METHODS

We retrospectively studied 586 children with a diagnosis of CPP. Patients with familial CPP (n = 276) were selected for clinical and genetic analysis. Data from previous studies were grouped, encompassing sequencing of MKRN3 and DLK1 genes in 204 patients. Large-scale parallel sequencing was performed in 48 individuals from 34 families.

RESULTS

The prevalence of familial CPP was estimated at 22%, with a similar frequency of maternal and paternal transmission. Pedigree analyses of families with maternal transmission suggested an autosomal dominant inheritance. Clinical and hormonal features, as well as treatment response to GnRHa, were similar among patients with different forms of transmission of familial CPP. MKRN3 loss-of-function mutations were the most prevalent cause of familial CPP, followed by DLK1 loss-of-function mutations, affecting, respectively, 22% and 4% of the studied families; both affected exclusively families with paternal transmission. Rare variants of uncertain significance were identified in CPP families with maternal transmission.

CONCLUSION

We demonstrated a similar prevalence of familial CPP with maternal and paternal transmission. MKRN3 and DLK1 loss-of-function mutations were the major causes of familial CPP with paternal transmission.

MKRN3 inhibits puberty onset via interaction with IGF2BP1 and regulation of hypothalamic plasticity

Makorin ring finger protein 3 (MKRN3) was identified as an inhibitor of puberty initiation with the report of loss-of-function mutations in association with central precocious puberty. Consistent with this inhibitory role, a prepubertal decrease in Mkrn3 expression was observed in the mouse hypothalamus. Here, we investigated the mechanisms of action of MKRN3 in the central regulation of puberty onset. We showed that MKRN3 deletion in hypothalamic neurons derived from human induced pluripotent stem cells was associated with significant changes in expression of genes controlling hypothalamic development and plasticity. Mkrn3 deletion in a mouse model led to early puberty onset in female mice. We found that Mkrn3 deletion increased the number of dendritic spines in the arcuate nucleus but did not alter the morphology of GnRH neurons during postnatal development. In addition, we identified neurokinin B (NKB) as an Mkrn3 target. Using proteomics, we identified insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) as another target of MKRN3. Interactome analysis revealed that IGF2BP1 interacted with MKRN3, along with several members of the polyadenylate-binding protein family. Our data show that one of the mechanisms by which MKRN3 inhibits pubertal initiation is through regulation of prepubertal hypothalamic development and plasticity, as well as through effects on NKB and IGF2BP1.

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.

Association of the Estrogen Receptor 1 Polymorphisms rs2046210 and rs9383590 with the Risk, Age at Onset and Prognosis of Breast Cancer

Estrogen receptor α (ERα), encoded by the ESR1 gene, is a key prognostic and predictive biomarker firmly established in routine diagnostics and as a therapeutic target of breast cancer, and it has a central function in breast cancer biology. Genetic variants at 6q25.1, containing the ESR1 gene, were found to be associated with breast cancer susceptibility. The rs2046210 and rs9383590 single nucleotide variants (SNVs) are located in the same putative enhancer region upstream of ESR1 and were separately identified as candidate causal variants responsible for these associations. Here, both SNVs were genotyped in a hospital-based case-control study of 409 female breast cancer patients and 422 female controls of a Central European (Austrian) study population. We analyzed the association of both SNVs with the risk, age at onset, clinically and molecularly relevant characteristics and prognosis of breast cancer. We also assessed the concordances between both SNVs and the associations of each SNV conditional on the other SNV. The minor alleles of both SNVs were found to be non-significantly associated with an increased breast cancer risk. Significant associations were found in specific subpopulations, particularly in patients with an age younger than 55 years. The minor homozygotes of rs2046210 and the minor homozygotes plus heterozygotes of rs9383590 exhibited a several-years-younger age at onset than the common homozygotes, which was more pronounced in ER-positive and luminal patients. Importantly, the observed associations of each SNV were not consistently nullified upon correction for the other SNV nor upon analyses in common homozygotes for the other SNV. We conclude that both SNVs remain independent candidate causal variants.

Causal relationships between birth weight, childhood obesity and age at menarche: A two-sample Mendelian randomization analysis

OBJECTIVES

Observational studies suggest birth weight and childhood obesity are closely associated with age at menarche. However, the relationships between them are currently inconsistent and it remains elusive whether such associations are causal. Therefore, the aim of the study was to investigate whether there existed causal relationships between birth weight, childhood obesity and age at menarche.

DESIGN, PATIENTS AND MEASUREMENTS

A two-sample Mendelian randomization (MR) study. The standard inverse variance weighted MR analyses were adopted to evaluate the causal effects of birth weight (n = 143,677), childhood body mass index (BMI) (n = 39,620) on age at menarche (n = 182,416) with summary statistics from large-scale genome-wide association studies (GWASs). Meanwhile, we validated our MR results with some sensitivity analyses including maximum likelihood, weighted-median and MR pleiotropy residual sum and outlier methods.

RESULTS

The present study showed that each one standard deviation (1-SD) lower birth weight was predicted to result in a 0.1479 years earlier of age at menarche (β = .1479, 95% confidence interval [CI] = 0.0422-0.2535; p = 0.0061). We also found that genetically predicted 1-SD increase in childhood BMI was causally associated with early age at menarche (β = -.3966, 95% CI = -0.5294 to -0.2639; p = 4.73E-09).

CONCLUSIONS

Our MR study suggests the causal effect of lower birth weight and higher childhood BMI on the increased risk of earlier menarche. It may be the opportune time to carry out weight control intervention in prenatal and early childhood development periods to prevent early menarche onset, thus decreasing the future adverse consequences.

Genome-wide association study of age at menarche in the Taiwan Biobank suggests NOL4 as a novel associated gene

Sexual maturation is a complex physiological process that involves multiple variables, such as genetic and environmental factors. Among females, age at menarche (AM) is a critical milestone for sexual maturation. This study aimed to identify genetic markers of AM using nationwide population cohort data in Taiwan. Females with self-reported AM between 10 and 16 years (N = 39,827) were eligible for the final analysis. To identify genetic signals related to AM, we conducted a genome-wide association study using a linear regression model and split-half meta-analysis method to verify our findings. The Functional Mapping and Annotation web-based platform was used for positional mapping and gene-based and gene-set analyses. The meta-analysis identified four significant loci, i.e., LIN28B (pooled P = 1.39 × 10^-21), NOL4 (pooled P = 8.94 × 10^-9), GPR45 (pooled P = 4.19 × 10^-11), and LOC105373831 (pooled P = 4.37 × 10^-8), that were associated with AM. MAGMA gene-based analysis revealed that LIN28B (P = 1.13 × 10^-8), NOL4 (P = 2.27 × 10^-7), RXRG (P = 4.34 × 10-⁷), ETV5 (P = 1.75 × 10^-6), and HACE1 (P = 1.82 × 10^-6) were significantly associated with AM, while the gene-set analysis identified a significantly enriched pathway involving mTOR signaling complex (FDR corrected P = 1.28 × 10^-2). The results replicated evidence for several genetic markers associated with AM in the Taiwanese female population. Our analysis identified a novel locus (rs7239368) in NOL4 associated with AM (β = 0.051 ± 0.009 years, pooled P = 8.94 × 10^-9), whereas additional research is needed to validate its molecular role in sexual maturation.

Effects of half-dose spiomet treatment in girls with early puberty and accelerated bone maturation: a multicenter, randomized, placebo-controlled study protocol

BACKGROUND

A "mismatch" sequence of less prenatal weight gain and more postnatal weight gain may lead to ectopic lipid accumulation, and trigger the development of early adrenarche/pubarche and the activation of the gonadotropic axis resulting in early puberty and ending up in full-blown adolescent polycystic ovary syndrome (PCOS). In the present study, we assess whether a low-dose combination of generics that collectively reduce ectopic fat through different pathways can slow down the accelerated maturation in "mismatch" girls with early puberty.

METHODS

Randomized, placebo-controlled, multicenter, phase 2a, study in 64 girls [age, 8.0-9.3 years; birthweight (BW) for gestational age in lower tertile (-1.96< Z-score <-0.44), body mass index (BMI) in upper tertile (+0.44< Z-score < +1.96) and early progressive puberty (Tanner B2 at 7.7-9.0 years)]. Pharmacological intervention will be with a half-dose version of SPIOMET (mini-spiomet), a combination that reverts the PCOS phenotype in "mismatch" adolescents; mini-spiomet will contain spironolactone (25 mg/day, to raise brown adipose tissue activity), pioglitazone (3.75 mg/day, to raise adiponectin and insulin sensitivity), and metformin (425 mg/day, to raise AMPK activity and GDF15). Recruitment: 1 year; double-blind treatment: 1 year; open follow-up: 1 year; analyses and reporting: 1 year.

INTERVENTIONS

randomization (1:1) for placebo vs mini-spiomet.

PRIMARY OUTCOME

annualized bone age advancement (0-1 year) by BoneXpert; secondary outcomes: insulin, IGF-I, high-molecular-weight adiponectin (HMW-adip), sex hormone binding globulin (SHBG), ultra-sensitive C-reactive protein (usCRP), androgens, luteinizing hormone (LH), follicle-stimulating hormone (FSH), oestradiol, growth-and-differentiation factor 15 (GDF15), C-X-C motif chemokine ligand-14 (CXCL14), safety parameters, and quantification of hepato-visceral fat.

DISCUSSION

The present study, if successful, may provide a first proof of the concept that the rapid maturation of girls with an upward mismatch between pre- and post-natal weight gain can be slowed down with a fixed low-dose combination of old and safe generics jointly targeting a reduction of ectopic fat without necessarily lowering body weight.

TRIAL REGISTRATION

EudraCT 2021-006766-21. Registered on May 30, 2022.

Identification and immune features of cuproptosis-related molecular clusters in polycystic ovary syndrome

Polycystic ovary syndrome (PCOS), a common reproductive endocrine disease, has clinically heterogeneous characteristics. Recently, cuproptosis causes several diseases by killing cells. Hence, we aimed to explore cuproptosis-related molecular clusters in PCOS and construct a prediction model. Based on the GSE5090, GSE43264, GSE98421, and GSE124226 datasets, an analysis of cuproptosis regulators and immune features in PCOS was conducted. In 25 cases of PCOS, the molecular clusters of cuproptosis-related genes and the immune cell infiltration associated with PCOS were investigated. Weighted gene co-expression network analysis was used to identify differentially expressed genes within clusters. Next, we compared the performance of the random forest model, support vector machine model, generalized linear model, and eXtreme Gradient Boosting for deciding the optimum machine model. Validation of the predictive effectiveness was accomplished through nomogram, calibration curve, decision curve analysis, and using other two datasets. PCOS and non-PCOS controls differed in the dysregulation of cuproptosis-related genes and the activation of immunoreaction. Two cuproptosis-related molecular clusters associated with PCOS were identified. Significant heterogeneity was noted in immunity between the two clusters based on the analysis of immune infiltration. The immune-related pathways related to cluster-specific differentially expressed genes in Cluster1 were revealed by functional analysis. With a relatively low residual error and root mean square error and a higher area under the curve (1.000), the support vector machine model demonstrated optimal discriminative performance. An ultimate 5-gene-based support vector machine model was noted to perform satisfactorily in the other two validation datasets (area under the curve = 1.000 for both). Moreover, the nomogram, calibration curve, and decision curve analysis showed that PCOS subtypes can be accurately predicted. Our study results helped demonstrate a comprehensive understanding of the complex relationship between cuproptosis and PCOS and establish a promising prediction model for assessing the risk of cuproptosis in patients with PCOS.