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.     

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 August 2012 – 30 September 2012

This article documents the addition of 83 microsatellite marker loci and 96 pairs of single‐nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Bembidion lampros, Inimicus japonicus, Lymnaea stagnalis, Panopea abbreviata, Pentadesma butyracea, Sycoscapter hirticola and Thanatephorus cucumeris (anamorph: Rhizoctonia solani). These loci were cross‐tested on the following species: Pentadesma grandifolia and Pentadesma reyndersii. This article also documents the addition of 96 sequencing primer pairs and 88 allele‐specific primers or probes for Plutella xylostella.     

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.     

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.     

Small‐ and large‐scale heterogeneity in genetic variation across the collard flycatcher genome: implications for estimating genetic diversity in nonmodel organisms

Population genetic studies in nonmodel organisms are often hampered by a lack of reference genomes that are essential for whole‐genome resequencing. In the light of this, genotyping methods have been developed to effectively eliminate the need for a reference genome, such as genotyping by sequencing or restriction site‐associated DNA sequencing (RAD‐seq). However, what remains relatively poorly studied is how accurately these methods capture both average and variation in genetic diversity across an organism's genome. In this issue of Molecular Ecology Resources, Dutoit et al. (2016) use whole‐genome resequencing data from the collard flycatcher to assess what factors drive heterogeneity in nucleotide diversity across the genome. Using these data, they then simulate how well different sequencing designs, including RAD sequencing, could capture most of the variation in genetic diversity. They conclude that for evolutionary and conservation‐related studies focused on the estimating genomic diversity, researchers should emphasize the number of loci analysed over the number of individuals sequenced.     

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

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.     

Using genotyping data to assign markers to their chromosome type and to infer the sex of individuals: a Bayesian model‐based classifier

The recent democratization of next‐generation‐sequencing‐based approaches towards nonmodel species has made it cost‐effective to produce large genotyping data sets for a wider range of species. However, when no detailed genome assembly is available, poor knowledge about the organization of the markers within the genome might hamper the optimal use of this abundant information. At the most basic level of genomic organization, the type of chromosome (autosomes, sex chromosomes, mitochondria or chloroplast in plants) may remain unknown for most markers which might be limiting or even misleading in some applications, particularly in population genetics. Conversely, the characterization of sex‐linked markers allows molecular sexing of the individuals. In this study, we propose a Bayesian model‐based classifier named detsex, to assign markers to their chromosome type and/or to perform sexing of individuals based on genotyping data. The performance of detsex is further evaluated by a comprehensive simulation study and by the analysis of real data sets from various origins (microsatellite and SNP data derived from genotyping assay designs and NGS experiments). Irrespective of the origin of the markers or the size of the data set, detsex was proved efficient (i) to identify the sex‐linked markers, (ii) to perform molecular sexing of the individuals and (iii) to perform basic quality check of the genotyping data sets. The underlying structure of the model also allows to consider each of these potential applications either separately or jointly.     

Permanent genetic resources added to Molecular Ecology Resources Database 1 October 2011-30 November 2011

This article documents the addition of 139 microsatellite marker loci and 90 pairs of single‐nucleotide polymorphism sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Aglaoctenus lagotis, Costus pulverulentus, Costus scaber, Culex pipiens, Dascyllus marginatus, Lupinus nanus Benth, Phloeomyzus passerini, Podarcis muralis, Rhododendron rubropilosum Hayata var. taiwanalpinum and Zoarces viviparus. These loci were cross‐tested on the following species: Culex quinquefasciatus, Rhododendron pseudochrysanthum Hay. ssp. morii (Hay.) Yamazaki and R. pseudochrysanthum Hayata. This article also documents the addition of 48 sequencing primer pairs and 90 allele‐specific primers for Engraulis encrasicolus.     

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.     

A universal and reliable assay for molecular sex identification of three‐spined sticklebacks (Gasterosteus aculeatus)

In heterogametic species, biological differences between the two sexes are ubiquitous, and hence, errors in sex identification can be a significant source of noise and bias in studies where sex‐related sources of variation are of interest or need to be controlled for. We developed and validated a universal multimarker assay for reliable sex identification of three‐spined sticklebacks (Gasterosteus aculeatus). The assay makes use of genotype scores from three sex‐linked loci and utilizes Bayesian probabilistic inference to identify sex of the genotyped individuals. The results, validated with 286 phenotypically sexed individuals from six populations of sticklebacks representing all major genetic lineages (cf. Pacific, Atlantic and Japan Sea), indicate that in contrast to commonly used single‐marker‐based sex identification assays, the developed multimarker assay should be 100% accurate. As the markers in the assay can be scored from agarose gels, it provides a quick and cost‐efficient tool for universal sex identification of three‐spined sticklebacks. The general principle of combining information from multiple markers to improve the reliability of sex identification is transferable and can be utilized to develop and validate similar assays for other species.     

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.     

Novel genetic sex markers reveal high frequency of sex reversal in wild populations of the agile frog (Rana dalmatina) associated with anthropogenic land use

Populations of ectothermic vertebrates are vulnerable to environmental pollution and climate change because certain chemicals and extreme temperatures can cause sex reversal during early ontogeny (i.e. genetically female individuals develop male phenotype or vice versa), which may distort population sex ratios. However, we have troublingly little information on sex reversals in natural populations, due to unavailability of genetic sex markers. Here, we developed a genetic sexing method based on sex‐linked single nucleotide polymorphism loci to study the prevalence and fitness consequences of sex reversal in agile frogs (Rana dalmatina). Out of 125 juveniles raised in laboratory without exposure to sex‐reversing stimuli, 6 showed male phenotype but female genotype according to our markers. These individuals exhibited several signs of poor physiological condition, suggesting stress‐induced sex reversal and inferior fitness prospects. Among 162 adults from 11 wild populations in North‐Central Hungary, 20% of phenotypic males had female genotype according to our markers. These individuals occurred more frequently in areas of anthropogenic land use; this association was attributable to agriculture and less strongly to urban land use. Female‐to‐male sex‐reversed adults had similar body mass as normal males. We recorded no events of male‐to‐female sex reversal either in the laboratory or in the wild. These results support recent suspicions that sex reversal is widespread in nature, and suggest that human‐induced environmental changes may contribute to its pervasiveness. Furthermore, our findings indicate that sex reversal is associated with stress and poor health in early life, but sex‐reversed individuals surviving to adulthood may participate in breeding.     

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 June 2011-31 July 2011

This article documents the addition of 112 microsatellite marker loci and 24 pairs of single nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Agelaius phoeniceus, Austrolittorina cincta, Circus cyaneus, Circus macrourus, Circus pygargus, Cryptocoryne × purpurea Ridl. nothovar. purpurea, Mya arenaria, Patagioenas squamosa, Prochilodus mariae, Scylla serrata and Scytalopus speluncae. These loci were cross-tested on the following species: Cryptocoryne × purpurea nothovar. purpurea, Cryptocoryne affinis, Cryptocoryne ciliata, Cryptocoryne cordata var. cordata, Cryptocoryne elliptica, Cryptocoryne griffithii, Cryptocoryne minima, Cryptocoryne nurii and Cryptocoryne schulzei. This article also documents the addition of 24 sequencing primer pairs and 24 allele-specific primers or probes for Aphis glycines.     

A resource of genome‐wide single‐nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel

Reduced representation genome sequencing such as restriction‐site‐associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single‐nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome‐wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome‐wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19 703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long‐term effective population size was estimated to range between 132 000 and 1 320 000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82 425 loci and 376 918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.     

Estimation of population allele frequencies from next‐generation sequencing data: pool‐versus individual‐based genotyping

Molecular markers produced by next‐generation sequencing (NGS) technologies are revolutionizing genetic research. However, the costs of analysing large numbers of individual genomes remain prohibitive for most population genetics studies. Here, we present results based on mathematical derivations showing that, under many realistic experimental designs, NGS of DNA pools from diploid individuals allows to estimate the allele frequencies at single nucleotide polymorphisms (SNPs) with at least the same accuracy as individual‐based analyses, for considerably lower library construction and sequencing efforts. These findings remain true when taking into account the possibility of substantially unequal contributions of each individual to the final pool of sequence reads. We propose the intuitive notion of effective pool size to account for unequal pooling and derive a Bayesian hierarchical model to estimate this parameter directly from the data. We provide a user‐friendly application assessing the accuracy of allele frequency estimation from both pool‐ and individual‐based NGS population data under various sampling, sequencing depth and experimental error designs. We illustrate our findings with theoretical examples and real data sets corresponding to SNP loci obtained using restriction site–associated DNA (RAD) sequencing in pool‐ and individual‐based experiments carried out on the same population of the pine processionary moth (Thaumetopoea pityocampa). NGS of DNA pools might not be optimal for all types of studies but provides a cost‐effective approach for estimating allele frequencies for very large numbers of SNPs. It thus allows comparison of genome‐wide patterns of genetic variation for large numbers of individuals in multiple populations.     

Searching for sex‐reversals to explain population demography and the evolution of sex chromosomes

Sex determination can be purely genetic (as in mammals and birds), purely environmental (as in many reptiles), or genetic but reversible by environmental factors during a sensitive period in life, as in many fish and amphibians ( Wallace et al. 1999 ; Baroiller et al. 2009a ; Stelkens & Wedekind 2010 ). Such environmental sex reversal (ESR) can be induced, for example, by temperature changes or by exposure to hormone‐active substances. ESR has long been recognized as a means to produce more profitable single‐sex cultures in fish farms ( Cnaani & Levavi‐Sivan 2009 ), but we know very little about its prevalence in the wild. Obviously, induced feminization or masculinization may immediately distort population sex ratios, and distorted sex ratios are indeed reported from some amphibian and fish populations ( Olsen et al. 2006 ; Alho et al. 2008 ; Brykov et al. 2008 ). However, sex ratios can also be skewed by, for example, segregation distorters or sex‐specific mortality. Demonstrating ESR in the wild therefore requires the identification of sex‐linked genetic markers (in the absence of heteromorphic sex chromosomes) followed by comparison of genotypes and phenotypes, or experimental crosses with individuals who seem sex reversed, followed by sexing of offspring after rearing under non‐ESR conditions and at low mortality. In this issue, Alho et al. (2010) investigate the role of ESR in the common frog (Rana temporaria) and a population that has a distorted adult sex ratio. They developed new sex‐linked microsatellite markers and tested wild‐caught male and female adults for potential mismatches between phenotype and genotype. They found a significant proportion of phenotypic males with a female genotype. This suggests environmental masculinization, here with a prevalence of 9%. The authors then tested whether XX males naturally reproduce with XX females. They collected egg clutches and found that some had indeed a primary sex ratio of 100% daughters. Other clutches seemed to result from multi‐male fertilizations of which at least one male had the female genotype. These results suggest that sex‐reversed individuals affect the sex ratio in the following generation. But how relevant is ESR if its prevalence is rather low, and what are the implications of successful reproduction of sex‐reversed individuals in the wild?     

Interpreting the flock algorithm: a reply to Anderson & Barry (2015)

Anderson & Barry (Molecular Ecology Resources, 2015, 10, 1020–1030) compared a reprogrammed version of flock (Duchesne & Turgeon , Molecular Ecology Resources, 2009, 9, 1333–1344), flockture , to a particular model of structure (Pritchard , Genetics, 2000, 155, 945–959) that they propose is equivalent to flock , a non‐MCMC, non‐Bayesian algorithm. They conclude that structure performs better than flockture at clustering individuals from simulated populations with very low level of differentiation (F_ST c. 0.008) based on 15 microsatellites or 96 SNPs. We rather consider that both algorithms failed, with proportions of correct allocations lower than 50%. The authors also noted the slightly better performance of flockture with SNPs at intermediate F_ST values (c. 0.02–0.04) but did not comment. Finally, we disagree with the way the processing time of each program was compared. When compared on the basis of a run leading to a clustering solution, the main output of any clustering algorithm, flock , is, as users can readily experience, much faster. In all, we feel that flock performs at least as well as structure as a clustering algorithm. Moreover, flock has two major assets: high speed and clear, well validated, rules to estimate K, the number of populations. It thus provides a valuable addition to the set of tools at the disposal of the many researchers dealing with real empirical data sets.     

Permanent genetic resources added to Molecular Ecology Resources Database 1 December 2011-31 January 2012

This article documents the addition of 473 microsatellite marker loci and 71 pairs of single‐nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Barteria fistulosa, Bombus morio, Galaxias platei, Hematodinium perezi, Macrocentrus cingulum Brischke (a.k.a. M. abdominalis Fab., M. grandii Goidanich or M. gifuensis Ashmead), Micropogonias furnieri, Nerita melanotragus, Nilaparvata lugens Stål, Sciaenops ocellatus, Scomber scombrus, Spodoptera frugiperda and Turdus lherminieri. These loci were cross‐tested on the following species: Barteria dewevrei, Barteria nigritana, Barteria solida, Cynoscion acoupa, Cynoscion jamaicensis, Cynoscion leiarchus, Cynoscion nebulosus, Cynoscion striatus, Cynoscion virescens, Macrodon ancylodon, Menticirrhus americanus, Nilaparvata muiri and Umbrina canosai. This article also documents the addition of 116 sequencing primer pairs for Dicentrarchus labrax.