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

Polygenic scores for dizygotic twinning: insights into the genetic architecture of female fertility

PURPOSE

Natural dizygotic twinning (DZT) results from hyper-ovulation and is an indicator of female fertility. However, some traits linked to DZ twinning are also associated with infertility. We examined the relationship between DZT and female (in)fertility using recent GWAS findings.

METHODS

We investigated the genetic architecture of DZT and compared polygenic scores (PGS) for DZT between mothers of naturally conceived DZ twin pregnancies and mothers who required fertility treatments (MAR) in the Netherlands Twin Register (NTR) and the Norwegian Mother, Father, and Child Cohort Study (MoBa). We also calculated genetic correlations between DZT and seven fertility related traits.

RESULTS

DZT has a low polygenicity, with only 0.20% of SNPs estimated to have a nonzero effect. The DZT PGS explains 1.6% of variance in DZT liability, and we observe an odds ratio of 2.29 between the first and the tenth PGS deciles. The DZT PGS distinguishes between mothers of naturally conceived pregnancies and mothers who received MAR and is associated with a shorter time to pregnancy in mothers of singletons. The lowest PGSs were observed for mothers who received hormonal ovulation induction, indicating maternal fertility issues. DZT showed genetic correlations with anovulatory infertility (rg =  - 0.698) and PCOS (rg =  - 0.278), and endometriosis (rg = 0.279).

CONCLUSIONS

Female fertility appears to exists on a genetic spectrum, with anovulation/infertility at one end and DZT at the other. Results suggest that the DZT PGS can be of added value to evaluate female fertility and be incorporated in clinical practice in the future.

Twin Births in 42 Sub-Saharan African Countries from 1986 to 2016: Frequency, Trends and Factors of Variation

Since the 1970s, twin birth rates have increased sharply in developed countries. In Africa, where the rate is the highest globally, its evolution and variation are poorly understood. This article aims to estimate the twinning rate in sub-Saharan African (SSA) countries throughout 1986-2016 and analyze its spatial, temporal, and ethnic variations. It also seeks to identify social and demographic factors associated with a high probability of twin births and outline a forecast of the twinning rate. We used data from 174 Demographic and Health Surveys from 42 countries. We supplemented them with the UN World Population Prospects (WPP). The twinning rate was calculated by reporting the number of twin births per thousand total births. We used logistic regression to analyze the factors associated with twin births. We projected the twinning rate based on WPP. The overall SSA twinning rate is 17.4 per 1000, but it has changed very little over time, and we expect it will grow a little between 2015 and 2050, increasing at most from 17.4 per 1000 to 18.4 per 1000. We also show significant differences in the twinning rate in SSA according to mother ethnicity. Most ethnic groups with high twinning rates belong to the large Bantu ethnic family. SSA remains the 'land of twins', with the twinning rate changing slowly. However, specific health policies must target twin births in SSA to address the public health challenges they present.

A Point Mutation in an Otherwise Dispensable Upstream Fshb Enhancer Moderately Impairs Fertility in Female Mice

Follicle-stimulating hormone (FSH) is necessary for fertility in both sexes as a regulator of gametogenesis and hormone synthesis. Humans with loss-of-function mutations within the gene encoding the FSH beta subunit (FSHB) are infertile. Similarly, female Fshb knock-out mice are infertile and fail to ovulate, and males are subfertile. We recently reported the discovery and characterization of an upstream enhancer of FSHB located 26 Kb upstream of the transcriptional start site in humans (-17 Kb in mouse) that also amplifies activin and gonadotropin-releasing hormone induction of FSHB. Notably, the upstream enhancer contains a polymorphic, fertility-associated site in humans, rs11031006 (G/A), and the minor allele (A) increased enhancer activity in vitro as compared to the major allele (G), likely by increasing the affinity of an SF1 binding element. To investigate the role of the novel enhancer and rs11031006 variant in vivo, we created mouse models to assess deletion of the upstream enhancer and the impact of the G>A point mutation at the rs11031006-equivalent base. A full characterization of the -17 Kb enhancer deletion model revealed no apparent differences in fertility or serum FSH/LH levels, nor did a larger deletion that also included an additional putative regulatory element. In contrast, female mice homozygous for the mutated A allele at the rs11031006-equivalent position had fewer litters over a 120-day fertility assay, abnormal estrous cycling at 10 months, and reduced pituitary Lhb transcript abundance. Overall, while the mouse -17 Kb Fshb enhancer is dispensable for fertility, the rs11031006-equivalent G>A mutation results in subfertility in females.

Transcriptomic and Metabolomic Analyses in Monozygotic and Dizygotic Twins

Monozygotic (MZ) and dizygotic (DZ) twins are studied to understand genetic and environmental influences on complex traits, however the mechanisms behind twinning are not completely understood. (Epi)genomic studies identified SNPs associated with DZ twinning and DNA methylation sites with MZ twinning. To find molecular biomarkers of twinning, we compared transcriptomics and metabolomics data from MZ and DZ twins. We analyzed 42,663 RNA transcripts in 1453 MZ twins and 1294 DZ twins from the Netherlands Twin Register (NTR), followed by sex-stratified analyses. The top 5% transcripts with lowest p-values were analyzed for replication in 217 MZ and 158 DZ twins from the older Finnish Twin cohort (FTC). In the NTR, one transcript (PURG) was significantly differentially expressed between MZ and DZ twins; but this did not replicate in FTC. Pathway analyses highlighted the WNT-pathway, previously associated with MZ twinning, and the TGF-B and SMAD pathway, previously associated with DZ twinning. Meta-analysis of 169 serum metabolites in 2797 MZ and 2040 DZ twins from the NTR, FTC and FinnTwin12, showed no metabolomic differences. Overall, we did not find replicable transcript-level expression differences in blood between MZ and DZ twins, but highlighted the TGF-B/SMAD pathway as a potential transcriptional biomarker for DZ twinning.

Polygenic Scores for Dizygotic Twinning: Insights into the Genetic Architecture of Female Fertility

Natural dizygotic twinning (DZT) results from hyper-ovulation and is considered an indicator of female fertility. DZT has low polygenicity, with only 0.20% of SNPs estimated to have a nonzero effect. A polygenic score (PGS) for DZT explains 1.6% of variance in DZT liability and we observe an odds ratio of 2.29 between the 1st and the 10th PGS decile. The PGS is higher in mothers of naturally conceived twins compared to mothers who received fertility (MAR) treatments in both the Netherlands and Norway. The largest differences were observed for mothers who received hormonal ovulation induction, indicating maternal fertility issues. A higher PGS was also linked to shorter time to pregnancy. DZT showed significant negative genetic correlations with anovulatory infertility (rg = -0.698) and PCOS (rg = -0.278), and a positive genetic correlation with endometriosis (rg = 0.279). These findings suggest DZT PGS is an important, yet under-recognized fertility marker.

Evolution of Litter Size: Proximate and Ultimate Mechanisms

Relative reproductive success and failure are the ultimate determinants of Darwinian fitness. As such, reproductive traits and variations therein have an immediate and considerable impact on the evolutionary trajectory of lineages. Historically, significant attention has been paid to the ecological and evolutionary processes (ultimate factors) that shape the diversity and canalization of reproductive traits within groups to better our understanding of organismal diversity and population or species resilience. In contrast, the physiological systems that mediate variation within and among species (i.e., the proximate factors) in reproductive traits remain a significant black box. To date, there is comparatively little information about how proximate mechanisms constrain or promote evolutionary potential in reproductive traits. In this mini-review, we focus on litter size in Eutherian mammals as a trait with relatively well-defined diversity (litter sizes are well-described both within and across species) and for which some genetic determinants have been identified. We discuss both the ultimate and potential proximate determinants of litter size with special attention to the breadth of physiological traits that may act as "toggle" switches for evolution of litter size. We close with a brief discussion of the role that physiological plasticity may play in the evolution of litter size and lay out several forward-looking areas for future research.

Proteome-wide Mendelian randomization and functional studies uncover therapeutic targets for polycystic ovarian syndrome

Polycystic ovarian syndrome (PCOS) is an endocrine syndrome that affects a large portion of women worldwide. This proteogenomic and functional study aimed to uncover candidate therapeutic targets for PCOS. We comprehensively investigated the causal association between circulating proteins and PCOS using two-sample Mendelian randomization analysis. Cis-protein quantitative trait loci were derived from six genome-wide association studies (GWASs) on plasma proteome. Genetic associations with PCOS were obtained from a large-scale GWAS meta-analysis, FinnGen cohort, and UK Biobank. Colocalization analyses were performed to prioritize the causal role of candidate proteins. Protein-protein interaction (PPI) and druggability evaluation assessed the druggability of candidate proteins. We evaluated the enrichment of tier 1 and 2 candidate proteins in individuals with PCOS and a mouse model and explored the potential application of the identified drug target. Genetically predicted levels of 65 proteins exhibited associations with PCOS risk, with 30 proteins showing elevated levels and 35 proteins showing decreased levels linked to higher susceptibility. PPI analyses revealed that FSHB, POSTN, CCN2, and CXCL11 interacted with targets of current PCOS medications. Eighty medications targeting 20 proteins showed their potential for repurposing as therapeutic targets for PCOS. EGLN1 levels were elevated in granulosa cells and the plasma of individuals with PCOS and in the plasma and ovaries of dehydroepiandrosterone (DHEA)-induced PCOS mouse model. As an EGLN1 inhibitor, administration of roxadustat in the PCOS mouse model elucidated the EGLN1-HIF1α-ferroptosis axis in inducing PCOS and validated its therapeutic effect in PCOS. Our study identifies candidate proteins causally associated with PCOS risk and suggests that targeting EGLN1 provides a promising treatment strategy.

Genetics of female and male reproductive traits and their relationship with health, longevity and consequences for offspring

Substantial shifts in reproductive behaviors have recently taken place in many high-income countries including earlier age at menarche, advanced age at childbearing, rising childlessness and a lower number of children. As reproduction shifts to later ages, genetic factors may become increasingly important. Although monogenic genetic effects are known, the genetics underlying human reproductive traits are complex, with both causal effects and statistical bias often confounded by socioeconomic factors. Here, we review genome-wide association studies (GWASs) of 44 reproductive traits of both female and male individuals from 2007 to early 2024, examining reproductive behavior, reproductive lifespan and aging, infertility and hormonal concentration. Using the GWAS Catalog as a basis, from 159 relevant studies, we isolate 37 genes that harbor association signals for four or more reproductive traits, more than half of which are linked to rare Mendelian disorders, including ten genes linked to reproductive-related disorders: FSHB, MCM8, DNAH2, WNT4, ESR1, IGSF1, THRB, BRWD1, CYP19A1 and PTPRF. We also review the relationship of reproductive genetics to related health and behavioral traits, aging and longevity and the effect of parental age on offspring outcomes as well as reflecting on limitations, open questions and challenges in this fast-moving field.

Clustering Identifies Subtypes With Different Phenotypic Characteristics in Women With Polycystic Ovary Syndrome

CONTEXT

Hierarchical clustering (HC) identifies subtypes of polycystic ovary syndrome (PCOS).

OBJECTIVE

This work aimed to identify clinically significant subtypes in a PCOS cohort diagnosed with the Rotterdam criteria and to further characterize the distinct subtypes.

METHODS

Clustering was performed using the variables body mass index (BMI), luteinizing hormone (LH), follicle-stimulating hormone, dehydroepiandrosterone sulfate, sex hormone-binding globulin (SHBG), testosterone, insulin, and glucose. Subtype characterization was performed by analyzing the variables estradiol, androstenedione, dehydroepiandrosterone, cortisol, anti-Müllerian hormone (AMH), total follicle count (TFC), lipid profile, and blood pressure. Study participants were girls and women who attended our university hospital for reproductive endocrinology screening between February 1993 and February 2021. In total, 2502 female participants of European ancestry, aged 13 to 45 years with PCOS (according to the Rotterdam criteria), were included. A subset of these (n = 1067) fulfilled the National Institutes of Health criteria (ovulatory dysfunction and hyperandrogenism). Main outcome measures included the identification of distinct PCOS subtypes using cluster analysis. Additional clinical variables associated with these subtypes were assessed.

RESULTS

Metabolic, reproductive, and background PCOS subtypes were identified. In addition to high LH and SHBG levels, the reproductive subtype had the highest TFC and levels of AMH (all P < .001). In addition to high BMI and insulin levels, the metabolic subtype had higher low-density lipoprotein levels and higher systolic and diastolic blood pressure (all P < .001). The background subtype had lower androstenedione levels and features of the other 2 subtypes.

CONCLUSION

Reproductive and metabolic traits not used for subtyping differed significantly in the subtypes. These findings suggest that the subtypes capture distinct PCOS causal pathways.

Maternal Characteristics in Natural and Medically Assisted Reproduction Dizygotic Twin Pregnancies

Previous studies have shown that mothers of naturally conceived dizygotic (DZ) twins tend to be taller, older, and smoke more than mothers of naturally conceived monozygotic (MZ) twin and mothers of singletons. Here, we investigate whether mothers of naturally conceived DZ twins differ from mothers who conceived their DZ twins after medically assisted reproduction (MAR) in eight maternal traits related to fertility based on observational survey data. We include data from 33,648 mothers from the Netherlands Twin Register (NTR) and 1660 mothers of twins from the Norwegian Mother, Father and Child Cohort Study (MoBA). We contrast mothers of naturally conceived DZ twins with mothers of MAR DZ twins. Next, we further segment the MAR group into mothers who underwent hormonal induction of ovulation but not in vitro fertilization (IVF) and those who IVF twins, comparing them both to each other and against the mothers of naturally conceived DZ twins. Mothers of naturally conceived DZ twins smoke more often, differ in body composition, have a higher maternal age and have more offspring before the twins than mothers of MZ twins. Compared to MAR DZ twin mothers, mothers of naturally conceived DZ twins have fewer miscarriages, lower maternal age and increased height, more offspring and are more often smokers. BMI before the twin pregnancy is similar in both natural and MAR DZ twin mothers. Mothers who received hormonal induction of ovulation (OI) have a lower maternal age, fewer miscarriages, and a higher number of offspring before their twin pregnancy than twin mothers who received IVF and/or intracytoplasmic sperm injection (ICSI) treatments. Our study shows that twin mothers are a heterogenous group and the differences between twin mothers should be taken into account in epidemiological and genetic research that includes twins.

Genetic regulation of ovulation rate and multiple births

Ovulation rate in many mammalian species is controlled to regulate the numbers of offspring and maximise reproductive success. Pathways that regulate ovulation rate still respond to genetic and environmental factors and show considerable variation within and between species. Genetic segregation, positional cloning, and association studies have discovered numerous mutations and genetic risk factors that contribute to this variation. Notable among the discoveries has been the role of mutations in bone morphogenetic protein 15 (BMP15 ), growth differentiation factor 9 (GDF9 ) and bone morphogenetic protein receptor type 1B (BMPR1B ) from the intra-ovarian signalling pathway contributing to the evidence that signalling from the oocyte is the key driver in follicle regulation rather than circulating gonadotrophin concentrations. Multiple variants in different domains of BMP15 and GDF9 result in partial or complete loss of function of the proteins providing insights into their functional roles and differential regulation contributing to species differences in ovulation rate. Early success encouraged many more studies in prolific strains of sheep, cattle and goats providing a valuable catalogue of genetic variants of large effect increasing ovulation rate and litter size. More recently, genetic association studies are beginning to identify genetic risk factors with smaller effects. Most genes implicated are from pathways with defined roles in regulation of the ovarian function. However, some genomic regions suggest regulation by novel genes. Continuing genetic and related functional studies will add further to our understanding of the detailed regulation of ovulation rate and litter size with implications for health and animal production systems.

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

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

Family History of Twinning and Fertility Traits in Nigerian Mothers of Dizygotic Twins

Familial twinning and fertility traits were investigated in Nigerian mothers of dizygotic (DZ) twins (MoDZT; N = 972) and controls (N = 525) who responded to our person-to-person interview, which included questions on pregnancy history and family history of DZ twinning. Controls were defined as women who are not twins themselves and do not have twins in their first-degree relatives. Over 95% of the participants were Yoruba. We found that Nigerian MoDZT had an average of 4.0 (±2.6) pairs of twins among their relatives, and of these, the prevalence of DZ twins was significantly higher than that of monozygotic (MZ) twins (45.9% vs. 25.8%). Controls had an average of 0.5 (±0.4) pairs, and over 95% of the controls had no twins in their relatives. These results suggest genetic influences on DZ twinning in Nigerians. MoDZT were significantly younger in their mean age at first child, and had higher parity than controls, suggesting increased fertility in MoDZT. As compared to mothers with a single set of twins, mothers (N = 130) with multiple sets had significantly more twins among their relatives (5.4 pairs vs. 3.7 pairs) and had their first twins at a younger age (28.4 vs. 30.7 years), indicating that mothers with multiple sets of twins might have higher genetic propensity for twinning associated with earlier age at twin pregnancy. Our findings argue for genomewide association studies for DZ twinning in Nigerians, and may help to develop intervention strategies to overcome infertility/subfertility problems.

New insights into the (epi)genetics of twinning

Spontaneous dizygotic (DZ) twins, i.e. twins conceived without the use of ARTs, run in families and their prevalence varies widely around the globe. In contrast, monozygotic (MZ) twins occur at a constant rate across time and geographical regions and, with some rare exceptions, do not cluster in families. The leading hypothesis for MZ twins, which arise when a zygote splits during preimplantation stages of development, is random occurrence. We have found the first series of genes underlying the liability of being the mother of DZ twins and have shown that being an MZ twin is strongly associated with a stable DNA methylation signature in child and adult somatic tissues. Because identical twins keep this molecular signature across the lifespan, this discovery opens up completely new possibilities for the retrospective diagnosis of whether a person is an MZ twin whose co-twin may have vanished in the early stages of pregnancy. Here, we summarize the gene finding results for mothers of DZ twins based on genetic association studies followed by meta-analysis, and further present the striking epigenetic results for MZ twins.

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 Skłodowska-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.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

Genetic and environmental influence on resting state networks in young male and female adults: a cartographer mapping study

We propose a unique, minimal assumption, approach based on variance analyses (compared with standard approaches) to investigate genetic influence on individual differences on the functional connectivity of the brain using 65 monozygotic and 65 dizygotic healthy young adult twin pairs' low-frequency oscillation resting state functional Magnetic Resonance Imaging (fMRI) data from the Human Connectome Project. Overall, we found high number of genetically-influenced functional (GIF) connections involving posterior to posterior brain regions (occipital/temporal/parietal) implicated in low-level processes such as vision, perception, motion, categorization, dorsal/ventral stream visuospatial, and long-term memory processes, as well as high number across midline brain regions (cingulate) implicated in attentional processes, and emotional responses to pain. We found low number of GIF connections involving anterior to anterior/posterior brain regions (frontofrontal > frontoparietal, frontotemporal, frontooccipital) implicated in high-level processes such as working memory, reasoning, emotional judgment, language, and action planning. We found very low number of GIF connections involving subcortical/noncortical networks such as basal ganglia, thalamus, brainstem, and cerebellum. In terms of sex-specific individual differences, individual differences in males were more genetically influenced while individual differences in females were more environmentally influenced in terms of the interplay of interactions of Task positive networks (brain regions involved in various task-oriented processes and attending to and interacting with environment), extended Default Mode Network (a central brain hub for various processes such as internal monitoring, rumination, and evaluation of self and others), primary sensorimotor systems (vision, audition, somatosensory, and motor systems), and subcortical/noncortical networks. There were >8.5-19.1 times more GIF connections in males than females. These preliminary (young adult cohort-specific) findings suggest that individual differences in the resting state brain may be more genetically influenced in males and more environmentally influenced in females; furthermore, standard approaches may suggest that it is more substantially nonadditive genetics, rather than additive genetics, which contribute to the differences in sex-specific individual differences based on this young adult (male and female) specific cohort. Finally, considering the preliminary cohort-specific results, based on standard approaches, environmental influences on individual differences may be substantially greater than that of genetics, for either sex, frontally and brain-wide. [Correction added on 10 May 2023, after first online publication: added: functional Magnetic Resonance Imaging. Added: individual differences in, twice. Added statement between furthermore … based on standard approaches.].

Identification of common genetic polymorphisms associated with down-regulated gonadotropin levels in an exome-wide association study

OBJECTIVE

To investigate whether common genetic polymorphisms are associated with gonadotropin levels after down-regulation with daily gonadotropin-releasing hormone agonist and whether the polymorphisms of candidate variants influence the ovarian response to exogenous gonadotropins.

DESIGN

Genetic association study.

SETTING

University-affiliated in vitro fertilization center.

PATIENTS

Subjects enrolled in an exploratory exome-wide association study (n = 862), a replication exome-wide association study (n = 86), and a classifier validation study (n = 148) were recruited from September 2016 to October 2018, September 2019 to September 2020, and January 2021 to December 2021, respectively. The included patients were aged ≤40 years and had a basal follicle-stimulating hormone (FSH) ≤12 IU/L.

INTERVENTIONS

All participants received a luteal phase down-regulation long protocol. Genome DNA was extracted from the peripheral blood leukocytes. For the exploratory and replication cohorts, exome sequencing was conducted on a HiSeq 2500 sequencing platform. The multiplex polymerase chain reaction amplification technique and next-generation sequencing also were performed in the exploratory and replication cohorts. For the samples of the validation cohort, Sanger sequencing was performed.

MAIN OUTCOME MEASURES

The primary endpoint was the gonadotropin levels after down-regulation, and the secondary endpoints were hormone levels and follicle diameters during stimulation, the total dose of FSH, duration of FSH stimulation, number of oocytes retrieved, and clinical pregnancy rate.

RESULTS

In the exploratory cohort, we identified that FSHB rs6169 (P=2.71 × 10-24) and its single-nucleotide polymorphisms in high linkage disequilibrium were associated with the down-regulated FSH level. The same locus was confirmed in the replication cohort. Women carrying the C allele of FSHB rs6169 exhibited higher average estradiol level during stimulation (P=6.82 × 10-5), shorter duration of stimulation, and less amount of exogenous FSH (Pduration=0.0002; Pdose=0.0024). In the independent validation set, adding rs6169 genotypes into the prediction model for FSH level after down-regulation enhanced the area under the curve from 0.560 to 0.712 in a logistic regression model, and increased prediction accuracy by 41.05% when a support vector machine classifier was applied.

CONCLUSION

The C allele of FSHB rs6169 is a susceptibility site for the relatively high level of FSH after down-regulation, which may be associated with increased ovarian FSH sensitivity.

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.

Maximizing the value of twin studies in health and behaviour

In the classical twin design, researchers compare trait resemblance in cohorts of identical and non-identical twins to understand how genetic and environmental factors correlate with resemblance in behaviour and other phenotypes. The twin design is also a valuable tool for studying causality, intergenerational transmission, and gene-environment correlation and interaction. Here we review recent developments in twin studies, recent results from twin studies of new phenotypes and recent insights into twinning. We ask whether the results of existing twin studies are representative of the general population and of global diversity, and we conclude that stronger efforts to increase representativeness are needed. We provide an updated overview of twin concordance and discordance for major diseases and mental disorders, which conveys a crucial message: genetic influences are not as deterministic as many believe. This has important implications for public understanding of genetic risk prediction tools, as the accuracy of genetic predictions can never exceed identical twin concordance rates.

Zwillingsschwangerschaften nach natürlicher Konzeption: wie stark ist die erbliche Komponente?

Die Faszination, die von Zwillingen, Zwillingsschwangerschaften und deren Ätiologie ausgeht, ist über die Jahrhunderte unverändert groß geblieben. Beim Menschen sind Zwillinge nach natürlicher Konzeption ein relativ häufiges Ereignis. Aufgrund des steigenden maternalen Alters und der reproduktionsmedizinischen Technologien hat die Anzahl an Zwillingen in der industrialisierten Gesellschaft seit den 1980er-Jahren stetig zugenommen. Dank des besseren Verständnisses und der Fortschritte in der Pränataldiagnostik, klinischen Genetik, Molekulargenetik und Entwicklungsbiologie ist es mittlerweile ersichtlicher geworden, warum und wie Zwillingsschwangerschaften beim Menschen entstehen. Man unterscheidet zwei Arten von Zwillingen: monozygote („eineiige“) und dizygote („zweieiige“) Zwillinge, jede dieser Formen mit unterschiedlichen Inzidenzen, Häufungen und Ursachen. Die Tendenz, dizygote Zwillinge nach natürlicher Konzeption zu empfangen, ist ein komplexes Geschehen, das erheblichen geografischen, ethnischen, saisonalen und familiären Einflüssen unterliegt und eine klare genetische Komponente aufweist. Die Rate an monozygoten Zwillingen ist weltweit relativ konstant, wobei die Ursache noch nicht gänzlich geklärt ist. Wie man jetzt weiß, spielen auch bei dieser Entität erbliche Ursachen eine Rolle.

Deconstructing a Syndrome: Genomic Insights Into PCOS Causal Mechanisms and Classification

Polycystic ovary syndrome (PCOS) is among the most common disorders in women of reproductive age, affecting up to 15% worldwide, depending on the diagnostic criteria. PCOS is characterized by a constellation of interrelated reproductive abnormalities, including disordered gonadotropin secretion, increased androgen production, chronic anovulation, and polycystic ovarian morphology. It is frequently associated with insulin resistance and obesity. These reproductive and metabolic derangements cause major morbidities across the lifespan, including anovulatory infertility and type 2 diabetes (T2D). Despite decades of investigative effort, the etiology of PCOS remains unknown. Familial clustering of PCOS cases has indicated a genetic contribution to PCOS. There are rare Mendelian forms of PCOS associated with extreme phenotypes, but PCOS typically follows a non-Mendelian pattern of inheritance consistent with a complex genetic architecture, analogous to T2D and obesity, that reflects the interaction of susceptibility genes and environmental factors. Genomic studies of PCOS have provided important insights into disease pathways and have indicated that current diagnostic criteria do not capture underlying differences in biology associated with different forms of PCOS. We provide a state-of-the-science review of genetic analyses of PCOS, including an overview of genomic methodologies aimed at a general audience of non-geneticists and clinicians. Applications in PCOS will be discussed, including strengths and limitations of each study. The contributions of environmental factors, including developmental origins, will be reviewed. Insights into the pathogenesis and genetic architecture of PCOS will be summarized. Future directions for PCOS genetic studies will be outlined.

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.

A GWAS in Idiopathic/Unexplained Infertile Men Detects a Genomic Region Determining Follicle-Stimulating Hormone Levels

CONTEXT

Approximately 70% of infertile men are diagnosed with idiopathic (abnormal semen parameters) or unexplained (normozoospermia) infertility, with the common feature of lacking etiologic factors. Follicle-stimulating hormone (FSH) is essential for initiation and maintenance of spermatogenesis. Certain single-nucleotide variations (SNVs; formerly single-nucleotide polymorphisms [SNPs]) (ie, FSHB c.-211G > T, FSHR c.2039A > G) are associated with FSH, testicular volume, and spermatogenesis. It is unknown to what extent other variants are associated with FSH levels and therewith resemble causative factors for infertility.

OBJECTIVE

We aimed to identify further genetic determinants modulating FSH levels in a cohort of men presenting with idiopathic or unexplained infertility.

METHODS

We retrospectively (2010-2018) selected 1900 men with idiopathic/unexplained infertility. In the discovery study (n = 760), a genome-wide association study (GWAS) was performed (Infinium PsychArrays) in association with FSH values (Illumina GenomeStudio, v2.0). Minor allele frequencies (MAFs) were analyzed for the discovery and an independent normozoospermic cohort. In the validation study (n = 1140), TaqMan SNV polymerase chain reaction was conducted for rs11031005 and rs10835638 in association with andrological parameters.

RESULTS

Imputation revealed 9 SNVs in high linkage disequilibrium, with genome-wide significance (P < 4.28e-07) at the FSHB locus 11p.14.1 being associated with FSH. The 9 SNVs accounted for up to a 4.65% variance in FSH level. In the oligozoospermic subgroup, this was increased up to 6.95% and the MAF was enhanced compared to an independent cohort of normozoospermic men. By validation, a significant association for rs11031005/rs10835638 with FSH (P = 4.71e-06/5.55e-07) and FSH/luteinizing hormone ratio (P = 2.08e-12/6.4e-12) was evident.

CONCLUSIONS

This GWAS delineates the polymorphic FSHB genomic region as the main determinant of FSH levels in men with unexplained or idiopathic infertility. Given the essential role of FSH, molecular detection of one of the identified SNVs that causes lowered FSH and therewith decreases spermatogenesis could resolve the idiopathic/unexplained origin by this etiologic factor.

Mothers with higher twinning propensity had lower fertility in pre-industrial Europe

Historically, mothers producing twins gave birth, on average, more often than non-twinners. This observation has been interpreted as twinners having higher intrinsic fertility - a tendency to conceive easily irrespective of age and other factors - which has shaped both hypotheses about why twinning persists and varies across populations, and the design of medical studies on female fertility. Here we show in >20k pre-industrial European mothers that this interpretation results from an ecological fallacy: twinners had more births not due to higher intrinsic fertility, but because mothers that gave birth more accumulated more opportunities to produce twins. Controlling for variation in the exposure to the risk of twinning reveals that mothers with higher twinning propensity - a physiological predisposition to producing twins - had fewer births, and when twin mortality was high, fewer offspring reaching adulthood. Twinning rates may thus be driven by variation in its mortality costs, rather than variation in intrinsic fertility.

Differentially Expressed Circular RNA Profile Signatures Identified in Prolificacy Trait of Yunshang Black Goat Ovary at Estrus Cycle

CircRNAs acting as miRNA sponges play important roles in the growth process of animal individuals. The prolificacy trait of goats is involved in many pathways, however, the variation of circRNA expression profiles in the different phases of the estrus cycle at high and low fecundity groups is still unknown. Here, we analyzed the circRNA profiles of ovarian tissues among high and low fecundity groups in the follicular phase (HF vs LF), high and low fecundity groups in the luteal phase (HL vs LL), and high and low fecundity in the whole estrus cycle (HF vs HL and LF vs LL) using RNA-seq. A total of 283 (114 upregulated and 169 downregulated), 559 (299 upregulated and 260 downregulated), 449 (254 upregulated and 195 downregulated), and 314 (210 upregulated and 104 downregulated) differentially expressed (DE) circRNAs were screened in HF vs LF, HF vs HL, HL vs LL, and LF vs LL groups, respectively. Enrichment analysis suggested that the targeting of DE circRNAs was mainly enriched in oocyte meiosis, the GnRH signaling pathway, and estrogen signaling pathway. After integrating our previous study of miRNA-seq, there were 56 miRNAs that could target to 192 DE circRNAs, including the miR-133 family (including miR-133a-3p and miR-133b), miR-129-3p, and miR-21, which also had important influence on the prolificacy trait of goats. Then, 18 circRNAs with coding potential were obtained by four software predictions, and 9 circRNAs were validated by RT-qPCR. Together, circRNAs play a key role in the prolificacy trait and the transformation of the follicular phase to the luteal phase in the estrus cycle of goats.

The Heritability of Twinning in Seven Large Historic Pedigrees

It is widely recognized that dizygotic twinning (DZT) runs in families, but estimates of heritability from twin and family data are remarkably scarce and vary considerably. Here, we traced seven large, sometimes historical, multigeneration pedigrees from West Africans, fin de siècle French Jews, Canadians (two pedigrees), and the French royal family, in which twin births were recorded. We estimated heritability of twinning (of all types) as zygosity information was not available, diluting the true DZT heritability by a third or so. The estimates in the range 8-20% are remarkably consistent across time (8-19 generations) and ethnicities and also consistent with twin and family estimates.

A multi-level investigation of the genetic relationship between endometriosis and ovarian cancer histotypes

Endometriosis is associated with increased risk of epithelial ovarian cancers (EOCs). Using data from large endometriosis and EOC genome-wide association meta-analyses, we estimate the genetic correlation and evaluate the causal relationship between genetic liability to endometriosis and EOC histotypes, and identify shared susceptibility loci. We estimate a significant genetic correlation (rg) between endometriosis and clear cell (rg = 0.71), endometrioid (rg = 0.48), and high-grade serous (rg = 0.19) ovarian cancer, associations supported by Mendelian randomization analyses. Bivariate meta-analysis identified 28 loci associated with both endometriosis and EOC, including 19 with evidence for a shared underlying association signal. Differences in the shared risk suggest different underlying pathways may contribute to the relationship between endometriosis and the different histotypes. Functional annotation using transcriptomic and epigenomic profiles of relevant tissues/cells highlights several target genes. This comprehensive analysis reveals profound genetic overlap between endometriosis and EOC histotypes with valuable genomic targets for understanding the biological mechanisms linking the diseases.

A major QTL at the LHCGR/FSHR locus for multiple birth in Holstein cattle

BACKGROUND

Twin and multiple births are rare in cattle and have a negative impact on the performance and health of cows and calves. Therefore, selection against multiple birth would be desirable in dairy cattle breeds such as Holstein. We applied different methods to decipher the genetic architecture of this trait using de-regressed breeding values for maternal multiple birth of ~ 2500 Holstein individuals to perform genome-wide association analyses using ~ 600 K imputed single nucleotide polymorphisms (SNPs).

RESULTS

In the population studied, we found no significant genetic trend over time of the estimated breeding values for multiple birth, which indicates that this trait has not been selected for in the past. In addition to several suggestive non-significant quantitative trait loci (QTL) on different chromosomes, we identified a major QTL on chromosome 11 for maternal multiple birth that explains ~ 16% of the total genetic variance. Using a haplotype-based approach, this QTL was fine-mapped to a 70-kb window on chromosome 11 between 31.00 and 31.07 Mb that harbors two functional candidate genes (LHCGR and FSHR). Analysis of whole-genome sequence data by linkage-disequilibrium estimation revealed a regulatory variant in the 5'-region of LHCGR as a possible candidate causal variant for the identified major QTL. Furthermore, the identified haplotype showed significant effects on stillbirth and days to first service.

CONCLUSIONS

QTL detection and subsequent identification of causal variants in livestock species remain challenging in spite of the availability of large-scale genotype and phenotype data. Here, we report for the first time a major QTL for multiple birth in Holstein cattle and provide evidence for a linked variant in the non-coding region of a functional candidate gene. This discovery, which is a first step towards the understanding of the genetic architecture of this polygenic trait, opens the path for future selection against this undesirable trait, and thus contributes to increased animal health and welfare.

Distal Enhancer Potentiates Activin- and GnRH-Induced Transcription of FSHB

FSH is critical for fertility. Transcription of FSHB, the gene encoding the beta subunit, is rate-limiting in FSH production and is regulated by both GnRH and activin. Activin signals through SMAD transcription factors. Although the mechanisms and importance of activin signaling in mouse Fshb transcription are well-established, activin regulation of human FSHB is less well understood. We previously reported a novel enhancer of FSHB that contains a fertility-associated single nucleotide polymorphism (rs10031006) and requires a region resembling a full (8 base-pair) SMAD binding element (SBE). Here, we investigated the role of the putative SBE within the enhancer in activin and GnRH regulation of FSHB. In mouse gonadotrope-derived LβT2 cells, the upstream enhancer potentiated activin induction of both the human and mouse FSHB proximal promoters and conferred activin responsiveness to a minimal promoter. Activin induction of the enhancer required the SBE and was blocked by the inhibitory SMAD7, confirming involvement of the classical SMAD signaling pathway. GnRH induction of FSHB was also potentiated by the enhancer and dependent on the SBE, consistent with known activin/GnRH synergy regulating FSHB transcription. In DNA pull-down, the enhancer SBE bound SMAD4, and chromatin immunoprecipitation demonstrated SMAD4 enrichment at the enhancer in native chromatin. Combined activin/GnRH treatment elevated levels of the active transcriptional histone marker, histone 3 lysine 27 acetylation, at the enhancer. Overall, this study indicates that the enhancer is directly targeted by activin signaling and identifies a novel, evolutionarily conserved mechanism by which activin and GnRH can regulate FSHB transcription.

Searching for female reproductive aging and longevity biomarkers

Female reproductive aging is, in a way, a biological phenomenon that develops along canonical molecular pathways; however, it has particular features. Recent studies revealed complexity of the interconnections between reproductive aging and aging of other systems, and even suggested a cause-effect uncertainty between them. It was also shown that reproductive aging can impact aging processes in an organism at the level of cells, tissues, organs, and systems. Women at the end of their reproductive lives are characterized by the accelerated incidence of age-related diseases. Timing of the onset of menarche and menopause and variability in the duration of reproductive life carry a latent social risk: not having enough information about the reproductive potential, women keep on postponing childbirth. Identification and use of the most accurate and sensitive aging biomarkers enable the prediction of menopause timing and quantification of the true biological and reproductive ages of an organism. We discuss current views on reproductive aging and peculiarities of using available biomarkers of aging. We also consider latest advances in the search for potential genetic markers of reproductive aging. Finally, we posit the importance of determining the female biological age and highlight potential research directions in this area.

Genetic and in silico analysis show a role of SMAD3 on recurrent pregnancy loss

Recurrent pregnancy loss (RPL) is one of the most common reproductive failures affecting 1-5% of couples. Smad3 is an effector of signalling of the Transforming Growth Factors-β superfamily (TGF-β), regulating the transcription of several target genes of these cytokines. The objective of this study was to evaluate the influence of a variant on SMAD3 (rs17293443) in RPL. A case-control study was carried out with 149 women who experienced RPL and 159 controls, as well as bioinformatics tools to determine the role of this variant in this condition. Our study showed an allelic (p = 0.023) and genotypic (p < 0.01) association of this variant with the RPL. Our functional in silico predictions suggest that this variant causes a change in SMAD3 expression levels. Alterations in the expression of this gene can directly compromise the Smad3-dependent signalling pathway that is fundamental for key processes for gestation such as steroid hormone regulation and implantation, as demonstrated by ontologies analyses performed and the literature. Our findings regarding the involvement of Smad3 on RPL are a novelty in this field, and they seem to be promising to the clinical management of this condition.

Do Males Affect Twinning Events? A Review of Current Findings/Twin Research Reviews: Monozygotic Twins Discordant for Parkinson's Disease; Fetal Loss in Twin Pregnancies Following Prenatal Diagnosis; Uterine Rupture and Repair in an Early Twin Pregnancy; Twin Study of Affectionate Communication/Human Interest: Conjoined Twins in a Triplet Set; Identical Twin Nurses Deliver Identical Twins; Identical Twins Discordant for COVID-19 Recovery Course; Identical Twins Pass Away from COVID-19; Archeological Finds of Oldest Identical Twins

Research into the origins of twinning has focused mostly on contributions from the female side of the family. A review of current findings suggests that possible male contributions to twinning events have been overlooked. This section is followed by brief reviews of twin research concerning monozygotic twins discordant for Parkinson's disease, fetal loss in twin pregnancies following prenatal diagnosis, uterine rupture and repair in an early twin pregnancy and a twin study of affectionate communication. The concluding portion of this article presents human interest stories involving twins that are both informative and poignant, namely conjoined twins in a triplet set, identical twin nurses who delivered identical twins, identical twins discordant for COVID-19 recovery course, identical twins who passed away from COVID-19 and archeological finds of the oldest identical twins.

Colocalization analysis of polycystic ovary syndrome to identify potential disease-mediating genes and proteins

Polycystic ovary syndrome (PCOS) is a common complex disease in women with a strong genetic component and downstream consequences for reproductive, metabolic and psychological health. There are currently 19 known PCOS risk loci, primarily identified in women of Han Chinese or European ancestry, and 14 of these risk loci were identified or replicated in a genome-wide association study of PCOS performed in up to 10,074 cases and 103,164 controls of European descent. However, for most of these loci the gene responsible for the association is unknown. We therefore use a Bayesian colocalization approach (Coloc) to highlight genes in PCOS-associated regions that may have a role in mediating the disease risk. We evaluated the posterior probabilities of evidence consistent with shared causal variants between 14 PCOS genetic risk loci and intermediate cellular phenotypes in one protein (N = 3301) and two expression quantitative trait locus datasets (N = 31,684 and N = 80–491). Through these analyses, we identified seven proteins or genes with evidence of a possibly shared causal variant for almost 30% of known PCOS signals, including follicle stimulating hormone and ERBB3, IKZF4, RPS26, SUOX, ZFP36L2, and C8orf49. Several of these potential effector proteins and genes have been implicated in the hypothalamic–pituitary–gonadal signalling pathway and provide an avenue for functional follow-up in order to demonstrate a causal role in PCOS pathophysiology.

Genetic Regulation of Physiological Reproductive Lifespan and Female Fertility

There is substantial genetic variation for common traits associated with reproductive lifespan and for common diseases influencing female fertility. Progress in high-throughput sequencing and genome-wide association studies (GWAS) have transformed our understanding of common genetic risk factors for complex traits and diseases influencing reproductive lifespan and fertility. The data emerging from GWAS demonstrate the utility of genetics to explain epidemiological observations, revealing shared biological pathways linking puberty timing, fertility, reproductive ageing and health outcomes. The observations also identify unique genetic risk factors specific to different reproductive diseases impacting on female fertility. Sequencing in patients with primary ovarian insufficiency (POI) have identified mutations in a large number of genes while GWAS have revealed shared genetic risk factors for POI and ovarian ageing. Studies on age at menopause implicate DNA damage/repair genes with implications for follicle health and ageing. In addition to the discovery of individual genes and pathways, the increasingly powerful studies on common genetic risk factors help interpret the underlying relationships and direction of causation in the regulation of reproductive lifespan, fertility and related traits.

FSHB Transcription is Regulated by a Novel 5' Distal Enhancer With a Fertility-Associated Single Nucleotide Polymorphism

The pituitary gonadotropins, follicle-stimulating hormone (FSH) and luteinizing hormone, signal the gonads to regulate male and female fertility. FSH is critical for female fertility as it regulates oocyte maturation, ovulation, and hormone synthesis. Multiple genome-wide association studies (GWAS) link a 130 Kb locus at 11p14.1, which encompasses the FSH beta-subunit (FSHB) gene, with fertility-related traits that include polycystic ovary syndrome, age of natural menopause, and dizygotic twinning. The most statistically significant single nucleotide polymorphism from several GWAS studies (rs11031006) resides within a highly conserved 450 bp region 26 Kb upstream of the human FSHB gene. Given that sequence conservation suggests an important biological function, we hypothesized that the region could regulate FSHB transcription. In luciferase assays, the conserved region enhanced FSHB transcription and gel shifts identified a binding site for Steroidogenic factor 1 (SF1) contributing to its function. Analysis of mouse pituitary single-cell ATAC-seq demonstrated open chromatin at the conserved region exclusive to a gonadotrope cell-type cluster. Additionally, enhancer-associated histone markers were identified by immunoprecipitation of chromatin from mouse whole pituitary and an immortalized mouse gonadotrope-derived LβT2 cell line at the conserved region. Furthermore, we found that the rs11031006 minor allele upregulated FSHB transcription via increased SF1 binding to the enhancer. All together, these results identify a novel upstream regulator of FSHB transcription and indicate that rs11031006 can modulate FSH levels.

Community perceptions on causes of high dizygotic twinning rate in Igbo-Ora, South-west Nigeria: A qualitative study

BACKGROUND

Dizygotic (DZ, non-identical) twinning rates vary widely across different regions in the world. With a DZ twinning rate of 45 per 1000 live births, Igbo-Ora Community in South-west Nigeria has the highest dizygotic (DZ) twinning rate in the world. Although several postulations exist on the causes of high DZ twinning rates in Igbo-Ora, no study has yet been conclusive on a definite causative agent.

OBJECTIVE

Using qualitative methods, this study explored the perceptions and beliefs of Igbo-Ora residents about the causes of high DZ twinning rates.

METHODS

Focus group discussion sessions and key informant interviews were organized among fathers and mothers of twins, those without twins, and health care providers. Key informant interviews were also held with persons considered to be custodians of culture who may have knowledge relevant to twinning such as traditional rulers, and traditional birth attendants; as well as health care providers, mothers and fathers of twins, and adult twins.

RESULTS

The results showed three factors featuring as the leading perceived causes of twinning in the community. These included twinning being an act of God, hereditary, and being due to certain foods consumed in the community. Contrary to reports that the consumption of a species of yam (Dioscorea rotundata) may be responsible for the DZ twinning in this Community; yam was not prioritized by the respondents as associated with twinning. In contrast, participants repeatedly mentioned the consumption of "ilasa" a soup prepared with okra leaves (Abelmoschus esculenta) with water that is obtained from the community, and "amala" a local delicacy produced from cassava (Manihot esculenta) as the most likely dietary factors responsible for twinning in the community.

CONCLUSION

Since the same foods are consumed in neighboring communities that have lower rates of twinning, we conjecture that nutritional and other environmental factors may produce epigenetic modifications that influence high DZ twinning rates in Igbo-Ora community. We conclude that more directed scientific studies based on these findings are required to further elucidate the etiology of the high rate of DZ twinning in Igbo-Ora.

Epigenetic consequences of hormonal interactions between opposite-sex twin fetuses

Previous studies reported inconsistent evidence about some phenotypic traits of females in human opposite-sex twins (opposite-sex females [OSF]) being distinct from females in same-sex twins (SSF). Comparatively, less evidence showed significant differences between males in OS twins (opposite-sex males [OSM]) and males in same-sex twins (SSM). The twin testosterone transfer hypothesis suggests that prenatal exposure of testosterone in utero may be a possible explanation for the differential traits in OSF; however, the underlying mechanism is unknown. Here, we investigated the potential epigenetic effects of hormone interactions and their correlation to the observed phenotypic traits. In the study, DNA methylomic data from 54 newborn twins and histone modification data (H3K4me3, H3K4me1, H3K27me3, and H3K27ac) from 14 newborn twins, including same-sex females (SSF), OS twins, and same-sex males (SSM) were generated. We found that OSF were clearly distinguishable from SSF by DNA methylome, while OSM were distinguishable from SSM by H3K4me1 and H3K4me3. To be more specific, compared to SSF, OSF showed a stronger correlation to males (OSM and SSM) in genome-wide DNA methylation. Further, the DNA methylomic differences between OSF and SSF were linked to the process involving cognitive functions and nervous system regulation. The differential H3K4me3 between OSM and SSM was linked to immune responses. These findings provide epigenetic evidence for the twin testosterone transfer hypothesis and offer novel insights on how prenatal hormone exposure in utero may be linked to the reported differential traits of OS twins.

Genetic analysis of endometriosis and depression identifies shared loci and implicates causal links with gastric mucosa abnormality

Evidence from observational studies indicates that endometriosis and depression often co-occur. However, conflicting evidence exists, and the etiology as well as biological mechanisms underlying their comorbidity remain unknown. Utilizing genome-wide association study (GWAS) data, we comprehensively assessed the relationship between endometriosis and depression. Single nucleotide polymorphism effect concordance analysis (SECA) found a significant genetic overlap between endometriosis and depression (P_{Fsig-permuted} = 9.99 × 10^-4). Linkage disequilibrium score regression (LDSC) analysis estimated a positive and highly significant genetic correlation between the two traits (r_G = 0.27, P = 8.85 × 10^-27). A meta-analysis of endometriosis and depression GWAS (sample size = 709,111), identified 20 independent genome-wide significant loci (P < 5 × 10^-8), of which eight are novel. Mendelian randomization analysis (MR) suggests a causal effect of depression on endometriosis. Combining gene-based association results across endometriosis and depression GWAS, we identified 22 genes with a genome-wide significant Fisher's combined P value (FCP_gene < 2.75 × 10^-6). Genes with a nominal gene-based association (P_gene < 0.05) were significantly enriched across endometriosis and depression (P_{binomial-test} = 2.90 × 10^-4). Also, genes overlapping the two traits at P_gene < 0.1 (P_{binomial-test} = 1.31 × 10^-5) were significantly enriched for the biological pathways 'cell-cell adhesion', 'inositol phosphate metabolism', 'Hippo-Merlin signaling dysregulation' and 'gastric mucosa abnormality'. These results reveal a shared genetic etiology for endometriosis and depression. Indeed, additional analyses found evidence of a causal association between each of endometriosis and depression and at least one abnormal condition of gastric mucosa. Our study confirms the comorbidity of endometriosis and depression, implicates links with gastric mucosa abnormalities in their causal pathways and reveals potential therapeutic targets for further investigation.

Migraine, Human Genetics and a Passion for Science

This note reflects on my collaborations with Nick Martin and the GenEpi group over the past 20 years. Over the past two decades, our work together has focused on gene mapping and understanding the genetic architecture of a wide range of traits with particular foci on migraine and common baldness. Our migraine research has included latent class and twin analyses cumulating in genome-wide association analyses which had identified 44 (34 new) risk variants for migraine. Leveraging these results through polygenic risk score analyses identified subgroups of patients likely to respond to triptans (an acute migraine drug), providing the first step toward precision medicine in migraine [Kogelman et al. (2019) Neurology Genetics, 5, e364].

An age-dependent ovulatory strategy explains the evolution of dizygotic twinning in humans

Dizygotic twinning, the simultaneous birth of siblings when multiple ova are released, is an evolutionary paradox. Twin-bearing mothers often have elevated fitness, but despite twinning being heritable, twin births occur only at low frequencies in human populations. We resolve this paradox by showing that twinning and non-twinning are not competing strategies; instead, dizygotic twinning is the outcome of an adaptive conditional ovulatory strategy of switching from single to double ovulation with increasing age. This conditional strategy, when coupled with the well-known decline in fertility as women age, maximizes reproductive success and explains the increase and subsequent decrease in the twinning rate with maternal age that is observed across human populations. We show that the most successful ovulatory strategy would be to always double ovulate as an insurance against early fetal loss, but to never bear twins. This finding supports the hypothesis that twinning is a by-product of selection for double ovulation rather than selection for twinning.

The Genetics of Human DZ Twinning

In the course of twin studies whose main focus was elucidation of genetic and environmental factors on behavioral traits, many twin researchers became aware of the strong tendency for dizygotic (DZ) twinning to run in families. Over four decades, Nick Martin and others initiated hormone and ultrasound studies, performed segregation and pedigree analyses, tested candidate genes, carried out linkage projects in sister pairs and formed large collaborations to illuminate the genetics of DZ twinning by genome-wide association studies and meta-analysis. This article summarizes the early work on hormone and genetic studies and describes the meta-analyses that have at last met with success in finding the first genes that predispose to DZ twinning, which also appear to influence many other female reproductive traits.

Demographic and evolutionary trends in ovarian function and aging

BACKGROUND

The human female reproductive lifespan is regulated by the dynamics of ovarian function, which in turn is influenced by several factors: from the basic molecular biological mechanisms governing folliculogenesis, to environmental and lifestyle factors affecting the ovarian reserve between conception and menopause. From a broader point of view, global and regional demographic trends play an additional important role in shaping the female reproductive lifespan, and finally, influences on an evolutionary scale have led to the reproductive senescence that precedes somatic senescence in humans.

OBJECTIVE AND RATIONALE

The narrative review covers reproductive medicine, by integrating the molecular mechanisms of ovarian function and aging with short-term demographic and long-term evolutionary trends.

SEARCH METHODS

PubMed and Google Scholar searches were performed with relevant keywords (menopause, folliculogenesis, reproductive aging, reproductive lifespan and life history theory). The reviewed articles and their references were restricted to those written in English.

OUTCOMES

We discuss and summarize the rapidly accumulating information from large-scale population-based and single-reproductive-cell genomic studies, their constraints and advantages in the context of female reproductive aging as well as their possible evolutionary significance on the life history trajectory from foetal-stage folliculogenesis until cessation of ovarian function in menopause. The relevant environmental and lifestyle factors and demographic trends are also discussed in the framework of predominant evolutionary hypotheses explaining the origin and maintenance of menopause.

WIDER IMPLICATIONS

The high speed at which new data are generated has so far raised more questions than it has provided solid answers and has been paralleled by a lack of satisfactory interpretations of the findings in the context of human life history theory. Therefore, the recent flood of data could offer an unprecedented tool for future research to possibly confirm or rewrite human evolutionary reproductive history, at the same time providing novel grounds for patient counselling and family planning strategies.

Shared Molecular Genetic Mechanisms Underlie Endometriosis and Migraine Comorbidity

Observational epidemiological studies indicate that endometriosis and migraine co-occur within individuals more than expected by chance. However, the aetiology and biological mechanisms underlying their comorbidity remain unknown. Here we examined the relationship between endometriosis and migraine using genome-wide association study (GWAS) data. Single nucleotide polymorphism (SNP) effect concordance analysis found a significant concordance of SNP risk effects across endometriosis and migraine GWAS. Linkage disequilibrium score regression analysis found a positive and highly significant genetic correlation (rG = 0.38, P = 2.30 × 10-25) between endometriosis and migraine. A meta-analysis of endometriosis and migraine GWAS data did not reveal novel genome-wide significant SNPs, and Mendelian randomisation analysis found no evidence for a causal relationship between the two traits. However, gene-based analyses identified two novel loci for migraine. Also, we found significant enrichment of genes nominally associated (Pgene < 0.05) with both traits (Pbinomial-test = 9.83 × 10-6). Combining gene-based p-values across endometriosis and migraine, three genes, two (TRIM32 and SLC35G6) of which are at novel loci, were genome-wide significant. Genes having Pgene < 0.1 for both endometriosis and migraine (Pbinomial-test = 1.85 ×10-°3) were significantly enriched for biological pathways, including interleukin-1 receptor binding, focal adhesion-PI3K-Akt-mTOR-signaling, MAPK and TNF-α signalling. Our findings further confirm the comorbidity of endometriosis and migraine and indicate a non-causal relationship between the two traits, with shared genetically-controlled biological mechanisms underlying the co-occurrence of the two disorders.

The Netherlands Twin Register: Longitudinal Research Based on Twin and Twin-Family Designs

The Netherlands Twin Register (NTR) is a national register in which twins, multiples and their parents, siblings, spouses and other family members participate. Here we describe the NTR resources that were created from more than 30 years of data collections; the development and maintenance of the newly developed database systems, and the possibilities these resources create for future research. Since the early 1980s, the NTR has enrolled around 120,000 twins and a roughly equal number of their relatives. The majority of twin families have participated in survey studies, and subsamples took part in biomaterial collection (e.g., DNA) and dedicated projects, for example, for neuropsychological, biomarker and behavioral traits. The recruitment into the NTR is all inclusive without any restrictions on enrollment. These resources - the longitudinal phenotyping, the extended pedigree structures and the multigeneration genotyping - allow for future twin-family research that will contribute to gene discovery, causality modeling, and studies of genetic and cultural inheritance.

Avera Twin Register Growing Through Online Consenting and Survey Collection

The aim of the Avera Twin Register (ATR) is to establish a prospective longitudinal repository of twins, multiples, siblings and family members’ biological samples to study environmental and genetic influences on health and disease. Also, it is our intention to contribute to international genome-wide association study (GWAS) twin consortia when appropriate sample size is achieved within the ATR. The ATR is young compared with existing registers and continues to collect a longitudinal repository of biological specimens, survey data and health information. Data and biological specimens were originally collected via face-to-face appointments or the postal department and consisted of paper-informed consents and questionnaires. Enrollment of the ATR began on May 18, 2016 and is located in Sioux Falls, South Dakota, a rural and frontier area in the Central United States with a regional population of approximately 880,000. The original target area for the ATR was South Dakota and the four surrounding states: Minnesota, Iowa, North Dakota and Nebraska. The ATR has found a need to expand that area based on twin and multiple siblings who live in various areas surrounding these states. A description of the state of the ATR today and its transition to online data collection and informed consent will be presented. The ATR collects longitudinal data on lifestyle, including diet and activity levels, aging, plus complex traits and diseases. All twins and multiples participating in the ATR are genotyped on the Illumina Global Screening Array and receive zygosity results.

Modelling for the inheritance of multiple births and fertility in endangered equids: Determining risk factors and genetic parameters in donkeys (Equus asinus)

Multiple births or twinning in equids are dangerous, undesirable situations that compromise the life of the dam and resulting offspring. However, embryo vitrification and freezing techniques take advantage of individuals whose multiple ovulations allow flushing more fertilised embryos from the oviduct to be collected, increasing the productivity and profitability of reproductive techniques. Embryo preservation is especially important in highly endangered populations such as certain donkey (Equus asinus) breeds; for which conventional reproductive techniques have previously been deemed inefficient. For instance, becoming an effective alternative to artificial insemination with frozen semen to preserve the individuals' genetic material. The objective of this study was to examine the historical foaling records of Andalusian donkeys to estimate prevalence, risk factors, phenotypic and genetic parameters for multiple births, assessing the cumulative foal number born per animal, maximum foal number per birth and multiple birth number per animal. We designed a Bayesian General Animal Mixed Model with single records considering the 'fixed' effects of birth year, birth season, birth month, sex, farm, location, and husbandry system. Age was considered and included as a linear and quadratic covariate. Gibbs sampling reported heritability estimates ranging from 0.18 ± 0.101 to 0.24 ± 0.078. Genetic and phenotypic correlations ranged from 0.496 ± 0.298 to 0.846 ± 0.152 and 0.206 ± 0.063 to 0.607 ± 0.054, respectively. Predicted breeding values obtained enable the potential selection against/for these traits, offering a new perspective for donkey breeding and conservation.

Genetic Similarity Assessment of Twin-Family Populations by Custom-Designed Genotyping Array

Twin registries often take part in large collaborative projects and are major contributors to genome-wide association (GWA) meta-analysis studies. In this article, we describe genotyping of twin-family populations from Australia, the Midwestern USA (Avera Twin Register), the Netherlands (Netherlands Twin Register), as well as a sample of mothers of twins from Nigeria to assess the extent, if any, of genetic differences between them. Genotyping in all cohorts was done using a custom-designed Illumina Global Screening Array (GSA), optimized to improve imputation quality for population-specific GWA studies. We investigated the degree of genetic similarity between the populations using several measures of population variation with genotype data generated from the GSA. Visualization of principal component analysis (PCA) revealed that the Australian, Dutch and Midwestern American populations exhibit negligible interpopulation stratification when compared to each other, to a reference European population and to globally distant populations. Estimations of fixation indices (FST values) between the Australian, Midwestern American and Netherlands populations suggest minimal genetic differentiation compared to the estimates between each population and a genetically distinct cohort (i.e., samples from Nigeria genotyped on GSA). Thus, results from this study demonstrate that genotype data from the Australian, Dutch and Midwestern American twin-family populations can be reasonably combined for joint-genetic analysis.

Appraising the role of previously reported risk factors in epithelial ovarian cancer risk: A Mendelian randomization analysis

BACKGROUND

Various risk factors have been associated with epithelial ovarian cancer risk in observational epidemiological studies. However, the causal nature of the risk factors reported, and thus their suitability as effective intervention targets, is unclear given the susceptibility of conventional observational designs to residual confounding and reverse causation. Mendelian randomization (MR) uses genetic variants as proxies for risk factors to strengthen causal inference in observational studies. We used MR to evaluate the association of 12 previously reported risk factors (reproductive, anthropometric, clinical, lifestyle, and molecular factors) with risk of invasive epithelial ovarian cancer, invasive epithelial ovarian cancer histotypes, and low malignant potential tumours.

METHODS AND FINDINGS

Genetic instruments to proxy 12 risk factors were constructed by identifying single nucleotide polymorphisms (SNPs) that were robustly (P < 5 × 10-8) and independently associated with each respective risk factor in previously reported genome-wide association studies. These risk factors included genetic liability to 3 factors (endometriosis, polycystic ovary syndrome, type 2 diabetes) scaled to reflect a 50% higher odds liability to disease. We obtained summary statistics for the association of these SNPs with risk of overall and histotype-specific invasive epithelial ovarian cancer (22,406 cases; 40,941 controls) and low malignant potential tumours (3,103 cases; 40,941 controls) from the Ovarian Cancer Association Consortium (OCAC). The OCAC dataset comprises 63 genotyping project/case-control sets with participants of European ancestry recruited from 14 countries (US, Australia, Belarus, Germany, Belgium, Denmark, Finland, Norway, Canada, Poland, UK, Spain, Netherlands, and Sweden). SNPs were combined into multi-allelic inverse-variance-weighted fixed or random effects models to generate effect estimates and 95% confidence intervals (CIs). Three complementary sensitivity analyses were performed to examine violations of MR assumptions: MR-Egger regression and weighted median and mode estimators. A Bonferroni-corrected P value threshold was used to establish strong evidence (P < 0.0042) and suggestive evidence (0.0042 < P < 0.05) for associations. In MR analyses, there was strong or suggestive evidence that 2 of the 12 risk factors were associated with invasive epithelial ovarian cancer and 8 of the 12 were associated with 1 or more invasive epithelial ovarian cancer histotypes. There was strong evidence that genetic liability to endometriosis was associated with an increased risk of invasive epithelial ovarian cancer (odds ratio [OR] per 50% higher odds liability: 1.10, 95% CI 1.06-1.15; P = 6.94 × 10-7) and suggestive evidence that lifetime smoking exposure was associated with an increased risk of invasive epithelial ovarian cancer (OR per unit increase in smoking score: 1.36, 95% CI 1.04-1.78; P = 0.02). In analyses examining histotypes and low malignant potential tumours, the strongest associations found were between height and clear cell carcinoma (OR per SD increase: 1.36, 95% CI 1.15-1.61; P = 0.0003); age at natural menopause and endometrioid carcinoma (OR per year later onset: 1.09, 95% CI 1.02-1.16; P = 0.007); and genetic liability to polycystic ovary syndrome and endometrioid carcinoma (OR per 50% higher odds liability: 0.89, 95% CI 0.82-0.96; P = 0.002). There was little evidence for an association of genetic liability to type 2 diabetes, parity, or circulating levels of 25-hydroxyvitamin D and sex hormone binding globulin with ovarian cancer or its subtypes. The primary limitations of this analysis include the modest statistical power for analyses of risk factors in relation to some less common ovarian cancer histotypes (low grade serous, mucinous, and clear cell carcinomas), the inability to directly examine the association of some ovarian cancer risk factors that did not have robust genetic variants available to serve as proxies (e.g., oral contraceptive use, hormone replacement therapy), and the assumption of linear relationships between risk factors and ovarian cancer risk.

CONCLUSIONS

Our comprehensive examination of possible aetiological drivers of ovarian carcinogenesis using germline genetic variants to proxy risk factors supports a role for few of these factors in invasive epithelial ovarian cancer overall and suggests distinct aetiologies across histotypes. The identification of novel risk factors remains an important priority for the prevention of epithelial ovarian cancer.

Gonadotropins and their receptors: coevolution, genetic variants, receptor imaging, and functional antagonists

Gonadotropins belong to the family of dimeric glycoprotein hormones and regulate gonadal physiology mediated by G protein-coupled, seven-transmembrane receptors. These glycoprotein hormones are widely used in the clinic to promote ovarian follicle development and for treating some cases of male infertility. We traced the coevolution of dimeric gonadotropin hormones and their receptors, together with thyrotropin and its receptor. We updated recent findings on human genetic variants of these genes and their association with dizygotic twining, polycystic ovarian syndrome, primary ovarian insufficiency, male-limited precocious puberty, and infertility. In addition to the known physiological roles of gonadotropin-receptor signaling in gonadal tissues, we also discussed emerging understanding of extragonadal functions of gonadotropins in bones and adipose tissues, together with recent advances in in vivo imaging of gonadotropin receptors in live animals. Recent development of gonadotropin receptor agonists and antagonists were summarized with an emphasis on the development of functional antagonists for FSH receptors to alleviate osteoporosis and obesity associated with menopause.