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.     

Let's talk about sex: A rigorous statistical framework to assign the sex of individuals from reduced‐representation sequencing data

Molecular markers have been used to identify the sex of sampled individuals for several decades, but the time‐consuming development phase prevented their application in many systems. Recently, a growing number of papers have applied reduced‐representation sequencing (RRS) protocols to the identification of sex‐specific markers without the use of test crosses or prior genomic information. While such an approach has great advantages in terms of versatility and ease of use, the “shotgun sequencing” nature of RRS data sets leads to a high amount of missing data, which results in statistical challenges to the confident assignment of sex to individuals. In this issue of Molecular Ecology Resources, Stovall et al. (Molecular Ecology Resources, 18, 2018) provide a statistical framework to answer two questions: (1) how many individuals of one sex only must possess a genotype for this locus to be considered significantly sex‐specific? and (2) How many sex‐specific loci must an individual of unknown sex possess (in a given data set) to be confidently assigned a sex? The statistical pipeline introduced, and applied to samples of New Zealand fur seal (Arctocephalus forsteri) to identify 90 sex‐specific loci, should be broadly applicable to a large number of species and constitutes a nice addition to the molecular ecology toolkit in the genomics era.     

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.     

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.     

Climate change,sex reversal and lability of sex‐determining systems

Sex reversal at high temperatures during embryonic development (e.g., ZZ females) provides the opportunity for new genotypic crosses (e.g., ZZ male × ZZ female). This raises the alarming possibility that climatic warming could lead to the loss of an entire chromosome—one member of the sex chromosome pair (the Y or W)—and the transition of populations to environmental sex determination (ESD). Here we examine the evolutionary dynamics of sex‐determining systems exposed to climatic warming using theoretical models. We found that the loss of sex chromosomes is not an inevitable consequence of sex reversal. A large frequency of ZZ sex reversal (50% reversal from male to female) typically divides the outcome between loss of the ZW genotype and the stable persistence of ZZ males, ZW females and ZZ females. The amount of warming associated with sex chromosome loss depended on several features of wild populations—environmental fluctuation, immigration, heritable variation in temperature sensitivity and differential fecundity of sex‐reversed individuals. Chromosome loss was partially or completely buffered when sex‐reversed individuals suffered a reproductive fitness cost, when immigration occurred or when heritable variation for temperature sensitivity existed. Thus, under certain circumstances, sex chromosomes may persist cryptically in systems where the environment is the predominant influence on sex.     

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.     

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.     

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.     

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.     

Inference of sex‐specific expansion patterns in human populations from Y‐chromosome polymorphism

Studying the current distribution of genetic diversity in humans has important implications for our understanding of the history of our species. We analyzed a set of linked STR and SNP loci from the paternally inherited Y chromosome to infer the past demography of 55 African and Eurasian populations, using both the parametric and nonparametric coalescent‐based methods implemented in the BEAST application. We inferred expansion events in most sedentary farmer populations, while we found constant effective population sizes for both nomadic hunter‐gatherers and seminomadic herders. Our results differed, on several aspects, from previous results on mtDNA and autosomal markers. First, we found more recent expansion patterns in Eurasia than in Africa. This discrepancy, substantially stronger than the ones found with the other kind of markers, may result from a lower effective population size for men, which might have made male‐transmitted markers more sensitive to the out‐of‐Africa bottleneck. Second, we found expansion signals only for sedentary farmers but not for nomadic herders in Central Asia, while these signals were found for both kind of populations in this area when using mtDNA or autosomal markers. Expansion signals in this area may result from spatial expansion processes and may have been erased for the Y chromosome among the herders because of restricted male gene flow. Am J Phys Anthropol 157:217–225, 2015. © 2015 Wiley Periodicals, Inc.     

The sex chromosome system can influence the evolution of sex‐biased dispersal

Sex‐biased dispersal is a much‐discussed feature in literature on dispersal. Diverse hypotheses have been proposed to explain the evolution of sex‐biased dispersal, a difference in dispersal rate or dispersal distance between males and females. An early hypothesis has indicated that it may rely on the difference in sex chromosomes between males and females. However, this proposal was quickly rejected without a real assessment. We propose a new perspective on this hypothesis by investigating the evolution of sex‐biased dispersal when dispersal genes are sex‐linked, that is when they are located on the sex chromosomes. We show that individuals of the heterogametic sex disperse relatively more than do individuals of the homogametic sex when dispersal genes are sex‐linked rather than autosomal. Although such a sex‐biased dispersal towards the heterogametic sex is always observed in monogamous species, the mating system and the location of dispersal genes interact to modulate sex‐biased dispersal in monandry and polyandry. In the context of the multicausality of dispersal, we suggest that sex‐linked dispersal genes can influence the evolution of sex‐biased dispersal.     

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.     

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast‐Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.     

Genetic correlations and sex‐specific adaptation in changing environments

Females and males have conflicting evolutionary interests. Selection favors the evolution of different phenotypes within each sex, yet divergence between the sexes is constrained by the shared genetic basis of female and male traits. Current theory predicts that such “sexual antagonism” should be common: manifesting rapidly during the process of adaptation, and slow in its resolution. However, these predictions apply in temporally stable environments. Environmental change has been shown empirically to realign the direction of selection acting on shared traits and thereby alleviate signals of sexually antagonistic selection. Yet there remains no theory for how common sexual antagonism should be in changing environments. Here, we analyze models of sex‐specific evolutionary divergence under directional and cyclic environmental change, and consider the impact of genetic correlations on long‐run patterns of sex‐specific adaptation. We find that environmental change often aligns directional selection between the sexes, even when they have divergent phenotypic optima. Nevertheless, some forms of environmental change generate persistent sexually antagonistic selection that is difficult to resolve. Our results reinforce recent empirical observations that changing environmental conditions alleviate conflict between males and females. They also generate new predictions regarding the scope for sexually antagonistic selection and its resolution in changing environments.     

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.     

Virus epidemics can lead to a population‐wide spread of intragenomic parasites in a previously parasite‐free asexual population

Sexual reproduction is problematic to explain due to its costs, most notably the twofold cost of sex. Yet, sex has been suggested to be favourable in the presence of proliferating intragenomic parasites given that sexual recombination provides a mechanism to confine the accumulation of deleterious mutations. Kraaijeveld et al. compared recently the accumulation of transposons in sexually and asexually reproducing lines of the same species, the parasitoid wasp Leptopilina clavipes. They discovered that within asexually reproducing wasps, the number of gypsy‐like retrotransposons was increased fourfold, whereas other retrotransposons were not. Interestingly, gypsy‐like retrotransposons are closely related to retroviruses. Endogenous retroviruses are retroviruses that have integrated to the germ line cells and are inherited thereafter vertically. They can also replicate within the genome similarly to retrotransposons as well as form virus particles and infect previously uninfected cells. This highlights the possibility that endogenous retroviruses could play a role in the evolution of sexual reproduction. Here, we show with an individual‐based computational model that a virus epidemic within a previously parasite‐free asexual population may establish a new intragenomic parasite to the population. Moreover and in contrast to other transposons, the possibility of endogenous viruses to maintain a virus epidemic and simultaneously provide resistance to individuals carrying active endogenous viruses selects for the presence of active intragenomic parasites in the population despite their deleterious effects. Our results suggest that the viral nature of certain intragenomic parasites should be taken into account when sex and its benefits are being considered.     

Hybridization in Trypanosoma congolense does not challenge the predominant clonal evolution model. A comment on Tihon et al., 2017, Mol. Ecol.

Tihon et al. have just published in Mol. Ecol. a fine genomic study on Trypanosoma congolense, agent of Animal African Trypanosomiasis. They present very convincing evidence that T. congolense underwent several hybridization events between distinct genetic lines in Zambia. They claim that their data challenge our predominant clonal evolution model (PCE) of micropathogens. We point out the main tenets of our model and show that Tihon et al.'s claim is based on a misinterpretation of the PCE model. Actually, their data strongly support PCE in T. congolense at a microevolutionary level.     

Higher rates of sex evolve under K‐selection

The geographical distribution of sexual and related asexual species has been suggested to correlate with habitat stability; sexual species tend to be in stable habitats (K‐selection), whereas related asexual taxa tend to be in unstable habitats (r‐selection). We test whether this broad‐scale pattern can be re‐created at a microevolutionary scale by experimentally evolving populations of facultatively sexual rotifers under different ecological conditions. Consistent with the pattern in nature, we find that the rate of sex evolves to lower levels in the r‐selected than in K‐selection environments. We consider several different explanations for these results.     

An unusual ZZ/ZW sex chromosome system in Characidium fishes (Crenuchidae, Characiformes) with the presence of rDNA sites

Characidium fishes with a sex chromosome system form a monophyletic group. This work presents data of Characidium lanei from the South Atlantic basin (Brazil), including an unknown type of ZW sex chromosome system for the groups including the presence of rDNA sites on sex chromosomes.