Dmrt1 polymorphism covaries with sex‐determination patterns in Rana temporaria

Patterns of sex‐chromosome differentiation and gonadal development have been shown to vary among populations of Rana temporaria along a latitudinal transect in Sweden. Frogs from the northern‐boreal population of Ammarnäs displayed well‐differentiated X and Y haplotypes, early gonadal differentiation, and a perfect match between phenotypic and genotypic sex. In contrast, no differentiated Y haplotypes could be detected in the southern population of Tvedöra, where juveniles furthermore showed delayed gonadal differentiation. Here, we show that Dmrt1, a gene that plays a key role in sex determination and sexual development across all metazoans, displays significant sex differentiation in Tvedöra, with a Y‐specific haplotype distinct from Ammarnäs. The differential segment is not only much shorter in Tvedöra than in Ammarnäs, it is also less differentiated and associates with both delayed gonadal differentiation and imperfect match between phenotypic and genotypic sex. Whereas Tvedöra juveniles with a local Y haplotype tend to ultimately develop as males, those without it may nevertheless become functional XX males, but with strongly female‐biased progeny. Our findings suggest that the variance in patterns of sex determination documented in common frogs might result from a genetic polymorphism within a small genomic region that contains Dmrt1. They also substantiate the view that recurrent convergences of sex determination toward a limited set of chromosome pairs may result from the co‐option of small genomic regions that harbor key genes from the sex‐determination pathway.     

Geographic variation in sex‐chromosome differentiation in the common frog (Rana temporaria)

In sharp contrast with birds and mammals, sex‐determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex‐determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500‐km transect across Sweden. Genetic differentiation between sexes is strongest (F_ST = 0.152) in a northern‐boreal population, where male‐specific alleles and heterozygote excesses (F_IS = ?0.418 in males, +0.025 in females) testify to a male‐heterogametic system and lack of X‐Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (F_ST = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male‐specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X‐Y recombination, co‐option of an alternative sex‐chromosome pair, or a mixed sex‐determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the ‘sexual races’ described in common frogs in the 1930s.     

Within‐population polymorphism of sex‐determination systems in the common frog (Rana temporaria)

In sharp contrast with birds and mammals, the sex chromosomes of ectothermic vertebrates are often undifferentiated, for reasons that remain debated. A linkage map was recently published for Rana temporaria (Linnaeus, 1758) from Fennoscandia (Eastern European lineage), with a proposed sex‐determining role for linkage group 2 (LG2). We analysed linkage patterns in lowland and highland populations from Switzerland (Western European lineage), with special focus on LG2. Sibship analyses showed large differences from the Fennoscandian map in terms of recombination rates and loci order, pointing to large‐scale inversions or translocations. All linkage groups displayed extreme heterochiasmy (total map length was 12.2 cM in males, versus 869.8 cM in females). Sex determination was polymorphic within populations: a majority of families (with equal sex ratios) showed a strong correlation between offspring phenotypic sex and LG2 paternal haplotypes, whereas other families (some of which with female‐biased sex ratios) did not show any correlation. The factors determining sex in the latter could not be identified. This coexistence of several sex‐determination systems should induce frequent recombination of X and Y haplotypes, even in the absence of male recombination. Accordingly, we found no sex differences in allelic frequencies on LG2 markers among wild‐caught male and female adults, except in one high‐altitude population, where nonrecombinant Y haplotypes suggest sex to be entirely determined by LG2. Multifactorial sex determination certainly contributes to the lack of sex‐chromosome differentiation in amphibians.     

Novel evolutionary pathways of sex‐determining mechanisms

Evolutionary transitions between sex‐determining mechanisms (SDMs) are an enigma. Among vertebrates, individual sex (male or female) is primarily determined by either genes (genotypic sex determination, GSD) or embryonic incubation temperature (temperature‐dependent sex determination, TSD), and these mechanisms have undergone repeated evolutionary transitions. Despite this evolutionary lability, transitions from GSD (i.e. from male heterogamety, XX/XY, or female heterogamety, ZZ/ZW) to TSD are an evolutionary conundrum, as they appear to require crossing a fitness valley arising from the production of genotypes with reduced viability owing to being homogametic for degenerated sex chromosomes (YY or WW individuals). Moreover, it is unclear whether alternative (e.g. mixed) forms of sex determination can persist across evolutionary time. It has previously been suggested that transitions would be easy if temperature‐dependent sex reversal (e.g. XX male or XY female) was asymmetrical, occurring only in the homogametic sex. However, only recently has a mechanistic model of sex determination emerged that may allow such asymmetrical sex reversal. We demonstrate that selection for TSD in a realistic sex‐determining system can readily drive evolutionary transitions from GSD to TSD that do not require the production of YY or WW individuals. In XX/XY systems, sex reversal (female to male) occurs in a portion of the XX individuals only, leading to the loss of the Y allele (or chromosome) from the population as XX individuals mate with each other. The outcome is a population of XX individuals whose sex is determined by incubation temperature (TSD). Moreover, our model reveals a novel evolutionarily stable state representing a mixed‐mechanism system that has not been revealed by previous approaches. This study solves two long‐standing puzzles of the evolution of sex‐determining mechanisms by illuminating the evolutionary pathways and endpoints.     

Phylogeography,more than elevation,accounts for sex chromosome differentiation in Swiss populations of the common frog (Rana temporaria)*

Sex chromosomes in vertebrates range from highly heteromorphic (as in most birds and mammals) to strictly homomorphic (as in many fishes, amphibians, and nonavian reptiles). Reasons for these contrasted evolutionary trajectories remain unclear, but species such as common frogs with polymorphism in the extent of sex chromosome differentiation may potentially deliver important clues. By investigating 92 common frog populations from a wide range of elevations throughout Switzerland, we show that sex chromosome differentiation strongly correlates with alleles at the candidate sex-determining gene Dmrt1. Y-specific Dmrt1 haplotypes cluster into two main haplogroups, Y_A and Y_B, with a phylogeographic signal that parallels mtDNA haplotypes: Y_A populations, with mostly well-differentiated sex chromosomes, occur primarily south of the main alpine ridge that bisects Switzerland, whereas Y_B populations, with mostly undifferentiated (proto-)sex chromosomes, occur north of this ridge. Elevation has only a marginal effect, opposing previous suggestions of a major role for climate on sex chromosome differentiation. The Y-haplotype effect might result from differences in the penetrance of alleles at the sex-determining locus (such that sex reversal and ensuing X-Y recombination are more frequent in Y_B populations), and/or fixation of an inversion on Y_A (as supported by the empirical observation that Y_A haplotypes might not recombine in XY_A females).     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

TEMPERATURE‐DEPENDENT TURNOVERS IN SEX‐DETERMINATION MECHANISMS: A QUANTITATIVE MODEL

Sex determination is often seen as a dichotomous process: individual sex is assumed to be determined either by genetic (genotypic sex determination, GSD) or by environmental factors (environmental sex determination, ESD), most often temperature (temperature sex determination, TSD). We endorse an alternative view, which sees GSD and TSD as the ends of a continuum. Both effects interact a priori, because temperature can affect gene expression at any step along the sex‐determination cascade. We propose to define sex‐determination systems at the population‐ (rather than individual) level, via the proportion of variance in phenotypic sex stemming from genetic versus environmental factors, and we formalize this concept in a quantitative‐genetics framework. Sex is seen as a threshold trait underlain by a liability factor, and reaction norms allow modeling interactions between genotypic and temperature effects (seen as the necessary consequences of thermodynamic constraints on the underlying physiological processes). As this formalization shows, temperature changes (due to e.g., climatic changes or range expansions) are expected to provoke turnovers in sex‐ determination mechanisms, by inducing large‐scale sex reversal and thereby sex‐ratio selection for alternative sex‐determining genes. The frequency of turnovers and prevalence of homomorphic sex chromosomes in cold‐blooded vertebrates might thus directly relate to the temperature dependence in sex‐determination mechanisms.     

High‐density linkage maps fail to detect any genetic component to sex determination in a Rana temporaria family

Sex chromosome differentiation in Rana temporaria varies strikingly among populations or families: whereas some males display well‐differentiated Y haplotypes at microsatellite markers on linkage group 2 (LG₂), others are genetically undistinguishable from females. We analysed with RADseq markers one family from a Swiss lowland population with no differentiated sex chromosomes, and where sibship analyses had failed to detect any association between the phenotypic sex of progeny and parental haplotypes. Offspring were reared in a common tank in outdoor conditions and sexed at the froglet stage. We could map a total of 2177 SNPs (1123 in the mother, 1054 in the father), recovering in both adults 13 linkage groups (= chromosome pairs) that were strongly syntenic to Xenopus tropicalis despite > 200 My divergence. Sexes differed strikingly in the localization of crossovers, which were uniformly distributed in the female but limited to chromosome ends in the male. None of the 2177 markers showed significant association with offspring sex. Considering the very high power of our analysis, we conclude that sex determination was not genetic in this family; which factors determined sex remain to be investigated.     

Recombination is suppressed and variability reduced in a nascent Y chromosome

Several hypotheses have been elaborated to account for the evolutionary decay commonly observed in full-fledged Y chromosomes. Enhanced drift, background selection and selective sweeps, which are expected to result from reduced recombination, may all share responsibilities in the initial decay of proto-Y chromosomes, but little empirical information has been gathered so far. Here we take advantage of three markers that amplify on both of the morphologically undifferentiated sex chromosomes of the European tree frog (Hyla arborea) to show that recombination is suppressed in males (the heterogametic sex) but not in females. Accordingly, genetic variability is reduced on the Y, but in a way that can be accounted for by merely the number of chromosome copies per breeding pair, without the need to invoke background selection or selective sweeps.     

Inference of sex‐specific expansion patterns in human populations from Y‐chromosome polymorphism

Studying the current distribution of genetic diversity in humans has important implications for our understanding of the history of our species. We analyzed a set of linked STR and SNP loci from the paternally inherited Y chromosome to infer the past demography of 55 African and Eurasian populations, using both the parametric and nonparametric coalescent‐based methods implemented in the BEAST application. We inferred expansion events in most sedentary farmer populations, while we found constant effective population sizes for both nomadic hunter‐gatherers and seminomadic herders. Our results differed, on several aspects, from previous results on mtDNA and autosomal markers. First, we found more recent expansion patterns in Eurasia than in Africa. This discrepancy, substantially stronger than the ones found with the other kind of markers, may result from a lower effective population size for men, which might have made male‐transmitted markers more sensitive to the out‐of‐Africa bottleneck. Second, we found expansion signals only for sedentary farmers but not for nomadic herders in Central Asia, while these signals were found for both kind of populations in this area when using mtDNA or autosomal markers. Expansion signals in this area may result from spatial expansion processes and may have been erased for the Y chromosome among the herders because of restricted male gene flow. Am J Phys Anthropol 157:217–225, 2015. © 2015 Wiley Periodicals, Inc.     

No evidence that Y‐chromosome differentiation affects male fitness in a Swiss population of common frogs

The canonical model of sex‐chromosome evolution assigns a key role to sexually antagonistic (SA) genes on the arrest of recombination and ensuing degeneration of Y chromosomes. This assumption cannot be tested in organisms with highly differentiated sex chromosomes, such as mammals or birds, owing to the lack of polymorphism. Fixation of SA alleles, furthermore, might be the consequence rather than the cause of recombination arrest. Here we focus on a population of common frogs (Rana temporaria) where XY males with genetically differentiated Y chromosomes (nonrecombinant Y haplotypes) coexist with both XY° males with proto‐Y chromosomes (only differentiated from X chromosomes in the immediate vicinity of the candidate sex‐determining locus Dmrt1) and XX males with undifferentiated sex chromosomes (genetically identical to XX females). Our study finds no effect of sex‐chromosome differentiation on male phenotype, mating success or fathering success. Our conclusions rejoin genomic studies that found no differences in gene expression between XY, XY° and XX males. Sexual dimorphism in common frogs might result more from the differential expression of autosomal genes than from sex‐linked SA genes. Among‐male variance in sex‐chromosome differentiation seems better explained by a polymorphism in the penetrance of alleles at the sex locus, resulting in variable levels of sex reversal (and thus of X‐Y recombination in XY females), independent of sex‐linked SA genes.     

Female‐only sex‐linked amplified fragment length polymorphism markers support ZW/ZZ sex determination in the giant freshwater prawn Macrobrachium rosenbergii

Sex determination mechanisms in many crustacean species are complex and poorly documented. In the giant freshwater prawn, Macrobrachium rosenbergii, a ZW/ZZ sex determination system was previously proposed based on sex ratio data obtained by crosses of sex‐reversed females (neomales). To provide molecular evidence for the proposed system, novel sex‐linked molecular markers were isolated in this species. Amplified fragment length polymorphism (AFLP) using 64 primer combinations was employed to screen prawn genomes for DNA markers linked with sex loci. Approximately 8400 legible fragments were produced, 13 of which were uniquely identified in female prawns with no indication of corresponding male‐specific markers. These AFLP fragments were reamplified, cloned and sequenced, producing two reliable female‐specific sequence characterized amplified region (SCAR) markers. Additional individuals from two unrelated geographic populations were used to verify these findings, confirming female‐specific amplification of single bands. Detection of internal polymorphic sites was conducted by designing new primer pairs based on these internal fragments. The internal SCAR fragments also displayed specificity in females, indicating high levels of variation between female and male specimens. The distinctive feature of female‐linked SCAR markers can be applied for rapid detection of prawn gender. These sex‐specific SCAR markers and sex‐associated AFLP candidates unique to female specimens support a sex determination system consistent with female heterogamety (ZW) and male homogamety (ZZ).     

Genomic characterization of sex‐identification markers in Sebastes carnatus and Sebastes chrysomelas rockfishes

Fish have evolved a variety of sex‐determining (SD) systems including male heterogamy (XY), female heterogamy (ZW) and environmental SD. Little is known about SD mechanisms of Sebastes rockfishes, a highly speciose genus of importance to evolutionary and conservation biology. Here, we characterize the sex determination system in the sympatrically distributed sister species Sebastes chrysomelas and Sebastes carnatus. To identify sex‐specific genotypic markers, double digest restriction site – associated DNA sequencing (ddRAD‐seq) of genomic DNA from 40 sexed individuals of both species was performed. Loci were filtered for presence in all of the individuals of one sex, absence in the other sex and no heterozygosity. Of the 74 965 loci present in all males, 33 male‐specific loci met the criteria in at least one species and 17 in both. Conversely, no female‐specific loci were detected, together providing evidence of an XY sex determination system in both species. When aligned to a draft reference genome from Sebastes aleutianus, 26 sex‐specific loci were interspersed among 1168 loci that were identical between sexes. The nascent Y chromosome averaged 5% divergence from the X chromosome and mapped to reference Sebastes genome scaffolds totalling 6.9Mbp in length. These scaffolds aligned to a single chromosome in three model fish genomes. Read coverage differences were also detected between sex‐specific and autosomal loci. A PCR‐RFLP assay validated the bioinformatic results and correctly identified sex of five additional individuals of known sex. A sex‐determining gene in other teleosts gonadal soma‐derived factor (gsdf) was present in the model fish chromosomes that spanned our sex‐specific markers.     

Conflict over condition‐dependent sex allocation can lead to mixed sex‐determination systems

Theory suggests that genetic conflicts drive turnovers between sex‐determining mechanisms, yet these studies only apply to cases where sex allocation is independent of environment or condition. Here, we model parent–offspring conflict in the presence of condition‐dependent sex allocation, where the environment has sex‐specific fitness consequences. Additionally, one sex is assumed to be more costly to produce than the other, which leads offspring to favor a sex ratio less biased toward the cheaper sex in comparison to the sex ratio favored by mothers. The scope for parent–offspring conflict depends on the relative frequency of both environments: when one environment is less common than the other, parent–offspring conflict can be reduced or even entirely absent, despite a biased population sex ratio. The model shows that conflict‐driven invasions of condition‐independent sex factors (e.g., sex chromosomes) result either in the loss of condition‐dependent sex allocation, or, interestingly, lead to stable mixtures of condition‐dependent and condition‐independent sex factors. The latter outcome corresponds to empirical observations in which sex chromosomes are present in organisms with environment‐dependent sex determination. Finally, conflict can also favor errors in environmental perception, potentially resulting in the loss of condition‐dependent sex allocation without genetic changes to sex‐determining loci.     

The evolution of sex‐determining mechanisms: lessons from temperature‐sensitive mutations in sex determination genes in Caenorhabditis elegans

Sexual reproduction is one of the most taxonomically conserved traits, yet sex‐determining mechanisms (SDMs) are quite diverse. For instance, there are numerous forms of environmental sex determination (ESD), in which an organism’s sex is determined not by genotype, but by environmental factors during development. Important questions remain regarding transitions between SDMs, in part because the organisms exhibiting unique mechanisms often make difficult study organisms. One potential solution is to utilize mutant strains in model organisms better suited to answering these questions. We have characterized two such strains of the model nematode Caenorhabditis elegans. These strains harbour temperature‐sensitive mutations in key sex‐determining genes. We show that they display a sex ratio reaction norm in response to rearing temperature similar to other organisms with ESD. Next, we show that these mutations also cause deleterious pleiotropic effects on overall fitness. Finally, we show that these mutations are fundamentally different at the genetic sequence level. These strains will be a useful complement to naturally occurring taxa with ESD in future research examining the molecular basis of and the selective forces driving evolutionary transitions between sex determination mechanisms.     

Biased sex‐ratio and sex‐biased heterozygote disadvantage affect the maintenance of a genetic polymorphism and the properties of hybrid zones

The evolution of biodiversity is a major issue of modern biology, and it is becoming increasingly topical as the ongoing erosion of diversity puts serious threats on human well‐being. An elementary mechanism that allows maintaining diversity is the interplay between dispersal and heterozygote selective disadvantage, which can lead to self‐sustainable spatial genetic structures and is central to the stability of hybrid zones. Theoretical studies supporting the importance of this mechanism assume a balanced sex‐ratio and a heterozygote disadvantage equally affecting both sexes, despite the multiplicity of empirical evidence that (i) adult sex‐ratio is usually biased towards either male or female and that (ii) heterozygote disadvantage often affects a single sex. We expanded the existing theory by weighting the strength of selection against heterozygote according to the biased in sex‐ratio and in heterozygote disadvantage. The range of conditions allowing for the maintenance of polymorphism can then either double or vanish. We discuss the implications of such finding for birds, mammals and insects diversity. Finally, we provide simple analytical predictions about the effect of those biased on the width and speed of hybrid zones and on the time for the spread of beneficial mutations through such zones.     

The sex chromosome system can influence the evolution of sex‐biased dispersal

Sex‐biased dispersal is a much‐discussed feature in literature on dispersal. Diverse hypotheses have been proposed to explain the evolution of sex‐biased dispersal, a difference in dispersal rate or dispersal distance between males and females. An early hypothesis has indicated that it may rely on the difference in sex chromosomes between males and females. However, this proposal was quickly rejected without a real assessment. We propose a new perspective on this hypothesis by investigating the evolution of sex‐biased dispersal when dispersal genes are sex‐linked, that is when they are located on the sex chromosomes. We show that individuals of the heterogametic sex disperse relatively more than do individuals of the homogametic sex when dispersal genes are sex‐linked rather than autosomal. Although such a sex‐biased dispersal towards the heterogametic sex is always observed in monogamous species, the mating system and the location of dispersal genes interact to modulate sex‐biased dispersal in monandry and polyandry. In the context of the multicausality of dispersal, we suggest that sex‐linked dispersal genes can influence the evolution of sex‐biased dispersal.     

Differential introgression of a female competitive trait in a hybrid zone between sex‐role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male–male competition in maintaining or eroding species barriers. Although female–female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex‐role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female–female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex‐role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.     

Climate change,sex reversal and lability of sex‐determining systems

Sex reversal at high temperatures during embryonic development (e.g., ZZ females) provides the opportunity for new genotypic crosses (e.g., ZZ male × ZZ female). This raises the alarming possibility that climatic warming could lead to the loss of an entire chromosome—one member of the sex chromosome pair (the Y or W)—and the transition of populations to environmental sex determination (ESD). Here we examine the evolutionary dynamics of sex‐determining systems exposed to climatic warming using theoretical models. We found that the loss of sex chromosomes is not an inevitable consequence of sex reversal. A large frequency of ZZ sex reversal (50% reversal from male to female) typically divides the outcome between loss of the ZW genotype and the stable persistence of ZZ males, ZW females and ZZ females. The amount of warming associated with sex chromosome loss depended on several features of wild populations—environmental fluctuation, immigration, heritable variation in temperature sensitivity and differential fecundity of sex‐reversed individuals. Chromosome loss was partially or completely buffered when sex‐reversed individuals suffered a reproductive fitness cost, when immigration occurred or when heritable variation for temperature sensitivity existed. Thus, under certain circumstances, sex chromosomes may persist cryptically in systems where the environment is the predominant influence on sex.