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.     

Using RAD‐seq to develop sex‐linked markers in a haplodiplontic alga

For many taxa, including isomorphic haplodiplontic macroalgae, determining sex and ploidy is challenging, thereby limiting the scope of some population demographic and genetic studies. Here, we used double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) to identify sex‐linked molecular markers in the widespread red alga Agarophyton vermiculophyllum. In the ddRAD‐seq library, we included 10 female gametophytes, 10 male gametophytes, and 16 tetrasporophytes from one native and one non‐native site (N = 40 gametophytes and N = 32 tetrasporophytes total). We identified seven putatively female‐linked and 19 putatively male‐linked sequences. Four female‐ and eight male‐linked markers amplified in all three life cycle stages. Using one female‐ and one male‐linked marker that were sex‐specific, we developed a duplex PCR and tested the efficacy of this assay on a subset of thalli sampled at two sites in the non‐native range. We confirmed ploidy based on the visual observation of reproductive structures and previous microsatellite genotyping at 10 polymorphic loci. For 32 vegetative thalli, we were able to assign sex and confirm ploidy in these previously genotyped thalli. These markers will be integral to ongoing studies of A. vermiculophyllum invasion. We discuss the utility of RAD‐seq over other approaches previously used, such as RAPDs (random amplified polymorphic DNA), for future work designing sex‐linked markers in other haplodiplontic macroalgae for which genomes are lacking.     

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.     

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.     

CONSERVED SEX CHROMOSOMES ACROSS ADAPTIVELY RADIATED ANOLIS LIZARDS

Vertebrates possess diverse sex‐determining systems, which differ in evolutionary stability among particular groups. It has been suggested that poikilotherms possess more frequent turnovers of sex chromosomes than homoiotherms, whose effective thermoregulation can prevent the emergence of the sex reversals induced by environmental temperature. Squamate reptiles used to be regarded as a group with an extensive variability in sex determination; however, we document how the rather old radiation of lizards from the genus Anolis, known for exceptional ecomorphological variability, was connected with stability in sex chromosomes. We found that 18 tested species, representing most of the phylogenetic diversity of the genus, share the gene content of their X chromosomes. Furthermore, we discovered homologous sex chromosomes in species of two genera (Sceloporus and Petrosaurus) from the family Phrynosomatidae, serving here as an outgroup to Anolis. We can conclude that the origin of sex chromosomes within iguanas largely predates the Anolis radiation and that the sex chromosomes of iguanas remained conserved for a significant part of their evolutionary history. Next to therian mammals and birds, Anolis lizards therefore represent another adaptively radiated amniote clade with conserved sex chromosomes. We argue that the evolutionary stability of sex‐determining systems may reflect an advanced stage of differentiation of sex chromosomes rather than thermoregulation strategy.     

Substantial differences in bias between single‐digest and double‐digest RAD‐seq libraries: A case study

The trade‐offs of using single‐digest vs. double‐digest restriction site‐associated DNA sequencing (RAD‐seq) protocols have been widely discussed. However, no direct empirical comparisons of the two methods have been conducted. Here, we sampled a single population of Gulf pipefish (Syngnathus scovelli) and genotyped 444 individuals using RAD‐seq. Sixty individuals were subjected to single‐digest RAD‐seq (sdRAD‐seq), and the remaining 384 individuals were genotyped using a double‐digest RAD‐seq (ddRAD‐seq) protocol. We analysed the resulting Illumina sequencing data and compared the two genotyping methods when reads were analysed either together or separately. Coverage statistics, observed heterozygosity, and allele frequencies differed significantly between the two protocols, as did the results of selection components analysis. We also performed an in silico digestion of the Gulf pipefish genome and modelled five major sources of bias: PCR duplicates, polymorphic restriction sites, shearing bias, asymmetric sampling (i.e., genotyping fewer individuals with sdRAD‐seq than with ddRAD‐seq) and higher major allele frequencies. This combination of approaches allowed us to determine that polymorphic restriction sites, an asymmetric sampling scheme, mean allele frequencies and to some extent PCR duplicates all contribute to different estimates of allele frequencies between samples genotyped using sdRAD‐seq versus ddRAD‐seq. Our finding that sdRAD‐seq and ddRAD‐seq can result in different allele frequencies has implications for comparisons across studies and techniques that endeavour to identify genomewide signatures of evolutionary processes in natural populations.     

Dynamic sex chromosomes in Old World chameleons (Squamata: Chamaeleonidae)

Much of our current state of knowledge concerning sex chromosome evolution is based on a handful of ‘exceptional’ taxa with heteromorphic sex chromosomes. However, classifying the sex chromosome systems of additional species lacking easily identifiable, heteromorphic sex chromosomes is indispensable if we wish to fully understand the genesis, degeneration and turnover of vertebrate sex chromosomes. Squamate reptiles (lizards and snakes) are a potential model clade for studying sex chromosome evolution as they exhibit a suite of sex‐determining modes yet most species lack heteromorphic sex chromosomes. Only three (of 203) chameleon species have identified sex chromosome systems (all with female heterogamety, ZZ/ZW). This study uses a recently developed method to identify sex‐specific genetic markers from restriction site‐associated DNA sequence (RADseq) data, which enables the identification of sex chromosome systems in species lacking heteromorphic sex chromosomes. We used RADseq and subsequent PCR validation to identify an XX/XY sex chromosome system in the veiled chameleon (Chamaeleo calyptratus), revealing a novel transition in sex chromosome systems within the Chamaeleonidae. The sex‐specific genetic markers identified here will be essential in research focused on sex‐specific, comparative, functional and developmental evolutionary questions, further promoting C. calyptratus’ utility as an emerging model organism.     

Identifying homomorphic sex chromosomes from wild‐caught adults with limited genomic resources

We demonstrate a genotyping‐by‐sequencing approach to identify homomorphic sex chromosomes and their homolog in a distantly related reference genome, based on noninvasive sampling of wild‐caught individuals, in the moor frog Rana arvalis. Double‐digest RADseq libraries were generated using buccal swabs from 30 males and 21 females from the same population. Search for sex‐limited markers from the unfiltered data set (411 446 RAD tags) was more successful than searches from a filtered data set (33 073 RAD tags) for markers showing sex differences in heterozygosity or in allele frequencies. Altogether, we obtained 292 putatively sex‐linked RAD loci, 98% of which point to male heterogamety. We could map 15 of them to the Xenopus tropicalis genome, all but one on chromosome pair 1, which seems regularly co‐opted for sex determination among amphibians. The most efficient mapping strategy was a three‐step hierarchical approach, where R. arvalis reads were first mapped to a low‐coverage genome of Rana temporaria (17 My divergence), then the R. temporaria scaffolds to the Nanorana parkeri genome (90 My divergence), and finally the N. parkeri scaffolds to the X. tropicalis genome (210 My). We validated our conclusions with PCR primers amplifying part of Dmrt1, a candidate sex determination gene mapping to chromosome 1: a sex‐diagnostic allele was present in all 30 males but in none of the 21 females. Our approach is likely to be productive in many situations where biological samples and/or genomic resources are limited.     

A mid‐Pleistocene rainforest corridor enabled synchronous invasions of the Atlantic Forest by Amazonian anole lizards

Shifts in the geographic distribution of habitats over time can promote dispersal and vicariance, thereby influencing large‐scale biogeographic patterns and ecological processes. An example is that of transient corridors of suitable habitat across disjunct but ecologically similar regions, which have been associated with climate change over time. Such connections likely played a role in the assembly of tropical communities, especially within the highly diverse Amazonian and Atlantic rainforests of South America. Although these forests are presently separated by open and dry ecosystems, paleoclimatic and phylogenetic evidence suggest that they have been transiently connected in the past. However, little is known about the timing, magnitude and the distribution of former forest connections. We employ sequence data at multiple loci from three codistributed arboreal lizards (Anolis punctatus, Anolis ortonii and Polychrus marmoratus) to infer the phylogenetic relationships among Amazonian and Atlantic Forest populations and to test alternative historical demographic scenarios of colonization and vicariance using coalescent simulations and approximate Bayesian computation (ABC). Data from the better‐sampled Anolis species support colonization of the Atlantic Forest from eastern Amazonia. Hierarchical ABC indicates that the three species colonized the Atlantic Forest synchronously during the mid‐Pleistocene. We find support of population bottlenecks associated with founder events in the two Anolis, but not in P. marmoratus, consistently with their distinct ecological tolerances. Our findings support that climatic fluctuations provided key opportunities for dispersal and forest colonization in eastern South America through the cessation of environmental barriers. Evidence of species‐specific histories strengthens assertions that biological attributes play a role in responses to shared environmental change.     

Character displacement in the midst of background evolution in island populations of Anolis lizards: A spatiotemporal perspective

Negative interactions between species can generate divergent selection that causes character displacement. However, other processes cause similar divergence. We use spatial and temporal replication across island populations of Anolis lizards to assess the importance of negative interactions in driving trait shifts. Previous work showed that the establishment of Anolis sagrei on islands drove resident Anolis carolinensis to perch higher and evolve larger toepads. To further test the interaction's causality and predictability, we resurveyed a subset of islands nine years later. Anolis sagrei had established on one island between surveys. We found that A. carolinensis on this island now perch higher and have larger toepads. However, toepad morphology change on this island was not distinct from shifts on six other islands whose Anolis community composition had not changed. Thus, the presence of A. sagrei only partly explains A. carolinensis trait variation across space and time. We also found that A. carolinensis on islands with previously established A. sagrei now perch higher than a decade ago, and that current A. carolinensis perch height is correlated with A. sagrei density. Our results suggest that character displacement likely interacts with other evolutionary processes in this system, and that temporal data are key to detecting such interactions.     

Population structure of Atlantic mackerel inferred from RAD‐seq‐derived SNP markers: effects of sequence clustering parameters and hierarchical SNP selection

Restriction‐site‐associated DNA sequencing (RAD‐seq) and related methods are revolutionizing the field of population genomics in nonmodel organisms as they allow generating an unprecedented number of single nucleotide polymorphisms (SNPs) even when no genomic information is available. Yet, RAD‐seq data analyses rely on assumptions on nature and number of nucleotide variants present in a single locus, the choice of which may lead to an under‐ or overestimated number of SNPs and/or to incorrectly called genotypes. Using the Atlantic mackerel (Scomber scombrus L.) and a close relative, the Atlantic chub mackerel (Scomber colias), as case study, here we explore the sensitivity of population structure inferences to two crucial aspects in RAD‐seq data analysis: the maximum number of mismatches allowed to merge reads into a locus and the relatedness of the individuals used for genotype calling and SNP selection. Our study resolves the population structure of the Atlantic mackerel, but, most importantly, provides insights into the effects of alternative RAD‐seq data analysis strategies on population structure inferences that are directly applicable to other species.     

Genetic variation in Miscanthus × giganteus and the importance of estimating genetic distance thresholds for differentiating clones

Miscanthus × giganteus (Mxg) is an important bioenergy feedstock crop, however, genetic diversity among legacy cultivars may be severely constrained. Only one introduction from Japan to Denmark of this sterile, triploid, vegetatively propagated crop was recorded in the 1930s. We sought to determine if the Mxg cultivars in North America were all synonyms, and if they were derived from the European introduction. We used 64 nuclear and five chloroplast simple sequence repeat (SSR) markers to estimate genetic similarity for 27 Mxg accessions from North America, and compared them with six accessions from Europe, including the species’ type‐specimen. A subset of accessions was also evaluated by restriction‐site associated DNA sequencing (RAD‐seq). In addition, we assessed the potential of new crosses to increase Mxg genetic diversity by comparing eight new triploid Mxg progeny grown from seed, along with samples of the parental species M. sacchariflorus and M. sinensis. Estimates of genotyping error rates were essential for distinguishing between experimental error and true genotypic differences among accessions. Given differences in estimated error rates and costs per marker for SSRs and RAD‐seq, the former is currently more cost‐effective for determining if two accessions are genetically identical. We concluded that all of the Mxg legacy cultivars were derived via vegetative propagation from a single genet. In contrast with the Mxg legacy cultivars, genetic similarity to the type‐specimen of eight new triploid Mxg progeny ranged from 0.46 to 0.56. Though genetic diversity among the Mxg legacy cultivars is critically low, new crosses can provide much‐needed variation to growers.     

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.     

Identification of the sex‐determining region in flathead grey mullet (Mugil cephalus)

Elucidation of the sex‐determination mechanism in flathead grey mullet (Mugil cephalus) is required to exploit its economic potential by production of genetically determined monosex populations and application of hormonal treatment to parents rather than to the marketed progeny. Our objective was to construct a first‐generation linkage map of the M. cephalus in order to identify the sex‐determining region and sex‐determination system. Deep‐sequencing data of a single male was assembled and aligned to the genome of Nile tilapia (Oreochromis niloticus). A total 245 M. cephalus microsatellite markers were designed, spanning the syntenic tilapia genome assembly at intervals of 10 Mb. In the mapping family of full‐sib progeny, 156 segregating markers were used to construct a first‐generation linkage map of 24 linkage groups (LGs), corresponding to the number of chromosomes. The linkage map spanned approximately 1200 cM with an average inter‐marker distance of 10.6 cM. Markers segregating on LG9 in two independent mapping families showed nearly complete concordance with gender (R² = 0.95). The sex determining locus was fine mapped within an interval of 8.6 cM on LG9. The sex of offspring was determined only by the alleles transmitted from the father, thus indicating an XY sex‐determination system.     

A cryptic heterogametic transition revealed by sex-linked DNA markers in Palearctic green toads

In sharp contrast to birds and mammals, most cold‐blooded vertebrates have homomorphic (morphologically undifferentiated) sex chromosomes. This might result either from recurrent X‐Y recombination (occurring e.g. during occasional events of sex reversal) or from frequent turnovers (during which sex‐determining genes are overthrown by new autosomal mutations). Evidence for turnovers is indeed mounting in fish, but very few have so far been documented in amphibians, possibly because of practical difficulties in identifying sex chromosomes. Female heterogamety (ZW) has long been established in Bufo bufo, based on sex reversal and crossing experiments. Here, we investigate a sex‐linked marker identified from a laboratory cross between Palearctic green toads (Bufo viridis subgroup). The F₁ offspring produced by a female Bufo balearicus and a male Bufo siculus were phenotypically sexed, displaying an even sex ratio. A sex‐specific marker detected in highly reproducible AFLP genotypes was cloned. Sequencing revealed a noncoding, microsatellite‐containing fragment. Reamplification and genotyping of families of this and a reciprocal cross showed B. siculus to be male heterogametic (XY) and suggested the same system for B. balearicus. Our results thus reveal a cryptic heterogametic transition within bufonid frogs and help explain patterns of hybrid fitness within the B. viridis subgroup. Turnovers of genetic sex‐determination systems may be more frequent in amphibians than previously thought and thus contribute to the prevalence of homomorphic sex chromosomes in this group.     

A mobile sex‐determining region,male‐specific haplotypes and rearing environment influence age at maturity in Chinook salmon

Variation in age at maturity is an important contributor to life history and demographic variation within and among species. The optimal age at maturity can vary by sex, and the ability of each sex to evolve towards its fitness optimum depends on the genetic architecture of maturation. Using GWAS of RAD sequencing data, we show that age at maturity in Chinook salmon exhibits sex‐specific genetic architecture, with age at maturity in males influenced by large (up to 20 Mb) male‐specific haplotypes. These regions showed no such effect in females. We also provide evidence for translocation of the sex‐determining gene between two different chromosomes. This has important implications for sexually antagonistic selection, particularly that sex linkage of adaptive genes may differ within and among populations based on chromosomal location of the sex‐determining gene. Our findings will facilitate research into the genetic causes of shifting demography in Chinook salmon as well as a better understanding of sex determination in this species and Pacific salmon in general.     

Climate‐mediated adaptation after mainland colonization of an ancestrally subtropical island lizard,Anolis carolinensis

Climate‐mediated evolution plays an integral role in species migration and range expansion. Gaining a clearer understanding of how climate affects demographic history and adaptation provides fundamental insight into the generation of intra‐ and interspecific diversity. In this study, we used the natural colonization of the green anole (Anolis carolinensis) from the island of Cuba to mainland North America to investigate the role of evolution at the niche, phenotypic and genetic levels after long‐term establishment in a novel environment. The North American green anole occupies a broader range of thermal habitats than its Cuban sister species. We documented niche expansion in the mainland green anole, mediated primarily through adaptation to winter temperatures. Common garden experiments strongly suggest a genetic component to differences in thermal performance found between populations in different temperature regimes. Analysis of geographic variation in population structure based on 53 486 single nucleotide variants from RAD loci revealed increased genetic isolation between populations in different vs. similar thermal environments. Selection scans for environment‐allele correlations reveal 19 genomic loci of known function that may have played a role in the physiological adaptation of A. carolinensis to temperate environments on the mainland.     

EVOLUTION OF SEX DETERMINATION SYSTEMS WITH HETEROGAMETIC MALES AND FEMALES IN SILENE

The plant genus Silene has become a model for evolutionary studies of sex chromosomes and sex‐determining mechanisms. A recent study performed in Silene colpophylla showed that dioecy and the sex chromosomes in this species evolved independently from those in Silene latifolia, the most widely studied dioecious Silene species. The results of this study show that the sex‐determining system in Silene otites, a species related to S. colpophylla, is based on female heterogamety, a sex determination system that is unique among the Silene species studied to date. Our phylogenetic data support the placing of S. otites and S. colpophylla in the subsection Otites and the analysis of ancestral states suggests that the most recent common ancestor of S. otites and S. colpophylla was most probably dioecious. These observations imply that a switch from XX/XY sex determination to a ZZ/ZW system (or vice versa) occurred in the subsection Otites. This is the first report of two different types of heterogamety within one plant genus of this mostly nondioecious plant family.