Trio, a novel high fecundity allele: I. Transcriptome analysis of granulosa cells from carriers and noncarriers of a major gene for bovine ovulation rate

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.

Associations between polymorphisms in the myostatin gene with calving difficulty and carcass merit in cattle

A fully functional myostatin gene inhibits muscle fiber growth. The objective of the present study was to quantify the association between 21 known myostatin mutations with both calving and carcass traits in 12 cattle breeds. The myostatin genotypes of 32,770 dam-progeny combinations were used in the association analysis of calving dystocia, with the genotypes of 129,803 animals used in the mixed model association analyses of carcass weight, conformation, and fat score. The mixed model included additive genetic, maternal, and permanent environmental effects where appropriate. The mutant genotypes of nt821, Q204X, and F94L were all associated (P < 0.01) with more calving difficulty when present in either the dam or the progeny. The nt821 deletion had the greatest association with calving difficulty when the homozygous deletion was present in either the calf (0.37 points greater calving difficulty score relative to calves carrying no copies of the deletion based on a one to four scale) or the dam (1.30 points greater calving difficulty score relative to dams carrying no copies of the deletion), although the association between the calf's nt821 genotype and calving difficulty differed depending on the nt821 genotype of the dam. With the exception of nt748_78, nt414, and nt374_51, all other seven segregating myostatin variants were associated (range of allele substitution effect size relative to animals with no copies of the mutant allele) with carcass weight (2.36 kg lighter to 15.56 kg heavier), all 10 segregating variants with conformation (0.15 units less conformed to 2.24 units more conformed assessed on a scale of 1 to 15), and all segregating variants other than E226X with carcass fat (0.23 units less carcass fat cover to 3.85 units more carcass fat cover assessed on a scale of 1 to 15). Of these, the F94L, Q204X, and nt821 mutations generally had the greatest association with all three carcass traits, giving rise to heavier and more conformed carcasses. Despite the antagonistic genetic relationship between calving difficulty and carcass traits, the nt374_51, F94L, and E226X mutations were all associated with improved carcass merit while having minimal expected consequences on calving difficulty. Thus, animals carrying these mutation(s) may have favorable genetic merit for calving difficulty and carcass merit. Furthermore, depending on the dam genotype, a bull with two copies of the nt821 mutation can produce progeny with improved carcass merit while minimizing calving problems.

Identifying genetic variants and pathways influencing daughter averages for twinning in North American Holstein cattle and evaluating the potential for genomic selection

Multiple birth in dairy cattle is a detrimental trait both economically for producers and for animal health. Genetics of twinning is complex and has led to several quantitative trait loci regions being associated with increased twinning. To identify variants associated with this trait, calving records from 2 time periods were used to estimate daughter averages for twinning for Holstein bulls. Multiple analyses were conducted and compared including GWAS, genomic prediction, and gene set enrichment analysis for pathway detection. Although pathway analysis did not yield many congruent pathways of interest between data sets, it did indicate two of interest. Both pathways have ties to the strong candidate region on BTA11 from the genome-wide association analysis across data sets. This region does not overlap with previously identified quantitative trait loci regions for twinning or ovulation rate in cattle. The strongest associated SNPs were upstream from 2 candidate genes LHCGR and FSHR, which are involved in folliculogenesis. Genomic prediction showed a moderate correlation accuracy (0.43) when predicting genomic breeding values for bulls with estimates from calving records from 2010 to 2016. Future analysis of the region on BTA11 and the relation of the candidate genes could improve this accuracy.

Effects of FecB Mutation on Estrus, Ovulation, and Endocrine Characteristics in Small Tail Han Sheep

The Booroola fecundity gene (FecB) has a mutation that was found to increase the ovulation rate and litter size in Booroola Merino sheep. This mutation is also associated with the fecundity of small-tail han (STH) sheep, an important maternal breed used to produce hybrid offspring for mutton production in China. Previous research showed that the FecB gene affects reproduction in STH sheep, based on litter size records. However, the effects of this gene on estrus, ovulation, and endocrine characteristics in these sheep remain unclear. Here, we analyzed the traits mentioned earlier and compared them among the three FecB genotypes of STH ewes using estrus synchronization. Overall, 53 pluriparous ewes were selected from among 890 STH ewes and subjected to FecB genotyping for experiments to characterize estrous and ovulation rates. FecB heterozygous (+B) ewes presented an earlier onset of estrus (42.9 ± 2.2 h) and a shorter estrous cycle (17.2 ± 0.2 days) (P ≤ 0.05). The ovulation rates increased with the increasing copy number of the B allele (P ≤ 0.01). Ovulation time showed no significant differences among the three FecB genotypes. The serum concentrations of follicle-stimulating hormone (FSH), luteinizing hormone, estrogen (E2), and progesterone (P4) were measured in 19 of the ewes. Serum concentrations of E2 and FSH dramatically varied around the time of behavioral estrus. In FecB mutant homozygous (BB) ewes, E2 concentration had two peaks, which were higher (P ≤ 0.05) than those of ++ genotypes. FSH concentration of BB ewes was higher (P ≤ 0.05) than that of the ++ ewes just after estrus. The expression of the estrogen receptor 1 (ESR1) gene in the +B genotype was higher than in the other genotypes. Based on the data for the reproductive performance of STH ewes with the three FecB genotypes, our study suggests that the development of follicles in ewes with the B allele is dependent on the response to FSH regulated by E2 in the early stage. +B ewes, exhibiting moderate ovulation and litter size and a shorter estrous cycle, can be highly recommended in sheep crossbreeding systems for commercial mutton production. Moreover, this study provides useful information to conserve better and use the genetic resources of STH sheep in China.

Selection of fewer dominant follicles in Trio carriers given GnRH antagonist and luteinizing hormone action replaced by nonpulsatile human chorionic gonadotropin†

Studying selection of multiple dominant follicles (DFs) in monovulatory species can advance our understanding of mechanisms regulating selection of single or multiple DFs. Carriers of the bovine high fecundity Trio allele select multiple DFs, whereas half-sib noncarriers select a single DF. This study compared follicle selection during endogenous gonadotropin pulses versus during ablation of pulses with Acyline (GnRH antagonist) and luteinizing hormone (LH) action replaced with nonpulsatile human chorionic gonadotropin (hCG) treatment in Trio carriers (n = 28) versus noncarriers (n = 32). On Day 1.5 (Day 0 = ovulation), heifers were randomized: (1) Control, untreated; (2) Acyline, two i.m. doses (Days 1.5 and D3) of 3 μg/kg; (3) hCG, single i.m. dose of 50 IU hCG on Day 1.5 followed by daily doses of 100 IU; and (4) Acyline + hCG. Treatments with nonpulsatile hCG were designed to replace LH action in heifers treated with Acyline. Acyline treatment resulted in cessation of follicle growth on Day 3 with smaller (P < 0.0001) maximum follicle diameter in Trio carriers (6.6 ± 0.2 mm) than noncarriers (8.7 ± 0.4 mm). Replacement of LH action (hCG) reestablished follicle diameter deviation and maximum diameter of DFs in both genotypes (8.9 ± 0.3 mm and 13.1 ± 0.5 mm; P < 0.0001). Circulating follicle stimulating hormone (FSH) was greater in Acyline-treated than in controls. Finally, Acyline + hCG decreased (P < 0.0001) the number of DFs from 2.7 ± 0.2 to 1.3 ± 0.2 in Trio carriers, with most heifers having only one DF. This demonstrates the necessity for LH in acquisition of dominance in Trio carriers (~6.5 mm) and noncarriers (~8.5 mm) and provides evidence for a role of GnRH-induced FSH/LH pulses in selection of multiple DFs in Trio carriers and possibly other physiologic situations with increased ovulation rate.

Anti-Müllerian hormone production in the ovary: a comparative study in bovine and porcine granulosa cells†

In this study, we aimed to determine the origin of the difference, in terms of anti-Müllerian hormone production, existing between the bovine and porcine ovaries. We first confirmed by quantitative real-time-Polymerase-Chain Reaction, ELISA assay and immunohistochemistry that anti-Müllerian hormone mRNA and protein production are very low in porcine ovarian growing follicles compared to bovine ones. We then have transfected porcine and bovine granulosa cells with vectors containing the luciferase gene driven by the porcine or the bovine anti-Müllerian hormone promoter. These transfection experiments showed that the porcine anti-Müllerian hormone promoter is less active and less responsive to bone morphogenetic protein stimulations than the bovine promoter in both porcine and bovine cells. Moreover, bovine but not porcine granulosa cells were responsive to bone morphogenetic protein stimulation after transfection of a plasmidic construction including a strong response element to the bone morphogenetic proteins (12 repetitions of the GCCG sequence) upstream of the luciferase reporter gene. We also showed that SMAD6, an inhibitor of the SMAD1-5-8 pathway, is strongly expressed in porcine compared to the bovine granulosa cells. Overall, these results suggest that the low expression of anti-Müllerian hormone in porcine growing follicles is due to both a lack of activity/sensitivity of the porcine anti-Müllerian hormone promoter, and to the lack of responsiveness of porcine granulosa cells to bone morphogenetic protein signaling, potentially due to an overexpression of SMAD6 compared to bovine granulosa cells. We propose that the low levels of anti-Müllerian hormone in the pig would explain the poly-ovulatory phenotype in this species.

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.

A Comparative Analysis of Oocyte Development in Mammals

Sexual reproduction requires the fertilization of a female gamete after it has undergone optimal development. Various aspects of oocyte development and many molecular actors in this process are shared among mammals, but phylogeny and experimental data reveal species specificities. In this chapter, we will present these common and distinctive features with a focus on three points: the shaping of the oocyte transcriptome from evolutionarily conserved and rapidly evolving genes, the control of folliculogenesis and ovulation rate by oocyte-secreted Growth and Differentiation Factor 9 and Bone Morphogenetic Protein 15, and the importance of lipid metabolism.

Trio a novel bovine high-fecundity allele: II. Hormonal profile and follicular dynamics underlying the high ovulation rate

The newly discovered Trio high-fecundity allele produces multiple ovulations in cattle. This study evaluated (1) size and growth rates of follicles in Trio carriers during a synchronized follicular wave, induced by follicle aspiration; (2) follicle-stimulating hormone (FSH) patterns associated with the follicular wave; (3) size of corpora lutea (CL) and circulating progesterone; and (4) intrafollicular estradiol concentrations prior to normal deviation. Trio carriers had mean dominant follicles that were significantly smaller in diameter and volume than noncarriers. Onset of diameter deviation occurred at ∼3 days after the last follicle aspiration in both genotypes despite Trio carriers having much smaller individual follicles. Follicles of Trio carriers grew at a slower rate than noncarrier follicles (∼65% in mm/day or ∼30% in mm3/day) resulting in much smaller individual dominant follicles (∼25% volume). However, total dominant follicle volume, calculated as the sum of all dominant follicles in each animal, was similar in carriers and noncarriers of Trio throughout the entire follicular wave. Circulating FSH was greater in Trio carriers during the 24 h encompassing deviation. Trio carriers had significantly more ovulations than noncarriers, and individual CL volume was smaller, although total luteal tissue volume and circulating P4 were not different. Thus, increased ovulation rate in Trio carriers relates to smaller individual follicles (one-third the volume) near the time of deviation due to slower follicle growth rate, although time of deviation is similar, with increased circulating FSH near deviation leading to selection of multiple dominant follicles in Trio carriers with similar total follicle volume.

Trio, a novel bovine high fecundity allele: III. Acquisition of dominance and ovulatory capacity at a smaller follicle size

The acquisition of dominance and ovulatory capacity was evaluated in follicles from cows that were carriers or half-sibling noncarriers of the Trio allele. Follicle size at acquisition of follicular dominance was determined by evaluating whether follicles ovulate after GnRH challenge (ovulatory capacity-experiment 1) and by determination of intrafollicular concentrations of estradiol and free insulin like growth factor 1 (IGF1) and relative mRNA expression of cytochrome P450 family 19 subfamily A member 1 (CYP19A1), luteinizing hormone/choriogonadotropin receptor (LHCGR), and pappalysin 1 (PAPPA, previously known as pregnancy-associated plasma protein A, pappalysin 1) in granulosa cells from follicles of different sizes (experiment 2). Ovulatory capacity developed in follicles at 8.3 mm (50% ovulatory capacity) in noncarriers but at smaller sizes (5.5 mm) in Trio carriers. Similarly, in experiment 2, follicles of Trio carriers acquired a dominant phenotype, as determined by intrafollicular estradiol and CYP19A1, LHCGR, and PAPPA mRNA expression in granulosa cells, at significantly smaller sizes but at a similar time after wave emergence. Overall, dominance/ovulatory capacity was acquired when follicles of Trio carriers were ∼30% the size (volume basis) of follicles in noncarriers. In addition, follicles in Trio carriers appear to acquire dominance in a hierarchal manner, as demonstrated by the progressively greater number of follicles with a dominant phenotype between days 2 and 4 after wave emergence. Thus, results from this study provide further support for a physiological model in which selection of multiple follicles in Trio allele carriers is characterized by acquisition of dominance at a smaller follicle size but at a similar time in the follicular wave with multiple follicles acquiring dominance in a hierarchal sequence.

Complex genetics of female fertility

Variation in reproductive lifespan and female fertility have implications for health, population size and ageing. Fertility declines well before general signs of menopause and is also adversely affected by common reproductive diseases, including polycystic ovarian syndrome (PCOS) and endometriosis. Understanding the factors that regulate the timing of puberty and menopause, and the relationships with fertility are important for individuals and for policy. Substantial genetic variation exists for common traits associated with reproductive lifespan and for common diseases influencing female fertility. Genetic studies have identified mutations in genes contributing to disorders of reproduction, and in the last ten years, genome-wide association studies (GWAS) have transformed our understanding of common genetic contributions to these complex traits and diseases. These studies have made great progress towards understanding the genetic factors contributing to variation in traits and diseases influencing female 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. Many variants implicate DNA damage/repair genes in variation in the age at menopause 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.

The local regulation of folliculogenesis by members of the transforming growth factor superfamily and its relevance for advanced breeding programmes

Regulation of the growth and maturation of the ovarian follicle is critical for normal reproductive function. Alterations in this growth can lead to pathological conditions, such as cystic follicles, reduced oocyte quality, or an abnormal endocrine environment leading to poor fertility. Alterations in follicular growth also influence the number of follicles ovulating and thus can change litter size. Both endocrine factors, such as follicle stimulating hormone and luteinizing hormone, as well as local factors, are known to regulate follicular growth and development. This review will focus on the role of local factors in regulation of ovarian follicular growth in ruminants, with a focus on members of the transforming growth factor superfamily. The potential role of these factors in regulating proliferation, apoptosis, steroidogenesis and responsiveness to gonadotrophins will be considered.

Proteomic analysis of follicular fluid in carriers and non-carriers of the Trio allele for high ovulation rate in cattle

This study was conducted to characterise differences in follicular fluid proteins between carriers and non-carriers of a bovine allele for high ovulation rate. A total of four non-carrier and five carrier females were used in an initial study with four and six additional non-carriers and carriers respectively used in a validation study. Emergence of the follicular wave was synchronised and the ovaries containing the dominant follicle(s) were extracted by ovariectomy for follicular fluid collection. A hexapeptide ligand library was used to overcome the masking effect of high-abundance proteins and to increase detection of low-abundance proteins in tandem mass spectrometry. After correcting for multiple comparisons, only two proteins, glia-derived nexin precursor (SERPINE2) and inhibin β B chain precursor (INHBB), were significantly differentially expressed (false-discovery rate <0.05). In a replicate study of analogous design differential expression was confirmed (P < 0.05). Joint analysis of results from the two studies indicated that three additional proteins were consistently differentially expressed between genotypes. For three of these five, previous studies have indicated that expression is increased by transforming growth factor-β–bone morphogenetic protein signalling; their reduction in follicular fluid from carrier animals is consistent with the ~9-fold overexpression of SMAD family member 6 (SMAD6) in carriers that is inhibitory to this pathway.