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.

Genetic contribution of reproductive traits to risk of uterine leiomyomata: a large-scale, genome-wide, cross-trait analysis.

BACKGROUND: Although phenotypic associations between female reproductive characteristics and uterine leiomyomata have long been observed in epidemiologic investigations, the shared genetic architecture underlying these complex phenotypes remains unclear. OBJECTIVE: We aimed to investigate the shared genetic basis, pleiotropic effects, and potential causal relationships underlying reproductive traits (age at menarche, age at natural menopause, and age at first birth) and uterine leiomyomata. STUDY DESIGN: With the use of large-scale, genome-wide association studies conducted among women of European ancestry for age at menarche (n=329,345), age at natural menopause (n=201,323), age at first birth (n=418,758), and uterine leiomyomata (ncases/ncontrols=35,474/267,505), we performed a comprehensive, genome-wide, cross-trait analysis to examine systematically the common genetic influences between reproductive traits and uterine leiomyomata. RESULTS: Significant global genetic correlations were identified between uterine leiomyomata and age at menarche (rg, -0.17; P=3.65×10-10), age at natural menopause (rg, 0.23; P=3.26×10-07), and age at first birth (rg, -0.16; P=1.96×10-06). Thirteen genomic regions were further revealed as contributing significant local correlations (P<.05/2353) to age at natural menopause and uterine leiomyomata. A cross-trait meta-analysis identified 23 shared loci, 3 of which were novel. A transcriptome-wide association study found 15 shared genes that target tissues of the digestive, exo- or endocrine, nervous, and cardiovascular systems. Mendelian randomization suggested causal relationships between a genetically predicted older age at menarche (odds ratio, 0.88; 95% confidence interval, 0.85-0.92; P=1.50×10-10) or older age at first birth (odds ratio, 0.95; 95% confidence interval, 0.90-0.99; P=.02) and a reduced risk for uterine leiomyomata and between a genetically predicted older age at natural menopause and an increased risk for uterine leiomyomata (odds ratio, 1.08; 95% confidence interval, 1.06-1.09; P=2.30×10-27). No causal association in the reverse direction was found. CONCLUSION: Our work highlights that there are substantial shared genetic influences and putative causal links that underlie reproductive traits and uterine leiomyomata. The findings suggest that early identification of female reproductive risk factors may facilitate the initiation of strategies to modify potential uterine leiomyomata risk.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Investigating the shared genetic architecture of uterine leiomyoma and breast cancer: A genome-wide cross-trait analysis.

Little is known regarding the shared genetic architecture or causality underlying the phenotypic association observed for uterine leiomyoma (UL) and breast cancer (BC). Leveraging summary statistics from the hitherto largest genome-wide association study (GWAS) conducted in each trait, we investigated the genetic overlap and causal associations of UL with BC overall, as well as with its subtypes defined by the status of estrogen receptor (ER). We observed a positive genetic correlation between UL and BC overall (rg = 0.09, p = 6.00 × 10-3), which was consistent in ER+ subtype (rg = 0.06, p = 0.01) but not in ER- subtype (rg = 0.06, p = 0.08). Partitioning the whole genome into 1,703 independent regions, local genetic correlation was identified at 22q13.1 for UL with BC overall and with ER+ subtype. Significant genetic correlation was further discovered in 9 out of 14 functional categories, with the highest estimates observed in coding, H3K9ac, and repressed regions. Cross-trait meta-analysis identified 9 novel loci shared between UL and BC. Mendelian randomization demonstrated a significantly increased risk of BC overall (OR = 1.09, 95% CI = 1.01-1.18) and ER+ subtype (OR = 1.09, 95% CI = 1.01-1.17) for genetic liability to UL. No reverse causality was found. Our comprehensive genome-wide cross-trait analysis demonstrates a shared genetic basis, pleiotropic loci, as well as a putative causal relationship between UL and BC, highlighting an intrinsic link underlying these two complex female diseases.

Mid-pass whole genome sequencing enables biomedical genetic studies of diverse populations.

BACKGROUND: Historically, geneticists have relied on genotyping arrays and imputation to study human genetic variation. However, an underrepresentation of diverse populations has resulted in arrays that poorly capture global genetic variation, and a lack of reference panels. This has contributed to deepening global health disparities. Whole genome sequencing (WGS) better captures genetic variation but remains prohibitively expensive. Thus, we explored WGS at "mid-pass" 1-7x coverage. RESULTS: Here, we developed and benchmarked methods for mid-pass sequencing. When applied to a population without an existing genomic reference panel, 4x mid-pass performed consistently well across ethnicities, with high recall (98%) and precision (97.5%). CONCLUSION: Compared to array data imputed into 1000 Genomes, mid-pass performed better across all metrics and identified novel population-specific variants with potential disease relevance. We hope our work will reduce financial barriers for geneticists from underrepresented populations to characterize their genomes prior to biomedical genetic applications.

Development of a Screening Platform to Identify Small Molecules That Modify ELP1 Pre-mRNA Splicing in Familial Dysautonomia.

Familial dysautonomia (FD) is an autonomic and sensory neuropathy caused by a mutation in the splice donor site of intron 20 of the ELP1 gene. Variable skipping of exon 20 leads to a tissue-specific reduction in the level of ELP1 protein. We have shown that the plant cytokinin kinetin is able to increase cellular ELP1 protein levels in vivo and in vitro through correction of ELP1 splicing. Studies in FD patients determined that kinetin is not a practical therapy due to low potency and rapid elimination. To identify molecules with improved potency and efficacy, we developed a cell-based luciferase splicing assay by inserting renilla (Rluc) and firefly (Fluc) luciferase reporters into our previously well-characterized ELP1 minigene construct. Evaluation of the Fluc/Rluc signal ratio enables a fast and accurate way to measure exon 20 inclusion. Further, we developed a secondary assay that measures ELP1 splicing in FD patient-derived fibroblasts. Here we demonstrate the quality and reproducibility of our screening method. Development and implementation of this screening platform has allowed us to efficiently screen for new compounds that robustly and specifically enhance ELP1 pre-mRNA splicing.

Clinical, pathologic, cytogenetic, and molecular profiling in self-identified black women with uterine leiomyomata.

Black women are disproportionately affected by uterine leiomyomata (UL), or fibroids, compared to other racial groups, having a greater lifetime risk of developing UL and an earlier age of diagnosis. In order to elucidate molecular and genetic mechanisms responsible for the increased prevalence and morbidity associated with UL in black women, clinical, pathologic, cytogenetic, and select molecular profiling (MED12 mutation analysis) of 75 self-reported black women undergoing surgical treatment for UL was performed. Our observations are broadly representative of previous cytogenetic studies of UL: karyotypically abnormal tumors were detected in 30.7% of women and 17.4% of analyzed tumors. No notable association was observed between race and increased occurrence of cytogenetic abnormalities that might contribute to any population-specific morbidity or prevalence rate. Our data on MED12 mutation analyses (73.2% of tumors harbored a MED12 mutation) provide additional support for a significant role of MED12 in tumorigenesis. Although the effect of MED12-mediated tumorigenesis appears significant irrespective of race, other genetic events such as the distribution of karyotypic abnormalities appear differently in black women. This case series indicates that presently recognized genetic and molecular characteristics of UL do not appear to explain the increased prevalence and morbidity of UL in black women.

Erratum to: A multi-stage genome-wide association study of uterine fibroids in African Americans.

The article "A multi-stage genome-wide association study of uterine fibroids in African Americans", written by Jacklyn N. Hellwege, was originally published Online First without open access. After publication in volume 136, issue 10, page 1363-1373 the author decided to opt for Open Choice and to make the article an open access publication. Therefore, the copyright of the article has been changed to

A multi-stage genome-wide association study of uterine fibroids in African Americans.

Uterine fibroids are benign tumors of the uterus affecting up to 77% of women by menopause. They are the leading indication for hysterectomy, and account for $34 billion annually in the United States. Race/ethnicity and age are the strongest known risk factors. African American (AA) women have higher prevalence, earlier onset, and larger and more numerous fibroids than European American women. We conducted a multi-stage genome-wide association study (GWAS) of fibroid risk among AA women followed by in silico genetically predicted gene expression profiling of top hits. In Stage 1, cases and controls were confirmed by pelvic imaging, genotyped and imputed to 1000 Genomes. Stage 2 used self-reported fibroid and GWAS data from 23andMe, Inc. and the Black Women's Health Study. Associations with fibroid risk were modeled using logistic regression adjusted for principal components, followed by meta-analysis of results. We observed a significant association among 3399 AA cases and 4764 AA controls at rs739187 (risk-allele frequency = 0.27) in CYTH4 (OR (95% confidence interval) = 1.23 (1.16-1.30), p value = 7.82 × 10-9). Evaluation of the genetic association results with MetaXcan identified lower predicted gene expression of CYTH4 in thyroid tissue as significantly associated with fibroid risk (p value = 5.86 × 10-8). In this first multi-stage GWAS for fibroids among AA women, we identified a novel risk locus for fibroids within CYTH4 that impacts gene expression in thyroid and has potential biological relevance for fibroids.

Genetic Association Studies in Uterine Fibroids: Risk Alleles Presage the Path to Personalized Therapies.

Uterine leiomyoma (UL) is the most common tumor of the female reproductive system. Epidemiological analyses, including familial aggregation, twin studies, and racial discrepancies in disease prevalence and morbidity, indicated genetic factors influence risk for developing UL. Genome-wide association studies (GWASs) are a powerful method for identifying genetic variants that are associated with elevated risk for a common, complex disease. To date, three genome-wide scans for UL have been performed: a GWAS in Japanese women, a genome-wide linkage and association study in women of European decent, and an admixture-based analysis in African American women. Results from each of the three genome-wide scans performed have had varying success in identifying unique loci associated with predisposition to developing UL. Here, we address the evidence for a genetic basis for UL risk, discuss genetic association studies and their results, and identify challenges and future directions for UL GWAS analyses.