Conserved ZZ/ZW sex chromosomes in Caribbean croaking geckos (Aristelliger: Sphaerodactylidae)

Current understanding of sex chromosome evolution is largely dependent on species with highly degenerated, heteromorphic sex chromosomes, but by studying species with recently evolved or morphologically indistinct sex chromosomes we can greatly increase our understanding of sex chromosome origins, degeneration and turnover. Here, we examine sex chromosome evolution and stability in the gecko genus Aristelliger. We used RADseq to identify sex‐specific markers and show that four Aristelliger species, spanning the phylogenetic breadth of the genus, share a conserved ZZ/ZW system syntenic with avian chromosome 2. These conserved sex chromosomes contrast with many other gecko sex chromosome systems by showing a degree of stability among a group known for its dynamic sex‐determining mechanisms. Cytogenetic data from A. expectatus revealed homomorphic sex chromosomes with an accumulation of repetitive elements on the W chromosome. Taken together, the large number of female‐specific A. praesignis RAD markers and the accumulation of repetitive DNA on the A. expectatus W karyotype suggest that the Z and W chromosomes are highly differentiated despite their overall morphological similarity. We discuss this paradoxical situation and suggest that it may, in fact, be common in many animal species.     

ANOLIS SEX CHROMOSOMES ARE DERIVED FROM A SINGLE ANCESTRAL PAIR

To explain the frequency and distribution of heteromorphic sex chromosomes in the lizard genus Anolis, we compared the relative roles of sex chromosome conservation versus turnover of sex‐determining mechanisms. We used model‐based comparative methods to reconstruct karyotype evolution and the presence of heteromorphic sex chromosomes onto a newly generated Anolis phylogeny. We found that heteromorphic sex chromosomes evolved multiple times in the genus. Fluorescent in situ hybridization (FISH) of repetitive DNA showed variable rates of Y chromosome degeneration among Anolis species and identified previously undetected, homomorphic sex chromosomes in two species. We confirmed homology of sex chromosomes in the genus by performing FISH of an X‐linked bacterial artificial chromosome (BAC) and quantitative PCR of X‐linked genes in multiple Anolis species sampled across the phylogeny. Taken together, these results are consistent with long‐term conservation of sex chromosomes in the group. Our results pave the way to address additional questions related to Anolis sex chromosome evolution and describe a conceptual framework that can be used to evaluate the origins and evolution of heteromorphic sex chromosomes in other clades.     

Long-term sex reversal by oestradiol in amniotes with heteromorphic sex chromosomes

Oestradiol application during embryonic development reverses the sex of male embryos and results in normal female differentiation in reptiles lacking heteromorphic sex chromosomes, but fails to do so in birds and mammals with heteromorphic sex chromosomes. It is not clear whether the evolution of heteromorphic sex chromosomes in amniotes is accompanied by insensitivity to oestradiol, or if the association between oestradiol insensitivity and heteromorphic sex chromosomes can be attributable to phylogenetic constraints in these taxa. Turtles provide an ideal system to examine the potential relationship between oestradiol insensitivity and sex chromosome heteromorphy, since there are species with heteromorphic sex chromosomes that are closely related to species lacking heteromorphic sex chromosomes. We investigated this relationship by examining the long-term effects of oestradiol-17beta application on sex determination in Staurotypus triporcatus and Staurotypus salvinii, two turtle species with male heterogamety. After raising the turtles in the lab for 3 years, we found follicular and Müllerian duct morphology in oestradiol-treated turtles that was identical to that of untreated females. The lasting sex reversal suggests that the evolutionary transition between systems lacking heteromorphic sex chromosomes and those with heteromorphic sex chromosomes is not constrained by a fundamental mechanistic difference.     

Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing

A major barrier to evolutionary studies of sex determination and sex chromosomes has been a lack of information on the types of sex‐determining mechanisms that occur among different species. This is particularly problematic in groups where most species lack visually heteromorphic sex chromosomes, such as fish, amphibians and reptiles, because cytogenetic analyses will fail to identify the sex chromosomes in these species. We describe the use of restriction site‐associated DNA (RAD) sequencing, or RAD‐seq, to identify sex‐specific molecular markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the lizard Anolis carolinensis. We performed RAD‐seq on seven male and ten female A. carolinensis and found one male‐specific molecular marker. Anolis carolinensis has previously been shown to possess male heterogamety and the recently published A. carolinensis genome facilitated the characterization of the sex‐specific RAD‐seq marker. We validated the male specificity of the new marker using PCR on additional individuals and also found that it is conserved in some other Anolis species. We discuss the utility of using RAD‐seq to identify sex‐determining mechanisms in other species with cryptic or homomorphic sex chromosomes and the implications for the evolution of male heterogamety in Anolis.     

Dynamic gene order on the <Emphasis Type="Italic">Silene latifolia</Emphasis> Y chromosome

Dioecious Silene latifolia evolved heteromorphic sex chromosomes within the last ten million years, making it a species of choice for studies of the early stages of sex chromosome evolution in plants. About a dozen genes have been isolated from its sex chromosomes and basic genetic and deletion maps exist for the X and Y chromosomes. However, discrepancies between Y chromosome maps led to the proposal that individual Y chromosomes may differ in gene order. Here, we use an alternative approach, with fluorescence in situ hybridization (FISH), to locate individual genes on S. latifolia sex chromosomes. We demonstrate that gene order on the Y chromosome differs between plants from two populations. We suggest that dynamic gene order may be a general property of Y chromosomes in species with XY systems, in view of recent work demonstrating that the gene order on the Y chromosomes of humans and chimpanzees are dramatically different.     

Sex chromosomes in flowering plants

Sex chromosomes in dioecious and polygamous plants evolved as a mechanism for ensuring outcrossing to increase genetic variation in the offspring. Sex specificity has evolved in 75% of plant families by male sterile or female sterile mutations, but well-defined heteromorphic sex chromosomes are known in only four plant families. A pivotal event in sex chromosome evolution, suppression of recombination at the sex determination locus and its neighboring regions, might be lacking in most dioecious species. However, once recombination is suppressed around the sex determination region, an incipient Y chromosome starts to differentiate by accumulating deleterious mutations, transposable element insertions, chromosomal rearrangements, and selection for male-specific alleles. Some plant species have recently evolved homomorphic sex chromosomes near the inception of this evolutionary process, while a few other species have sufficiently diverged heteromorphic sex chromosomes. Comparative analysis of carefully selected plant species together with some fish species promises new insights into the origins of sex chromosomes and the selective forces driving their evolution.     

Genetic sex assignment in wild populations using genotyping‐by‐sequencing data: A statistical threshold approach

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping‐by‐sequencing (GBS) and other reduced‐representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex‐specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex‐specific markers in the locus discovery phase. We then used (i) a sex‐specific locus threshold (SSLT) to identify significantly male‐specific loci within our data set; and (ii) a significant sex‐assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male‐specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross‐validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex‐specific loci, we developed a panel of three sex‐specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex‐linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex‐specific markers across a range of species and RRS data sets.     

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.     

Molecular cytogenetic evidence of rearrangements on the Y chromosome of the threespine stickleback fish

To identify the processes shaping vertebrate sex chromosomes during the early stages of their evolution, it is necessary to study systems in which genetic sex determination was recently acquired. Previous cytogenetic studies suggested that threespine stickleback fish (Gasterosteus aculeatus) do not have a heteromorphic sex chromosome pair, although recent genetic studies found evidence of an XY genetic sex-determination system. Using fluorescence in situ hybridization (FISH), we report that the threespine stickleback Y chromosome is heteromorphic and has suffered both inversions and deletion. Using the FISH data, we reconstruct the rearrangements that have led to the current physical state of the threespine stickleback Y chromosome. These data demonstrate that the threespine Y is more degenerate than previously thought, suggesting that the process of sex chromosome evolution can occur rapidly following acquisition of a sex-determining region.     

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.     

Characterization of sex chromosomes in rainbow trout and coho salmon using fluorescence in situ hybridization (FISH)

With the aim of characterizing the sex chromosomes of rainbow trout (Oncorhynchus mykiss) and to identify the sex chromosomes of coho salmon (O. kisutch), we used molecular markers OmyP9, 5S rDNA, and a growth hormone gene fragment (GH2), as FISH probes. Metaphase chromosomes were obtained from lymphocyte cultures from farm specimens of rainbow trout and coho salmon. Rainbow trout sex marker OmyP9 hybridizes on the sex chromosomes of rainbow trout, while in coho salmon, fluorescent signals were localized in the medial region of the long arm of one subtelocentric chromosome pair. This hybridization pattern together with the hybridization of a GH2 intron probe on a chromosome pair having the same morphology, suggests that a subtelocentric pair could be the sex chromosomes in this species. We confirm that in rainbow trout, one of the two loci for 5S rDNA genes is on the X chromosome. In males of this species that lack a heteromorphic sex pair (XX males), the 5S rDNA probe hybridized to both subtelocentrics This finding is discussed in relation to the hypothesis of intraspecific polymorphism of sex chromosomes in rainbow trout.     

Evolution of ZZ/ZW and XX/XY sex-determination systems in the closely related medaka species, <Emphasis Type="Italic">Oryzias hubbsi</Emphasis> and <Emphasis Type="Italic">O. dancena</Emphasis>

A DM-domain gene on the Y chromosome was identified as the sex-determining gene in the medaka, Oryzias latipes, and named DMY (also known as dmrt1bY). However, this gene is absent in most Oryzias fishes, suggesting that closely related species have another sex-determining gene. In fact, it has been demonstrated that the Y chromosome in O. dancena is not homologous to that in O. latipes, whereas both species have an XX/XY sex-determination system. Through a progeny test of sex-reversed fish and a linkage analysis of isolated sex-linked DNA markers, we show that O. hubbsi, which is one of the most closely related species to O. dancena, has a ZZ/ZW system. In addition, genetic and fluorescence in situ hybridization mapping of the sex-linked markers revealed that sex chromosomes in O. hubbsi and O. dancena are not homologous, indicating different origins of these ZW and XY sex chromosomes. Furthermore, we found that O. hubbsi has morphologically heteromorphic sex chromosomes, in which the W chromosome has 4,6-diamidino-2-phenylindole (DAPI)-positive heterochromatin blocks and is larger than the Z chromosome, although such differentiated sex chromosomes have not been observed in other Oryzias species. These findings suggest that a variety of sex-determining mechanisms and sex chromosomes have evolved in Oryzias.     

Accumulation of Y-specific satellite DNAs during the evolution of <Emphasis Type="Italic">Rumex acetosa</Emphasis> sex chromosomes

The study of the molecular structure of young heteromorphic sex chromosomes of plants has shed light on the evolutionary forces that control the differentiation of the X and Y during the earlier stages of their evolution. We have used the model plant Rumex acetosa, a dioecious species with multiple sex chromosomes, 2n = 12 + XX female and 2n = 12 + XY₁Y₂ male, to analyse the significance of repetitive DNA accumulation during the differentiation of the Y. A bulk segregant analysis (BSA) approach allowed us to identify and isolate random amplified polymorphic DNA (RAPD) markers linked to the sex chromosomes. From a total of 86 RAPD markers in the parents, 6 markers were found to be linked to the Ys and 1 to the X. Two of the Y-linked markers represent two AT-rich satellite DNAs (satDNAs), named RAYSII and RAYSIII, that share about 80% homology, as well as with RAYSI, another satDNA of R. acetosa. Fluorescent in situ hybridisation demonstrated that RAYSII is specific for Y₁, whilst RAYSIII is located in different clusters along Y₁ and Y₂. The two satDNAs were only detected in the genome of the dioecious species with XX/XY₁Y₂ multiple sex chromosome systems in the subgenus Acetosa, but were absent from other dioecious species with an XX/XY system of the subgenera Acetosa or Acetosella, as well as in gynodioecious or hermaphrodite species of the subgenera Acetosa, Rumex and Platypodium. Phylogenetic analysis with different cloned monomers of RAYSII and RAYSIII from both R. acetosa and R. papillaris indicate that these two satDNAs are completely separated from each other, and from RAYSI, in both species. The three Y-specific satDNAs, however, evolved from an ancestral satDNA with repeating units of 120 bp, through intermediate satDNAs of 360 bp. The data therefore support the idea that Y-chromosome differentiation and heterochromatinisation in the Rumex species having a multiple sex chromosome system have occurred by different amplification events from a common ancestral satDNA. Since dioecious species with multiple XX/XY₁Y₂ sex chromosome systems of the section Acetosa appear to have evolved from dioecious species with an XX/XY system, the amplification of tandemly repetitive elements in the Ys of the section Acetosa is a recent evolutionary process that has contributed to an increase in the size and differentiation of the already non-recombining Y chromosomes.     

Increased differentiation and reduced gene flow in sex chromosomes relative to autosomes between lineages of the brown creeper Certhia americana

The properties of sex chromosomes, including patterns of inheritance, reduced levels of recombination, and hemizygosity in one of the sexes may result in the faster fixation of new mutations via drift and natural selection. Due to these patterns and processes, the two rules of speciation to describe the genetics of postzygotic isolation, Haldane's rule and the large‐X effect, both explicitly include quicker evolution on sex chromosomes relative to autosomes. Because sex‐linked mutations may be the first to become fixed in the speciation process, and appear to be due to stronger genetic drift (in birds), we may identify pronounced genetic differentiation in sex chromosomes in taxa experiencing recent speciation and diverging mainly via genetic drift. Here, we use nine sex‐linked and 21 autosomal genetic markers to investigate differential divergence and introgression between marker types in Certhia americana. We identified increased levels of genetic differentiation and reduced levels of gene flow on sex chromosomes relative to autosomes. This pattern is similar to those observed in other recently‐divergent avian species, providing another case study of the earlier role of sex chromosomes in divergence, relative to autosomes. Additionally, we identify three markers that may be under selection between Certhia americana lineages.     

Contrasting patterns of X‐chromosome divergence underlie multiple sex‐ratio polymorphisms in stalk‐eyed flies

Sex‐linked segregation distorters cause offspring sex ratios to differ from equality. Theory predicts that such selfish alleles may either go to fixation and cause extinction, reach a stable polymorphism or initiate an evolutionary arms race with genetic modifiers. The extent to which a sex ratio distorter follows any of these trajectories in nature is poorly known. Here, we used X‐linked sequence and simple tandem repeat data for three sympatric species of stalk‐eyed flies (Teleopsis whitei and two cryptic species of T. dalmanni) to infer the evolution of distorting X chromosomes. By screening large numbers of field and recently laboratory‐bred flies, we found no evidence of males with strongly female‐biased sex ratio phenotypes (SR) in one species but high frequencies of SR males in the other two species. In the two species with SR males, we find contrasting patterns of X‐chromosome evolution. T. dalmanni‐1 shows chromosome‐wide differences between sex‐ratio (X^SR) and standard (X^ST) X chromosomes consistent with a relatively old sex‐ratio haplotype based on evidence including genetic divergence, an inversion polymorphism and reduced recombination among X^SR chromosomes relative to X^ST chromosomes. In contrast, we found no evidence of genetic divergence on the X between males with female‐biased and nonbiased sex ratios in T. whitei. Taken with previous studies that found evidence of genetic suppression of sex ratio distortion in this clade, our results illustrate that sex ratio modification in these flies is undergoing recurrent evolution with diverse genomic consequences.     

Different origins of bird and reptile sex chromosomes inferred from comparative mapping of chicken Z-linked genes

Recent progress of chicken genome projects has revealed that bird ZW and mammalian XY sex chromosomes were derived from different autosomal pairs of the common ancestor; however, the evolutionary relationship between bird and reptilian sex chromosomes is still unclear. The Chinese soft-shelled turtle (Pelodiscus sinensis) exhibits genetic sex determination, but no distinguishable (heteromorphic) sex chromosomes have been identified. In order to investigate this further, we performed molecular cytogenetic analyses of this species, and thereby identified ZZ/ZW-type micro-sex chromosomes. In addition, we cloned reptile homologues of chicken Z-linked genes from three reptilian species, the Chinese soft-shelled turtle and the Japanese four-striped rat snake (Elaphe quadrivirgata), which have heteromorphic sex chromosomes, and the Siam crocodile (Crocodylus siamensis), which exhibits temperature-dependent sex determination and lacks sex chromosomes. We then mapped them to chromosomes of each species using FISH. The linkage of the genes has been highly conserved in all species: the chicken Z chromosome corresponded to the turtle chromosome 6q, snake chromosome 2p and crocodile chromosome 3. The order of the genes was identical among the three species. The absence of homology between the bird Z chromosome and the snake and turtle Z sex chromosomes suggests that the origin of the sex chromosomes and the causative genes of sex determination are different between birds and reptiles.     

Fissions,fusions, and translocations shaped the karyotype and multiple sex chromosome constitution of the northeast‐Asian wood white butterfly,Leptidea amurensis

Previous studies have shown a dynamic karyotype evolution and the presence of complex sex chromosome systems in three cryptic Leptidea species from the Western Palearctic. To further explore the chromosomal particularities of Leptidea butterflies, we examined the karyotype of an Eastern Palearctic species, Leptidea amurensis. We found a high number of chromosomes that differed between the sexes and slightly varied in females (i.e. 2n = 118–119 in females and 2n = 122 in males). The analysis of female meiotic chromosomes revealed multiple sex chromosomes with three W and six Z chromosomes. The curious sex chromosome constitution [i.e. W_1–3/Z_1–6 (females) and Z_1–6/Z_1–6 (males)] and the observed heterozygotes for a chromosomal fusion are together responsible for the sex‐specific and intraspecific variability in chromosome numbers. However, in contrast to the Western Palearctic Leptidea species, the single chromosomal fusion and static distribution of cytogenetic markers (18S rDNA and H3 histone genes) suggest that the karyotype of L. amurensis is stable. The data obtained for four Leptidea species suggest that the multiple sex chromosome system, although different among species, is a common feature of the genus Leptidea. Furthermore, inter‐ and intraspecific variations in chromosome numbers and the complex meiotic pairing of these multiple sex chromosomes indicate the role of chromosomal fissions, fusions, and translocations in the karyotype evolution of Leptidea butterflies.     

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.