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.     

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.     

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

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.     

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.     

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.     

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.     

Isolation and characterization of 145 polymorphic microsatellite loci for the common frog (Rana temporaria)

We describe primers and polymerase chain reaction conditions to amplify 145 di‐, tri‐ and tetranucleotide microsatellite loci from the common frog (Rana temporaria), a species commonly used as a model in ecological and evolutionary research. Primers were tested on 46 individuals from two Fennoscandian populations and yielded an average of six to nine alleles per locus (range = 1–30) depending on the population. Average observed heterozygosities in the two populations were 0.16 (range = 0–0.91) and 0.36 (range = 0–1).     

Assessment of the potential of temporin peptides from the frog Rana temporaria (Ranidae) as anti‐diabetic agents

Temporin A (FLPLIGRVLSGIL‐NH₂), temporin F (FLPLIGKVLSGIL‐NH₂), and temporin G (FFPVIGRILNGIL‐NH₂), first identified in skin secretions of the frog Rana temporaria, produced concentration‐dependent stimulation of insulin release from BRIN‐BD11 rat clonal β‐cells at concentrations ≥1 nM, without cytotoxicity at concentrations up to 3 μM. Temporin A was the most effective. The mechanism of insulinotropic action did not involve an increase in intracellular Ca²⁺ concentrations. Temporins B, C, E, H, and K were either inactive or only weakly active. Temporins A, F, and G also produced a concentration‐dependent stimulation of insulin release from 1.1B4 human‐derived pancreatic β‐cells, with temporin G being the most potent and effective, and from isolated mouse islets. The data indicate that cationicity, hydrophobicity, and the angle subtended by the charged residues in the temporin molecule are important determinants for in vitro insulinotropic activity. Temporin A and F (1 μM), but not temporin G, protected BRIN‐BD11 cells against cytokine‐induced apoptosis (P < 0.001) and augmented (P < 0.001) proliferation of the cells to a similar extent as glucagon‐like peptide‐1. Intraperitoneal injection of temporin G (75 nmol/kg body weight) together with a glucose load (18 mmol/kg body weight) in C57BL6 mice improved glucose tolerance with a concomitant increase in insulin secretion whereas temporin A and F administration was without significant effect on plasma glucose levels. The study suggests that combination therapy involving agents developed from the temporin A and G sequences may find application in Type 2 diabetes treatment.     

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

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high‐ and low‐altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high‐ and low‐altitude sites using three common temperature treatments. Local adaptation was assessed using Q_ST–F_ST analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high‐altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.     

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.     

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.     

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.     

Microsatellite variation and fine-scale population structure in the wood frog (Rana sylvatica)

We investigated genetic population structure in wood frogs (Rana sylvatica) from a series of Prairie Pothole wetlands in the northern Great Plains. Amphibians are often thought to exist in demographic metapopulations, which require some movement between populations, yet genetic studies have revealed strong subdivision among populations, even at relatively fine scales (several km). Wood frogs are highly philopatric and studies of dispersal suggest that they may exhibit subdivision on a scale of approximately 1-2 km. We used microsatellites to examine population structure among 11 breeding assemblages separated by as little as 50 m up to approximately 5.5 km, plus one population separated from the others by 20 km. We found evidence for differentiation at the largest distances we examined and among a few neighbouring ponds, but most populations were strikingly similar in allele frequencies, suggesting high gene flow among all but the most distant populations. We hypothesize that the few significant differences among neighbouring populations at the finest scale may be a transient effect of extinction-recolonization founder events, driven by periodic drying of wetlands in this hydrologically dynamic landscape.     

The paternal‐sex‐ratio (PSR) chromosome in natural populations of Nasonia (Hymenoptera: Chalcidoidea)

Selfish genetic elements may be important in promoting evolutionary change. Paternal sex ratio (PSR) is a selfish B chromosome that causes all‐male families in the haplodiploid parasitic wasp Nasonia vitripennis, by inducing paternal genome loss in fertilized eggs. The natural distribution and frequency of this chromosome in North American populations of N. vitripennis was investigated using a combination of phenotypic and molecular assays. Sampling throughout North America failed to recover PSR except from populations in the Great Basin area of western North America. Extensive sampling of Great Basin populations revealed PSR in frequencies ranging from 0 to 6% at different collection sites, and extended its distribution to Idaho and Wyoming. Intensive sampling in upstate New York did not detect the chromosome. Frequencies of the maternal‐sex ratio distorter (MSR), son killer (SK) and virgin females ranged from 0 to 12%. Paternal sex ratio may be restricted to the Great Basin because its spread is hampered by geographical barriers, or because populations in other areas are not conducive to PSR maintenance. However, it cannot be ruled out that PSR occurs in other regions at very low frequencies. The apparent limited distribution and low frequency of PSR suggest that it will have relatively little impact on genome evolution in Nasonia.     

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.     

Sexual differentiation in the developing mouse brain: contributions of sex chromosome genes

Neural sexual differentiation begins during embryogenesis and continues after birth for a variable amount of time depending on the species and brain region. Because gonadal hormones were the first factors identified in neural sexual differentiation, their role in this process has eclipsed investigation of other factors. Here, we use a mouse with a spontaneous translocation that produces four different unique sets of sex chromosomes. Each genotype has one normal X‐chromosome and a unique second sex chromosome creating the following genotypes: XY^*x, XX, XY^*, XX^Y*. This Y^* mouse line is used by several laboratories to study two human aneuploid conditions: Turner and Klinefelter syndromes. As sex chromosome number affects behavior and brain morphology, we surveyed brain gene expression at embryonic days 11.5 and 18.5 to isolate X‐chromosome dose effects in the developing brain as possible mechanistic changes underlying the phenotypes. We compared gene expression differences between gonadal males and females as well as individuals with one vs. two X‐chromosomes. We present data showing, in addition to genes reported to escape X‐inactivation, a number of autosomal genes are differentially expressed between the sexes and in mice with different numbers of X‐chromosomes. Based on our results, we can now identify the genes present in the region around the chromosomal break point that produces the Y^* model. Our results also indicate an interaction between gonadal development and sex chromosome number that could further elucidate the role of sex chromosome genes and hormones in the sexual differentiation of behavior.     

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.