- Open Access
Gene changes may minimize masculinizing and defeminizing influences of exposure to male cotwins in female callitrichine primates
Biology of Sex Differencesvolume 7, Article number: 28 (2016)
Sexual differentiation in female mammals can be altered by the proximity of male littermates in utero, a phenomenon known as the intrauterine position effect (IUP). Among simian primates, callitrichines (marmosets and tamarins) are likely candidates for IUP, since they exhibit obligate dizygotic twinning and fetuses share extensive vascularization in utero. In this paper, we determined whether female reproductive parameters are altered by gestating with a male twin and evaluated changes in genes associated with anti-Müllerian and steroid hormones in twinning callitrichine primates.
We assessed the impact of gestation with male cotwins on reproductive performance and survivorship in female marmosets (Callithrix) and lion tamarins (Leontopithecus), contrasting measures for females gestated with one or more littermates (M+) or no male littermates (0M). We compared targeted coding regions for genes involved in steroidal and anti-Müllerian hormone mediation of sexual differentiation for representatives of twinning callitrichines (Callithrix, Saguinus, and Leontopithecus) with closely related New World primates that produce single births (Saimiri and Callimico).
IUP effects in females were absent in female callitrichine primates: age at first ovulation, average litter size, and the proportion of stillborn infants, and lifetime survivorship did not differ between M+ and 0M females. We documented multiple nonsynonymous substitutions in genes associated with steroid synthesis, transport, and cellular action (SRD5A2, CYP19A1, SHBG, and AR) and with anti-Müllerian hormone (AMH and AMHR2) in callitrichines. In the only callitrichine to produce single infants (Callimico), two genes contained nonsynonymous substitutions relative to twinning callitrichines (CYP19A1 and AMRHR2); these substitutions were identical with nontwinning Saimiri and humans, suggesting a reversion to an ancestral sequence.
In spite of a shared placental vasculature with opposite-sex twins throughout embryonic and fetal development, female callitrichine primates gestated with a male cotwin exhibit no decrement in reproductive performance relative to females gestated with female cotwins. Hence, IUP effects on female reproduction in callitrichines are modest. We have identified mutations in candidate genes relevant for steroid hormone signaling and metabolism, and especially in AMH-related genes, that are likely to alter protein structure and function in the callitrichines. These mutations may confer protection for females from the masculinizing and defeminizing influences of gestating with a male cotwin.
Sexual differentiation in mammals represents a complex temporal cascade of genetic, endocrine, and neuroendocrine processes. Sexual differentiation toward male phenotypes is associated with the presence of the Sry gene, whose protein products result in the development of testes from undifferentiated gonadal tissue . The testes are endocrinologically active early in fetal life, secreting testosterone and other androgenic metabolites, resulting in the masculinization of secondary sexual characteristics, including the brain. The autosomal gene Sox9 also plays an important role in sexual differentiation, primarily through the upregulation of anti-Müllerian hormone (AMH) production from testicular Sertoli cells [2, 3]. AMH production results in the regression of the Müllerian ducts and their derivatives, including the fallopian tubes, uterus, and upper vagina, suppressing the development of female-typical accessory reproductive structures . A host of environmental variables work in concert with these genetic and endocrine factors to canalize the process of sexual differentiation, and the sex and location of littermates in utero have been identified as important variables that can modulate the process of differentiation. The intrauterine position effect (IUP) is a well-established phenomenon in litter-bearing species, with the majority of research conducted in rodents. Female fetuses that are gestated in closer proximity to male fetuses become partially masculinized, presumably as a consequence of the vascular diffusion of endocrine secretions of male origin via the fetoplacental circulatory system [5, 6].
Both androgens and AMH are active in mediating IUP effects. Androgen-related masculinization of females via proximity to male littermates is supported by both morphological and behavioral evidence. Female rodents that occupy uterine positions between two male fetuses (2M) tend to have masculinized external genitalia and delayed vaginal opening and also exhibit significant delays in puberty, offspring production, and lifetime reproductive performance, relative to females that were not gestated adjacent to male littermates (0M; ), although these effects are not universal across species (e.g., [7–9]). Female rodents that gestate between two males also display higher rates of male-typical behavior than females that gestate between two females, including more territorial behavior, enhanced aggression and male-typical sexual behavior in response to low doses of androgens, enhanced spatial memory, and a greater tendency to disperse from natal groups (see review in ). The involvement of AMH in IUP-related defeminization is demonstrated by the freemartin syndrome  in cattle and other domestic ungulates. Reproductive sterility in females is often a consequence of gestating with a male cotwin, and interfetal transfer of male-derived AMH results in the regression of female accessory reproductive system [11–13]. Thus, there are at least two independent endocrine mechanisms whereby fetal males can alter the timing and trajectories of sexual differentiation in female littermates.
IUP effects have not been evaluated extensively in the simian primates (monkeys and hominoids), since single births are normative for most primates. However, the callitrichine primates (marmoset and tamarin monkeys, family Cebidae) provide a unique test case for IUP effects in anthropoid primates. First, all callitrichine species except one (Callimico) have obligate “litters,” producing modal litter sizes of two offspring [14–16], and candidate genes for twinning have been identified . Since callitrichine twins are dizygotic, roughly 50 % of births involve mixed-sex male and female cotwins. Callitrichines also exhibit a host of embryonic, fetal, and fetoplacental features that facilitate the transfer of signaling molecules between cotwins, including hormones involved in sexual differentiation. Following implantation, embryos have a period of delayed development , when embryos develop a fused placental chorion and develop extensive vascular anastomoses by day 29 of embryonic development and persist throughout fetal development [18–20]. The efficacy of intertwin transfer via these shared vascular connections is confirmed by the observation that an estimated 95 % of pregnancies that result in twins are chimeric in hematopoietic tissues [21, 22] and possibly other cell lines . Like freemartin cattle, female callitrichines gestated with male cotwins possess cells with both XX and XY chromosomal elements [24–27] as a consequence of the transfer of embryonic stem cells between twins. However, unlike cattle, chimeric female callitrichine primates are not sterile.
Therefore, given the extensive vascular communication between callitrichine twins, one would expect that masculinizing signals from male fetuses are highly likely to be diffused to female cotwins, potentially leading to masculinization and hence reduced reproductive capacity in females. However, only two studies have explicitly addressed this issue in callitrichines, with differing outcomes. Ardito et al.  contrasted female reproductive performance in marmosets that possessed both XX and XY lymphocytes (and hence gestated with a male cotwin) relative to females with only XX lymphocytes. A higher proportion of females with XX/XY karyotypes were among the colony’s best reproductive females (measured by number and survivorship of offspring) than XX-only females. In contrast, Rutherford et al.  demonstrated that female marmosets that gestated with a male cotwin had a greater lifetime pregnancy loss than females that gestated with a female cotwin, but rates of pregnancy loss were not adjusted for lifetime reproductive tenure.
The present study had two aims. First, we quantitatively assessed the potential reproductive and survivorship costs to females that shared a uterine environment with at least one male in three species of callitrichines, using multiple measures of reproductive performance. We addressed this question by analyzing long-term captive breeding records in female marmosets (Callithrix spp.) and golden lion tamarins (Leontopithecus rosalia) born in litters with no known male cotwins (0M) and those born in litters with at least one male cotwin (M+). Second, we conducted a comparative analysis of coding regions for candidate genes implicated in sexual differentiation [30, 31] in several genera of twinning callitrichine primates. These included genes that code for aromatase (CYP19A1), 5α-reductase (SRD5A2), the nuclear androgen receptor (AR), the glycoprotein that serves as a carrier molecule for sex steroids (SHBG), anti-Müllerian hormone (AMH), and the type II anti-Müllerian hormone receptor (AMRHR2). We contrasted these sequences with those from the closely related squirrel monkey (Saimiri, family Cebidae) and the callitrichine Callimico, both of which produce single offspring. If reproduction in female callitrichines is negatively impacted by gestation with a male cotwin, we predicted few large impact nonsynonymous substitutions (NSs) in gene sequences in the twinning callitrichines relative to related species that produce single offspring. However, if female reproductive development and function in twinning callitrichines are not impacted by gestation with a male cotwin, we would expect significant differences in nucleotide and predicted amino acid sequences for these genes in twinning callitrichines in ways that suggest protective mechanisms against masculinization of females that gestate with male cotwins.
Reproductive performance and offspring survivorship in females born to mixed-sex (M+) vs. same-sex (0M) litters were evaluated for two genera of callitrichine primates: marmosets (Callithrix spp.) and golden lion tamarins (Leontopithecus rosalia). Records for marmosets were derived from the long-term breeding program at the University of Nebraska at Omaha (UNO) and included data on two closely related species: C. kuhlii and C. geoffroyi. We derived our data for lion tamarins from the international studbook, a digital database on captive breeding in this species. For both marmosets and lion tamarins, females were included in our analyses if breeding records contained information regarding sex composition of litters from which they derived, dates of pairing in a reproductive context, and dates of birth for offspring produced by these females. For each species, the number of females used in the analyses, their distribution as a function of sex composition of the litter into which they were born, and the average number of litters produced throughout their lifetime is indicated in Table 1. We did not genotype 0M females for XX/XY chimerism, and since fetal loss can occur in callitrichines , some 0M females may have had exposure to male cotwins at some point in embryonic and fetal development. However, litter-size reduction verified by ultrasound in marmosets has a low incidence (7.2 % of pregnancies; ), thus potentially misclassified 0M females represent a very small proportion of the overall sample.
Reproductive and life-history measures
Our predictor variable for female reproduction was derived from the sex composition of the litter into which breeding females were born. We classified breeding females into all-female litters (0M) or females that had undergone gestation with at least one male fetus (M+). We tested four reproductive parameters to differentiate the effects of litter sex ratios on prenatal and postnatal female reproductive performance. The first three measures evaluate the impact of exposure to male siblings in utero on the development of ovarian function, fetoplacental factors affecting litter size, and the adequacy of gestational ability. Latency to first parturition: length of time (in days) between pairing of the female with an unrelated, reproductively intact adult male and first parturition. We excluded pregnancies in females that resulted in preterm birth of offspring, and age at pairing with a breeding male was not available for 13 female Leontopithecus (7 0M and 6 M+). Age at first ovulation: for a small number of female Callithrix (0M: n = 6; M+: n = 10), we had archived urine samples that allowed us to definitively determine the age at first ovulation based on urinary pregnanediol concentrations. Details on endocrine methods and criteria for first ovulation can be found in . Litter size: we calculated the average number of full-term fetuses (both stillborn and alive) produced per parturition event. Females that had given birth to at least one litter were included in this analysis, and litters of all sizes (range = one to five offspring) were included in these calculations. Proportion of stillborn infants: we calculated the proportion of the stillborn infants (i.e., offspring that were physically full-term but who were dead upon discovery and showed no signs of maternal abuse or physical damage) from the total number of full-term fetuses produced per parturition event. Female longevity: we collected data on age at death for females from both genera and the age of currently living females for the use in a censored survivorship analysis (see below).
For latency to first parturition, age at first ovulation, litter size, and proportion stillborn infants, we contrasted means for 0M breeding females with M+ females using Welch’s t test. We used the Cox proportional hazards model to test for differences in lifetime survivorship of 0M and M+ females.
We identified a cluster of genes that have been demonstrated to be critical in the process of normative sexual differentiation in mammals. Four genes are involved in steroid metabolism, receptor function, and transport [31, 35]: CYP19A1, a locus coding for the aromatase enzyme that converts C19 androgens into C18 estrogens; SRD5A2, which codes for the enzyme 5α-reductase that serves to convert testosterone into dihydrotestosterone; AR, which codes for nuclear androgen receptors; and SHBG, which codes for sex hormone-binding globulin, a transport glycoprotein that binds to sex steroids and renders them biologically inert. Two genes associated with the early embryonic differentiation of internal genitalia and associated structures were also analyzed: AMH, which codes for anti-Müllerian hormone (AMH), a protein that inhibits the embryonic development of female-typical Müllerian ducts that are destined to become uterine structures, and AMRHR2, which codes for the type II AMH-receptor that mediates Müllerian duct regression in the presence of AMH.
Publically available exomic coding sequences for each gene were aligned firstly across a wide variety of primate species, including human, gibbon, rhesus monkey, common marmoset (C. jacchus), squirrel monkey (Saimiri boliviensis), tarsier, bushbaby, and mouse lemur; with outgroup species constituting cattle and sheep. All published sequence data were accessed from UCSC Genome Browser, NCBI, or Ensembl. Sequences were aligned in MEGA 6.0 using ClustalW. Regions of the candidate genes that contained Callithrix-specific NSs relative to comparison species were selected for further sequence analysis. Specific regions identified by these analyses for each candidate gene and details of PCR amplification are noted in Additional file 1.
We generated de novo nucleotide sequences via traditional Sanger sequencing for the candidate genes in five additional callitrichine species, including C. kuhlii, C. penicillata, L. rosalia, Saguinus midas, and C. goeldii. The institutional source of tissue or extracted DNA for each species is listed in Additional file 2. We sequenced selected coding regions for the candidate genes from at least two unrelated individuals for each species, one male and one female. Our primary comparisons of interest were (1) the contrast between predicted amino acid sequences for the candidate genes between the twinning callitrichines and the closely related singleton-producing squirrel monkey (Saimiri; both members of the Family Cebidae) and (2) between the twinning callitrichine species and Callimico, the only callitrichine primate that produces single offspring. Genomic DNA was extracted using the DNeasy Blood and Tissue Kit (Qiagen) following the manufacturer’s protocol. DNA yield from the extractions was quantified using Nanodrop® (Thermo Scientific), and DNA quality was measured using 0.8 % agarose gel electrophoresis. The same sets of PCR primers were used to amplify coding regions in our candidate genes for all species, and the details of PCR protocols are found in Additional file 2. PCR products were purified using QIAquick PCR Purification Kit (Qiagen) and sequenced directly using an Applied Biosystems (ABI) 3730 48-capillary electrophoresis DNA analyzer, in the High-Throughput DNA Sequencing and Genotyping Core Facility at the University of Nebraska Medical Center. Sequences generated de novo for this project are available at NCBI [GenBank: 1882226].
For all New World primates, amino acid substitutions were classified as NSs (and hence potentially relevant for altered protein function) if the nucleotide sequence predicted variable amino acids at each position among species. We further classified NSs as conservative or radical according changes in polarity, charge, and volume: NSs with a change in one or more categories were classified as radical, while NSs with no changes in the three categories were classified as conservative . To estimate the impact of the identified NSs in each gene, we submitted comparative sequences of predicted amino acids in Saimiri and the callitrichine primates to an analytical algorithm (PROVEAN; ) designed to assess the potential significance of the callitrichine NSs for protein structure, function, and putative phenotypic changes. PROVEAN employs an alignment algorithm to predict the functional impact of protein sequence variations (e.g., single and multiple amino acid substitutions, deletions, and insertions) and can contrast variation between any organisms (in our case, variation between Saimiri and Callithrix) and is not limited to contrasts with human sequences (e.g., SIFT, Polyphen: [38, 39]).
Exposure to a male fetus during embryonic and fetal development had only modest effects on subsequent reproductive performance in female marmosets and lion tamarins (Table 2). Age at first ovulation determined from urinary pregnanediol profiles did not differ between 0M and M+ Callithrix females. Overall, female marmosets gave birth to larger litters (2.52 ± 0.15 offspring, mean ± SEM) than female lion tamarins (1.84 ± 0.003 offspring), but average litter sizes did not differ between 0M and M+ females in either species. Pregnancy outcomes were also unaffected by litter sex ratios, as the proportion of stillborn infants across the females’ reproductive lifetime did not differ between 0M and M+ females. For both genera, the latency to first parturition after pairing with an appropriate male breeding partner was affected by litter sex ratio. In Leontopithecus, M+ females took longer to produce their first litter after pairing than did 0M females (d = 0.37).
Female lifetime survivorship
Overall lifetime survivorship of females in both taxa was unaffected by the presence or absence of a male cotwin. Figure 1 portrays survivorship curves for females, and the analyses (both censored and uncensored lifetime measures) revealed no significant differences in the pattern of survivorship between 0M and M+ females for either species (Leontopithecus: χ 2(1) = 2.38, p = 0.12; Callithrix: χ 2(1) = 0.11, p = 0.74).
In the targeted exonic regions of the six candidate genes, we documented 53 NSs in one or more genera of callitrichine primates, relative to Saimiri. Substitutions in each gene are presented in Fig. 2 and a summary of these NSs in Table 3. Further details on specific protein positions, amino acid substitutions, and predicted functional consequences across all species sampled can be found in Additional file 3. Of the 53 NSs, amino acid substitutions at 26 positions were documented across all callitrichine genera. Each of the target genes contained at least two callitrichine-wide substitutions (range 2–8). We utilized PROVEAN  to identify NSs in callitrichines, relative to Saimiri, that were predicted to lead to significant alteration in protein structure/function. We also used SIFT  and PolyPhen  to identify likely deleterious NSs in callitrichines, relative to protein sequences in humans.
Each of the four genes associated with steroid metabolism, transport, and nuclear activation expressed at least three callitrichine-specific substitutions, and many amino acid positions varied across the callitrichines, relative to Saimiri. All amino acid substitutions in CYP19A1, SHBG, and AR were identified by PROVEAN as neutral substitutions with regard to predicted protein structure and function. However, P64S in SRD5A2 was identified as a substitution that would lead to altered protein structure/function (and was also targeted by SIFT as a deleterious change). We noted one substitution in CYP19A1 (T268A) that was present in all callitrichines except Callimico, and the predicted residue at this location for Callimico was T, the predicted residue for Saimiri and humans, but this NS was predicted to lead to a neutral change in protein structure/function.
Amino acid substitutions were more numerous in the AMH-related target genes than in the steroid-related candidate genes. For AMH, eight of 14 amino acid positions contained substitutions across all callitrichines. PROVEAN identified three substitutions predicted to significantly alter protein structure/function. These included a deletion (G40del) in Callithrix and two substitutions (T177M,R in Callithrix and Callimico, respectively; R473P in all callitrichines). Although there were only two callitrichine-wide substitutions in AMHR2 (positions 190 and 293), a deletion at positions 268–272 in Callithrix and Callimico and H289L in Callithrix and Leontopithecus was identified as producing significant alterations in protein structure/function. As with CYP19A1, a substitution (R176G) in AMHR2 was present in all callitrichines except Callimico, which expressed the same residue as Saimiri and human. Relative to steroid-relevant genes, then, NSs were more numerous for AMH-related genes, and these changes were more likely to lead to altered protein structure and/or function.
Given the features of dizygotic twinning and early shared intrauterine vascularization in most callitrichine primates, they would appear to be obvious candidates for the masculinization of females gestated with male cotwins, and consequently, there should be evidence of protection of female fetuses from masculinizing and defeminizing signals from male cotwins. With the exception of latency to produce the first litter, our data reveal that females incur little or no cost to reproduction or survivorship as a result of a shared uterine environment with a male cotwin. These results hint strongly at one or more protective mechanisms in callitrichines that might minimize masculinization and defeminisation in females from mixed-sex litters. In a suite of candidate genes that regulate two important endocrine cascades that mediate sexual differentiation (AMH and androgens), we documented a number of nonsynonymous substitutions in the coding regions for callitrichines that were different from human and macaque monkeys. Importantly, these NSs also differed from the squirrel monkey (Saimiri), which provides an important phylogenetically relevant contrast, since Saimiri belongs to the same Family (Cebidae) as callitrichines but produces only single offspring.
Latency to produce the first litter after pairing with a male was longer in M+ than 0M females in Leontopithecus but not in Callithrix, although the absence of a significant effect in Callithrix may be a consequence of a relatively small sample size and hence low statistical power. There appear to be species differences in the effects of gestating close to males on the onset and success of female reproduction. In mice, while many features of female morphology and behavior are masculinized among 2M females , 2M females do not appear to have delayed conception or reduced offspring survivorship, relative to 0M females [40, 41]. In pigs, however, females born into male-biased litters have higher reproductive failure during their first pregnancy and have reduced reproductive success over their first four litters, relative to females born in female-biased litters . The extent to which delayed conception in M+ callitrichine females represents a common trait among species in this taxon, or is unique to the reproductive biology and life history of Leontopithecus, will require analyses characterized by both larger sample sizes and additional species.
Our candidate genes were selected on two premises: (1) both AMH and androgens are critical in the process of sexual differentiation in mammals  and (2) mutations in these genes are associated with abnormalities in sexual differentiation, particularly undermasculinization and defeminization, in human populations . In our six candidate genes, we identified a total of 53 NSs leading to a change in predicted amino acids that differentiated one or more species of callitrichines from Saimiri, with 21 NSs that were common to all callitrichine genera sampled. AMH and SRD5A2 each had six callitrichine-specific substitutions. SHBG had the fewest amino acid substitution of the candidate genes, and only two of the four substitutions were present among all callitrichines. Across genes, the majority of substitutions represented a radical change in the physicochemical properties of the residue .
Multiple single nucleotide polymorphisms (SNPs) in androgen-relevant genes have been documented that alter the process of masculinization [43, 44]. Analyses of SHBG in female-masculinized spotted hyenas (Crocuta crocuta), relative to nonmasculinized hyenid species, reveal multiple substitutions that reduce production of SHBG, thereby increasing bioactive sex steroids in this species . To the extent that protection against masculinization for female callitrichines is mediated partly by modified SHBG, we would expect the NSs to increase either sex steroid affinity or production of the binding protein. SNPs that lead to reductions in efficacy or loss-of-function for enzymes involved in steroid biosynthesis have been documented for both CYP19A1 [44–46] and SRD5A2 [47–49], and these changes produce an undermasculinized phenotype in males. Our data revealed a callitrichine-specific NSs in SRD5A2 at protein position 49 (A49S), and polymorphisms at this position in human alter aromatase activity in bioassay and produce hypospadic phenotypes in males [50–53]. Further, PROVEAN identified a callitrichine-specific amino acid substitution (P64S) in SRD5A2 that was predicted to significantly alter protein function in 5α-reductase, relative to Saimiri. Finally, polymorphisms in the coding region for AR are also known to produce undermasculinization in human [54, 55]. We identified nine NSs in AR across callitrichines, with six common among all genera sampled. However, none of these substitutions were predicted to significantly alter AR structure or function by PROVEAN.
Dysfunction in the AMH signaling system can derive from nucleotide substitutions in both the ligand and the type II receptor genes [4, 56–58]. Our comparative analyses of callitrichines vs. Saimiri revealed multiple NSs in these genes, with significantly altered protein function as a consequence of three substitutions in AMH and two in AMHR2. The clinical condition of persistent Müllerian duct syndrome, in which males exhibit normal external genitalia but possess fallopian tubes and uteri, is associated with mutations in the coding regions for both AMH ligand and receptor genes [59, 60] in regions of exons in close proximity to NSs we identified in callitrichine primates. These clinical findings suggest that AMH system polymorphisms reduce the canalizing of early reproductive development toward a male phenotype, and the NSs we identified represent important candidates for female protection from AMH-mediated masculinization in callitrichines.
There are a growing number of reports of female callitrichines that present with masculinized features, but these reports suggest that the masculinizing effects are limited to external genitalia and are not associated with atypical development of primary gonads or accessory reproductive structures. Multiple callitrichine breeding colonies have reported the occurrence of females born with ambiguous masculinized genitalia, including hypertrophied clitori and congenitally fused labia with a single urethral opening, in contrast to normative urethra and vaginal structures. These atypical females have been reported in both Callithrix and Leontopithecus. In cases where the sex of the female’s cotwin was reported, it was universally a male, and karyotyping or genotyping of females presenting with atypical genitalia revealed the presence of Y chromosomes in lymphocytes or Y chromosome-associated zinc-finger protein (Y-ZFP) [25, 26]. Endocrine analyses of females indicated normative ovarian function including progesterone profiles consistent with ovulatory cycles [25, 61, 62] and low levels of testosterone  and its metabolic precursor androstenedione . Ultrasound or post-mortem analyses revealed female-typical structures in the fallopian tubes, uterine horns and uterus, cervix, and vagina, ovaries that displayed primary and secondary follicles and corpora lutea, and no signs of testicular tissue [24, 25, 61, 63]. In one case in which the fused labia were surgically separated, the female became pregnant and delivered a normal twin litter after pairing with a male . Collectively, these data suggest that while there may be masculinizing effects of a male cotwin on the structure of external genitalia in females, the development of the accessory reproductive system in females is relatively insensitive to the presence of males in utero. Given that the differentiation of accessory reproductive structures is strongly regulated by AMH, this pathway is likely to be more important in protecting female offspring from the defeminizing consequences of exposure to a male cotwin. Our data show that more NSs, and a greater proportion NSs with significant predicted changes in protein structure and function are found in AMH-related genes in callitrichine species, suggesting that this system is has undergone selection in twinning species and hence may serve an important role in female protection against defeminnization in utero.
A proposal has been forwarded regarding a shift from prenatal to postnatal staging of development in callitrichines, including the realm of sexual differentiation . Abbott and Hearn [65, 66] suggested that the majority of differentiation of behavioral sex differences, and some aspects of genital morphology, may occur postnatally as an evolutionary adaptation to minimize the impact of male-derived endocrine signals for masculinization and defeminization. This notion is supported by normative endocrine data, since testosterone levels are elevated in males but not females from 15 to 100 days postpartum [67–69]. Immediate postnatal treatment of females with exogenous testosterone partially masculinizes female external genitalia [66, 70]. Neonatal castration produces adult males with slightly demasculinized genitalia, and testosterone treatment of these castrates as adults results in increased penile and scrotal size . Behavioral sexual differentiation in callitrichines appears to be organized in part by postnatal endocrine factors, since neonatal castration in males is associated with lower rates of sexual and aggressive behavior as adults [72, 73] and neonatal androgen treatment in females enhances male-typical sexual behavior and reduces female-typical behavior . The observation that untreated females gestated with male cotwins showed no male-typical behavior and perfectly competent female sexual behavior is evidence for both the postnatal staging hypothesis and for the presence of mechanisms to minimize masculinization of females in utero. However, marmoset genitalia are typically well-differentiated at birth, and immediate postnatal androgen manipulations only slightly modify genital morphology . Thus, delayed timing in behavioral sexual differentiation, along with NSs-associated altered function in hormones, enzymes, and receptors involved in sexual differentiation of accessory reproductive structures, may work together to minimize masculinization and defeminization of females in twinning callitrichine primates.
Callimico, the only callitrichine that produces single-offspring litters, represents an interesting “control” species for our analyses. Recent molecular phylogenies have placed this genus as a sister taxa of Callithrix [74, 75], suggesting that Callimico does not represent an ancestral callitrichine from which twinning species subsequently evolved but rather a derived species that lost the ability to produce twin litters. Callimico shares a number of NSs with the other callitrichines, relative to Saimiri. However, we noted two NSs in Callimico that led to a reversion to the Saimiri and human amino acid residue in AMHR2 and CYP19A1. To the extent that NSs associated with female protection are neutral for single-litter offspring, there is no reason to expect reversion. However, in the event that NSs associated with sexual differentiation in twinning species are deleterious for single-litter species (e.g., reduction in masculinization) then reversions are to be expected. The two NSs we identified thus represent good candidates for further evaluation for critical roles in sexual differentiation and female protection from masculinization/defeminzation.
Our sequence results, while suggesting important differences in the genes involved in the steroid- and AMH-mediated components of sexual differentiation, do not definitively demonstrate that the NSs identified in callitrichines represent the key “protective” factors that minimize female masculinization and defeminization in utero, nor do they imply that they are the only components of the complex signaling cascade the yields male-female differences. Future studies that highlight the full coding and promoter regions of the genes we selected for analysis and additional candidate genes involved in sexual differentiation will provide a more complete list of candidate alleles and NSs that may protect female fetuses in twinning vs. nontwinning simian primates. Further, knowledge of the pattern of expression of our candidate genes and their protein products, both in a developmental context and with regard to tissue specificity of expression between males and females, will assist in understanding the full story of prenatal and postnatal sexual differentiation.
In sum, our data show that the reproductive and longevity consequences of sharing a common uterine environment and extensive vascular connections with a male cotwin for female callitrichines are minimal, with no effects on female longevity, litter size, number of stillborn infants, or age at the onset of first ovulation. These findings are in marked contrast to the extreme masculinization and sterility in the freemartin condition in domestic cattle and other species , with which callitrichines share the common traits of XY chimerism in females and extensive shared vascularization among opposite-sex cotwins. The notion that selection may have favored mechanisms that protect female callitrichines against masculinization and defeminization in utero is supported by the large number of NSs in a suite of genes that are critical for morphological, neuroendocrine, and behavioral sexual differentiation. Our finding that more amino acid substitutions were documented in AMH-related genes than in genes associated with steroid synthesis, transport, and cellular action is consistent with clinical cases and colony-wide surveys of masculinized external genitalia in females that gestate with male cotwins but normative accessory reproductive structures and sex-typical behavior in these females.
AMH, anti-Müllerian hormone; AMHR2, anti-Müllerian hormone receptor type 2; AR, androgen receptor; CYP19A1, aromatase; IUP, intrauterine position effect; NS, nonsynonymous nucleotide substitutions; SHBG, sex hormone-binding globulin; SRD5A2, steroid-5-alpha reductase
Quinn A, Koopman P. The molecular genetics of sex determination and sex reversal in mammals. Semin Reprod Med. 2012;30:351–63.
Lasala C, Schteingart HF, Arouche N, Bedecarrás P, Grinspon RP, Picard J-Y, et al. SOX9 and SF1 are involved in cyclic AMP-mediated upregulation of anti-Müllerian gene expression in the testicular prepubertal Sertoli cell line SMAT1. Am J Physiol-Endocrinol Metab. 2011;301:E539–47.
Rey R, Lukas-Croisier C, Lasala C, Bedecarrás P. AMH/MIS: what we know already about the gene, the protein and its regulation. Mol Cell Endocrinol. 2003;211:21–31.
Josso N, Picard JY, Rey R, Di Clemente N. Testicular anti-Mullerian hormone: history, genetics, regulation and clinical applications. Pediatr Endocrinol Rev PER. 2006;3:347–58.
Ryan BC, Vandenbergh JG. Intrauterine position effects. Neurosci Biobehav Rev. 2002;26:665–78.
vom Saal FS. Variation in phenotype due to random intrauterine positioning of male and female fetuses in rodents. J Reprod Fertil. 1981;62:633–50.
Hacklaender K, Arnold W. Litter sex ratio affects lifetime reproductive success of free-living female Alpine marmots Marmota marmota†. Mammal Rev. 2012;42:310–3.
Korsoff P, Bogl LH, Korhonen P, Kangas AJ, Soininen P, Ala-Korpela M, et al. A comparison of anthropometric, metabolic, and reproductive characteristics of young adult women from opposite-sex and same-sex twin pairs. Front Endocrinol. 2014;5:28.
Monclús R, Blumstein DT. Litter sex composition affects life-history traits in yellow-bellied marmots. J Anim Ecol. 2012;81:80–6.
Lillie FR. Sex-determination and sex-differentiation in mammals. Proc Natl Acad Sci U S A. 1917;3:464–70.
Capel B, Coveney D. Frank Lillie’s freemartin: illuminating the pathway to 21st century reproductive endocrinology. J Exp Zool Comp Exp Biol. 2004;301:853–6.
Mishina Y, Behringer RR. The in vivo function of Müllerian-inhibinting substance during mammalian sexual development. Adv Dev Biol. 1996;4:1–25.
Padula AM. The freemartin syndrome: an update. Anim Reprod Sci. 2005;87:93–109.
Harris RA, Tardif SD, Vinar T, Wildman DE, Rutherford JN, Rogers J, et al. Evolutionary genetics and implications of small size and twinning in callitrichine primates. Proc Natl Acad Sci. 2014;111:1467–72.
Ross CN, Fite JE, Jensen H, French JA. Demographic review of a captive colony of callitrichids (Callithrix kuhlii). Am J Primatol. 2007;69:234–40.
Smucny DA, Abbott DH, Mansfield KG, Schultz-Darken NJ, Yamamoto ME, Alencar AI, et al. Reproductive output, maternal age, and survivorship in captive common marmoset females (Callithrix jacchus). Am J Primatol. 2004;64:107–21.
Chambers PL, Hearn JP. Embryonic, foetal and placental development in the common marmoset monkey (Callithrix jacchus). J Zool. 1985;207:545–61.
Phillips IR. The embryology of the common marmoset (Callithrix jacchus). Adv Anat Embryol Cell Biol. 1975;52:3–47.
Rutherford JN, Tardif SD. Placental efficiency and intrauterine resource allocation strategies in the common marmoset pregnancy. Am J Phys Anthr. 2008;137:60–8.
Wislocki GB. Observations on twinning in marmosets. Am J Anat. 1939;64:445–83.
Benirschke K, Anderson JM, Brownhill LE. Marrow chimerism in marmosets. Science. 1962;138:513–5.
Gengozian N, Batson JS, Eide P. Hematologic and cytogenetic evidence for hematopoietic chimerism in the marmoset, Tamarinus nigricollis. Cytogenet Genome Res. 1964;3:384–93.
Ross CN, French JA, Orti G. Germ-line chimerism and paternal care in marmosets (Callithrix kuhlii). Proc Natl Acad Sci U S A. 2007;104:6278–82.
Goldschmidt B, Moraes IA, Souza LM, Paulino FS, Pissinatti A, Marsico FF, et al. Occurrence of virilixation signals in a female marmoset Leontopithecus chrysomelas (Callitrichidae; Primates) with 46, XX/46, XY chimerism. Isr J Vet Med. 2005;60:86.
Sanchez-Morgado JM, Haworth R, Morris TH. XY female marmoset (Callithrix jacchus). Comp Med. 2003;53:539–44.
Smith AS, Birnie AK, French JA. Prenatal androgens affect development and behavior in primates. In Building Babies: Primate Development in Proximate and Ultimate Perspective. Edited by Clancy KBH, Hinde, K, Rutherford, JN. New York: Springer; 2013, 103–132.
Sweeney CG, Curran E, Westmoreland SV, Mansfield KG, Vallender EJ. Quantitative molecular assessment of chimerism across tissues in marmosets and tamarins. BMC Genomics. 2012;13:98.
Ardito G, Lamberti L, Bigatti P, Crovella S, Oberto G. No correlation between chimerism and fertility in Callithrix jacchus (Callithricidae, Primates). Int J Anthropol. 1995;10:15–9.
Rutherford JN, Colon DGL, Ross CN, Tardif SD, et al. Developmental origins of pregnancy loss in the adult female common marmoset monkey (Callithrix jacchus). PLoS One. 2014;9:e96845.
Ahmed SF, Hughes IA. The genetics of male undermasculinization. Clin Endocrinol (Oxf). 2002;56:1–18.
Warne GL, Kanumakala S. Molecular endocrinology of sex differentiation. Semin Reprod Med. 2002;29:169–80.
Jaquish CE, Tardif SD, Toal RL, Carson RL. Patterns of prenatal survival in the common marmoset (Callithrix jacchus). J Med Primatol. 1996;25:57–63.
Tardif S, Ross C, Smucny D. Building marmoset babies: trade-offs and cutting bait. In Building Babies: Primate Development in Proximate and Ultimate Perspective. Edited by Clancy KBH, Hinde, K, Rutherford, JN. New York: Springer; 2013, 169–186.
Mustoe AC, Jensen HA, French JA. Describing ovarian cycles, pregnancy characteristics, and the use of contraception in female white-faced marmosets, Callithrix geoffroyi. Am J Primatol. 2012;74:1044–53.
Hiort O. Clinical and molecular aspects of androgen insensitivity. Endocr Dev. 2013;24:33–40.
Zhang J. Rates of conservative and radical nonsynonymous nucleotide substitutions in mammalian nuclear genes. J Mol Evol. 2000;50:56–68.
Choi Y, Chan AP. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015;31:2745–2747.
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;7–20.
Sim N-L, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012;40:W452–7.
Zielinski WJ, Vandenbergh JG. Effect of intrauterine position and social density on age of first reproduction in wild-type female house mice (Mus musculus). J Comp Psychol. 1991;105:134.
Zielinski WJ, vom Saal FS, Vandenbergh JG. The effect of intrauterine position on the survival, reproduction and home range size of female house mice (Mus musculus). Behav Ecol Sociobiol. 1992;30:185–91.
Drickamer LC, Arthur RD, Rosenthal TL. Conception failure in swine: importance of the sex ratio of a female’s birth litter and tests of other factors. J Anim Sci. 1997;75:2192–6.
Hammond GL. Diverse roles for sex hormone-binding globulin in reproduction. Biol Reprod. 2011;85:431–41.
Wu T-S, Hammond GL. Naturally occurring mutants inform SHBG structure and function. Mol Endocrinol. 2014;28:1026–38.
Hammond GL, Miguel-Queralt S, Yalcinkaya TM, Underhill C, Place NJ, Glickman SE, et al. Phylogenetic comparisons implicate sex hormone-binding globulin in “masculinization” of the female spotted hyena (Crocuta crocuta). Endocrinology. 2012;153:1435–43.
Ma CX, Adjei AA, Salavaggione OE, Coronel J, Pelleymounter L, Wang L, et al. Human aromatase: gene resequencing and functional genomics. Cancer Res. 2005;65:11071–82.
Rochira V, Carani C. Aromatase deficiency in men: a clinical perspective. Nat Rev Endocrinol. 2009;5:559–68.
Stratakis CA. An aroma of complexity: how the unique genetics of aromatase (CYP19A1) explain diverse phenotypes from hens and hyenas to human gynecomastia, and testicular and other tumors. J Clin Endocrinol Metab. 2013;98:4676–81.
Berra M, Williams EL, Muroni B, Creighton SM, Honour JW, Rumsby G, et al. Recognition of 5α-reductase-2 deficiency in an adult female 46XY DSD clinic. Eur J Endocrinol. 2011;164:1019–25.
Di Marco C, Bulotta AL, Varetti C, Dosa L, Michelucci A, Baldinotti F, et al. Ambiguous external genitalia due to defect of 5-α-reductase in seven Iraqi patients: prevalence of a novel mutation. Gene. 2013;526:490–3.
Tsai M-C, Chou Y-Y, Lin S-J, Tsai L-P. A novel SRD5A2 mutation in a Taiwanese newborn with ambiguous genitalia. Kaohsiung J Med Sci. 2012;28:231–5.
Makridakis NM, di Salle E, Reichardt JK. Biochemical and pharmacogenetic dissection of human steroid 5α-reductase type II. Pharmacogenet Genomics. 2000;10:407–13.
Silver RI, Russell DW. 5 α-Reductase type 2 mutations are present in some boys with isolated hypospadias. J Urol. 1999;162:1142–5.
Lagarde WH, Blackwelder AJ, Minges JT, Hnat AT, French FS, Wilson EM. Androgen receptor exon 1 mutation causes androgen insensitivity by creating phosphorylation site and inhibiting melanoma antigen-A11 activation of NH2-and carboxyl-terminal interaction-dependent transactivation. J Biol Chem. 2012;287:10905–15.
Matias PM, Donner P, Coelho R, Thomaz M, Peixoto C, Macedo S, et al. Structural evidence for ligand specificity in the binding domain of the human androgen receptor Implications for pathogenic gene mutations. J Biol Chem. 2000;275:26164–71.
Belville C, Van Vlijmen H, Ehrenfels C, Pepinsky B, Rezaie AR, Picard J-Y, et al. Mutations of the anti-Mullerian hormone gene in patients with persistent Mullerian duct syndrome: biosynthesis, secretion, and processing of the abnormal proteins and analysis using a three-dimensional model. Mol Endocrinol. 2004;18:708–21.
Josso N, Belville C, di Clemente N, Picard J-Y. AMH and AMH receptor defects in persistent Müllerian duct syndrome. Hum Reprod Update. 2005;11:351–6.
Nishi MY, Domenice S, Maciel-Guerra AT, Zaba Neto A, da Silva MACP, Costa EMF, et al. Analysis of anti-Müllerian hormone (AMH) and its receptor (AMHR2) genes in patients with persistent Müllerian duct syndrome. Arq Bras Endocrinol Metabol. 2012;56:473–8.
Abduljabbar M, Taheini K, Picard J-Y, Cate RL, Josso N. Mutations of the AMH type II receptor in two extended families with persistent Müllerian duct syndrome: lack of phenotype/genotype correlation. Horm Res Paediatr. 2012;77:291–7.
Wongprasert H, Somanunt S, De Filippo R, Picard JY, Pitukcheewanont P. A novel mutation of anti-Mullerian hormone gene in persistent Mullerian duct syndrome presented with bilateral cryptorchidism: a case report. J Pediatr Urol. 2013;9:e147–9.
Isachenko EF, Nayudu PL, Isachenko VV, Nawroth F, Michelmann HW. Congenitally caused fused labia in the common marmoset (Callithrix jacchus). J Med Primatol. 2002;31:350–5.
Wedi E, Nayudu PL, Michelmann HW. A case report of spontaneous opening of congenitally fused labia in a female common marmoset (Callithrix jacchus) followed by pregnancy and birth of twins. J Med Primatol. 2011;40:351–3.
Niimi K, Oguchi A, Nishio K, Okano Y, Takahashi E. Congenital malformation of the vaginal orifice, imperforate vagina, in the common marmoset (Callithrix jacchus). J Vet Med Sci. 2015;77:345.
Missler M, Wolff JR, Rothe H, Heger W, Merker HJ, Treiber A, et al. Developmental biology of the common marmoset: proposal for a“postnatal staging”. J Med Primatol. 1992;21:285–98.
Abbott DH. Differentiation of sexual behaviour in female marmoset monkeys: effects of neonatal testosterone or a male co-twin. Prog Brain Res. 1984;61:349–58.
Abbott DH, Hearn JP. The effects of neonatal exposure to testosterone on the development of behaviour in female marmoset monkeys. Sex, hormones, and behaviour. 1979;62:299–313.
Abbott DH, Hearn JP. Physical, hormonal and behavioural aspects of sexual development in the marmoset monkey, Callithrix jacchus. J Reprod Fertil. 1978;53:155–66.
Birnie AK, Smith AS, Nali C, French JA. Social and developmental influences on urinary androgen levels in young male white-faced marmosets (Callithrix geoffroyi). Am J Primatol. 2011;73:378–85.
Dixson AF. Plasma testosterone concentrations during postnatal development in the male common marmoset. Folia Primatol Basel. 1986;47:166–70.
Abbott DH, Hearn JP. The effects of neonatal exposure to testosterone on the development of behaviour in female marmoset monkeys. Sex Horm Behav. 1979;62:299–313.
Dixson AF. Sexual and aggressive behaviour of adult male marmosets (Callithrix jacchus) castrated neonatally, prepubertally, or in adulthood. Physiol Behav. 1993;54:301–7.
Dixson AF. Effects of testosterone propionate upon the sexual and aggressive behavior of adult male marmosets (Callithrix jacchus) castrated as neonates. Horm Behav. 1993;27:216–30.
Epple G, Alveario MC, Belcher AM. Copulatory behavior of adult tamarins (Saguinus fuscicollis) castrated as neonates or juveniles: effect of testosterone treatment. Horm Behav. 1990;24:470–83.
Perelman P, Johnson WE, Roos C, Seuánez HN, Horvath JE, Moreira MA, et al. A molecular phylogeny of living primates. PLoS Genet. 2011;7, e1001342.
Wildman DE, Jameson NM, Opazo JC, Soojin VY. A fully resolved genus level phylogeny of neotropical primates (Platyrrhini). Mol Phylogenet Evol. 2009;53:694–702.
We thank J.H. Taylor for comments and critiques on previous versions of this paper and the institutions that kindly provided the tissues or DNA samples for our study. This study is supported in part by a grant to JAF (NIH HD42884).
JAF designed research. DR, ACM, and JC sequenced and analyzed the genetic data. BF, JC, LR, and JM analyzed the breeding records. JAF, BF, JC, DR, ACM, and LR wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Regions of candidate genes sequenced de novo for this study and PCR primers used to amplify selected regions. (DOCX 19 kb)
Tissue or DNA source for gene sequencing. (DOCX 67 kb)
Summary of amino acid residue substitutions and predicted change in structure/function (PROVEAN, SIFT, and PolyPhen) for the six candidate genes for human, Saimiri, and four callitrichine primates. (XLSX 19 kb)