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.     

Dmrt1 polymorphism and sex‐chromosome differentiation in Rana temporaria

Sex‐determination mechanisms vary both within and among populations of common frogs, opening opportunities to investigate the molecular pathways and ultimate causes shaping their evolution. We investigated the association between sex‐chromosome differentiation (as assayed from microsatellites) and polymorphism at the candidate sex‐determining gene Dmrt1 in two Alpine populations. Both populations harboured a diversity of X‐linked and Y‐linked Dmrt1 haplotypes. Some males had fixed male‐specific alleles at all markers (“differentiated” Y chromosomes), others only at Dmrt1 (“proto‐” Y chromosomes), while still others were genetically indistinguishable from females (undifferentiated X chromosomes). Besides these XX males, we also found rare XY females. The several Dmrt1 Y haplotypes differed in the probability of association with a differentiated Y chromosome, which we interpret as a result of differences in the masculinizing effects of alleles at the sex‐determining locus. From our results, the polymorphism in sex‐chromosome differentiation and its association with Dmrt1, previously inferred from Swedish populations, are not just idiosyncratic features of peripheral populations, but also characterize highly diverged populations in the central range. This implies that an apparently unstable pattern has been maintained over long evolutionary times.     

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.     

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.     

A reciprocal translocation radically reshapes sex‐linked inheritance in the common frog

X and Y chromosomes can diverge when rearrangements block recombination between them. Here we present the first genomic view of a reciprocal translocation that causes two physically unconnected pairs of chromosomes to be coinherited as sex chromosomes. In a population of the common frog (Rana temporaria), both pairs of X and Y chromosomes show extensive sequence differentiation, but not degeneration of the Y chromosomes. A new method based on gene trees shows both chromosomes are sex‐linked. Furthermore, the gene trees from the two Y chromosomes have identical topologies, showing they have been coinherited since the reciprocal translocation occurred. Reciprocal translocations can thus reshape sex linkage on a much greater scale compared with inversions, the type of rearrangement that is much better known in sex chromosome evolution, and they can greatly amplify the power of sexually antagonistic selection to drive genomic rearrangement. Two more populations show evidence of other rearrangements, suggesting that this species has unprecedented structural polymorphism in its sex chromosomes.     

Sex reversal and primary sex ratios in the common frog (Rana temporaria)

Sex reversal has been suggested to have profound implications for the evolution of sex chromosomes and population dynamics in ectotherms. Occasional sex reversal of genetic males has been hypothesized to prevent the evolutionary decay of nonrecombining Y chromosomes caused by the accumulation of deleterious mutations. At the same time, sex reversals can have a negative effect on population growth rate. Here, we studied phenotypic and genotypic sex in the common frog (Rana temporaria) in a subarctic environment, where strongly female‐biased sex ratios have raised the possibility of frequent sex reversals. We developed two novel sex‐linked microsatellite markers for the species and used them with a third, existing marker and a Bayesian modelling approach to study the occurrence of sex reversal and to determine primary sex ratios in egg clutches. Our results show that a significant proportion (0.09, 95% credible interval: 0.04–0.18) of adults that were genetically female expressed the male phenotype, but there was no evidence of sex reversal of genetic males that is required for counteracting the degeneration of Y chromosome. The primary sex ratios were mostly equal, but three clutches consisted only of genetic females and three others had a significant female bias. Reproduction of the sex‐reversed genetic females appears to create all‐female clutches potentially skewing the population level adult sex‐ratio consistent with field observations. However, based on a simulation model, such a bias is expected to be small and transient and thus does not fully explain the observed female‐bias in the field.     

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).     

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.     

History vs. current demography: explaining the genetic population structure of the common frog (Rana temporaria)

The amount of genetic variability at neutral marker loci is expected to decrease, and the degree of genetic differentiation among populations to increase, as a negative function of effective population size. We assessed the patterns of genetic variability and differentiation at seven microsatellite loci in the common frog (Rana temporaria) in a hierarchical sampling scheme involving three regions (208-885 km apart), three subregions within regions and nine populations (5-20 km apart) within subregions, and related the variability and differentiation estimates to variation in local population size estimates. Genetic variability within local populations decreased significantly with increasing latitude, as well as with decreasing population size and regional site occupancy (proportion of censured localities occupied). The positive relationship between population size and genetic variability estimates was evident also when the effect of latitude (cf. colonization history) was accounted for. Significant genetic differentiation was found at all hierarchical levels, and the degree of population differentiation tended to increase with increasing latitude. Isolation by distance was evident especially at the regional sampling level, and its strength increased significantly towards the north in concordance with decreasing census and marker-based neighbourhood size estimates. These results are in line with the conjecture that the influence of current demographic factors can override the influence of historical factors on species population genetic structure. Further, the observed reductions in genetic variability and increased degree of population differentiation towards the north are in line with theoretical and empirical treatments suggesting that effective population sizes decline towards the periphery of a species' range.     

Latitudinal countergradient variation in the common frog (Rana temporaria) development rates--evidence for local adaptation

Adaptive genetic differentiation along a climatic gradient as a response to natural selection is not necessarily expressed at phenotypic level if environmental effects on population mean phenotypes oppose the genotypic effects. This form of cryptic evolution--called countergradient variation--has seldom been explicitly demonstrated for terrestrial vertebrates. We investigated the patterns of phenotypic and genotypic differentiation in developmental rates of common frogs (Rana temporaria) along a ca. 1600 km latitudinal gradient across Scandinavia. Developmental rates in the field were not latitudinally ordered, but displayed large variation even among different ponds within a given latitudinal area. In contrast, development rates assessed in the laboratory increased strongly and linearly with increasing latitude, suggesting a genetic capacity for faster development in the northern than the southern larvae. Experiments further revealed that environmental effects (temperature and food) could easily override the genetic effects on developmental rates, providing a possible mechanistic explanation as to why the genetic differentiation was not seen in the samples collected from the wild. Our results suggest that the higher developmental rates of the northern larvae are likely to be related to selection stemming from seasonal time constrains, rather than from selection dictated by low ambient temperatures per se. All in all, the results provide a demonstration of environmental effects concealing substantial latitudinally ordered genetic differentiation understandable in terms of adaptation to clinal variation in time constrains.     

A sex‐linked SCAR marker in Bryonia dioica (Cucurbitaceae), a dioecious species with XY sex‐determination and homomorphic sex chromosomes

Genetic crosses between the dioecious Bryonia dioica (Cucurbitaceae) and the monoecious B. alba in 1903 provided the first clear evidence for Mendelian inheritance of dioecy and made B. dioica the first organism for which XY sex‐determination was experimentally proven. Applying molecular tools to this system, we developed a sex‐linked sequence‐characterized amplified region (SCAR) marker for B. dioica and sequenced it for individuals representing the full geographic range of the species from Scotland to North Africa. For comparison, we also sequenced this marker for representatives of the dioecious B. cretica, B. multiflora and B. syriaca, and monoecious B. alba. In no case did any individual, male or female, yield more than two haplotypes. In northern Europe, we found strong linkage between our marker and sex, with all Y‐sequences being identical to each other. In southern Europe, however, the linkage between our marker and sex was weak, with recombination detected within both the X‐ and the Y‐homologues. Population genetic analyses suggest that the SCAR marker experienced different evolutionary pressures in northern and southern Europe. These findings fit with phylogenetic evidence that the XY system in Bryonia is labile and suggest that the genus may be a good system in which to study the early steps of sex chromosome evolution.     

Matriline effects on metamorphic traits in a natural system in the European common frog (Rana temporaria)

Successful reproduction is an important determinant of the fitness of an individual and of the dynamics of populations. Offspring of the European common frog (Rana temporaria) exhibit a high degree of variability in metamorphic traits. However, environmental factors alone cannot explain this phenotypic variability, and the influence of genetic factors remains to be determined. Here, we tested whether the maternal genotype influences developmental time, body size, and body condition of offspring in a forest pond in Germany. We collected fertilized eggs from all 57 clutches deposited in the pond. We used multilocus genotypes based on seven microsatellite loci to assign metamorphosed offspring to mothers and to determine the number of fathers for a single matriline. We tested the influence of genetic effects in the same environment by comparing variability of metamorphic traits within and between full‐sib offspring grouped to matrilines and tested whether multiple paternity increases the variability of metamorphic traits in a single matriline. The variability in size and body condition was higher within matrilines than between them, which indicates that these traits are more strongly influenced by environmental effects, which are counteracting underlying genetic effects. The developmental time varied considerably between matrilines and variability increased with the effective number of fathers, suggesting an additive genetic effect of multiple paternity. Our results show that metamorphic traits are shaped by environmental as well as genetic effects.     

First-generation linkage map for the common frog Rana temporaria reveals sex-linkage group

The common frog (Rana temporaria) has become a model species in the fields of ecology and evolutionary biology. However, lack of genomic resources has been limiting utility of this species for detailed evolutionary genetic studies. Using a set of 107 informative microsatellite markers genotyped in a large full-sib family (800 F1 offspring), we created the first linkage map for this species. This partial map-distributed over 15 linkage groups-has a total length of 1698.8 cM. In line with the fact that males are the heterogametic sex in this species and a reduction of recombination is expected, we observed a lower recombination rate in the males (map length: 1371.5 cM) as compared with females (2089.8 cM). Furthermore, three loci previously documented to be sex-linked (that is, carrying male-specific alleles) in adults from the wild mapped to the same linkage group. The linkage map described in this study is one of the densest ones available for amphibians. The discovery of a sex linkage group in Rana temporaria, as well as other regions with strongly reduced male recombination rates, should help to uncover the genetic underpinnings of the sex-determination system in this species. As the number of linkage groups found (n=15) is quite close to the actual number of chromosomes (n=13), the map should provide a useful resource for further evolutionary, ecological and conservation genetic work in this and other closely related species.     

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).     

Sex‐chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea)

Occasional XY recombination is a proposed explanation for the sex‐chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW‐European Hyla arborea populations identified male‐specific alleles at sex‐linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex‐linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex‐chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY‐recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).     

Molecular characterization of major histocompatibility complex class II alleles in the common frog, Rana temporaria

Major histocompatibility complex (MHC) class II genes, which play a major role in the immune system response, are some of the most polymorphic genes in vertebrates. We developed polymerase chain reaction primers for part of the second exon of an expressed MHC class II gene in the common frog, Rana temporaria. We genotyped this locus in five frog populations in southeast England and detected eight alleles in 215 individuals. Five or six alleles were detected in each population with a maximum of two alleles per individual, indicating that only a single locus was amplified. We also inferred the possible existence of a null allele. There were 23 variable nucleotide sites (out of 136) and 13 variable amino acid sites (out of 44), many of which corresponded to amino acids involved in antigen recognition. We detected a significant excess of nonsynonymous substitutions at antigen binding sites, indicating that this gene is under positive selection. The level of variation we found was similar to that in other amphibian MHC class II loci, such as those in Bombina bombina, Xenopus laevis and Ambystoma tigrinum.     

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.     

A microsatellite‐based linkage map for song sparrows (Melospiza melodia)

Although linkage maps are important tools in evolutionary biology, their availability for wild populations is limited. The population of song sparrows (Melospiza melodia) on Mandarte Island, Canada, is among the more intensively studied wild animal populations. Its long‐term pedigree data, together with extensive genetic sampling, have allowed the study of a range of questions in evolutionary biology and ecology. However, the availability of genetic markers has been limited. We here describe 191 new microsatellite loci, including 160 high‐quality polymorphic autosomal, 7 Z‐linked and 1 W‐linked markers. We used these markers to construct a linkage map for song sparrows with a total sex‐averaged map length of 1731 cM and covering 35 linkage groups, and hence, these markers cover most of the 38–40 chromosomes. Female and male map lengths did not differ significantly. We then bioinformatically mapped these loci to the zebra finch (Taeniopygia guttata) genome and found that linkage groups were conserved between song sparrows and zebra finches. Compared to the zebra finch, marker order within small linkage groups was well conserved, whereas the larger linkage groups showed some intrachromosomal rearrangements. Finally, we show that as expected, recombination frequency between linked loci explained the majority of variation in gametic phase disequilibrium. Yet, there was substantial overlap in gametic phase disequilibrium between pairs of linked and unlinked loci. Given that the microsatellites described here lie on 35 of the 38–40 chromosomes, these markers will be useful for studies in this species, as well as for comparative genomics studies with other species.