Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

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.

Update Breast Cancer 2023 Part 1 - Early Stage Breast Cancer

With abemaciclib (monarchE study) and olaparib (OlympiA study) gaining approval in the adjuvant treatment setting, a significant change in the standard of care for patients with early stage breast cancer has been established for some time now. Accordingly, some diverse developments are slowly being transferred from the metastatic to the adjuvant treatment setting. Recently, there have also been positive reports of the NATALEE study. Other clinical studies are currently investigating substances that are already established in the metastatic setting. These include, for example, the DESTINY Breast05 study with trastuzumab deruxtecan and the SASCIA study with sacituzumab govitecan. In this review paper, we summarize and place in context the latest developments over the past months.

Causality between COVID‐19 and female reproductive function: A Mendelian randomization study

Coronavirus disease 2019 (COVID‐19) has experienced a global pandemic, and currently, the emergence of its variants has posed challenges in terms of prevention and treatment. Nonetheless, the effect of COVID‐19 infection on female reproductive function is unclear. This study aimed to systematically evaluate for the first time the causal effect of COVID‐19 on female reproductive function. Genetic correlations were assessed using linkage disequilibrium score regression. Mendelian randomization (MR) analysis was performed using summary statistics of two variables, including COVID‐19 severity and eight female reproductive traits. The three degrees of severity had genetically significant associations with sex hormone‐binding globulin (SHBG) concentrations (rg = –0.153, p = 0.004; rg = –0.187, p < 0.001; rg = –0.180, p = 0.003). Additionally, MR showed that SHBG (β = –0.020, p = 0.040) and total testosterone levels (β = –0.061, p = 0.009) followed a decreasing trend, as the COVID‐19 infection higher. No significant genetic association was found between COVID‐19 infection and total estradiol concentrations, menstruation, and female infertility. Simultaneously, MR found no causal relationships between COVID‐19 infection and total estradiol concentrations, menstruation, and female infertility (all p > 0.05). In conclusion, COVID‐19 was causally associated with lower SHBG and total testosterone concentrations, offering invaluable insights that will help guide clinical decision‐making.

Multi-trait genome-wide association study identifies a novel endometrial cancer risk locus that associates with testosterone levels

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

Father absence and pubertal timing in Korean boys and girls

Background and objectives

Pubertal timing is a key life history trait with long-term health consequences in both sexes. Evolutionary theory has guided extensive research on developmental influences, in particular growing up without a father, on earlier menarche. Far less is known whether a similar relationship exists for boys, especially beyond western contexts. We used longitudinal data from the nationally representative sample of Korean adolescents, which provided us with a unique opportunity for studying male puberty using a hitherto underutilized biomarker: age at first nocturnal ejaculation.

Methodology

We pre-registered and tested a prediction that growing up in father-absent households is associated with earlier puberty in both sexes. Large sample size (>6000) allowed testing the effect of father absence, which remains relatively uncommon in Korea, while adjusting for potential confounders using Cox proportional-hazard models.

Results

Self-reported age at first nocturnal ejaculation was on average 13.8 years, falling within the range known from other societies. Unlike previous findings mostly for white girls, we did not find evidence that Korean girls in father-absent households had a younger age at menarche. Boys in father-absent households reported having their first nocturnal ejaculation 3 months earlier on average, and the difference was evident before age 14.

Conclusion and implications

The association between father absence and pubertal timing appears sex- and age-dependent, and these differences may further interact with cultural norms regarding gender roles. Our study also highlights the utility of the recalled age of first ejaculation for male puberty research, which has lagged in both evolutionary biology and medicine.

The genetics of autism and steroid-related traits in prenatal and postnatal life

Background

Autism likelihood is a largely heritable trait. Autism prevalence has a skewed sex ratio, with males being diagnosed more often than females. Steroid hormones play a mediating role in this, as indicated by studies of both prenatal biology and postnatal medical conditions in autistic men and women. It is currently unclear if the genetics of steroid regulation or production interact with the genetic liability for autism.

Methods

To address this, two studies were conducted using publicly available datasets, which focused respectively on rare genetic variants linked to autism and neurodevelopmental conditions (study 1) and common genetic variants (study 2) for autism. In Study 1 an enrichment analysis was conducted, between autism-related genes (SFARI database) and genes that are differentially expressed (FDR<0.1) between male and female placentas, in 1st trimester chorionic villi samples of viable pregnancies (n=39). In Study 2 summary statistics of genome wide association studies (GWAS) were used to investigate the genetic correlation between autism and bioactive testosterone, estradiol and postnatal PlGF levels, as well as steroid-related conditions such as polycystic ovaries syndrome (PCOS), age of menarche, and androgenic alopecia. Genetic correlation was calculated based on LD Score regression and results were corrected for multiple testing with FDR.

Results

In Study 1, there was significant enrichment of X-linked autism genes in male-biased placental genes, independently of gene length (n=5 genes, p<0.001). In Study 2, common genetic variance associated with autism did not correlate to the genetics for the postnatal levels of testosterone, estradiol or PlGF, but was associated with the genotypes associated with early age of menarche in females (b=-0.109, FDR-q=0.004) and protection from androgenic alopecia for males (b=-0.135, FDR-q=0.007).

Conclusion

The rare genetic variants associated with autism appear to interact with placental sex differences, while the common genetic variants associated with autism appear to be involved in the regulation of steroid-related traits. These lines of evidence indicate that the likelihood for autism is partly linked to factors mediating physiological sex differences throughout development.

Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.

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.

Women's reproductive traits and major depressive disorder: A two-sample Mendelian randomization study

BACKGROUND

Evidence suggested strong associations between women's reproductive factors and major depressive disorder (MDD), but their causalities are unclear.

METHODS

Using female-specific SNPs as genetic instruments obtained from large-scale genome-wide association studies for women's reproductive traits, we designed two-sample univariable and multivariable Mendelian randomization (MR) analysis to evaluate the causal effects of women's reproductive traits on MDD. For both univariable MR (UVMR) and multivariable MR (MVMR), the inverse variance weighting estimates were reported as main results. MR-Egger, weighted median, and generalized summary-data-based MR (GSMR) methods for UVMR, and MVMR-Egger and MVMR-robust methods for MVMR were used as sensitivity analyses. Negative control analyses, MVMR of age at first birth (AFB) and age at first sexual intercourse (AFS) on MDD, and sex-combined genetic variants for AFB and AFS were performed to enhance the robustness of our study.

RESULTS

There was substantial evidence for associations of genetically predicted later age at menarche (AAM) (odds ratio (OR) = 0.97, 95 % confidence interval (CI) = 0.94-0.99, P = 0.007), AFB (OR = 0.91, 95 % CI = 0.86-0.97, P = 0.002) and AFS (OR = 0.70, 95 % CI = 0.60-0.80, P < 0.001) with lower MDD risk in UVMR. After adjustment of BMI and educational attainment using MVMR, we found consistently significant causal effects of AAM (OR = 0.95, 95 % CI = 0.92-0.99, P = 0.006), AFB (OR = 0.88, 95 % CI = 0.84-0.91, P < 0.001) and AFS (OR = 0.71, 95 % CI = 0.64-0.79, P < 0.001) on MDD.

CONCLUSIONS

Our results provide compelling evidence that early AAM, AFB, and AFS are risk factors for MDD. Promoting the cognition of reproductive health care for women may reduce the risk of MDD.

Genetic prediction of age at menarche, age at natural menopause and type 2 diabetes: A Mendelian randomization study

BACKGROUND AND AIMS

The relationship between reproductive factors and type 2 diabetes (T2D) is controversial; therefore, we explored the causal relationship of age at menarche (AAM), age at natural menopause (ANM), with the risk of T2D and glycemic traits using two-sample Mendelian randomization.

METHODS AND RESULTS

We used publicly available data at the summary level of genome-wide association studies, where AAM (N = 329,345), ANM (N = 69,360), T2D (N = 464,389). The inverse variance weighting (IVW) method was employed as the primary method. To demonstrate the robustness of the results, we also conducted various sensitivity analysis methods including the MR-Egger regression, the weighted median (WM) and the MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) test. After excluding IVs associated with confounders, we found a causal association between later AAM and reduced risk of T2D (OR 0.81 [95% CI 0.75, 0.87]; P = 2.20 × 10^-8), lower levels of FI (β -0.04 [95% CI -0.06, -0.01]; P = 2.19 × 10^-3), FPG (β -0.03 [95% CI -0.05, -0.007]; P = 9.67 × 10^-5) and HOMA-IR (β -0.04 [95% CI -0.06, -0.01]; P = 4,95 × 10^-3). As for ANM, we only found a causal effect with HOMA-IR (β -0.01 [95% CI -0.02, -0.005]; P = 1.77 × 10^-3), but not with T2D.

CONCLUSIONS

Our MR study showed a causal relationship between later AAM and lower risk of developing T2D, lower FI, FPG and HOMA-IR levels. This may provide new insights into the prevention of T2D in women.

Association between early life adversity and allostatic load in girls with precocious puberty

CONTEXT

The mechanisms underlying the elevated long-term health risk in girls with precocious puberty remain unclear, but might result from physiological wear and tear associated with greater exposure to early life adversity.

OBJECTIVE

This study aims to explore early life adversity in girls with precocious puberty and its association with allostatic load.

METHODS

Early life adversity and hair cortisol concentration were measured among 213 girls with precocious puberty (8.21 ± 1.07). Allostatic load score is constructed by using 13 physiological biomarkers representing four systems and hair cortisol concentration. Multivariate linear regression models have estimated the associations between cumulative early life adversity exposure with total and system-specific allostatic load scores. Associations between cumulative early life adversity and the risk of high allostatic load (3 + high-risk biomarkers) were tested using binary logistics regression.

RESULTS

More than two-thirds (67.6%) of girls with central precocious puberty reported two or more early life adversity exposure. Compared to those with no early life adversity exposure, girls who reported early life adversity score ≥ 2 had significantly higher total allostatic load score (β: 1.20-1.64, P < 0.001). Metabolic system was more sensitive to cumulative early life adversity exposure, each form of early life adversity exposure was associated with 0.48-unit increases in metabolic allostatic load score (95%CI: 0.06, 0.90, P = 0.026). Girls reported early life adversity score ≥ 3 were three times more likely to have a high allostatic load compared with those without early life adversity exposure in both unadjusted and adjusted models (OR_adjusted=3.83, 95%CI: 1.17, 12.55, P = 0.001).

CONCLUSION

Multisystem physiological dysregulation is observed in girls with central precocious puberty, which might result from cumulative wear-and-tear associated with early life adversity.

Women's reproductive traits and cerebral small-vessel disease: A two-sample Mendelian randomization study

Background

Observational studies have suggested that women's reproductive factors (age at menarche (AAM), age at first birth (AFB), age at first sexual intercourse (AFS), age at natural menopause (ANM), and pregnancy loss) may influence the risk of cerebral small-vessel disease (CSVD) although the causality remains unclear.

Methods

We conducted two-sample univariable Mendelian randomization (UVMR) and multivariable MR (MVMR) to simultaneously investigate the causal relationships between five women's reproductive traits and CSVD clinical [intracerebral hemorrhage (ICH) by location or small-vessel ischemic stroke (SVS)] and subclinical measures [white matter hyperintensities (WMH), fractional anisotropy (FA), and mean diffusivity (MD)], utilizing data from large-scale genome-wide association studies of European ancestry. For both UVMR and MVMR, the inverse-variance-weighted (IVW) estimates were reported as the main results. The MR-Egger, weighted median, generalized summary-data-based MR (GSMR), and MR-pleiotropy residual sum and outlier (MR-PRESSO) methods for UVMR and MVMR-Egger, and the MVMR-robust methods for MVMR were used as sensitivity analyses. Sex-combined instruments for AFS and AFB were used to assess the impact of sex instrumental heterogeneity. Positive control analysis was implemented to measure the efficacy of selected genetic instruments.

Results

We found no evidence to support causal associations between genetic liability for women's reproductive factors and the risk of CSVD in UVMR (all P-values &gt; 0.05). Using MVMR, the results were consistent with the findings of UVMR after accounting for body mass index and educational attainment (all P-values &gt; 0.05). Sensitivity analyses also provided consistent results. The putative positive causality was observed between AAM, ANM, and ovarian cancer, ensuring the efficacy of selected genetic instruments.

Conclusion

Our findings do not convincingly support a causal effect of women's reproductive factors on CSVD. Future studies are warranted to investigate specific estrogen-related physiological changes in women, which may inform current researchers on the causal mechanisms involved in cerebral small-vessel disease progression.

Sex-Specific Reproductive Factors Augment Cardiovascular Disease Risk in Women: A Mendelian Randomization Study

Background Observational studies suggest that reproductive factors are associated with cardiovascular disease, but these are liable to influence by residual confounding. This study explores the causal relevance of reproductive factors on cardiovascular disease in women using Mendelian randomization. Methods and Results Uncorrelated (r2<0.001), genome-wide significant (P<5×10-8) single-nucleotide polymorphisms were extracted from sex-specific genome-wide association studies of age at first birth, number of live births, age at menarche, and age at menopause. Inverse-variance weighted Mendelian randomization was used for primary analyses on outcomes of atrial fibrillation, coronary artery disease, heart failure, ischemic stroke, and stroke. Earlier genetically predicted age at first birth increased risk of coronary artery disease (odds ratio [OR] per year, 1.49 [95% CI, 1.28-1.74], P=3.72×10-7) heart failure (OR, 1.27 [95% CI, 1.06-1.53], P=0.009), and stroke (OR, 1.25 [95% CI, 1.00-1.56], P=0.048), with partial mediation through body mass index, type 2 diabetes, blood pressure, and cholesterol traits. Higher genetically predicted number of live births increased risk of atrial fibrillation (OR for <2, versus 2, versus >2 live births, 2.91 [95% CI, 1.16-7.29], P=0.023), heart failure (OR, 1.90 [95% CI, 1.28-2.82], P=0.001), ischemic stroke (OR, 1.86 [95% CI, 1.03-3.37], P=0.039), and stroke (OR, 2.07 [95% CI, 1.22-3.52], P=0.007). Earlier genetically predicted age at menarche increased risk of coronary artery disease (OR per year, 1.10 [95% CI, 1.06-1.14], P=1.68×10-6) and heart failure (OR, 1.12 [95% CI, 1.07-1.17], P=5.06×10-7); both associations were at least partly mediated by body mass index. Conclusions These results support a causal role of a number of reproductive factors on cardiovascular disease in women and identify multiple modifiable mediators amenable to clinical intervention.

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.

A Review of Major Danish Biobanks: Advantages and Possibilities of Health Research in Denmark

Biobank research may lead to an improved understanding of disease etiology and advance personalized medicine. Denmark (population ~5.9 million) provides a unique setting for population-based health research. The country is a rich source of biobanks and the universal, tax-funded healthcare system delivers routinely collected data to numerous registries and databases. By virtue of the civil registration number (assigned uniquely to all Danish citizens), biological specimens stored in biobanks can be combined with clinical and demographic data from these population-based health registries and databases. In this review, we aim to provide an understanding of advantages and possibilities of biobank research in Denmark. As knowledge about the Danish setting is needed to grasp the full potential, we first introduce the Danish healthcare system, the Civil Registration System, the population-based registries, and the interface with biobanks. We then describe the biobank infrastructures, comprising the Danish National Biobank Initiative, the Bio- and Genome Bank Denmark, and the Danish National Genome Center. Further, we briefly provide an overview of fourteen selected biobanks, including: The Danish Newborn Screening Biobank; The Danish National Birth Cohort; The Danish Twin Registry Biobank; Diet, Cancer and Health; Diet, Cancer and Health - Next generations; Danish Centre for Strategic Research in Type 2 Diabetes; Vejle Diabetes Biobank; The Copenhagen Hospital Biobank; The Copenhagen City Heart Study; The Copenhagen General Population Study; The Danish Cancer Biobank; The Danish Rheumatological Biobank; The Danish Blood Donor Study; and The Danish Pathology Databank. Last, we inform on practical aspects, such as data access, and discuss future implications.

Reassessing the causal role of obesity in breast cancer susceptibility: a comprehensive multivariable Mendelian randomization investigating the distribution and timing of exposure

BACKGROUND

Previous Mendelian randomization (MR) studies on obesity and risk of breast cancer adopted a small number of instrumental variables and focused mainly on the crude total effect. We aim to investigate the independent causal effect of obesity on breast cancer susceptibility, considering the distribution of fat, covering both early and late life.

METHODS

Using an enlarged set of female-specific genetic variants associated with adult general [body mass index (BMI)] and abdominal obesity [waist-to-hip ratio (WHR) with and without adjustment for BMI, WHR and WHRadjBMI] as well as using sex-combined genetic variants of childhood obesity (childhood BMI), we performed a two-sample univariable MR to re-evaluate the total effect of each obesity-related exposure on overall breast cancer (Ncase = 133 384, Ncontrol = 113 789). We further looked into its oestrogen receptor (ER)-defined subtypes (NER+ = 69 501, NER- = 21 468, Ncontrol = 105 974). Multivariable MR was applied to estimate the independent causal effect of each obesity-related exposure on breast cancer taking into account confounders as well as to investigate the independent effect of adult and childhood obesity considering their inter-correlation.

RESULTS

In univariable MR, the protective effects of both adult BMI [odds ratio (OR) = 0.89, 95% CI = 0.83-0.96, P = 2.06 × 10-3] and childhood BMI (OR = 0.78, 95% CI = 0.70-0.87, P = 4.58 × 10-6) were observed for breast cancer overall. Comparable effects were found in ER+ and ER- subtypes. Similarly, genetically predicted adult WHR was also associated with a decreased risk of breast cancer overall (OR = 0.87, 95% CI = 0.80-0.96, P = 3.77 × 10-3), restricting to ER+ subtype (OR = 0.88, 95% CI = 0.80-0.98, P = 1.84 × 10-2). Conditional on childhood BMI, the effect of adult general obesity on breast cancer overall attenuated to null (BMI: OR = 1.00, 95% CI = 0.90-1.10, P = 0.96), whereas the effect of adult abdominal obesity attenuated to some extent (WHR: OR = 0.90, 95% CI = 0.82-0.98, P = 1.49 × 10-2; WHRadjBMI: OR = 0.92, 95% CI = 0.86-0.99, P = 1.98 × 10-2). On the contrary, an independent protective effect of childhood BMI was observed in breast cancer overall, irrespective of adult measures (adjusted for adult BMI: OR = 0.84, 95% CI = 0.77-0.93, P = 3.93 × 10-4; adjusted for adult WHR: OR = 0.84, 95% CI = 0.76-0.91, P = 6.57 × 10-5; adjusted for adult WHRadjBMI: OR = 0.80, 95% CI = 0.74-0.87, P = 1.24 × 10-7).

CONCLUSION

Although successfully replicating the inverse causal relationship between adult obesity-related exposures and risk of breast cancer, our study demonstrated such effects to be largely (adult BMI) or partly (adult WHR or WHRadjBMI) attributed to childhood obesity. Our findings highlighted an independent role of childhood obesity in affecting the risk of breast cancer as well as the importance of taking into account the complex interplay underlying correlated exposures.

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.

Genetic architecture of self-limited delayed puberty and congenital hypogonadotropic hypogonadism

Distinguishing between self limited delayed puberty (SLDP) and congenital hypogonadotropic hypogonadism (CHH) may be tricky as they share clinical and biochemical characteristics. and appear to lie within the same clinical spectrum. However, one is classically transient (SDLP) while the second is typically a lifetime condition (CHH). The natural history and long-term outcomes of these two conditions differ significantly and thus command distinctive approaches and management. Because the first presentation of SDLP and CHH is very similar (delayed puberty with low LH and FSH and low sex hormones), the scientific community is scrambling to identify diagnostic tests that can allow a correct differential diagnosis among these two conditions, without having to rely on the presence or absence of phenotypic red flags for CHH that clinicians anyway seem to find hard to process. Despite the heterogeneity of genetic defects so far reported in DP, genetic analysis through next-generation sequencing technology (NGS) had the potential to contribute to the differential diagnostic process between SLDP and CHH. In this review we will provide an up-to-date overview of the genetic architecture of these two conditions and debate the benefits and the bias of performing genetic analysis seeking to effectively differentiate between these two conditions.

Analysis of transcriptional changes in the immune system associated with pubertal development in a longitudinal cohort of children with asthma

Puberty is an important developmental period marked by hormonal, metabolic and immune changes. Puberty also marks a shift in sex differences in susceptibility to asthma. Yet, little is known about the gene expression changes in immune cells that occur during pubertal development. Here we assess pubertal development and leukocyte gene expression in a longitudinal cohort of 251 children with asthma. We identify substantial gene expression changes associated with age and pubertal development. Gene expression changes between pre- and post-menarcheal females suggest a shift from predominantly innate to adaptive immunity. We show that genetic effects on gene expression change dynamically during pubertal development. Gene expression changes during puberty are correlated with gene expression changes associated with asthma and may explain sex differences in prevalence. Our results show that molecular data used to study the genetics of early onset diseases should consider pubertal development as an important factor that modifies the transcriptome.

The Circadian Clock, Nutritional Signals and Reproduction: A Close Relationship

The circadian rhythm, which is necessary for reproduction, is controlled by clock genes. In the mouse uterus, the oscillation of the circadian clock gene has been observed. The transcription of the core clock gene period (Per) and cryptochrome (Cry) is activated by the heterodimer of the transcription factor circadian locomotor output cycles kaput (Clock) and brain and muscle Arnt-like protein-1 (Bmal1). By binding to E-box sequences in the promoters of Per1/2 and Cry1/2 genes, the CLOCK-BMAL1 heterodimer promotes the transcription of these genes. Per1/2 and Cry1/2 form a complex with the Clock/Bmal1 heterodimer and inactivate its transcriptional activities. Endometrial BMAL1 expression levels are lower in human recurrent-miscarriage sufferers. Additionally, it was shown that the presence of BMAL1-depleted decidual cells prevents trophoblast invasion, highlighting the importance of the endometrial clock throughout pregnancy. It is widely known that hormone synthesis is disturbed and sterility develops in Bmal1-deficient mice. Recently, we discovered that animals with uterus-specific Bmal1 loss also had poor placental development, and these mice also had intrauterine fetal death. Furthermore, it was shown that time-restricted feeding controlled the uterine clock's circadian rhythm. The uterine clock system may be a possibility for pregnancy complications, according to these results. We summarize the most recent research on the close connection between the circadian clock and reproduction in this review.

AMH concentrations in infancy and mid-childhood predict ovarian activity in adolescence: A long-term longitudinal study of healthy girls

BACKGROUND

Anti-Müllerian hormone (AMH) is produced by granulosa cells in small growing ovarian follicles. In adult women, serum concentrations of AMH reflect the ovarian reserve of resting primordial follicles, and low AMH is associated with risk of early menopause. In contrast, patients with polycystic ovary syndrome (PCOS) have elevated AMH. The primary aim of this study was to evaluate the individual tracking of serum AMH concentrations, as well as whether AMH in early childhood reflects ovarian activity in adolescence.

METHODS

In this large longitudinal study of healthy girls were examined from infancy to adolescence (1997-2019) including physical examination, assessment of serum concentrations of reproductive hormones (in infancy, median age 0.3 yrs; mid-childhood, 7.2 yrs; puberty, 11.3 yrs; and adolescence, 15.9 yrs), transabdominal ultrasound (TAUS, puberty and adolescence) and magnetic resonance imaging (MRI, puberty) of the ovaries.

FINDINGS

Each girl maintained her relative AMH concentration (expressed as standard deviation (SD) scores) over time; mean variation of individual age adjusted AMH concentrations was 0.56 ± 0.31 SD.Serum concentrations of AMH in adolescence correlated with AMH in infancy and childhood; infancy: r = 0.347; mid-childhood: r = 0.637; puberty: r = 0.675, all p < 0.001.AMH correlated negatively with FSH concentrations in all age groups (infancy: r = -0.645, p < 0.001; mid-childhood: r = -0.222, p < 0.001; puberty: r = -0.354, p < 0.001; adolescence: n = 275, r = -0.175, p = 0.004).Serum AMH concentrations in mid-childhood correlated with the number of follicles in puberty (TAUS and MRI) as well as in adolescence (TAUS); e.g. total number of follicles: TAUS puberty (r = 0.607), MRI puberty (r = 0.379), TAUS adolescence (r = 0.414), all p < 0.001.AMH concentration in infancy as well as in mid-childhood predicted low AMH (<10 pmol/L) in adolescence; AMH infancy <7.5 pmol/L as predictor of low AMH in adolescence: sensitivity 0.71, specificity 0.70, AUC 0.759; AMH mid-childhood < 8.4 pmol/L as predictor of low AMH in adolescence: sensitivity 0.88, specificity 0.87, AUC 0.949.Girls with high serum AMH concentration in mid-childhood (AMH >30.0 pmol/L vs. other girls) had higher adolescent LH (median 4.53 vs. 3.29 U/L p = 0.041), LH/FSH ratio (1.00 vs 0.67, p = 0.019), testosterone (1.05 vs 0.81 nmol/L, p = 0.005), total number of follicles (23 vs. 19, p = 0.004), and higher prevalence of irregular cycles (10/15 = 67% vs. 28/113 = 25%, p = 0.002).

INTERPRETATION

The present findings suggest remarkably stable ovarian activity from small growing follicles in healthy girls, supporting AMH in early life as a useful clinical tool to predict future ovarian activity.

FUNDING

The work was supported by The Center on Endocrine Disruptors (CeHoS) under The Danish Environmental Protection Agency and The Ministry of Environment and Food (grant number: MST-621-00 065), the EU (QLK4-CT1999-01422; QLK4-2001-00269), the Novo Nordisk Foundation and The Danish Ministry of Science Technology and Innovation (2107-05-0006). A.S.B. is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - 464240267. KM receives honoraria from Novo Nordisk A/S for teaching at the Danish annual postgraduate course of pituitary diseases.

Mammalian puberty: a fly perspective

Mammalian puberty and Drosophila metamorphosis, despite their evolutionary distance, exhibit similar design principles and conservation of molecular components. In this Viewpoint, we review recent advances in this area and the similarities between both processes in terms of the signaling pathways and neuroendocrine circuits involved. We argue that the detection and uptake of peripheral fat by Drosophila prothoracic endocrine cells induces endomembrane remodeling and ribosomal maturation, leading to the acquisition of high biosynthetic and secretory capacity. The absence of this fat-neuroendocrine interorgan communication leads to giant, obese, non-pupating larvae. Importantly, human leptin is capable of signaling the pupariation process in Drosophila, and its expression prevents obesity and triggers maturation in mutants that do not pupate. This implies that insect metamorphosis can be used to address issues related to the biology of leptin and puberty.

Precocious puberty under stressful conditions: new understanding and insights from the lessons learnt from international adoptions and the COVID-19 pandemic

Neuro-biological variations in the timing of sexual maturation within a species are part of an evolved strategy that depend on internal and external environmental conditions. An increased incidence of central precocious puberty (CPP) has been described in both adopted and "covid-19 pandemic" children. Until recently, it was hypothesised that the triggers for CPP in internationally adopted children were likely to be better nutrition, greater environmental stability, and improved psychological wellbeing. However, following data collected during and after the coronavirus (COVID-19) global pandemic, other possibilities must be considered. In a society with high levels of child wellbeing, the threat to life presented by an unknown and potentially serious disease and the stressful environment created by lockdowns and other public health measures could trigger earlier pubertal maturation as an evolutionary response to favour early reproduction. The main driver for increased rates of precocious and rapidly progressive puberty during the pandemic could have been the environment of "fear and stress" in schools and households. In many children, CPP may have been triggered by the psychological effects of living without normal social contact, using PPE, being near adults concerned about financial and other issues and the fear of getting ill. The features and time of progression of CPP in children during the pandemic are similar to those observed in adopted children. This review considers the mechanisms regulating puberty with a focus on neurobiological and evolutionary mechanisms, and analyses precocious puberty both during the pandemic and in internationally adopted children searching for common yet unconsidered factors in an attempt to identify the factors which may have acted as triggers. In particular, we focus on stress as a potential factor in the early activation of the hypothalamic-pituitary-gonadal axis and its correlation with rapid sexual maturation.

Bring the life stages into the domain of basic and clinical pharmacology

Completely distinct physiological conditions and immune responses exist among different human life stages. Age is not always consistent with the life stage. We proposed to incorporate the concept of the life stages into basic and clinical pharmacology, including clinical trials, drug labels, and drug usage in clinical practice. Life-stage-based medical treatment is the application of medicine according to life stages such as prepuberty, reproductive, and aging. A large number of diseases are life-stage-dependent. Many medications and therapy have shown various age effects but not been recognized as life-stage-dependent. The same dosage and drug applications used in different life stages lead to divergent outcomes. Incorporating life stages in medicine and drug usage will enhance the efficacy and precision of the medication in disease treatment.

Exploring adolescence as a key life history stage in bioarchaeology

Adolescence is a unique period in the life history of an individual. It is characterized by a myriad of changes that bioarchaeologists are only just coming to appreciate, related to sexual maturation, linear growth, immunological transformation, and emotional and cognitive development. New methods allow us to measure this age of transition through the stages of the adolescent growth, as a proxy for the physical development associated with sexual maturation (puberty). This review outlines ways bioarchaeologists may draw on research developments from the fields of human biology, evolutionary theory and neurobiology to advance a more holistic approach to the study of adolescence in the past. It considers current theoretical and analytical approaches to highlight the research potential of this critical stage of life history. This synthesis integrates the most recent research in the medical sciences concerned with body and brain development, and outlines the biological processes involved with sexual and physical maturation of the adolescent. The goal of this review is to help inform potentially rewarding areas of research that bioarchaeologists can contribute to and draw from, as well as the challenges and limitations, theoretical and methodological questions, and ways in which we can develop the study of adolescence in the discipline going forward.

Dissecting the loci underlying maturation timing in Atlantic salmon using haplotype and multi-SNP based association methods

Characterizing the role of different mutational effect sizes in the evolution of fitness-related traits has been a major goal in evolutionary biology for a century. Such characterization in a diversity of systems, both model and non-model, will help to understand the genetic processes underlying fitness variation. However, well-characterized genetic architectures of such traits in wild populations remain uncommon. In this study, we used haplotype-based and multi-SNP Bayesian association methods with sequencing data for 313 individuals from wild populations to test the mutational composition of known candidate regions for sea age at maturation in Atlantic salmon (Salmo salar). We detected an association at five loci out of 116 candidates previously identified in an aquaculture strain with maturation timing in wild Atlantic salmon. We found that at four of these five loci, variation explained by the locus was predominantly driven by a single SNP suggesting the genetic architecture of this trait includes multiple loci with simple, non-clustered alleles and a locus with potentially more complex alleles. This highlights the diversity of genetic architectures that can exist for fitness-related traits. Furthermore, this study provides a useful multi-SNP framework for future work using sequencing data to characterize genetic variation underlying phenotypes in wild populations.

Women's reproductive traits and ischemic stroke: a two-sample Mendelian randomization study

OBJECTIVE

We conducted a Mendelian randomization (MR) study to disentangle causal associations between women's reproductive behaviors and ischemic stroke (IS) and investigate the roles of two modifiable risk factors (body mass index (BMI) and educational attainment (EA)) in these associations.

METHODS

Using summary-level data from large-scale genome-wide association studies, we performed univariable MR to examine whether there is genetic evidence that women's reproductive traits are causally associated with IS and its subtypes. Multivariable MR and MR mediation analysis were used to investigate whether BMI and EA are common mechanisms or mediators for these associations. A set of sensitivity analyses were conducted to test valid MR assumptions.

RESULTS

We observed consistent and statistically significant associations across female and sex-combined analyses for earlier age at first birth (AFB) and age at first sexual intercourse (AFS) with a higher risk of IS and large-artery atherosclerotic stroke (LAS) risk in the primary analysis. The odds ratios of IS per 1 SD increase in genetically predicted early AFB and AFS were 0.93 (95% CI, 0.86-0.99; p = 0.046) and 0.83 (95% CI, 0.70-0.97, p = 0.020), respectively. Further analyses indicated that BMI played a shared role in AFS and IS/LAS while EA played a shared role in AFS/AFB and IS/LAS as well as a mediator in the path from AFS to IS/LAS.

INTERPRETATION

These findings may inform prevention strategies and interventions directed toward relative women's reproductive behaviors and IS. Future studies are warranted to explore other factors related to EA which are responsible for these causalities.

Effects of Elevated Maternal Adiposity on Offspring Reproductive Health: A Perspective From Epidemiologic Studies

One in seven couples in developed countries suffers from infertility. Maternal overweight or obesity have detrimental and lasting effects on offspring cardiometabolic health, and although substantially more data are needed, hormonal imbalances in utero resulting from excessive maternal adiposity could also disrupt reproductive programming and affect the future reproductive health of offspring. Therefore, this mini-review evaluates the human epidemiologic evidence that maternal overweight/obesity could be associated with poor reproductive health outcomes in offspring. We searched PubMed for relevant studies using terms such as "maternal obesity" and "reproductive development." While the human epidemiologic literature is limited, studies have thus far observed that maternal obesity is associated with disrupted external genital development and several other markers of reproductive health across the lifespan. Specifically, maternal obesity is associated with higher risks of hypospadias and cryptorchidism in males and disrupted anogenital distance both in males and females. Maternal obesity has also been linked to earlier age at menarche in daughters, and precocious puberty in both sons and daughters. Finally, daughters of women with overweight or obesity have higher risks of developing polycystic ovarian syndrome, which has implications for fertility. This body of research suggests that in utero exposure to maternal obesity could disrupt reproductive system development, but substantially more evidence is needed, as almost no human epidemiologic studies have evaluated the long-term consequences of maternal obesity with regard to offspring fertility/fecundity.

Meta-analysis of genome-wide association studies of hoarding symptoms in 27,537 individuals

Hoarding Disorder (HD) is a mental disorder characterized by persistent difficulties discarding or parting with possessions, often resulting in cluttered living spaces, distress, and impairment. Its etiology is largely unknown, but twin studies suggest that it is moderately heritable. In this study, we pooled phenotypic and genomic data from seven international cohorts (N = 27,537 individuals) and conducted a genome wide association study (GWAS) meta-analysis of parent- or self-reported hoarding symptoms (HS). We followed up the results with gene-based and gene-set analyses, as well as leave-one-out HS polygenic risk score (PRS) analyses. To examine a possible genetic association between hoarding symptoms and other phenotypes we conducted cross-trait PRS analyses. Though we did not report any genome-wide significant SNPs, we report heritability estimates for the twin-cohorts between 26-48%, and a SNP-heritability of 11% for an unrelated sub-cohort. Cross-trait PRS analyses showed that the genetic risk for schizophrenia and autism spectrum disorder were significantly associated with hoarding symptoms. We also found suggestive evidence for an association with educational attainment. There were no significant associations with other phenotypes previously linked to HD, such as obsessive-compulsive disorder, depression, anxiety, or attention-deficit hyperactivity disorder. To conclude, we found that HS are heritable, confirming and extending previous twin studies but we had limited power to detect any genome-wide significant loci. Much larger samples will be needed to further extend these findings and reach a "gene discovery zone". To move the field forward, future research should not only include genetic analyses of quantitative hoarding traits in larger samples, but also in samples of individuals meeting strict diagnostic criteria for HD, and more ethnically diverse samples.

Using genetic instruments to estimate the causal effect of hormonal reproductive factors on osteoarthritis

Objectives

Hormonal reproductive factors have been considered to play an important role in the etiology of osteoarthritis (OA). We performed Mendelian randomization (MR) to examine whether a causal effect existed between them.

Methods

MR was performed by using publicly released genome-wide association study (GWAS) summary statistics to estimate the causal associations of three relevant exposures, including age at menarche (AAM), age at natural menopause (ANM) and age at first birth (AFB), with the risk of OA. We employed several MR methods, including inverse-variance weighted (IVW), MR-Egger regression, weighted median and weighted mode, to estimate the causality. We performed a sensitivity analysis by manually pruning pleiotropic variants associated with the known confounder body mass index (BMI).

Results

The instrumental variables that achieved genome-wide significance, including 349 AAM single nucleotide polymorphisms (SNPs), 121 AAM SNPs, 54 ANM SNPs, and 10 AFB SNPs, were incorporated into the operation. IVW analysis indicated that each additional year in AFB was associated with a decreasing risk of hip and/or knee OA and overall OA (hip and/or knee OA: OR = 0.79, 95% CI: 0.64-0.93, P = 1.33 × 10^-3; overall OA: OR = 0.80, 95% CI: 0.68-0.92, P = 1.80 × 10^-4). In addition, our results suggested that AAM exerted a causal effect on knee OA in an unfavorable manner (OR = 0.86, 95% CI: 0.76-0.95, P = 1.58 × 10^-3). After accounting for the effect of BMI, the causal effect association between AFB and hip and/or knee OA was also examined (IVW: OR = 0.78, 95% CI: 0.66-0.92, P = 3.22 × 10^-3).

Conclusion

Our findings add a growing body of evidence surrounding the unfavorable effects of early AFB on OA risk, suggesting the essential for relevant health problem management in susceptible populations.

Sex Hormone Candidate Gene Polymorphisms Are Associated with Endometriosis

The present study was designed to examine whether sex hormone polymorphisms proven by GWAS are associated with endometriosis risk. Unrelated female participants totaling 1376 in number (395 endometriosis patients and 981 controls) were recruited into the study. Nine single-nucleotide polymorphisms (SNPs) which GWAS correlated with circulating levels of sex hormones were genotyped using a TaqMan allelic discrimination assay. FSH-lowering, and LH- and testosterone-heightening polymorphisms of the FSHB promoter (allelic variants A rs11031002 and C rs11031005) exhibit a protective effect for endometriosis (OR = 0.60-0.68). By contrast, the TT haplotype loci that were GWAS correlated with higher FSH levels and lower LH and testosterone concentrations determined an increased risk for endometriosis (OR = 2.03). Endometriosis-involved epistatic interactions were found between eight loci of sex hormone genes (without rs148982377 ZNF789) within twelve genetic simulation models. In silico examination established that 8 disorder-related loci and 80 proxy SNPs are genome variants affecting the expression, splicing, epigenetic and amino acid conformation of the 34 genes which enrich the organic anion transport and secondary carrier transporter pathways. In conclusion, the present study showed that sex hormone polymorphisms proven by GWAS are associated with endometriosis risk and involved in the molecular pathophysiology of the disease due to their functionality.

Genetics of irritable bowel syndrome: shifting gear via biobank-scale studies

The pathophysiology of irritable bowel syndrome (IBS) is multifactorial and probably involves genetic predisposition and the effect of environmental factors. Unlike other gastrointestinal diseases with a heritable component, genetic research in IBS has been scarce and mostly characterized by small underpowered studies, leading to inconclusive results. The availability of genomic and health-related data from large international cohorts and population-based biobanks offers unprecedented opportunities for long-awaited, well-powered genetic studies in IBS. This Review focuses on the latest advances that provide compelling evidence for the importance of genes involved in the digestion of carbohydrates, ion channel function, neurotransmitters and their receptors, neuronal pathways and the control of gut motility. These discoveries have generated novel information that might be further refined for the identification of predisposed individuals and selection of management strategies for patients. This Review presents a conceptual framework, the advantages and potential limitations of modern genetic research in IBS, and a summary of available evidence.

Mendelian randomization analysis of factors related to ovulation and reproductive function and endometrial cancer risk

BACKGROUND

Observational epidemiological studies suggest a link between several factors related to ovulation and reproductive function and endometrial cancer (EC) risk; however, it is not clear whether these relationships are causal, and whether the risk factors act independently of each other. The aim of this study was to investigate putative causal relationships between the number of live births, age at last live birth, and years ovulating and EC risk.  METHODS: We conducted a series of observational analyses to investigate various risk factors and EC risk in the UK Biobank (UKBB). Additionally, multivariate analysis was performed to elucidate the relationship between the number of live births, age at last live birth, and years ovulating and other related factors such as age at natural menopause, age at menarche, and body mass index (BMI). Secondly, we used Mendelian randomization (MR) to assess if these observed relationships were causal. Genome-wide significant single nucleotide polymorphisms (SNPs) were extracted from previous studies of woman's number of live births, age at menopause and menarche, and BMI. We conducted a genome-wide association analysis using the UKBB to identify SNPs associated with years ovulating, years using the contraceptive pill, and age at last live birth.

RESULTS

We found evidence for a causal effect of the number of live births (inverse variance weighted (IVW) odds ratio (OR): 0.537, p = 0.006), the number of years ovulating (IVW OR: 1.051, p = 0.014), in addition to the known risk factors BMI, age at menarche, and age at menopause on EC risk in the univariate MR analyses. Due to the close relationships between these factors, we followed up with multivariable MR (MVMR) analysis. Results from the MVMR analysis showed that number of live births had a causal effect on EC risk (OR: 0.783, p = 0.036) independent of BMI, age at menarche and age at menopause.

CONCLUSIONS

MVMR analysis showed that the number of live births causally reduced the risk of EC.

[Associations between early onset of puberty and obesity types in children: Based on both the cross-sectional study and cohort study]

OBJECTIVE

To explore and analyze the relationship between early onset of puberty and different types of obesity in children, by combining large sample cross-sectional survey data with long-term longitudinal cohort data, so as to provide clues for further clarifying the health hazards of early onset of puberty and obesity prevention and control.

METHODS

The research data were from the cross-sectional survey data of seven provinces(autonomous regions, municipalities) in China and the cohort data of adolescent development in Xiamen. The study first found the association between early onset of puberty and obesity by Logistic regression on the cross-sectional data, and then used Poisson regression to analyze the association between early puberty initiation and various types of obesity risk.

RESULTS

In the study, 43 137 and 1 266 children were included in the cross-sectional survey and cohort survey respectively. The cross-sectional study found that among the girls aged 10-13 years, compared with the girls of the same age who did not start puberty, the body mass index (BMI)-Z score of the girls in the puberty start group was 0.5-0.8 higher, and the waist circumference Z score was 0.4-0.7 higher, and the risk of various types of obesity was higher. At the same time, the early onset of puberty was positively correlated with simple obesity, central obesity and compound obesity, the OR (95%CI) were 1.86 (1.42-2.44), 1.95 (1.65-2.32) and 1.86 (1.41-2.45), respectively. No significant association was found in boys. According to the cohort data, in girls, the risk of simple obesity was 6.00 times [RR (95%CI): 6.00 (1.07-33.60)], the risk of central obesity was 3.30 times [RR (95%CI): 3.30 (1.22-8.92)], and the risk of compound obesity was 5.76 times [RR (95%CI): 5.76 (1.03-32.30)], compared with the group without early puberty initiation, while no association between early puberty initiation and obesity was found in boys.

CONCLUSION

Based on the cross-sectional survey and longitudinal cohort survey, it is confirmed that the early onset of puberty in girls may increase the risk of simple obesity, central obesity and compound obesity, while there is no significant correlation between puberty onset and obesity in boys.

Common Genetic Variation and Age of Onset of Anorexia Nervosa

Background

Genetics and biology may influence the age of onset of anorexia nervosa (AN). The aims of this study were to determine whether common genetic variation contributes to age of onset of AN and to investigate the genetic associations between age of onset of AN and age at menarche.

Methods

A secondary analysis of the Psychiatric Genomics Consortium genome-wide association study (GWAS) of AN was performed, which included 9335 cases and 31,981 screened controls, all from European ancestries. We conducted GWASs of age of onset, early-onset AN (<13 years), and typical-onset AN, and genetic correlation, genetic risk score, and Mendelian randomization analyses.

Results

Two loci were genome-wide significant in the typical-onset AN GWAS. Heritability estimates (single nucleotide polymorphism-h 2) were 0.01-0.04 for age of onset, 0.16-0.25 for early-onset AN, and 0.17-0.25 for typical-onset AN. Early- and typical-onset AN showed distinct genetic correlation patterns with putative risk factors for AN. Specifically, early-onset AN was significantly genetically correlated with younger age at menarche, and typical-onset AN was significantly negatively genetically correlated with anthropometric traits. Genetic risk scores for age of onset and early-onset AN estimated from independent GWASs significantly predicted age of onset. Mendelian randomization analysis suggested a causal link between younger age at menarche and early-onset AN.

Conclusions

Our results provide evidence consistent with a common variant genetic basis for age of onset and implicate biological pathways regulating menarche and reproduction.

Adverse Effects of Early Puberty Timing in Girls and Potential Solutions

Given the global secular declining trends of the age at puberty and its relevant mechanisms, as illustrated in the first part of this series, the present part will discuss the public health implications of early puberty and potential clinical and public health measures. Although the major effect of earlier maturation impacts adolescents' mental health and likelihood of engaging in risky behaviors, there are also effects in adulthood on cardiometabolic health, especially type 2 diabetes, and an increased risk of certain cancers, especially hormone-related cancers such as breast cancer. The paper ends with recommendations for clinical management, especially for girls who should receive further evaluation, as well as recommendations for the patient and her family and public health considerations.

Genetics of pubertal delay

The timing of pubertal development is strongly influenced by the genetic background, and clinical presentations of delayed puberty are often found within families with clear patterns of inheritance. The discovery of the underlying genetic regulators of such conditions, in recent years through next generation sequencing, has advanced the understanding of the pathogenesis of disorders of pubertal timing and the potential for genetic testing to assist diagnosis for patients with these conditions. This review covers the significant advances in the understanding of the biological mechanisms of delayed puberty that have occurred in the last two decades.

Genome-wide causal mediation analysis identifies genetic loci associated with uterine fibroids mediated by age at menarche

STUDY QUESTION

Could the direct contribution of genetic variants to the pathophysiology of uterine fibroids and the contribution mediated by age at menarche be different?

SUMMARY ANSWER

Age at menarche plays a mediation role in the genetic influence on uterine fibroids, and four causal genetic mechanisms underlying the age at menarche-mediated effects of common genetic loci on uterine fibroid development were identified.

WHAT IS KNOWN ALREADY

Uterine fibroids are common benign tumors developing from uterine smooth muscle. Genome-wide association studies (GWASs) have identified over 30 genetic loci associated with uterine fibroids in different ethnic populations. Several genetic variations in or nearby these identified loci were also associated with early age at menarche, one of the major risk factors of uterine fibroids. Although the results of GWASs reveal how genetic variations affect uterine fibroids, the genetic mechanism of uterine fibroids mediated by age at menarche remains elusive.

STUDY DESIGN, SIZE, DURATION

In this study, we conducted a genome-wide causal mediation analysis in two cohorts covering a total of 69 552 females of Han Chinese descent from the Taiwan Biobank (TWB). TWB is an ongoing community- and hospital-based cohort aiming to enroll 200 000 individuals from the general Taiwanese population between 30 and 70 years old. It has been enrolling Taiwanese study participants since 2012 and has extensive phenotypic data collected from 148 291 individuals as of May 2021.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We recruited individuals in two cohorts, with 13 899 females in TWB1 and 55 653 females in TWB2. The two sets of individuals are almost distinct, with only 730 individuals enrolled in both cohorts. Over 99% of the participants are Han Chinese. Approximately 21% of participants developed uterine fibroids. DNA samples from both cohorts were genotyped using two different customized chips (TWB1 and TWB2 arrays). After quality control and genotype imputation, 646 973 TWB1 single-nucleotide polymorphisms (SNPs) and 686 439 TWB2 SNPs were assessed in our analysis. There were 99 939 SNPs which overlapped between the TWB1 and TWB2 arrays, 547 034 TWB1 array-specific SNPs and 586 500 TWB2 array-specific SNPs. We performed GWASs for screening potential risk SNPs for age at menarche and for uterine fibroids. We subsequently identified causal mediation effects of risk SNPs on uterine fibroids mediated by age at menarche.

MAIN RESULTS AND THE ROLE OF CHANCE

In addition to known loci at LIN28B associated with age at menarche and loci at WNT4 associated with uterine fibroids, we identified 162 SNPs in 77 transcripts that were associated with menarche-mediated causal effects on uterine fibroids via four different causal genetic mechanisms: a both-harmful group with 52 SNPs, a both-protective group with 34 SNPs, a mediator-harmful group with 22 SNPs and a mediator-protective group with 54 SNPs. Among these SNPs, rs809302 in SLK significantly increased the risk of developing uterine fibroids by 3.92% through a mechanism other than age at menarche (P < 10-10), and rs371721345 in HLA-DOB was associated with a 2.70% decreased risk (P < 10-10) in the occurrence of uterine fibroids, mediated by age at menarche. These findings provide insights into the mechanism underlying the effect of genetic loci on uterine fibroids mediated by age at menarche.

LIMITATIONS, REASONS FOR CAUTION

A potential issue is that the present study relied upon self-reported age at menarche and uterine fibroid information. Due to the experimental design, the consistency between self-reports and medical records for uterine fibroids in Taiwan cannot be checked. Fortunately, the literature support that self-reporting even years later remains a practical means for collecting data on menarche and uterine fibroids. We found that the impact of under-reporting of uterine fibroids is less in our study. In addition, the rate of reporting a diagnosis of uterine fibroids was within the rates of medical diagnosis based on national health insurance data. Future work investigating the consistency between self-reports and medical records in Taiwan can remedy this issue.

WIDER IMPLICATIONS OF THE FINDINGS

This study is the first to investigate whether and to what extent age at menarche mediates the causal effects of genetic variants on uterine fibroids by using genome-wide causal mediation analysis. By treating age at menarche as a mediator, this report provides an insight into the genetic risk factors for developing uterine fibroids. Thus, this article represents a step forward in deciphering the role of intermediated risk factors in the genetic mechanism of disease.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the China Medical University, Taiwan (CMU110-ASIA-13 and CMU107-Z-04), the Ministry of Science and Technology, Taiwan (MOST 110-2314-B-039-058) and the International Joint Usage/Research Center, the Institute of Medical Science, the University of Tokyo, Japan (K2104). The authors have no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Genetic control of RNA splicing and its distinct role in complex trait variation

Most genetic variants identified from genome-wide association studies (GWAS) in humans are noncoding, indicating their role in gene regulation. Previous studies have shown considerable links of GWAS signals to expression quantitative trait loci (eQTLs) but the links to other genetic regulatory mechanisms, such as splicing QTLs (sQTLs), are underexplored. Here, we introduce an sQTL mapping method, testing for heterogeneity between isoform-eQTLeffects (THISTLE), with improved power over competing methods. Applying THISTLE together with a complementary sQTL mapping strategy to brain transcriptomic (n = 2,865) and genotype data, we identified 12,794 genes with cis-sQTLs at P < 5 × 10⁻⁸, approximately 61% of which were distinct from eQTLs. Integrating the sQTL data into GWAS for 12 brain-related complex traits (including diseases), we identified 244 genes associated with the traits through cis-sQTLs, approximately 61% of which could not be discovered using the corresponding eQTL data. Our study demonstrates the distinct role of most sQTLs in the genetic regulation of transcription and complex trait variation.

A powerful method for splicing quantitative trait loci (sQTL) mapping, THISTLE, is presented and applied to a collection of 2,865 brain samples. Integration with GWAS identifies 244 genes associated via cis-sQTLs, of which 61% were not identified using expression QTLs.

Mendelian randomization analyses of 23 known and suspected risk factors and biomarkers for breast cancer overall and by molecular subtypes

Many risk factors have been identified for breast cancer. The potential causality for some of them remains uncertain, and few studies have comprehensively investigated these associations by molecular subtypes. We performed a two-sample Mendelian randomization (MR) study to evaluate potential causal associations of 23 known and suspected risk factors and biomarkers with breast cancer risk overall and by molecular subtypes using data from the Breast Cancer Association Consortium. The inverse-variance weighted method was used to estimate odds ratios (OR) and 95% confidence interval (CI) for association of each trait with breast cancer risk. Significant associations with breast cancer risk were found for 15 traits, including age at menarche, age at menopause, body mass index, waist-to-hip ratio, height, physical activity, cigarette smoking, sleep duration, and morning-preference chronotype, and six blood biomarkers (estrogens, insulin-like growth factor-1, sex hormone-binding globulin [SHBG], telomere length, HDL-cholesterol and fasting insulin). Noticeably, an increased circulating SHBG was associated with a reduced risk of estrogen receptor (ER)-positive cancer (OR = 0.83, 95% CI: 0.73-0.94), but an elevated risk of ER-negative (OR = 1.12, 95% CI: 0.93-1.36) and triple negative cancer (OR = 1.19, 95% CI: 0.92-1.54) (P_{heterogeneity}  = 0.01). Fasting insulin was most strongly associated with an increased risk of HER2-negative cancer (OR = 1.94, 95% CI: 1.18-3.20), but a reduced risk of HER2-enriched cancer (OR = 0.46, 95% CI: 0.26-0.81) (P_{heterogeneity}  = 0.006). Results from sensitivity analyses using MR-Egger and MR-PRESSO were generally consistent. Our study provides strong evidence supporting potential causal associations of several risk factors for breast cancer and suggests potential heterogeneous associations of SHBG and fasting insulin levels with subtypes of breast cancer.

Causal Effect of Age at Menarche on the Risk for Depression: Results From a Two-Sample Multivariable Mendelian Randomization Study

A fair number of epidemiological studies suggest that age at menarche (AAM) is associated with depression, but the reported effect sizes are small, and there is evidence of residual confounding. Moreover, previous Mendelian randomization (MR) studies to avoid inferential problems inherent to epidemiological studies have provided mixed findings. To clarify the causal relationship between age at menarche and broadly defined depression risk, we used 360 genome-wide significantly AAM-related single-nucleotide polymorphisms (SNPs) as instrumental variable and data from the latest GWAS for the broadly defined depression risk on 807,553 individuals (246,363 cases and 561,190 controls). Multiple methods to account for heterogeneity of the instrumental variable (penalized weighted median, MR Lasso, and contamination mixture method), systematic and idiosyncratic pleiotropy (MR RAPS), and horizontal pleiotropy (MR PRESSO and multivariable MR using three methods) were used. Body mass index, education attainment, and total white blood count were considered pleiotropic phenotypes in the multivariable MR analysis. In the univariable [inverse-variance weighted (IVW): OR = 0.96, 95% confidence interval = 0.94–0.98, p = 0.0003] and multivariable MR analysis (IVW: OR = 0.96, 95% confidence interval = 0.94–0.99, p = 0.007), there was a significant causal effect of AAM on depression risk. Thus, the present study supports conclusions from previous epidemiological studies implicating AAM in depression without the pitfalls of residual confounding and reverse causation. Considering the adverse consequences of an earlier AAM on mental health, this finding should foster efforts to address risk factors that promote an earlier AAM.

TSABL: Trait Specific Annotation Based Locus predictor

BACKGROUND

The majority of Genome Wide Associate Study (GWAS) loci fall in the non-coding genome, making causal variants difficult to identify and study. We hypothesized that the regulatory features underlying causal variants are biologically specific, identifiable from data, and that the regulatory architecture that influences one trait is distinct compared to biologically unrelated traits.

RESULTS

To better characterize and identify these variants, we used publicly available GWAS loci and genomic annotations to build 17 Trait Specific Annotation Based Locus (TSABL) predictors to identify differences between GWAS loci associated with different phenotypic trait groups. We used a penalized binomial logistic regression model to select trait relevant annotations and tested all models on a holdout set of loci not used for training in any trait. We were able to successfully build models for autoimmune, electrocardiogram, lipid, platelet, red blood cell, and white blood cell trait groups. We used these models both to prioritize variants in existing loci and to identify new genomic regions of interest.

CONCLUSIONS

We found that TSABL models identified biologically relevant regulatory features, and anticipate their future use to enhance the design and interpretation of genetic studies.

Sex Hormones and Risk of Aneurysmal Subarachnoid Hemorrhage: A Mendelian Randomization Study

The risk of aneurysmal subarachnoid hemorrhage (aSAH) is increased in postmenopausal women compared with men of similar age, suggesting a role for sex hormones. We aimed to explore whether sex hormones, and age at menarche/menopause have a causal effect on aSAH risk by conducting a 2-sample MR study (Mendelian randomization).