Abstract
A fundamental tenet of inheritance in sexually reproducing organisms such as humans and laboratory mice is that gametes combine randomly at fertilization, thereby ensuring a balanced and statistically predictable representation of inherited variants in each generation. This principle is encapsulated in Mendel’s First Law. But exceptions are known. With transmission ratio distortion, particular alleles are preferentially transmitted to offspring. Preferential transmission usually occurs in one sex but not both, and is not known to require interactions between gametes at fertilization. A reanalysis of our published work in mice and of data in other published reports revealed instances where any of 12 mutant genes biases fertilization, with either too many or too few heterozygotes and homozygotes, depending on the mutant gene and on dietary conditions. Although such deviations are usually attributed to embryonic lethality of the underrepresented genotypes, the evidence is more consistent with genetically-determined preferences for specific combinations of egg and sperm at fertilization that result in genotype bias without embryo loss. This unexpected discovery of genetically-biased fertilization could yield insights about the molecular and cellular interactions between sperm and egg at fertilization, with implications for our understanding of inheritance, reproduction, population genetics, and medical genetics.
The Problem
Our understanding of inheritance in sexually reproducing organisms assumes, with good evidence, that the combination of egg and sperm that join at fertilization is largely independent of their genetic content. Equal transmission of alternative alleles through meiosis in heterozygotes ensures a balanced parental genetic contribution to offspring at each generation. Independent segregation and random union of gametes at fertilization are foundations of classical, quantitative, population, evolutionary, and medical genetics (East 1916; Holliday 1984; Crow 1991). Mendel’s First Law captures this principle, which is one of the few that applies generally in biology.
The most prominent exceptions to random segregation are the rare, naturally occurring examples of transmission ratio distortion (TRD) that have been described in fungi (Turner and Perkins 1979), corn (Rhoades and Vilkomerson 1942), flies (Morgan et al. 1925; Gershenson 1928; Sandler and Novitski 1957; Pimpinelli and Dimitri 1989; Larracuente and Presgraves 2012), mice (Dunn 1957; Montagutelli et al. 1996; Lyon 2003; Wu et al. 2005; Didion et al. 2015, 2016), humans (Chen et al. 2012; Huang et al. 2013), and other species (Scofield et al. 1982; Fishman and Willis 2005; Koide et al. 2008; Hoang et al. 2016). Biased sex ratios have also been reported (Hamilton 1967; Jaenike 2001; Tao et al. 2007a,b; Helleu et al. 2014, 2016; Rice 2014). These exceptions arise despite strong selective pressures to maintain normal segregation (Crow 1991) and sex ratios (Fisher 1930; Crow and Kimura 1970). Based on alternative functions in haploid gametes, one allele is preferentially transmitted to offspring at the expense of other alleles. In these cases, TRD drives allelic preference in one sex, regardless of the genetics of the mating partner. Reproductive performance is often not reduced because the normal number of gametes and zygotes is produced. TRD may arise during chromosome segregation in meiosis (meiotic drive), gametogenesis (gamete competition), or embryonic development (preferential lethality). These examples of TRD probe the nature of Mendel’s First Law by illuminating genetic, molecular, and cellular mechanisms that underlie meiosis, recombination, gametogenesis, and early development (East 1916; Crow 1979, 1988, 1991; Holliday 1984; Lyttle 1993; Herrmann et al. 1999; Pardo-Manuel de Villena and Sapienza 2001; Lyon 2003). Many of these “selfish genetic” systems are composed of several closely linked elements that not only lead to preferential transmission of the chromosome on which they are carried, but also confer sterility or lethality on homozygous carriers, thereby preventing fixation at the cost of reduced population fitness (Dunn 1957). Without these counterbalances, TRD systems would quickly replace their wild-type (WT) alleles in the population, their competitive advantage would be lost, and selective sweeps affecting variation at neighboring loci would be their only legacy.
With spontaneous, induced, and engineered single-gene mutations, one of the first tasks is to assess consequences on viability, fertility, and other phenotypes (White et al. 2013; Hrabe de Angelis et al. 2015; Dickinson et al. 2016; Meehan et al. 2017). Absence of mutant homozygotes is usually accepted as evidence for induced lethality and reduced numbers of heterozygotes is taken as evidence for a detrimental dosage effect (Figure 1). But sometimes the fit to Mendelian segregation is not explicitly examined. Litter size, which should be reduced proportionately to the number of missing genotypes, is often not reported. Backcrosses, which can provide information about parent-of-origin effects on gametogenesis and embryogenesis, are sometimes not included in study designs. As a result, whether particular cases of non-Mendelian segregation result from lethality or from other phenomena such as TRD remains ambiguous, despite claims in publications.
Regular and irregular outcomes of Mendelian segregation. (A) Conventional segregation in backcrosses and intercrosses for a single gene with two alleles, mutant (m) and wild-type (+). (B) Irregular segregation in an intercross with deficiency of all m/m homozygotes and half of +/m heterozygotes as an exemplar.
Consider an early controversy in mammalian genetics. Cuenot (1905), studying the absence of pure yellow segregants in crosses between mice heterozygous for the dominantly-acting agouti yellow (Ay) mutation, argued that gametes carrying the agouti yellow allele never join at fertilization. By contrast, Castle and Little, based on considerations of both segregation ratios and litter size, correctly concluded that homozygous yellow mice fail to complete development, with reduced litter size providing the critical evidence for lethality (Castle and Little 1910; Heaney et al. 2009). As Castle and Little showed, a full analysis is needed to establish with confidence the basis for unusual segregation.
Four steps are needed to transfer genetic and epigenetic information from one generation to the next through the germline. Meiosis converts the chromosome complement from diploid to haploid in the parental generation. Gametogenesis provides a cellular vehicle for the haploid genome. From a pair of haploid gametes, fertilization restores diploidy in the zygote. Finally, the germline is set aside early in embryonic development to renew these steps in the offspring generation. TRD has been reported for three of these steps: meiotic drive, gamete competition, and preferential embryo survival. In each case, dysfunction in one sex is sufficient for TRD that is largely independent of the genetics of mating partners.
The evidence reviewed here provides examples for TRD in the third step, fertilization, where genetic variants, acting in both sexes and in some cases depending on environmental (dietary) conditions, control the combination of gametes that join at fertilization to create zygotes. Our work on epigenetic inheritance in mice and a review of the mouse literature revealed strong evidence that mutations in any of 12 genes (Box 1), where either too many or too few heterozygotes were found among intercross progeny, together with absent or deficient homozygotes, were found without evidence for dead embryos or reduced litter size (Table 1 and Table 2; see also Supplemental Material, Table S1 and Table S2). Of these 12 genes, six involve single spontaneous or engineered mutations on an inbred genetic background. In another case, biased segregation was found only in crosses involving a pair of mutant genes (epistasis). Finally, five cases involved dietary folic acid supplementation of mice carrying a single-gene mutation affecting neural tube development, where segregation was biased on one diet but normal on the alternative diet, with similar litter sizes and rates of prenatal lethality. Normal segregation in backcrosses between heterozygotes and WT homozygotes argues that meiosis and gametogenesis function normally in each sex. These unusual results suggest that fertilization is genetically biased toward particular gametes based on their genetic content.
Genes that bias fertilization. Additional information can be found in the Mouse Genome Database (informatics.jax.org).
A1cf
Apobec1 complementation factor, chromosome (Chr) 19, 26.6 cM. A1cf is expressed primarily in the nucleus where it encodes the RNA-binding subunit for Apolipoprotein B mRNA editing enzyme, catalytic polypeptide 1 (APOBEC1), a cytidine deaminase that edits specific bases, sometimes in coding sequences but more usually in 3′-UTRs (Blanc and Davidson 2010; Blanc et al. 2014). It must have additional functions since deficiency leads to early embryonic lethality (Blanc et al. 2005), while APOBEC1-deficient mice are fully viable and fertile (Hirano et al. 1996; Nelson et al. 2012). Partial deficiency increases TGCT risk in a parent-of-origin manner (Carouge et al. 2016).
Ago2
Argonaute RNA-induced silencing complex (RISC) catalytic subunit 2, Chr 15, 33.9 cM. AGO2 is required for RNA-mediated gene silencing (RNAi) by RISC (Ender and Meister 2010). Guide RNAs (miRNAs and siRNAs) direct RISC to complementary RNAs that are targets for RISC-mediated gene silencing. AGO2 and serine/threonine protein phosphatase 2 catalytic subunit β isoform (PPP2CB) (see below) promote mitotic chromosome segregation in the Drosophila and Caenorhabditis elegans germline (Claycomb et al. 2009; Pek and Kai 2011a,b). Interestingly, AGO3 piRNA component Aubergine enhances transmission distortion for the segregation distortion (SD) in Drosophila (Gell and Reenan 2013), raising the possibility that AGO2 might have similar effects under appropriate circumstances. Kirsten RAS oncogene homologue (KRAS) signaling controls AGO2 sorting into exosomes (McKenzie et al. 2016) that transport RNAs, including tRNA fragments (Sharma et al. 2016) for intercellular signaling. RNAs transferred to, as well as produced in, sperm could have significant effects on gamete functions (Hosken and Hodgson 2014). Selectivity in exosome targeting could lead to functional differences among haploid gametes. Loss of siRNA but not miRNA AGO2 activity leads to meiotic catastrophe in MI oocytes (Galgano et al. 2008; Leibovich et al. 2010; Stein et al. 2015). In AGO2-deficient mice, miRNA levels are reduced substantially in oocytes (Kaneda et al. 2009) and in AGO2-deficient oocytes siRNA levels are reduced while retrotransposons and selected mRNA levels are increased (Watanabe et al. 2008). Homozygous-deficient mice display embryonic lethality with various defects in embryonic and extraembryonic organs and tissues. Partial deficiency increases TGCT risk in a parent-of-origin manner (Carouge et al. 2016).
Apob
Apolipoprotein B, Chr 12, 3.53 cM. APOB is widely expressed where it transports lipids such as cholesterol. APOB is encoded as a single, long mRNA. The shorter apoB-48 protein is produced after RNA editing of the apoB-100 transcript at residue 2180 (CAA→UAA), resulting in the creation of a stop codon and early translation termination. Homozygous deficiency leads to embryonic lethality, with embryo loss by embryonic day E9. Heterozygotes tend to have incomplete neural tube closure. Partial deficiency severely reduces fertility in males with sperm showing impaired motility and reduced ability to fertilize both in vivo and in vitro, arguing for a diploid rather than a haploid effect (Huang et al. 1995, 1996).
Apobec1
Apobec1, Chr 6, 57.7 cM. Apobec1 encodes a cytidine deaminase that edits C to U (read as T in coding regions) primarily in 3′-UTRs (Blanc and Davidson 2010). A1CF is the RNA-binding protein that targets specific mRNA sites for editing (Blanc et al. 2005). Apobec1 interacts with the Dnd1Ter mutation to increase TGCT risk in a conventional manner in the male germline and in a transgenerational manner in the female germline (Nelson et al. 2012).
Dnd1
Dead-end microRNA-mediated repression inhibitor 1, Chr 18, 19.5 cM. Dnd1 is an A1cf-related RNA-binding protein expressed in many tissues (Youngren et al. 2005). DND1 controls access of particular miRNAs to their mRNA targets in human TGCTs (Kedde et al. 2007). It is essential for germ cell differentiation (Bustamante-Marin et al. 2013) and also acts like A1cf and Ago2 in offspring to process inherited epigenetic changes from parents (Carouge et al. 2016). DND1 and NANOS2 load RNAs onto the CCR4-Not (CNOT) deadenylase complex for germ cell differentiation (Suzuki et al. 2016) and to suppress specific RNAs (Suzuki et al. 2016). The Ter mutation severely reduces fertility and is a potent modifier of TGCT susceptibility (Stevens 1973; Youngren et al. 2005), whereas the targeted deficiency results in biased fertilization (Zechel et al. 2013).
Ddx1
DEAD (Asp-Glu-Ala-Asp) box polypeptide 1, Chr 12, 6.4 cM. Ddx1 is expressed primarily in the nucleus of the fetal and adult testis and ovary, where it functions as an ATP-dependent RNA helicase to unwind RNA–RNA and RNA–DNA secondary structures for translation initiation, nuclear and mitochondrial splicing, ribosome and spliceosome assembly (including tRNAs), and pre-miRNA and polyA processing (Han et al. 2014; Popow et al. 2014). Reduced DDX1 activity promotes ovarian tumor growth (Han et al. 2014). In Drosophila, Ddx1 deficiency results in stress (starvation)-induced sterility in males and autophagy in egg chambers (Germain et al. 2015).
L3P
No information except Marean et al. (2011).
Lrp6
Low density lipoprotein receptor-related protein 6, Chr 6, 65.4 cM. Lrp6 is expressed at E11.5 and after in the female reproductive system and in fetal testes (Kofron et al. 2007). Lrp6 encodes a transmembrane cell surface protein involved in receptor-mediated endocytosis of lipoprotein and protein ligands. It can function alone or as a coreceptor with Frizzed for canonical Wnt/β-catenin signaling. LRP6 and other members of the Hedgehog and WNT pathways are expressed in human embryonic stem cells and testicular cancers (Dormeyer et al. 2008). Partial embryonic lethality, growth retardation, and various vertebral and skeletal abnormalities are found in mutant homozygotes, and subtle skeletal defects in mutant heterozygotes.
Ppp2cb
Serine/threonine-protein phosphatase 2A catalytic subunit beta isoform, Chr 8, 20.6 cM. Ppp2cb encodes the 2A catalytic subunit of the protein phosphatase 2A (PP2A) heterodimer and is expressed in the female and male reproductive systems from E15 through adulthood. In particular, it is highly expressed in mature spermatozoa and MII oocytes where it localizes at centromeres in meiosis and spindle poles in mitosis (Su et al. 2012a). It is a negative regulator of the mitogen-activated protein kinase (MAPK) pathway and plays a role in the DNA damage response, cell cycle control, apoptosis, and mRNA splicing. Inhibition of PPP2CB releases meiotic arrest and enables meiotic progression (Su et al. 2012a,b). Loss of PPP2CB in oocytes causes both failure of meiosis II exit and reduced fertility in females (Su et al. 2012a,b) and males (Pan et al. 2015). Related proteins play similar roles in meiotic control and chromosome segregation (Bizzari and Marston 2011; Chang et al. 2011; Chambon et al. 2013). PPP2CB deficiency affects sperm tails (Soler et al. 2009; Pan et al. 2015) and chromosome segregation in females (Su et al. 2012a). Other reports find viable and fertile Ppp2cb mutant homozygotes without obvious phenotype (Gu et al. 2012).
Pum1
Pumilio1 RNA-binding family member 1, Chr 4, 63.4 cM. Pum1 is a widely expressed cytoplasmic protein found in the ovary and testis throughout fetal development and adulthood. It encodes an RNA-binding protein that targets Pumilio Response Elements (PRE) in 3′-UTRs to recruit both CCR4-NOT deadenylase, other deadenylases, and miRNAs such as miR221 and miR222, which together repress the expression of genes such as p27 that maintain genome integrity (Kedde et al. 2010; Miles et al. 2012; Van Etten et al. 2012). Interestingly, DND1 blocks access of miR221 to its p27 mRNA target (Kedde et al. 2010), suggesting that DND1 acts competitively with PUM1 to control miRNA actions. In the testis, PUM1 acts as a post-transcriptional regulator of spermatogenesis by binding to the 3′-UTR of mRNAs coding for regulators of transformation-related protein 53 (TRP53) (Subasic et al. 2016) and also suppresses caspase- and TRP53-apoptosis in germ cells (Chen et al. 2012; Subasic et al. 2016). It is involved in embryonic stem cell renewal by facilitating the transition from pluripotency to differentiation (Leeb et al. 2014). PUM1-deficient males exhibit reduced weight of testes and seminiferous tubules, reduced number of sperm, and increased germ cell apoptosis and infertility (Chen et al. 2012). PUM1 also controls to the number of primordial ovarian follicles, meiosis, and reproductive competence in females (Mak et al. 2016). Finally, PUM1 contributes to the antiviral response (Narita et al. 2014), suggesting that it might play a more general role in the stress response to environmental and physiological conditions that often lead to transgenerationally-inherited epigenetic changes (Schaefer and Nadeau 2015).
Vangl2
Vang-like planar cell polarity protein 2, Chr 1, 79.5 cM. VANGL2 is found in the plasma membrane and cytoskeleton where it provides directional signals to cilia (Torban et al. 2012). It is expressed in female and male reproductive systems from E15 through adulthood. Both male homozygous mutants (Lp/Lp) and female heterozygous mutants (Lp/+) are sterile (Strong and Hollander 1949).
Zic2
Zinc finger protein of the cerebellum 2, Chr 14, 66.0 cM. ZIC2 is expressed in the female and male reproductive system from E15 through adulthood and represses transcription in the nucleus. ZIC2 promotes self-renewal of liver cancer stem cells by recruiting the nuclear remodeling factor complex to activate the octamer-binding transcription factor 4 (OCT4) pluripotency factor (Zhu et al. 2015). It also enhances transcription to promote differentiation of embryonic stem cells in Drosophila (Luo et al. 2015) and controls naïve vs. primed pluripotency state (Buecker et al. 2014). Deficiency results in neurulation defects and embryonic lethality in mice (Marean et al. 2011) and holoprosencephaly in humans (Brown et al. 1998).
Historically, the genetics of fertilization has been largely resistant to molecular studies (Furnes and Schimenti 2007; Wright and Bianchi 2016). Discovery of genes and gene–diet conditions that bias fertilization may be a breakthrough in understanding mechanisms of sperm–egg interactions at fertilization. Here, evidence for TRD resulting from nonrandom union of gametes in mice is reviewed, possible mechanisms proposed, genetic and reproductive considerations discussed, and genetic implications considered.
Evidence: Reanalysis of the Literature
Genetics
The principles and methods are outlined and then evidence summarized for simple and complicated single-gene effects, as well as an unusual two-gene effect.
As expected, most genetic variants segregate normally and show reduced litter size in proportion to genotype bias, both in our hands and in the literature (informatics.jax.org; www.komp.org), suggesting that biased fertilization is exceptional and results from specific rather than general dysfunctions. The evidence reported here involves single-gene mutations on an inbred strain background. Two cases of fertilization bias have been reported previously (Agulnik et al. 1993; Wedekind et al. 1996), although the responsible genetic factor is not known in either case.
Principles and methods:
Central to this evidence are three expected outcomes: (1) a biased genotype distribution and reduced litter size suggesting either selective embryonic lethality or reduced oocyte production; (2) biased genotype distributions in both backcrosses and intercrosses together with normal litter sizes suggesting a bias in meiotic segregation or gamete numbers or function; and (3) biased genotype distributions in intercrosses but not backcrosses together with normal litter sizes suggesting that gametogenesis is normal but fertilization is biased.
Three measures of TRD were used, depending on the particular question: (1) the frequency of m/+ heterozygotes among m/+ and +/+ segregants; (2) the frequency of the m allele among all offspring; (3) the ratio of m/+ heterozygotes to +/+ homozygotes, and (4) the ratio of m/m to +/+ homozygotes. To test for departures from expectations and to summarize evidence for nonrandom fertilization, emphasis was placed first on testing departures from Mendelian expectations (hypothesis H1), namely 1:1 in backcrosses and 1:2:1 in intercrosses (Figure 1), and then on measures of the nature and magnitude of these departures (Table 1 for genetic effects and Table 2 for gene–folic acid interactions, see also Table S1 and Table S2 for complete data, analytical methods, and results).
Effect size is an important but often neglected measure of phenotypic variation and was used here to provide a normalized quantitative measure of departures from expectations for genetic and gene–diet effects (Cohen 1988). According to accepted standards (Cohen 1988), effect sizes > 0.10 are classified as “small,” > 0.30 “medium,” and > 0.50 “strong.” These measures are independent of sample size. This combination of statistical tests and effect size estimates ensured that significant results and their phenotypic effects were highlighted.
Single genes, simple cases:
For Dnd1 and Ago2, segregation was highly unusual with significant deficiencies of m/+ and m/m genotypes among intercross progeny, - examples of “too few heterozygotes” (Table 1, see also Table S1). For backcrosses, genotype transmission ratios were close to 1:1 expectations and effect sizes were small, whereas for intercrosses these ratios were closer to 1 (=1:1) than 2 (=2:1) and effect sizes were medium (Table 1). Heterozygotes that are missing in intercrosses are found in expected numbers in backcrosses, suggesting that they do not carry intrinsic deficits. In addition, embryonic lethality does not account for the observed genotype ratios because litter sizes were similar among intercrosses and backcrosses (Table S1). For Dnd1, neither mutant homozygotes nor dead embryos were found at embryonic day E3.5 (Zechel et al. 2013). For Ago2, a 25% reduction in litter size would be consistent with the reported lethality of m/m mutant homozygotes (Morita et al. 2007), but not with the observed 50% deficiency of m/+ heterozygotes.
For A1cf, Ppp2cb, and Pum1, highly significant excesses of m/+ heterozygotes were found together with an absence of mutant homozygotes among intercross progeny, - examples of “too many heterozygotes” (Table 1, see also Table S1). For backcrosses, the ratios of m/+ heterozygotes to WT were close to 1:1 expectations and effect sizes were small, whereas for intercrosses, these ratios ranged from 3.1 to 9.7, instead of the expected ratio of 2, and effect sizes were large (Table 1). Despite these differences, average litter sizes were remarkably similar for backcrosses and intercrosses (Table S1). For Pum1, m/m embryos were not detected at E3.5 and litter size did not differ between intercrosses and backcrosses (Zhang et al. 2015) (Significant departures from expectations were also found in intercrosses after accounting for the modest departures in backcrosses, χ2 = 12.64, P < 0.001, see Table S1). Pum1 also showed modest departures in backcrosses, For A1cf, loss of homozygous embryos between E3.5 and 4.5 does not account for excess heterozygosity or for normal litter size (Blanc et al. 2005). Although litter size was not reported for Ppp2cb (Sasaki et al. 2007), the highly significant heterozygote excess (> 3:1) in intercrosses vs. backcrosses is striking and consistent with results for A1cf, Ddx1 (see below), and Pum1. Normal segregation in backcrosses with mutant heterozygotes shows that gametes are produced in comparable (1:1) numbers and functionality in each sex.
Ddx1: complicated single-gene effect:
An engineered deficiency of DEAD Box 1 helicase (Ddx1) and an induced epigenetic change in its WT allele provide strong evidence for biased fertilization (Hildebrandt et al. 2015; Hollick 2017). For the engineered mutation, m/m embryos were missing with no homozygotes detected at E3.5 in test crosses. In addition, a substantial deficit of WT segregants relative to m/+ heterozygotes was also observed in some crosses (Table 3 and Table S1), leading the authors to conclude that the engineered mutation induced a modified, perhaps paramutated allele (Hollick 2017) (designated “*”) that leads to lethality in both WT (+/+) and */* embryos resulting from crosses involving a */+ parent. However, review of litter sizes for test and control crosses revealed remarkably little evidence for embryo loss that would account for the putative loss of m/m, */*, and +/+ embryos. Strongly biased transmission without embryo loss argues that preferential fertilization is a more likely explanation.
Apobec1 and Dnd1: a complicated two-gene effect:
The last genetic example emerged in tests to determine whether Apobec1 and Dnd1 interact to modulate inherited risk for spontaneous testicular germ cell tumors (TGCTs; see below for additional information about TGCT origins, genetics, and biology) (Nelson et al. 2012). These genes independently affect risk in both a conventional and an inherited epigenetic manner (Nelson et al. 2012), but whether they interact in the genetic sense was uncertain. Surveys for TGCTs among intercross offspring revealed unexpected evidence for biased fertilization (Table 4). With maternal heterozygosity for Dnd1 and paternal heterozygosity for Apobec1, offspring occurred in the expected (1:1:1:1) Mendelian ratio for two independently segregating genes, but results for the reciprocal cross revealed a marked deficiency of all single- and double-mutant genotypes. If the number of WT (+/+, n = 64) progeny is accepted as the proper reference for the three other genotypic classes, where 1:1:1:1 = 64:64:64:64, then only 27% (51/192, where 192 = 3 × 64) of the expected single- and double-mutant segregants was found. Again, litter size did not differ between the reciprocal crosses (Nelson et al. 2012). Remarkably, both heterozygous and homozygous mutants for each gene are fully viable in separate crosses (Table S1) (Hirano et al. 1996; Blanc et al. 2005; Youngren et al. 2005; Zechel et al. 2013). Thus, in this but not the reciprocal cross, the majority of mutant segregants are missing, with evidence for full viability in other crosses and with no evidence for reduced litter size.
Gene–folate diet interactions
A related class of TRD involves dietary effects on single-gene models of neural tube defects (NTDs). The beneficial effect of dietary folate supplementation on a common birth defect is one of the greatest achievements in epidemiology (Smithells et al. 1980). NTDs such as anencephaly and spina bifida are the second most common congenital defect, occurring in ∼1 per 1000 live births and often leading to disability and mortality (Detrait et al. 2005; Copp et al. 2013). Mothers and fetuses frequently show reduced folate and elevated homocysteine levels (Smithells et al. 1980; Carmel and Jacobsen 2001). Dietary folate supplementation before and during pregnancy significantly reduces the occurrence and severity of cases, but many (∼50%) remain resistant to the beneficial effects of folate supplementation (Smithells et al. 1980; Blom et al. 2006; Copp et al. 2013). Reliable prediction of individual response to supplementation is currently impossible.
Several studies examined the effects of dietary supplementation with folic acid on development of the neural tube in mouse models (Juriloff and Harris 2000; Greene and Copp 2005; Zohn 2012). Some models respond favorably to folic acid (Barbera et al. 2002) and several respond to other nutrients such as methionine and inositol (Essien 1992; Greene and Copp 1997), but as in humans many are resistant (Juriloff and Harris 2000; Franke et al. 2003; Greene and Copp 2005; Zohn 2012). In the design for these studies, developmental response to various nutrients was tested among intercross progeny. In some cases, +/m heterozygotes are maintained on a test or control dietary formulation, bred together (+/m × +/m), and developmental consequences evaluated among +/+, +/m, and m/m offspring, where offspring are expected to occur in a 1:2:1 Mendelian ratio. In other cases, pregnant +/m females are treated only during critical developmental windows and then embryos or offspring are examined. The former is the only protocol that is instructive about possible dietary effects on gametogenesis and fertilization. Attributing effects to diet simply involves comparing results for the high- vs. low-dosage treatment groups for the identical genotypes on the same defined genetic background. Efficacy is indicated if the number (proportion) of affected m/m individuals is reduced.
During studies to identify mouse NTD models that are resistant to the benefits of folic acid treatments, two examples of fertilization bias were found (Nakouzi and Nadeau 2014). Reanalysis of published reports then revealed three additional NTD models that show biased fertilization in response to folic acid treatment (Table 2). In several cases (Apob, Lrp6, and Vangl2), significant deviations are found in the high folic acid (10 ppm; parts per million) test, but in other cases (L3P and Zic2) deviations were found in the low (2 ppm) test, with similar litter sizes for each mutant on the two diet protocols (Table S2). On one diet, transmission ratios were consistent with normal segregation and effect sizes were small, whereas on the alternative diet, transmission ratios were more divergent and effect sizes were medium (Table 2, see also Table S2). None showed too many heterozygotes. Departures from Mendelian expectations with dietary supplementation were not as strong as the genetic results, perhaps because optimal supplementation levels were not tested and perhaps because dietary consumption and metabolism differ among the various cohorts and are generally more difficult to control. Interestingly, for most models, the percentage of affected m/m was similar in the test and control protocols, suggesting that supplemental folic acid did not reduce the proportion of affected mutant homozygotes (Gray et al. 2010; Marean et al. 2011).
These results raise a provocative question: does folate correct a developmental defect in the neural tube, or do other explanations apply such as reducing the incidence of cases by biasing fertilization away from at-risk genotypes (Nakouzi and Nadeau 2014)? In humans, where NTD genetics is not as clearly understood as in mouse models, with genetic heterogeneity a significant issue and isolated cases common, distinguishing a protective effect vs. biased fertilization would be difficult.
Gene functions
Six of the seven TGCT genes encode proteins that are directly involved in various aspects of RNA biology: A1cf, RNA editing (Blanc and Davidson 2010; Snyder et al. 2017a,b); Apobec1, RNA editing (Blanc and Davidson 2010); Ago2, RNA interference (RNAi) (Ender and Meister 2010); Dnd1, miRNA control (Kedde et al. 2007); Ddx1, RNA helicase (Hildebrandt et al. 2015); and Pum1, translation repression (Zhang et al. 2015) (Box 1). The seventh gene, Ppp2cb, is a serine/threonine phosphatase (Sasaki et al. 2007) (Box 1). Four show inherited epigenetic effects on TGCT risk (A1cf, Apobec1, Ago2, and Dnd1) and one shows epigenetic effects on embryonic viability (Ddx1). Neither Ppp2cb nor Pum1 have been tested for TGCT or epigenetic effects. These proteins (except Ppp2cb) have specific RNA targets that are in turn effectors of molecular, developmental, and physiological functions. These functions could be shared or distinct in males and females depending on the nature of the targeted RNAs in each sex. Identifying these targets is critical for understanding the ways that these genes control gamete choice at fertilization.
By contrast, the five NTD genes appear to have diverse, seemingly unrelated functions without an obvious theme (Apob, L3P, Lrp6, Vangl2, and Zic2; Box 1). Perhaps RNA genes and Ppp2cb are directly involved in molecular and cellular mechanisms of gametogenesis and fertilization, whereas the various NTD genes sensitize folate and homocysteine metabolism to adverse interactions with pathways that directly affect gamete interactions at fertilization.
Hypotheses
Reversed meiosis
If meiotic divisions in females are ordered the way we have been taught, with the reductional division (MI) preceding the equational division (MII) (Figure 2), biased fertilization without conventional TRD would be difficult. During ovulation, primary oocytes resume meiosis and, over the next few hours, complete MI and arrest at metaphase in MII (Edwards and Gates 1959) Homologous (nonsister) chromatid pairs segregate at MI with one product going to the secondary oocyte and the other to the first polar body, which may divide again at MII (Zanders and Malik 2015; Dean 2016). Fertilization triggers completion of MII with sister chromatids segregating to the second polar body and the oocyte (Zanders and Malik 2015; Dean 2016).
Regular and reversed meiosis. Chromatids from alternative chromosomes are marked in red and black. MI-reductional (MI-Red) precedes MII-equational (MI-Eq) in regular meiosis, whereas in reversed meiosis, MI-Eq occurs first (Zanders and Malik 2015; Dean 2016). MII arrests at metaphase until sperm entry at fertilization. PB, polar body. Recombination is not included in these scenarios, although normally ∼40% of chromatids at the haploid phase have a crossover (Singer et al. 2004).
Recently, “reversed meiosis” was reported after induced ovulation in humans (Ottolini et al. 2015). Reversing the order of meiotic divisions during oogenesis, with the equational division occurring at MI rather than MII, leaves the secondary oocyte heterozygous for marker genes at fertilization (Figure 2). The genetics of fertilizing sperm could then bias MII segregation with one chromosome preferentially remaining in the oocyte and the other segregating to the second polar body. Although formal tests have yet not been reported in mice, the study by Agulnik et al. (1993) is consistent with reversed meiosis. Interestingly, genes such as Ago2 and Ppp2cb that control chromosome segregation in mitosis are strong candidates for determining the sequence of meiotic divisions (Pek and Kai 2011a,b; Huang et al. 2015), both of which are expressed in germ cells (Su et al. 2012a). Whether meiotic reversal is adaptive or is an anomaly resulting from reduced gene dosage and physiological stresses such as induced ovulation is unclear.
The alternative hypothesis that meiosis is conventional cannot account for fertilization bias in intercrosses but not backcrosses, because completion of MI prior to fertilization makes the reductional division independent of the mating partner and its fertilizing sperm. Bias would therefore be evident in both backcrosses and intercrosses.
Polyamines
The next question concerns the links between folate metabolism and gamete function. Although anomalies in folate metabolism are associated with reduced fertility (Forges et al. 2008; Mora-Esteves and Shin 2013; Aarabi et al. 2015; Ly et al. 2017), there is as yet no evidence for genotype-specific effects on haploid gametes. The challenging task then is to find an explanation connecting folate metabolism with TRD, given that transcription and translation are limited but not absent during spermatogenesis (Dadoune et al. 2004; Shima et al. 2004; Hammoud et al. 2009; Licatalosi 2016; Schuster et al. 2016; Zuo et al. 2016; Zhang et al. 2017).
A strong argument can be made that polyamine metabolism is the missing link. Most of this evidence involves spermatogenesis; the evidence for oogenesis is meager by comparison. Polyamines such as spermine, spermidine, putrescine, and cadaverine are low-molecular weight organic molecules that have at least two amino groups. They are present in all cells and most are associated with nucleic acids (Matthews 1993; Wallace et al. 2003; Persson 2009; Kahana 2016). They are involved in transcription, translation, histone modifications, autophagy, apoptosis, and many other molecular, cellular, and physiological activities (Hesse and Hoefgen 2003; Wallace et al. 2003; Alhonen et al. 2009; Lefevre et al. 2011; Igarashi and Kashiwagi 2015; Miller-Fleming et al. 2015).
Given their interdependent roles in essential molecular, epigenetic, and cellular functions, the polyamine and folate pathways are highly conserved and highly regulated, from yeast to plants and mammals, with sometimes overlooked roles in gametogenesis and fertility (Rubinstein and Breitbart 1991; Wolukau 2004; Roje 2006; Lefevre et al. 2011; Bauer et al. 2013; Aloisi et al. 2016; Chen et al. 2016). The substrates for polyamine synthesis are ornithine, which is derived from arginine and proline, and S-adenosylmethione (SAM), which is part of the homocysteine cycle in folate metabolism (Figure 3). SAM is the methyl donor for all methylation reactions for nucleic acids (DNA and RNAs, including tRNAs), proteins (including histones), lipids, and other molecules (Matthews 1993; Stover 2011; Salbaum and Kappen 2012; Gueant et al. 2013). Moreover, by utilizing acetyl CoA, polyamine catabolism magnifies the effects of folate deficiency (Figure 3). Acetyl CoA is the substrate for synthesizing betaine, which is the alternative methyl donor to synthesize SAM from homocysteine. Despite the essential role for methylation in many molecular and cellular functions, cells preserve polyamine metabolism at the expense of methylation, at least under in vitro and in vivo stress conditions (Kramer et al. 1987, 1988; Sun et al. 2002; Bistulfi et al. 2009, 2010).
Folate and polyamine metabolic pathways. [Adapted from references Carmel and Jacobsen (2001) and Wallace et al. (2003)]. Enzymes and metabolites that affect fertility, haploid gamete function, or embryonic stem cell (ESC) pluripotency are highlighted in red (see text for background and references). Genes that encode relevant enzymes are shown in italics with a bar showing their site of action. Blocked bars show proteins with inhibitory effects. Folate is the primary methyl donor; choline and betaine are alternative methyl donors. Amd1 links folate–homocysteine metabolism with polyamine metabolism by converting S-adenosylmethionine to decarboxylated S-adenosylmethionine, which with putrescine produces spermidine. Odc1 converts ornithine to putrescine and is the heavily regulated rate-limiting step in polyamine synthesis. Catabolism of spermine and spermidine consumes acetyl CoA, which also serves as a substrate for choline and betaine synthesis. 5-THF, 5-tetrahydrofolate; 5,10-THF, 5,10-methylenetetrahydrofolate; Amd1, S-adenosylmethione decarboxylase 1; Odc1, ornithine decarboxylase 1; THF, tetrahydrofolate.
Functionality of haploid gametes depends heavily on polyamine metabolism in many species. In Arabidopsis, the MAT3 SAM synthase gene is expressed in pollen, with mutants showing reduced pollen tube growth and seed set as well as changes in polyamine biosynthesis, tRNA levels, and histone methylation (Chen et al. 2016). In humans, polyamine deficiency results in infertility, which can be corrected with SAM or polyamine supplementation (Shohat et al. 1990; Rubinstein et al. 1995; Calandra et al. 1996; Lefevre et al. 2011). In cattle, spermine is essential for acrosomal function (Rubinstein et al. 1995). Overexpression of ornithine decarboxylase 1 (ODC1), the rate-limiting step in polyamine synthesis, also leads to infertility (Qian et al. 1985; Halmekyto et al. 1991; Kilpelainen et al. 2001; Tokuhiro et al. 2009). Several polyamine genes are expressed in haploid gametes (Figure 3) (Lefevre et al. 2011). ODC antizyme 3 (OAZ3) is a testis-specific inhibitor of ODC1 (Qian et al. 1985; Tosaka et al. 2000; Kilpelainen et al. 2001; Ike et al. 2002; Tokuhiro et al. 2009). Deficiency leads to sperm that cannot fertilize (Kilpelainen et al. 2001; Tokuhiro et al. 2009). Antizyme inhibitor 2 blocks the inhibitory effects of OAZ3 in haploid cells (Lopez-Contreras et al. 2009) as well as ODC1 overexpression (Suzuki et al. 1994). Spermidine/spermine N1-acetyltransferase and OAZ1 are differentially expressed with folate supplementation in the LRP6 mouse NTD model (Gray et al. 2010; Nakouzi and Nadeau 2014). OAZ1 is also differentially expressed in sperm from folate-deficient mice (Lambrot et al. 2013). Finally, the primary inputs for polyamine metabolism, namely the enzyme SAM decarboxylase 1 (AMD1), which catalyzes the conversion of SAM to decarboxylated SAM, and the amino acids arginine, ornithine, and proline, are critical for pluripotency control in embryonic stem cells (ESCs) (Figure 3) (Pierce et al. 1990; Nishimura et al. 2002; Casalino et al. 2011; Zhang et al. 2012; Zhao et al. 2012; D’Aniello et al. 2015; Lee et al. 2016). Given the close lineage relationships between germ cells and embryonic stem cells (ESCs) (Zwaka and Thomson 2005; Bustamante-Marin et al. 2013), similar effects in gametes would not be surprising. Polyamine metabolism is therefore a compelling candidate for biased fertilization, given its established impact on haploid gametes and its strong connection to folate metabolism.
Genetic and Reproductive Considerations
Sperm in the epididymis, gametes in the oviduct
The issue here is whether reproductive organs and sperm–egg recognition provide opportunities for gametic selection based on haploid effects (Holt and Fazeli 2015). Semen, as well as oviductal and uterine fluids, contain various proteins, nucleic acids, and small molecules that provide a chemically appropriate environment for sperm maturation and capacitation, for fertilization, and for implantation of genetically nonself embryos (Nieder and Macon 1987; Vinijsanun and Martin 1991; Wolfner 2002; Huang et al. 2004; Waberski et al. 2006; Holt and Fazeli 2010; Lane et al. 2014; Johnson et al. 2015; Jodar et al. 2016). Because mating usually occurs prior to ovulation, sperm can be in the oviduct hours before ovulation, often needing just minutes from ejaculation to arrive in the oviduct (Lewis and Wright 1935). Sperm and eggs induce gene expression changes in the oviduct that alter the biochemistry of oviductal fluids (Fazeli et al. 2004) and that in turn restrict sperm access (Holt and Fazeli 2015, 2016). Millions of sperm are released at ejaculation but usually < 100 reach the oviduct (Braden and Austin 1954). Not only do sperm compete for access to eggs (Clark et al. 1999; Clark 2002), the oviduct can select sperm based in part on their genetic content, including sex chromosomes (X or Y) (Alminana et al. 2014; Holt and Fazeli 2016) and on chromatin stability (Ardon et al. 2008).
Signaling before contact between sperm and egg has obvious advantages for haploid gametes. There are reports that small, transiently expressed peptides attract a minority of sperm that are only briefly responsive, presumably representing the 10% of sperm that are appropriately capacitated (Eisenbach and Ralt 1992; Eisenbach 1999). Thermotaxis (Bahat and Eisenbach 2010), chemotaxis (Holt and Fazeli 2015, 2016), and signaling molecules including olfactory receptors (Flegel et al. 2015), trace-amine-associated receptors (Zucchi et al. 2006; Chiellini et al. 2012), cannabinoid receptors (Agirregoitia et al. 2010), and calcium receptors (Dell’Aquila et al. 2006; De Santis et al. 2009; Mendoza et al. 2012), have been reported in gametes and gonads. These could play a role in sperm–egg attraction, but the evidence is limited both about mechanisms and especially about the impact of genetic variation on molecular interactions.
In Drosophila and other species, success of fertilizing sperm depends on the genetics of both the male and female mating partners, suggesting that ligand–receptor interactions are critical (Braden 1958; Chow et al. 2010). But how genetic variation affects affinity and signaling in these ligand–receptor pairs is largely unknown. Once sperm penetrate the glycoprotein coat surrounding the egg, the sperm and egg must recognize each other and fuse membranes. The zona pellucida hardens irreversibly to prevent polyspermy (Braden et al. 1954). Presumably fertilization bias must occur before membrane fusion and zona pellucida hardening, otherwise the conceptus would either persist as a viable embryo, or be lost with a corresponding reduction in litter size. Recognition between sperm and egg is based on ligand–receptor binding between Izumo1 (sperm) and Juno (egg); absence of either protein leads to infertility (Bianchi et al. 2014; Aydin et al. 2016; Bianchi and Wright 2016; Ohto et al. 2016). Although Juno, which is a member of the folate receptor family, no longer binds folate (Bianchi et al. 2014), residual functions might still depend on folate levels. Unexplained anomalies between other ligand–receptor pairs such as loss of one but not the other member of the pair resulting in infertility suggests that models of sperm–egg recognition remain incomplete (Cho et al. 1998). In at least one instance, union of sperm and egg brings together two distinct proteins that act together as a dimer to suppress mutagenesis in early embryos (Lord and Aitken 2015).
Several approaches have been employed to define mechanisms for sperm–egg recognition. For example, ENU mutagenesis has been used to find genes controlling gametogenesis and fertilization (Hrabe de Angelis et al. 2000; Clark et al. 2004; Sakuraba et al. 2005). Although many of these mutated genes affect germ cell biology and meiosis, genes affecting fertilization were not found (Furnes and Schimenti 2007). Fortuitously, our genetic studies of epigenetic TGCT risk and NTD dietary response discovered several such genes, thereby enabling new genetic approaches for studying mechanisms of gamete function at fertilization.
Haploid effects
Functional effects confined to haploid gametes are critical for fertilization bias. The four products of meiosis produced in males usually have an equal chance of fertilization, in part because syncytia provide small intercellular channels through which developing sperm share molecules, thereby minimizing the impact of haploid effects (Greenbaum et al. 2011). However, gametes carrying a t-haplotype, which are a classic example of TRD in the mouse (Dunn 1957; Lyon 2003), produce two critical elements that control sperm motility in a haploid-specific manner. Gametes that carry a t-haplotype produce molecules that pass through syncytial bridges to hyperactivate Rho signaling in both t- and +-bearing sperm, thereby compromising their motility. However, t-bearing sperm also produce a protective, haploid-specific variant of SMOK1 (sperm motility kinase) that does not pass through bridges and protects t-bearing but not +-bearing sperm, thereby providing a motility advantage to t-bearing sperm (Herrmann et al. 1999; Veron et al. 2009; Bauer et al. 2012). Because such effects are intrinsic to heterozygous males and would be found in both backcrosses and intercrosses, a corresponding effect would need to operate in females to be a sufficient explanation for biased fertilization.
Competition between diploid and haploid phases
Discoveries about biased fertilization are relevant to theories concerning sexual antagonism, namely the contrasting priorities between diploid organisms and their haploid gametes. Diploids strive for reproductive success vs. other diploids, whereas haploid gametes compete with each other for successful fertilization. Because gametic competition could reduce parental fertility, diploid cells should seek to reduce functional differences among gametes by limiting their transcription and translation. As long as gametes are functionally equivalent, diploids have the advantage over their haploid gametes. It is noteworthy then that many of the genes that bias fertilization affect aspects of RNA biology that impact translation (Box 1). Partial loss-of-function in mutant heterozygotes might enable phenotypic differences among gametes, leading to gamete competition and biased fertilization.
The window for bias
The length of the haploid phase in spermatogenesis and oogenesis places conditions on the mechanisms that might contribute to biased fertilization. In mammals, the haploid phase is long in males, with the MI division arrested during embryonic development and recommencing at puberty. Spermatogenesis then continues throughout reproductive life. The haploid phase begins with completion of meiosis in the testis and continues as they mature and capacitate, and while they pass through the epididymis, vas deferens, urethra, and uterus to the site of fertilization in the oviduct, a period that can last several days. By contrast, the haploid phase is remarkably short in females, lasting only from completion of MII, which fertilization triggers, until female and male pronuclei fuse. This brief window raises the likelihood that fertilization drives the bias in oocytes.
Centromeric- or gene-specific effects
Sometimes, TRD results from preferential segregation of centromeric elements that guide chromosome movement and segregation during karyokinesis (Lyttle 1989; Pimpinelli and Dimitri 1989; Pardo-Manuel de Villena and Sapienza 2001; Axelsson et al. 2010; Chmatal et al. 2014; Wei et al. 2017). Genes that are closely linked to the centromere would also show TRD with the degree of distortion depending on the recombination distance from the centromere and with genes located 50 cM or more showing normal 1:1 segregation. However, several genes that are located far from the centromere show substantial TRD and overall there is little evidence that TRD declines as a strong function of recombination distance from the centromere (Figure 4). However, the propensity of recombined chromatids to segregate preferentially to oocytes during reversed meiosis precludes a firm conclusion (Ottolini et al. 2015). Finally, the dispersed chromosomal locations for these genes (Figure 4, Table 1, and Table 2) argue against a selfish gene complex, as is found frequently with other TRD systems (Lyttle 1993; Lyon 2003). Together, these observations are consistent with gene-specific TRD rather than hitch-hiking resulting from close linkage to centromere-driven elements.
The difference (%) from expectations was calculated as the absolute value of the deviation between observed and expected frequency of m/+ heterozygotes among heterozygous and wild-type (WT) segregants. Difference provides a direct measure of the departure from expectations. Recombination distances are from the Mouse Genome Database (informatics.jax.org). The correlation between deviation and distance was modest (r = −0.54) and accounted for only 29% of the variation (r2), suggesting that centromeric drivers, if any, had modest effect on transmission ratios. Limiting the analysis to data points < 50 cM yielded r = −0.39 and r2 = 0.15.
Litter size
Although many factors such as the number of fertilized eggs, pre- and postimplantation mortality, and uterine capacity can affect litter size, various evidence shows that the primary determinant is the number of ovulated eggs. Selection for larger litter size increases ovulation rate (Falconer 1960; Joakimsen and Baker 1977; Bakker and Politiek 1978; Durrant et al. 1980), while selection for increased ovulation rate results in larger litters (Bradford 1969). Inbreeding reduces litter size because fewer eggs are ovulated (McCarthy 1967). Unilateral ovariectomy reduces litter size by 50% (Vinijsanun and Martin 1991), arguing against eggs held in reserve to compensate for failed fertilizations and embryo loss. In parallel, eggs mature within ovarian follicles that rupture to release eggs at ovulation, with the number of growing follicles determined by host genetics long before and independent of fertilization. Thus, the number of ovulated eggs available for fertilization appears to be the primary determinant of litter size. Whether any of these 12 mutant genes, together with regular or folate-supplemented diets, affect the number of ovulated eggs remains to be determined, although any such effects should be found in both backcrosses and intercrosses.
Lethality
Three lines of evidence argue against lethality as the explanation for departures from Mendelian expectations (Figure 1). First, litter size is not reduced in intercrosses vs. backcrosses for any genes that bias fertilization (Table 1, Table 2, and Table 3). Again, if mutant homozygotes show embryonic lethality, then litter size should be reduced by 25%, and if half the heterozygotes are also missing, then litter size should be further reduced to 50%. Second, in some cases [Dnd1 (Zechel et al. 2013), Ddx1], surveys at E3.5 failed to find mutant homozygotes or dead embryos (Hildebrandt et al. 2015), and Pum1 (Zhang et al. 2015). Third, loss of particular genotypes is not based on their inherently deleterious nature, given that mutant heterozygotes that are missing among intercross progeny are found in expected numbers among backcross progeny (Table 1), in reciprocal crosses (Table 4), or on alternative folic acid diets (Table 2). Thus, negative selection against embryos having particular genotypes does not seem to be involved.
Perspectives
Biased fertilization may be a previously unrecognized manifestation of genes that control TGCT susceptibility, epigenetic inheritance, and related germline abnormalities. The first evidence for fertilization bias was discovered during work on the control of inherited TGCT risk in the 129 family of inbred strains (Lam et al. 2007; Heaney et al. 2008; Nelson et al. 2012). These TGCTs are models for several classes of TGCTs in humans (Bustamante-Marin et al. 2013; Rijlaarsdam and Looijenga 2014). Like gametes, many originate from the germ cell lineage during fetal development (Kleinsmith and Pierce 1964; Stevens 1967). The germline origin of TGCTs and gametes suggests that they may share similar vulnerabilities to genetic and environmental perturbations. Perhaps no other lineage undergoes such dramatic transitions in developmental potential and with such profound implications for health, fertility, and perpetuation of the germline across generations (Jablonka and Lamb 1999). The germline is unipotent until fertilization, when it transitions to pluripotent embryonic cells that in turn differentiate into specialized somatic and germ cell lineages (Saitou and Yamaji 2012; Hackett and Surani 2014; Reik and Surani 2015). Various mechanisms preserve the genomic integrity and developmental capacity of the “mother of all stem cell lineages” (Donovan 1998) by repairing DNA defects (Khurana et al. 2010) maintaining cellular conditions (Gemble et al. 2015), suppressing transposon activity (Bourc’his and Bestor 2004; Bourc’his and Voinnet 2010; Senti and Brennecke 2010; Barau et al. 2016; Wylie et al. 2016a,b), and programming epigenetic state (Buecker et al. 2014; Oliveros-Etter et al. 2015). Failure of pluripotency control can lead to precocious differentiation of germ cells (Hargan-Calvopina et al. 2016; Sanchez et al. 2016), spontaneous transformation of germ cells during fetal development (Kleinsmith and Pierce 1964; Stevens 1967; Bustamante-Marin et al. 2013; Ferguson and Agoulnik 2013), infertility (Skakkebaek et al. 2001, 2016; Carouge et al. 2016), and other reproductive disorders (Ferguson and Agoulnik 2013). Dysfunction can also lead to inherited epigenetic changes that affect risk for TGCTs and urogenital abnormalities in offspring and later generations in the absence of genes that originally triggered these transgenerational effects (Lam et al. 2007; Heaney et al. 2008; Nelson et al. 2012; Carouge et al. 2016). Genetic anomalies together with dietary, hormonal, and other environmental influences may induce germline dysfunctions that manifest as conventional genetic effects, unconventional epigenetic inheritance, and now perhaps preference for particular gamete combinations at fertilization.
Acknowledgments
I thank Rosaline Godbout and her group for generously sharing litter size and segregation data from Hildebrandt et al. (2015) and for sharing her thoughts about their study and the interpretation proposed here. I also thank David Crews, Mary Ann Handel, and especially Hamish Spencer and Monika Ward for comments on a draft of this paper. Andy Clark, Aimee Dudley, John Eppig, Fernando Pardo de Villena, Jasper Rine, Jesse Riordan, Carmen Sapienza, and Patrick Stover provided helpful suggestions, insights, and speculations. I thank the reviewers for their many helpful comments. The following National Institutes of Health grants supported this work: Eunice Kennedy Shriver National Institute of Child Health and Human Development Pioneer Award DP1 HD075624, National Cancer Institute CA75056, and National Institute of Neurological Disorders and Stroke NS058979.
Footnotes
Supplemental material is available online at www.genetics.org/lookup/suppl/doi:10.1534/genetics.117.300109/-/DC1.
Communicating editor: J. C. Schimenti
- Received April 12, 2017.
- Accepted July 10, 2017.
- Copyright © 2017 by the Genetics Society of America