When sex chromosomes stop recombining, they start to accumulate differences. The sex-limited chromosome (Y or W) especially is expected to degenerate via the loss of nucleotide sequence and the accumulation of repetitive sequences. However, how early signs of degeneration can be detected in a new sex chromosome is still unclear. The sex-determining region of the octoploid strawberries is young, small, and dynamic. Using PacBio HiFi reads, we obtained a chromosome-scale assembly of a female (ZW) Fragaria chiloensis plant carrying the youngest and largest of the known sex-determining region on the W in strawberries. We fully characterized the previously incomplete sex-determining region, confirming its gene content, genomic location, and evolutionary history. Resolution of gaps in the previous characterization of the sex-determining region added 10 kb of sequence including a noncanonical long terminal repeat-retrotransposon; whereas the Z sequence revealed a Harbinger transposable element adjoining the sex-determining region insertion site. Limited genetic differentiation of the sex chromosomes coupled with structural variation may indicate an early stage of W degeneration. The sex chromosomes have a similar percentage of repeats but differ in their repeat distribution. Differences in the pattern of repeats (transposable element polymorphism) apparently precede sex chromosome differentiation, thus potentially contributing to recombination cessation as opposed to being a consequence of it.

Recombination often differs markedly between males and females. Here we present the first analysis of sex-specific recombination in Gasterosteus sticklebacks. Using whole-genome sequencing of 15 crosses between G. aculeatus and G. nipponicus, we localized 698 crossovers with a median resolution of 2.3 kb. We also used a bioinformatic approach to infer historical sex-averaged recombination patterns for both species. Recombination is greater in females than males on all chromosomes, and overall map length is 1.64 times longer in females. The locations of crossovers differ strikingly between sexes. Crossovers cluster toward chromosome ends in males, but are distributed more evenly across chromosomes in females. Suppression of recombination near the centromeres in males causes crossovers to cluster at the ends of long arms in acrocentric chromosomes, and greatly reduces crossing over on short arms. The effect of centromeres on recombination is much weaker in females. Genomic differentiation between G. aculeatus and G. nipponicus is strongly correlated with recombination rate, and patterns of differentiation along chromosomes are strongly influenced by male-specific telomere and centromere effects. We found no evidence for fine-scale correlations between recombination and local gene content in either sex. We discuss hypotheses for the origin of sexual dimorphism in recombination and its consequences for sexually antagonistic selection and sex chromosome evolution.

Sexual differentiation is fundamentally important for reproduction, yet the genetic triggers of this developmental process can vary, even between closely related species. Recent studies have uncovered, for example, variation in the genetic triggers for sexual differentiation within and between species of African clawed frogs (genus Xenopus). Here, we extend these discoveries by demonstrating that yet another sex determination system exists in Xenopus, specifically in the species Xenopus borealis This system evolved recently in an ancestor of X. borealis that had the same sex determination system as X. laevis, a system which itself is newly evolved. Strikingly, the genomic region carrying the sex determination factor in X. borealis is homologous to that of therian mammals, including humans. Our results offer insights into how the genetic underpinnings of conserved phenotypes evolve, and suggest an important role for cooption of genetic building blocks with conserved developmental roles.

Fish sex determination (SD) systems are varied, suggesting evolutionary changes including either multiple evolution origins of genetic SD from nongenetic systems (such as environmental SD) and/or turnover events replacing one genetic system by another. When genetic SD is found, cytological differentiation between the two members of the sex chromosome pair is often minor or undetectable. The turbot (Scophthalmus maximus), a valuable commercial flatfish, has a ZZ/ZW system and a major SD region on linkage group 5 (LG5), but there are also other minor genetic and environmental influences. We here report refined mapping of the turbot SD region, supported by comparative mapping with model fish species, to identify the turbot master SD gene. Six genes were located to the SD region, two of them associated with gonad development (sox2 and dnajc19). All showed a high association with sex within families (P = 0), but not at the population level, so they are probably partially sex-linked genes, but not SD gene itself. Analysis of crossovers in LG5 using two families confirmed a ZZ/ZW system in turbot and suggested a revised map position for the master gene. Genetic diversity and differentiation for 25 LG5 genetic markers showed no differences between males and females sampled from a wild population, suggesting a recent origin of the SD region in turbot. We also analyzed associations with markers of the most relevant sex-related linkage groups in brill (S. rhombus), a closely related species to turbot; the data suggest that an ancient XX/XY system in brill changed to a ZZ/ZW mechanism in turbot.

A master sex-determining gene, the Y chromosome-linked anti-Müllerian hormone (amhy) gene, has been described in two New World atheriniform species but little is known on the distribution, evolution, and function(s) of this gene in other Atheriniformes. Interestingly, amhy has been found to coexist with temperature-dependent sex determination (TSD), providing a unique opportunity to explore the interplay between genotypic and environmental sex determination. In this study, the search for an amhy homolog was extended to an Old World atheriniform, the cobaltcap silverside Hypoatherina tsurugae (Atherinidae). The full sequences, including the coding and noncoding regions, of the autosomal amh (amha) and a putative amhy were obtained. The deduced Amha and Amhy proteins comprised 511 and 340 amino acids (aa), respectively. PCR analysis with genomic DNA from wild adults and from laboratory-reared juveniles revealed a high, but not complete association of ∼95% between amhy and maleness. The spatiotemporal expression of amhy and amha during gonadal sex differentiation was analyzed by qRT-PCR and in situ hybridization (ISH). amhy transcription (in amhy-positive larvae) started before and peaked during histological differentiation of the gonads whereas amha was negligible during the same period in both genotypes. These results demonstrate that the amhy, although with some structural differences in relation to the amhy of some New World atheriniforms, is strongly associated with maleness and probably important for testicular development in this Old World atheriniform. Thus, amhy is a candidate sex determination gene in cobaltcap silverside and it will be key to scrutinize the mechanism of sex determination in this species.

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane's Rule.

The fern Ceratopteris richardii is an important model for studies of sex determination and gamete differentiation in homosporous plants. Here we use RNA-seq to de novo assemble a transcriptome and identify genes differentially expressed in young gametophytes as their sex is determined by the presence or absence of the male-inducing pheromone called antheridiogen. Of the 1,163 consensus differentially expressed genes identified, the vast majority (1,030) are up-regulated in gametophytes treated with antheridiogen. GO term enrichment analyses of these DEGs reveals that a large number of genes involved in epigenetic reprogramming of the gametophyte genome are up-regulated by the pheromone. Additional hormone response and development genes are also up-regulated by the pheromone. This C. richardii gametophyte transcriptome and gene expression dataset will prove useful for studies focusing on sex determination and differentiation in plants.

Neural specification is regulated by one or many transcription factors that control expression of effector genes that mediate function and determine neuronal type. Here we identify a novel role for one conserved proneural factor, the bHLH protein HLH-3, implicated in the specification of sex-specific ventral cord motor neurons in C. elegans Proneural genes act in early stages of neurogenesis in early progenitors, but here, we demonstrate a later role for hlh-3 First, we document that differentiation of the ventral cord type C motor neuron class (VC) within their neuron class, is dynamic in time and space. Expression of VC class-specific and subclass-specific identity genes is distinct through development and is dependent on the VC position along the A-P axis and their proximity to the vulva. Our characterization of the expression of VC class and VC subclass-specific differentiation markers in the absence of hlh-3 function reveals that VC fate specification, differentiation, and morphology requires hlh-3 function. Finally, we conclude that hlh-3 cell-autonomously specifies VC cell fate.

Many species that are extensively studied in the laboratory are less well characterized in their natural habitat, and laboratory strains represent only a small fraction of the variation in a species' genome. Here we investigate genomic variation in 3 natural North American populations of an agricultural pest and a model insect for many scientific disciplines, the tobacco hornworm (Manduca sexta). We show that hornworms from Arizona, Kansas, and North Carolina are genetically distinct, with Arizona being particularly differentiated from the other 2 populations using Illumina whole-genome resequencing. Peaks of differentiation exist across the genome, but here, we focus in on the most striking regions. In particular, we identify 2 likely segregating inversions found in the Arizona population. One inversion on the Z chromosome may enhance adaptive evolution of the sex chromosome. The larger, 8 Mb inversion on chromosome 12 contains a pseudogene which may be involved in the exploitation of a novel hostplant in Arizona, but functional genetic assays will be required to support this hypothesis. Nevertheless, our results reveal undiscovered natural variation and provide useful genomic data for both pest management and evolutionary genetics of this insect species.

The Caenorhabditis elegans somatic gonad differs greatly between the two sexes in its pattern of cell divisions, migration, and differentiation. Despite decades of study, the genetic pathways directing early gonadal development and establishing sexual dimorphism in the gonad remain largely unknown. To help define the genetic networks that regulate gonadal development, we employed cell-specific RNA-seq. We identified transcripts present in the somatic gonadal precursor cells and their daughter cells of each sex at the onset of sexual differentiation. We identified several hundred gonad-enriched transcripts, including the majority of known regulators of early gonadal development, and transgenic reporter analysis confirmed the effectiveness of this approach. Before the division of the somatic gonad precursors, few sex-biased gonadal transcripts were detectable; less than 6 hr later, after their division, we identified more than 250 sex-biased transcripts, of which about a third were enriched in the somatic gonad compared to the whole animal. This indicates that a robust sex-biased developmental program, some of it gonad-specific, initiates in the somatic gonadal precursor cells around the time of their first division. About 10% of male-biased transcripts had orthologs with male-biased expression in the early mouse gonad, suggesting possible conservation of gonad sex differentiation. Cell-specific analysis also identified approximately 70 previously unannotated mRNA isoforms that are enriched in the somatic gonad. Our data illustrate the power of cell-specific transcriptome analysis and suggest that early sex differentiation in the gonad is controlled by a relatively small suite of differentially expressed genes, even after dimorphism has become apparent.

Many salmonids have a male heterogametic (XX/XY) sex determination system, and they are supposed to have a conserved master sex-determining gene (sdY) that interacts at the protein level with Foxl2 leading to the blockage of the synergistic induction of Foxl2 and Nr5a1 of the cyp19a1a promoter. However, this hypothesis of a conserved master sex-determining role of sdY in salmonids is challenged by a few exceptions, one of them being the presence of naturally occurring "apparent" XY Chinook salmon, Oncorhynchus tshawytscha, females. Here, we show that some XY Chinook salmon females have a sdY gene (sdY-N183), with 1 missense mutation leading to a substitution of a conserved isoleucine to an asparagine (I183N). In contrast, Chinook salmon males have both a nonmutated sdY-I183 gene and the missense mutation sdY-N183 gene. The 3-dimensional model of SdY-I183N predicts that the I183N hydrophobic to hydrophilic amino acid change leads to a modification in the SdY β-sandwich structure. Using in vitro cell transfection assays, we found that SdY-I183N, like the wild-type SdY, is preferentially localized in the cytoplasm. However, compared to wild-type SdY, SdY-I183N is more prone to degradation, its nuclear translocation by Foxl2 is reduced, and SdY-I183N is unable to significantly repress the synergistic Foxl2/Nr5a1 induction of the cyp19a1a promoter. Altogether, our results suggest that the sdY-N183 gene of XY Chinook females is nonfunctional and that SdY-I183N is no longer able to promote testicular differentiation by impairing the synthesis of estrogens in the early differentiating gonads of wild Chinook salmon XY females.

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7-10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

The nematode Caenorhabditis elegans is an important model for studies of germ cell biology, including the meiotic cell cycle, gamete specification as sperm or oocyte, and gamete development. Fundamental to those studies is a genome-level knowledge of the germline transcriptome. Here, we use RNA-Seq to identify genes expressed in isolated XX gonads, which are approximately 95% germline and 5% somatic gonadal tissue. We generate data from mutants making either sperm [fem-3(q96)] or oocytes [fog-2(q71)], both grown at 22°. Our dataset identifies a total of 10,754 mRNAs in the polyadenylated transcriptome of XX gonads, with 2748 enriched in spermatogenic gonads, 1732 enriched in oogenic gonads, and the remaining 6274 not enriched in either. These spermatogenic, oogenic, and gender-neutral gene datasets compare well with those of previous studies, but double the number of genes identified. A comparison of the additional genes found in our study with in situ hybridization patterns in the Kohara database suggests that most are expressed in the germline. We also query our RNA-Seq data for differential exon usage and find 351 mRNAs with sex-enriched isoforms. We suggest that this new dataset will prove useful for studies focusing on C. elegans germ cell biology.

Recombination suppression on sex chromosomes often extends in a stepwise manner, generating evolutionary strata of differentiation between sex chromosomes. Sexual antagonism is a widely accepted explanation for evolutionary strata, postulating that sets of genes beneficial in only one sex are successively linked to the sex-determining locus. The anther-smut fungus Microbotryum lychnidis-dioicae has mating-type chromosomes with evolutionary strata, only some of which link mating-type genes. Male and female roles are non-existent in this fungus, but mating-type antagonistic selection can also generate evolutionary strata, although the life cycle of the fungus suggests it should be restricted to few traits. Here, we tested the hypothesis that mating-type antagonism may have triggered recombination suppression beyond mating-type genes in M. lychnidis-dioicae by searching for footprints of antagonistic selection in evolutionary strata not linking mating-type loci. We found that these evolutionary strata (i) were not enriched in genes upregulated in the haploid phase, where cells are of alternative mating types, (ii) carried no gene differentially expressed between mating types, and (iii) carried no genes displaying footprints of specialization in terms of protein sequences (dN/dS) between mating types after recommended filtering. Without filtering, eleven genes showed signs of positive selection in the strata not linking mating-type genes, which constituted an enrichment compared to autosomes, but their functions were not obviously involved in antagonistic selection. Thus, we found no strong evidence that antagonistic selection has contributed to extending recombination suppression beyond mating-type genes. Alternative hypotheses should therefore be explored to improve our understanding of the sex-related chromosome evolution.

Since the autosomal genome is shared between the sexes, sex-specific fitness optima present an evolutionary challenge. While sexually antagonistic selection might favor different alleles within females and males, segregation randomly reassorts alleles at autosomal loci between sexes each generation. This process of homogenization during transmission thus prevents between-sex allelic divergence generated by sexually antagonistic selection from accumulating across multiple generations. However, recent empirical studies have reported high male-female FST statistics. Here, we use a population genetic model to evaluate whether these observations could plausibly be produced by sexually antagonistic selection. To do this, we use both a single-locus model with nonrandom mate choice, and individual-based simulations to study the relationship between strength of selection, degree of between-sex divergence, and the associated genetic load. We show that selection must be exceptionally strong to create measurable divergence between the sexes and that the decrease in population fitness due to this process is correspondingly high. Individual-based simulations with selection genome-wide recapitulate these patterns and indicate that small sample sizes and sampling variance can easily generate substantial male-female divergence. We therefore conclude that caution should be taken when interpreting autosomal allelic differentiation between the sexes.

Coxsackievirus B3 (CVB3) contributes to the development of myocarditis, an inflammatory heart disease that predominates in males, and infection is a cause of unexpected death in young individuals. Although gonadal hormones contribute significantly to sex differences, sex chromosomes may also influence disease. Increasing evidence indicates that Chromosome Y (ChrY) genetic variants can impact biological functions unrelated to sexual differentiation. Using C57BL/6J (B6)-ChrY consomic mice, we show that genetic variation in ChrY has a direct effect on the survival of CVB3-infected animals. This effect is not due to potential Sry-mediated differences in prenatal testosterone exposure or to differences in adult testosterone levels. Furthermore, we show that ChrY polymorphism influences the percentage of natural killer T cells in B6-ChrY consomic strains but does not underlie CVB3-induced mortality. These data underscore the importance of investigating not only the hormonal regulation but also ChrY genetic regulation of cardiovascular disease and other male-dominant, sexually dimorphic diseases and phenotypes.

Exemplary systems for adaptive divergence are often characterized by their large degrees of phenotypic variation. This variation represents the outcome of generations of diversifying selection. However, adaptive radiations can also contain a hierarchy of differentiation nested within them where species display only subtle phenotypic differences that still have substantial effects on ecology, function, and ultimately fitness. Sexual dimorphisms are also common in species displaying adaptive divergence and can be the result of differential selection between sexes that produce ecological differences between sexes. Understanding the genetic basis of subtle variation (between certain species or sexes) is therefore important for understanding the process of adaptive divergence. Using cichlids from the dramatic adaptive radiation of Lake Malawi, we focus on understanding the genetic basis of two aspects of relatively subtle phenotypic variation. This included a morphometric comparison of the patterns of craniofacial divergence between two ecologically similar species in relation to the larger adaptive radiation of Malawi, and male-female morphological divergence between their F2 hybrids. We then genetically map craniofacial traits within the context of sex and locate several regions of the genome that contribute to variation in craniofacial shape that is relevant to sexual dimorphism within species and subtle divergence between closely related species, and possibly to craniofacial divergence in the Malawi radiation as a whole. To enhance our search for candidate genes we take advantage of population genomic data and a genetic map that is anchored to the cichlid genome to determine which genes within our QTL regions are associated with SNPs that are alternatively fixed between species. This study provides a holistic understanding of the genetic underpinnings of adaptive divergence in craniofacial shape.

Various master key regulators (MKRs) that control a binary switch of sex determination (SD) have been found in fish; these provide an excellent model for the study of vertebrate genetic SD. The SD region in flathead grey mullet has been previously mapped to a 1 Mbp region harboring 27 genes, of which one is follicle-stimulating hormone receptor (fshr). Although this gene is involved in gonad differentiation and function, it has not been considered as an MKR of SD. We systematically investigated polymorphism in mullet fshr using DNA shotgun sequences, and compared them between males and females. Capable of encoding nonconservative amino acid substitutions, c.1732G>A and c.1759T>G exhibited association with sex on a population level (N = 83; P ≤ 6.7 × 10-19). Hence, 1732 A and 1759 G represent a male-specific haplotype of the gene, designated as "fshry." Additional flanking SNPs showed a weaker degree of association with sex, delimiting the SD critical region to 143 nucleotides on exon 14. Lack of homozygotes for fshry, and the resulting divergence from Hardy-Weinberg equilibrium (N = 170; P ≤ 3.9 × 10-5), were compatible with a male heterogametic model (XY/XX). Capable of replacing a phenylalanine with valine, c.1759T>G alters a conserved position across the sixth transmembrane domain of vertebrate FSHRs. Amino acid substitutions in this position in vertebrates are frequently associated with constant receptor activation and consequently with FSH/FSHR signaling alteration; thus, indicating a potential role of fshr as an MKR of SD.

Genome size varies across the tree of life, with no clear correlation to organismal complexity or coding sequence, but with differences in non-coding regions. Phylogenetic methods have recently been incorporated to further disentangle this enigma, yet most of these studies have focused on widely diverged species. Few have compared patterns of genome size change in closely related species with known structural differences in the genome. As a consequence, the relationship between genome size and differences in chromosome number or inter-sexual differences attributed to XY systems are largely unstudied. We hypothesize that structural differences associated with chromosome number and X-Y chromosome differentiation, should result in differing rates and patterns of genome size change. In this study, we utilize the subgenera within the Drosophila to ask if patterns and rates of genome size change differ between closely related species with differences in chromosome numbers and states of the XY system. Genome sizes for males and females of 152 species are used to answer these questions (with 92 newly added or updated estimates). While we find no relationship between chromosome number and genome size or chromosome number and inter-sexual differences in genome size, we find evidence for differing patterns of genome size change between the subgenera, and increasing rates of change throughout time. Estimated shifts in rates of change in sex differences in genome size occur more often in Sophophora and correspond to known neo-sex events.