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.     

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.     

Dmrt1 polymorphism and sex‐chromosome differentiation in Rana temporaria

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

Dmrt1 polymorphism covaries with sex‐determination patterns in Rana temporaria

Patterns of sex‐chromosome differentiation and gonadal development have been shown to vary among populations of Rana temporaria along a latitudinal transect in Sweden. Frogs from the northern‐boreal population of Ammarnäs displayed well‐differentiated X and Y haplotypes, early gonadal differentiation, and a perfect match between phenotypic and genotypic sex. In contrast, no differentiated Y haplotypes could be detected in the southern population of Tvedöra, where juveniles furthermore showed delayed gonadal differentiation. Here, we show that Dmrt1, a gene that plays a key role in sex determination and sexual development across all metazoans, displays significant sex differentiation in Tvedöra, with a Y‐specific haplotype distinct from Ammarnäs. The differential segment is not only much shorter in Tvedöra than in Ammarnäs, it is also less differentiated and associates with both delayed gonadal differentiation and imperfect match between phenotypic and genotypic sex. Whereas Tvedöra juveniles with a local Y haplotype tend to ultimately develop as males, those without it may nevertheless become functional XX males, but with strongly female‐biased progeny. Our findings suggest that the variance in patterns of sex determination documented in common frogs might result from a genetic polymorphism within a small genomic region that contains Dmrt1. They also substantiate the view that recurrent convergences of sex determination toward a limited set of chromosome pairs may result from the co‐option of small genomic regions that harbor key genes from the sex‐determination pathway.     

No sex‐specific markers detected in bluegill sunfish Lepomis macrochirus by AFLP

The amplified fragment‐length polymorphism (AFLP) technique was used to identify sex‐specific markers in bluegill sunfish. A total of 12 835 loci were produced by 256 primer combinations, of which nine (0·73‰) exhibited presumed sex‐associated amplifications in the pooled samples; however, none of which revealed sex specificity upon individual evaluation.     

Novel evolutionary pathways of sex‐determining mechanisms

Evolutionary transitions between sex‐determining mechanisms (SDMs) are an enigma. Among vertebrates, individual sex (male or female) is primarily determined by either genes (genotypic sex determination, GSD) or embryonic incubation temperature (temperature‐dependent sex determination, TSD), and these mechanisms have undergone repeated evolutionary transitions. Despite this evolutionary lability, transitions from GSD (i.e. from male heterogamety, XX/XY, or female heterogamety, ZZ/ZW) to TSD are an evolutionary conundrum, as they appear to require crossing a fitness valley arising from the production of genotypes with reduced viability owing to being homogametic for degenerated sex chromosomes (YY or WW individuals). Moreover, it is unclear whether alternative (e.g. mixed) forms of sex determination can persist across evolutionary time. It has previously been suggested that transitions would be easy if temperature‐dependent sex reversal (e.g. XX male or XY female) was asymmetrical, occurring only in the homogametic sex. However, only recently has a mechanistic model of sex determination emerged that may allow such asymmetrical sex reversal. We demonstrate that selection for TSD in a realistic sex‐determining system can readily drive evolutionary transitions from GSD to TSD that do not require the production of YY or WW individuals. In XX/XY systems, sex reversal (female to male) occurs in a portion of the XX individuals only, leading to the loss of the Y allele (or chromosome) from the population as XX individuals mate with each other. The outcome is a population of XX individuals whose sex is determined by incubation temperature (TSD). Moreover, our model reveals a novel evolutionarily stable state representing a mixed‐mechanism system that has not been revealed by previous approaches. This study solves two long‐standing puzzles of the evolution of sex‐determining mechanisms by illuminating the evolutionary pathways and endpoints.     

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.     

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

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.     

Density‐dependent sex‐biased development of macroptery in a water strider

In wing‐polymorphic insects, wing morphs differ not only in dispersal capability but also in life history traits because of trade‐offs between flight capability and reproduction. When the fitness benefits and costs of producing wings differ between males and females, sex‐specific trade‐offs can result in sex differences in the frequency of long‐winged individuals. Furthermore, the social environment during development affects sex differences in wing development, but few empirical tests of this phenomenon have been performed to date. Here, I used the wing‐dimorphic water strider Tenagogerris euphrosyne to test how rearing density and sex ratio affect the sex‐specific development of long‐winged dispersing morphs (i.e., sex‐specific macroptery). I also used a full‐sib, split‐family breeding design to assess genetic effects on density‐dependent, sex‐specific macroptery. I reared water strider nymphs at either high or low densities and measured their wing development. I found that long‐winged morphs developed more frequently in males than in females when individuals were reared in a high‐density environment. However, the frequency of long‐winged morphs was not biased according to sex when individuals were reared in a low‐density environment. In addition, full‐sib males and females showed similar macroptery incidence rates at low nymphal density, whereas the macroptery incidence rates differed between full‐sib males and females at high nymphal density. Thus complex gene‐by‐environment‐by‐sex interactions may explain the density‐specific levels of sex bias in macroptery, although this interpretation should be treated with some caution. Overall, my study provides empirical evidence for density‐specific, sex‐biased wing development. My findings suggest that social factors as well as abiotic factors can be important in determining sex‐biased wing development in insects.     

Genomic footprints of an old affair: Single nucleotide polymorphism data reveal historical hybridization and the subsequent evolution of reproductive barriers in two recently diverged grasshoppers with partly overlapping distributions

Secondary contact in close relatives can result in hybridization and the admixture of previously isolated gene pools. However, after an initial period of hybridization, reproductive isolation can evolve through different processes and lead to the interruption of gene flow and the completion of the speciation process. Omocestus minutissimus and O. uhagonii are two closely related grasshoppers with partially overlapping distributions in the Central System mountains of the Iberian Peninsula. To analyse spatial patterns of historical and/or contemporary hybridization between these two taxa and understand how species boundaries are maintained in the region of secondary contact, we sampled sympatric and allopatric populations of the two species and obtained genome‐wide single nucleotide polymorphism data using a restriction site‐associated DNA sequencing approach. We used Bayesian clustering analyses to test the hypothesis of contemporary hybridization in sympatric populations and employed a suite of phylogenomic approaches and a coalescent‐based simulation framework to evaluate alternative hypothetical scenarios of interspecific gene flow. Our analyses rejected the hypothesis of contemporary hybridization but revealed past introgression in the area where the distributions of the two species overlap. Overall, these results point to a scenario of historical gene flow after secondary contact followed by the evolution of reproductive isolation that currently prevents hybridization among sympatric populations.     

Chimaerism detection in bovine twins,triplets and quadruplets using sex chromosome‐linked markers

The phenomenon of chimaerism occurs in the majority of cattle twin pregnancies. The objectives of this study were to develop a powerful diagnostic test for chimaerism in bovine male and female co‐twins using X and Y chromosome‐linked markers and to determine the extent of chimaerism in twins, triplets and quadruplets. We developed a multiplex PCR set of three polymorphic markers on chromosome X (DIK2865, DIK2283, AGLA257), where the presence of >1 and >2 alleles per marker is sufficient to prove chimaerism in males and females, respectively. In addition, a specific segment on chromosome Y (BOV97M) is included in the set to indicate chimaerism in females. Visualization of chimaeric alleles was best for DNA extracted from blood, fair for DNA from vaginal smears and failed for DNA extracted from hair. The power of chimaerism identification using this set of markers for DNA extracted from blood was calculated as 99% in males and virtually 100% in females. All females and males in heterosexual twins, triplets and quadruplets displayed evidence of a chimaeric allele in at least one and maximum of three of three X chromosome markers analysed. In addition, all females showed the presence of the BOV97M segment and were validated as chimaeric by the standard clinical diagnosis of impaired vaginal length. Quantitative PCR analysis of BOV97M copies in all twins vs. their sires showed a mean ratio of 45–68% in females and 39–49% in males, indicating a substantial symmetrical exchange of cells among all co‐twins. The proposed analysis of X and Y chromosome‐linked markers is advantageous to previous methods based on Y chromosome sequences only, because it detects chimaerism in both male and female co‐twins.     

Scale dependence of sex‐specific movement in a small‐bodied stream fish

相似文献   

Genetic mapping of horizontal stripes in Lake Victoria cichlid fishes: benefits and pitfalls of using RAD markers for dense linkage mapping

The genetic dissection of naturally occurring phenotypes sheds light on many fundamental and longstanding questions in speciation and adaptation and is a central research topic in evolutionary biology. Until recently, forward‐genetic approaches were virtually impossible to apply to nonmodel organisms, but the development of next‐generation sequencing techniques eases this difficulty. Here, we use the ddRAD‐seq method to map a colour trait with a known adaptive function in cichlid fishes, well‐known textbook examples for rapid rates of speciation and astonishing phenotypic diversification. A suite of phenotypic key innovations is related to speciation and adaptation in cichlids, among which body coloration features prominently. The focal trait of this study, horizontal stripes, evolved in parallel in several cichlid radiations and is associated with piscivorous foraging behaviour. We conducted interspecific crosses between Haplochromis sauvagei and H. nyererei and constructed a linkage map with 867 SNP markers distributed on 22 linkage groups and total size of 1130.63 cM. Lateral stripes are inherited as a Mendelian trait and map to a single genomic interval that harbours a paralog of a gene with known function in stripe patterning. Dorsolateral and mid‐lateral stripes were always coinherited and are thus under the same genetic control. Additionally, we directly quantify the genotyping error rates in RAD markers and offer guidelines for identifying and dealing with errors. Uncritical marker selection was found to severely impact linkage map construction. Fortunately, by applying appropriate quality control steps, a genotyping accuracy of >99.9% can be reached, thus allowing for efficient linkage mapping of evolutionarily relevant traits.     

Diversification of body shape in Sebastes rockfishes of the north‐east Pacific

Body shape variation is integrally related to many aspects of fish ecology, including locomotion and foraging, and can indicate the functional diversity of fish assemblages. Few studies have thoroughly characterized body shape in a diverse marine fish clade, or investigated both temporal and spatial patterns of variation in body shape disparity. Here, I use digital photographs to measure geometric body shape in 66 species of north‐east Pacific rockfish (Sebastes spp.), including a correction for error introduced by arching of specimens. Different components of interspecific shape variation show associations with fish size, depth habitat, trophic niche and phylogenetic relationships. Overall, the accumulation of body shape disparity appears to have been near‐constant over time, and shows little variation across the latitudinal range of rockfish.     

Assessing the genetic diversity of farmed and wild Rufiji tilapia (Oreochromis urolepis urolepis) populations using ddRAD sequencing

Rufiji tilapia (Oreochromis urolepis urolepis) is an endemic cichlid in Tanzania. In addition to its importance for biodiversity conservation, Rufiji tilapia is also attractive for farming due to its high growth rate, salinity tolerance, and the production of all‐male hybrids when crossed with Nile tilapia (Oreochromis niloticus). The aim of the current study was to assess the genetic diversity and population structure of both wild and farmed Rufiji tilapia populations in order to inform conservation and aquaculture practices. Double‐digest restriction‐site‐associated DNA (ddRAD) libraries were constructed from 195 animals originating from eight wild (Nyamisati, Utete, Mansi, Mindu, Wami, Ruaha, Kibasira, and Kilola) and two farmed (Bwawani and Chemchem) populations. The identified single nucleotide polymorphisms (SNPs; n = 2,182) were used to investigate the genetic variation within and among the studied populations. Genetic distance estimates (F_st) were low among populations from neighboring locations, with the exception of Utete and Chemchem populations (F_st = 0.34). Isolation‐by‐distance (IBD) analysis among the wild populations did not detect any significant correlation signal (r = .05; p‐value = .4) between the genetic distance and the sampling (Euclidean distance) locations. Population structure and putative ancestry were further investigated using both Bayesian (Structure) and multivariate approaches (discriminant analysis of principal components). Both analysis indicated the existence of three distinct genetic clusters. Two cross‐validation scenarios were conducted in order to test the efficiency of the SNP dataset for discriminating between farmed and wild animals or predicting the population of origin. Approximately 95% of the test dataset was correctly classified in the first scenario, while in the case of predicting for the population of origin 68% of the test dataset was correctly classified. Overall, our results provide novel insights regarding the population structure of Rufiji tilapia and a new database of informative SNP markers for both conservation management and aquaculture activities.     

Long‐term effects of calorie restriction on serum sex‐hormone concentrations in men

Calorie restriction (CR) slows aging and consistently reduces circulating sex hormones in laboratory animals. However, nothing is known regarding the long‐term effects of CR with adequate nutrition on serum sex‐hormone concentration in lean healthy humans. In this study, we measured body composition, and serum total testosterone, total 17‐β‐estradiol, sex hormone–binding globulin (SHBG), and dehydroepiandrosterone sulfate (DHEA‐S) concentrations in 24 men (mean age 51.5 ± 13 years), who had been practicing CR with adequate nutrition for an average of 7.4 ± 4.5 years, in 24 age‐ and body fat–matched endurance runners (EX), and 24 age‐matched sedentary controls eating Western diets (WD). We found that both the CR and EX volunteers had significantly lower body fat than the WD volunteers (total body fat, 8.7 ± 4.2%; 10.5 ± 4.4%; 23.2 ± 6.1%, respectively; P = 0.0001). Serum total testosterone and the free androgen index were significantly lower, and SHBG was higher in the CR group than in the EX and WD groups (P ≤ 0.001). Serum 17β‐estradiol and the estradiol:SHBG ratio were both significantly lower in the CR and EX groups than in the WD group (P ≤ 0.005). Serum DHEA‐S concentrations were not different between the three groups. These findings demonstrate that, as in long‐lived CR rodents, long‐term severe CR reduces serum total and free testosterone and increases SHBG concentrations in humans, independently of adiposity. More studies are needed to understand the role of this CR‐mediated reduction in sex hormones in modulating the pathogenesis of age‐associated chronic diseases such as cancer and the aging process itself.     

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.     

Biased sex‐ratio and sex‐biased heterozygote disadvantage affect the maintenance of a genetic polymorphism and the properties of hybrid zones

The evolution of biodiversity is a major issue of modern biology, and it is becoming increasingly topical as the ongoing erosion of diversity puts serious threats on human well‐being. An elementary mechanism that allows maintaining diversity is the interplay between dispersal and heterozygote selective disadvantage, which can lead to self‐sustainable spatial genetic structures and is central to the stability of hybrid zones. Theoretical studies supporting the importance of this mechanism assume a balanced sex‐ratio and a heterozygote disadvantage equally affecting both sexes, despite the multiplicity of empirical evidence that (i) adult sex‐ratio is usually biased towards either male or female and that (ii) heterozygote disadvantage often affects a single sex. We expanded the existing theory by weighting the strength of selection against heterozygote according to the biased in sex‐ratio and in heterozygote disadvantage. The range of conditions allowing for the maintenance of polymorphism can then either double or vanish. We discuss the implications of such finding for birds, mammals and insects diversity. Finally, we provide simple analytical predictions about the effect of those biased on the width and speed of hybrid zones and on the time for the spread of beneficial mutations through such zones.