Molecular and Physiological Effects of 17α-methyltestosterone on Sex Differentiation of Black Rockfish, Sebastes schlegelii

It is widely known that all-female fish production holds economic value for aquaculture. Sebastes schlegelii, a preeminent economic species, exhibits a sex dimorphism, with females surpassing males in growth. In this regard, achieving all-female black rockfish production could significantly enhance breeding profitability. In this study, we utilized the widely used male sex-regulating hormone, 17α-methyltestosterone (MT) at three different concentrations (20, 40, and 60 ppm), to produce pseudomales of S. schlegelii for subsequent all-female offspring breeding. Long-term MT administration severely inhibits the growth of S. schlegelii, while short term had no significant impact. Histological analysis confirmed sex reversal at all MT concentrations; however, both medium and higher MT concentrations impaired testis development. MT also influenced sex steroid hormone levels in pseudomales, suppressing E2 while increasing T and 11-KT levels. In addition, a transcriptome analysis revealed that MT down-regulated ovarian-related genes (cyp19a1a and foxl2) while up-regulating male-related genes (amh) in pseudomales. Furthermore, MT modulated the TGF-β signaling and steroid hormone biosynthesis pathways, indicating its crucial role in S. schlegelii sex differentiation. Therefore, the current study provides a method for achieving sexual reversal using MT in S. schlegelii and offers an initial insight into the underlying mechanism of sexual reversal in this species.

De Novo Assembly, Characterization, and Comparative Transcriptome Analysis of Mature Male and Female Gonads of Rabbitfish (Siganus oramin) (Bloch & Schneider, 1801)

The rabbitfish, Siganus oramin, is a commercially important table fish in southeastern China. However, there have been few studies on its gonad development and reproduction regulation. Comparative transcriptome analysis was first performed on adult male and female gonads of S. oramin. In total, 47,070 unigenes were successfully assembled and 22,737 unigenes were successfully annotated. Through comparative transcriptome analysis of male and female gonads, a total of 6722 differentially expressed genes were successfully identified, with 3528 upregulated genes and 3154 downregulated genes in the testes. In addition, 39 differentially expressed reproduction-related genes were identified. Finally, quantitative real-time PCR was used to validate the expression levels of several differentially expressed genes. These results provide important data for further studying the function of reproduction-related genes and the molecular mechanism regulating gonad development and reproduction in S. oramin.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

Sex matters in preclinical research

International Women's Day 2024 has a theme of inclusion. As publishers of preclinical research, we aim to show how inclusion of females in research advances scientific rigor and improves treatment reliability. Sexual reproduction is key to all life across the plant and animal kingdoms. Biological sex takes many forms that are morphologically differentiated during development: stamens versus pistils in plants; color and plumage in birds; fallopian tubes versus vas deferens in mammals; and differences in size, for instance, males are smaller in the fruit fly Drosophila melanogaster. Physical differences may be obvious, but many traits may be more obscure, including hormonal, physiological and metabolic factors. These traits have a big influence on disease and responses to treatment. Thus, we call for improved inclusion, analysis and reporting of sex as a biological variable in preclinical animal modeling research.

It's a trap?! Escape from an ancient, ancestral sex chromosome system and implication of Foxl2 as the putative primary sex-determining gene in a lizard (Anguimorpha; Shinisauridae)

Although sex determination is ubiquitous in vertebrates, mechanisms of sex determination vary from environmentally to genetically influenced. In vertebrates, genetic sex determination is typically accomplished with sex chromosomes. Groups like mammals maintain conserved sex chromosome systems, while sex chromosomes in most vertebrate clades are not conserved across similar evolutionary timescales. One group inferred to have an evolutionarily stable mode of sex determination is Anguimorpha, a clade of charismatic taxa including monitor lizards, Gila monsters, and crocodile lizards. The common ancestor of extant anguimorphs possessed a ZW system that has been retained across the clade. However, the sex chromosome system in the endangered, monotypic family of crocodile lizards (Shinisauridae) has remained elusive. Here, we analyze genomic data to demonstrate that Shinisaurus has replaced the ancestral anguimorph ZW system on LG7 with a novel ZW system on LG3. The linkage group, LG3, corresponds to chromosome 9 in chicken, and this is the first documented use of this syntenic block as a sex chromosome in amniotes. Additionally, this ~1 Mb region harbors approximately 10 genes, including a duplication of the sex-determining transcription factor, Foxl2, critical for the determination and maintenance of sexual differentiation in vertebrates, and thus a putative primary sex-determining gene for Shinisaurus.

It's a Trap?! Escape from an ancient, ancestral sex chromosome system and implication of Foxl2 as the putative primary sex determining gene in a lizard (Anguimorpha; Shinisauridae)

Although sex determination is ubiquitous in vertebrates, mechanisms of sex determination vary from environmentally- to genetically-influenced. In vertebrates, genetic sex determination is typically accomplished with sex chromosomes. Groups like mammals maintain conserved sex chromosome systems, while sex chromosomes in most vertebrate clades aren't conserved across similar evolutionary timescales. One group inferred to have an evolutionarily stable mode of sex determination is Anguimorpha, a clade of charismatic taxa including: monitor lizards, Gila monsters, and crocodile lizards. The common ancestor of extant anguimorphs possessed a ZW system that has been retained across the clade. However, the sex chromosome system in the endangered, monotypic family of crocodile lizards (Shinisauridae) has remained elusive. Here, we analyze genomic data to demonstrate that Shinisaurus has replaced the ancestral anguimorph ZW system on LG7 chromosome with a novel ZW system on LG3. The linkage group LG3 corresponds to chromosome 9 in chicken, and this is the first documented use of this syntenic block as a sex chromosome in amniotes. Additionally, this ~1Mb region harbors approximately 10 genes, including a duplication of the sex-determining transcription factor, Foxl2-critical for the determination and maintenance of sexual differentiation in vertebrates, and thus a putative primary sex determining gene for Shinisaurus.

Concerning the eXclusion in human genomics: the choice of sex chromosome representation in the human genome drastically affects the number of identified variants

Over the past 30 years, a community of scientists has pieced together every base pair of the human reference genome from telomere to telomere. Interestingly, most human genomics studies omit more than 5% of the genome from their analyses. Under "normal" circumstances, omitting any chromosome(s) from an analysis of the human genome would be a cause for concern, with the exception being sex chromosomes. Sex chromosomes in eutherians share an evolutionary origin as an ancestral pair of autosomes. In humans, they share 3 regions of high-sequence identity (∼98-100%), which, along with the unique transmission patterns of the sex chromosomes, introduce technical artifacts in genomic analyses. However, the human X chromosome bears numerous important genes, including more "immune response" genes than any other chromosome, which makes its exclusion irresponsible when sex differences across human diseases are widespread. To better characterize the possible effect of the inclusion/exclusion of the X chromosome on variants called, we conducted a pilot study on the Terra cloud platform to replicate a subset of standard genomic practices using both the CHM13 reference genome and the sex chromosome complement-aware reference genome. We compared the quality of variant calling, expression quantification, and allele-specific expression using these 2 reference genome versions across 50 human samples from the Genotype-Tissue Expression consortium annotated as females. We found that after correction, the whole X chromosome (100%) can generate reliable variant calls, allowing for the inclusion of the whole genome in human genomics analyses as a departure from the status quo of omitting the sex chromosomes from empirical and clinical genomics studies.

De novo assembly, characterization and comparative transcriptome analysis of gonads reveals sex-biased genes in Coreoperca whiteheadi

The wild Coreoperca whiteheadi is considered as the primordial species in sinipercine fish, which has valuable genetic information. Unfortunately, C. whiteheadi was listed as a near-threatened species because of the environmental pollution, over-exploitation and species invasion. Therefore, more genetic information is needed to have a better understanding of gonadal development in C. whiteheadi. Here, the first gonadal transcriptomes analysis of C. whiteheadi was conducted and 277.14 million clean reads were generated. A total of 96,753 unigenes were successfully annotated. By comparing ovary and testis transcriptomes, a total of 21,741 differentially expressed genes (DEGs) were identified, of which 12,057 were upregulated and 9684 were downregulated in testes. Among them, we also identified about 53 differentially expressed sex-biased genes. Subsequently, the expression of twenty-four DEGs were confirmed by real-time fluorescence quantitative PCR. Furthermore, the histological analysis was conducted on ovaries and testes of one-year-old C. whiteheadi. Our results provided basic support for further studies on the function of sex-biased genes and the molecular mechanism of sex determination and reproduction in C. whiteheadi.

A lizard is never late: Squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012 and 2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

Following the Pathway of W Chromosome Differentiation in Triportheus (Teleostei: Characiformes)

In this work, we trace the dynamics of satellite DNAs (SatDNAs) accumulation and elimination along the pathway of W chromosome differentiation using the well-known Triportheus fish model. Triportheus stands out due to a conserved ZZ/ZW sex chromosome system present in all examined species. While the Z chromosome is conserved in all species, the W chromosome is invariably smaller and exhibits differences in size and morphology. The presumed ancestral W chromosome is comparable to that of T. auritus, and contains 19 different SatDNA families. Here, by examining five additional Triportheus species, we showed that the majority of these repetitive sequences were eliminated as speciation was taking place. The W chromosomes continued degeneration, while the Z chromosomes of some species began to accumulate some TauSatDNAs. Additional species-specific SatDNAs that made up the heterochromatic region of both Z and W chromosomes were most likely amplified in each species. Therefore, the W chromosomes of the various Triportheus species have undergone significant evolutionary changes in a short period of time (15-25 Myr) after their divergence.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Sex-Linked Loci on the W Chromosome in the Multi-Ocellated Racerunner (Eremias multiocellata) Confirm Genetic Sex-Determination Stability in Lacertid Lizards

The multi-ocellated racerunner, Eremias multiocellata, was considered to have temperature-dependent sex determination (TSD), as its sex ratio can be influenced at different temperatures. However, such an observation contrasts with recent findings that suggest TSD is less common than previously thought. Here, a genotyping-by-sequencing (GBS) approach was employed to identify sex-linked markers in the E. multiocellata, for which the mechanism choice of TSD or GSD is still controversial. We preliminarily identified 119 sex-linked markers based on sex-associated sex-specific sequences, 97% of which indicated female heterogamety. After eliminating the false positives, 38 sex-linked markers were recognized, all of which showed the ZW/ZZ system. Then, eight of the novel markers were verified by PCR amplification from 15 populations of E. multiocellata, which support the GSD in E. multiocellata without geographic variation. To test the conservation of sex chromosome in Eremias, the eight markers were further cross-tested by PCR amplification in 10 individuals of the Mongolian racerunner (Eremias argus), two of which exhibited cross-utility. The novel sex-linked markers could be mapped on the W chromosome of the sand lizard (Lacerta agilis). Our finding that the sex-linked markers are shared in closely related species, along with a conserved synteny of the W chromosome, further supports the homology and conservation of sex chromosomes in the lacertid lizards.

Analysis in Proceratophrys boiei genome illuminates the satellite DNA content in a frog from the Brazilian Atlantic forest

Satellite DNAs (satDNAs) are one of the most abundant elements in genomes. Characterized as tandemly organized sequences that can be amplified into multiple copies, mainly in heterochromatic regions. The frog P. boiei (2n = 22, ZZ♂/ZW♀) is found in the Brazilian Atlantic forest and has an atypical pattern of heterochromatin distribution when compared to other anuran amphibians, with large pericentromeric blocks on all chromosomes. In addition, females of Proceratophrys boiei have a metacentric sex chromosome W showing heterochromatin in all chromosomal extension. In this work, we performed high-throughput genomic, bioinformatic, and cytogenetic analyses to characterize the satellite DNA content (satellitome) in P. boiei, mainly due to high amount of C-positive heterochromatin and the highly heterochromatic W sex chromosome. After all the analyses, it is remarkable that the satellitome of P. boiei is composed of a high number of satDNA families (226), making P. boiei the frog species with the highest number of satellites described so far. Consistent with the observation of large centromeric C-positive heterochromatin blocks, the genome of P. boiei is enriched with high copy number of repetitive DNAs, with total satDNA abundance comprising 16.87% of the genome. We successfully mapped via Fluorescence in situ hybridization the two most abundant repeats in the genome, PboSat01-176 and PboSat02-192, highlighting the presence of certain satDNAs sequences in strategic chromosomal regions (e.g., centromere and pericentromeric region), which leads to their participation in crucial processes for genomic organization and maintenance. Our study reveals a great diversity of satellite repeats that are driving genomic organization in this frog species. The characterization and approaches regarding satDNAs in this species of frog allowed the confirmation of some insights from satellite biology and a possible relationship with the evolution of sex chromosomes, especially in anuran amphibians, including P. boiei, for which data were not available.

Concerning the eXclusion in human genomics: The choice of sex chromosome representation in the human genome drastically affects number of identified variants

Over the past 30 years, a community of scientists have pieced together every base pair of the human reference genome from telomere-to-telomere. Interestingly, most human genomics studies omit more than 5% of the genome from their analyses. Under 'normal' circumstances, omitting any chromosome(s) from analysis of the human genome would be reason for concern-the exception being the sex chromosomes. Sex chromosomes in eutherians share an evolutionary origin as an ancestral pair of autosomes. In humans, they share three regions of high sequence identity (~98-100%), which-along with the unique transmission patterns of the sex chromosomes-introduce technical artifacts into genomic analyses. However, the human X chromosome bears numerous important genes-including more "immune response" genes than any other chromosome-which makes its exclusion irresponsible when sex differences across human diseases are widespread. To better characterize the effect that including/excluding the X chromosome may have on variants called, we conducted a pilot study on the Terra cloud platform to replicate a subset of standard genomic practices using both the CHM13 reference genome and sex chromosome complement-aware (SCC-aware) reference genome. We compared quality of variant calling, expression quantification, and allele-specific expression using these two reference genome versions across 50 human samples from the Genotype-Tissue-Expression consortium annotated as females. We found that after correction, the whole X chromosome (100%) can generate reliable variant calls-allowing for the inclusion of the whole genome in human genomics analyses as a departure from the status quo of omitting the sex chromosomes from empirical and clinical genomics studies.

The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis)

Fish sex determination can be affected by environmental temperature. This process relies on temperature-sensitive proteins such as heat shock proteins (HSPs). Our previous work found that heat shock cognate proteins (HSCs) may participate in high-temperature associated sex reversal of Chinese tongue sole (Cynoglossus semilaevis). However, the role of hsc genes in responding to high temperature and affecting sex determination/differentiation remains unclear. Here, by using C. semilaevis as model, we identified hsc70 and hsc70-like. hsc70 was abundant in the gonads with a testicular-higher expression at all gonadal development stages except for 6 months post fertilization (mpf). Intriguingly, hsc70-like showed higher expression in testes from 6 mpf on. Both long-term heat treatment during the temperature-sensitive sex-determining period and short-term heat stress at the end of this period caused different expression of hsc70/hsc70-like between sexes. The dual-luciferase assay results also suggested that these genes can respond to high temperature rapidly in vitro. Heat treatment of C. semilaevis testis cells overexpressed with hsc70/hsc70-like could affect the expression of sex-related genes sox9a and cyp19a1a. Our results indicated that hsc70 and hsc70-like were key regulators linking external high-temperature signals with sex differentiation in vivo and provide a new idea for understanding the mechanism by which high temperature affects sex determination/differentiation in teleosts.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

A Germline-Specific Regulator of Mitochondrial Fusion is Required for Maintenance and Differentiation of Germline Stem and Progenitor Cells

Maintenance and differentiation of germline stem and progenitor cells (GSPCs) is important for sexual reproduction. Here, the authors identify zebrafish pld6 as a novel germline-specific gene by cross-analyzing different RNA sequencing results, and find that pld6 knockout mutants develop exclusively into infertile males. In pld6 mutants, GSPCs fail to differentiate and undergo apoptosis, leading to masculinization and infertility. Mitochondrial fusion in pld6-depleted GSPCs is severely impaired, and the mutants exhibit defects in piRNA biogenesis and transposon suppression. Overall, this work uncovers zebrafish Pld6 as a novel germline-specific regulator of mitochondrial fusion, and highlights its essential role in the maintenance and differentiation of GSPCs as well as gonadal development and gametogenesis.

New insights into Xenopus sex chromosome genomics from the Marsabit clawed frog X. borealis

In many groups, sex chromosomes change frequently but the drivers of their rapid evolution are varied and often poorly characterized. With an aim of further understanding sex chromosome turnover, we investigated the polymorphic sex chromosomes of the Marsabit clawed frog, Xenopus borealis, using genomic data and a new chromosome-scale genome assembly. We confirmed previous findings that 54.1 Mb of chromosome 8L is sex-linked in animals from east Kenya and a laboratory strain, but most (or all) of this region is not sex-linked in natural populations from west Kenya. Previous work suggests possible degeneration of the Z chromosomes in the east population because many sex-linked transcripts of this female heterogametic population have female-biased expression, and we therefore expected this chromosome to not be present in the west population. In contrast, our simulations support a model where most or all of the sex-linked portion of the Z chromosome from the east acquired autosomal segregation in the west, and where much genetic variation specific to the large sex-linked portion of the W chromosome from the east is not present in the west. These recent changes are consistent with the hot-potato model, wherein sex chromosome turnover is favoured by natural selection if it purges a (minimally) degenerate sex-specific sex chromosome, but counterintuitively suggest natural selection failed to purge a Z chromosome that has signs of more advanced and possibly more ancient regulatory degeneration. These findings highlight complex evolutionary dynamics of young, rapidly evolving Xenopus sex chromosomes and set the stage for mechanistic work aimed at pinpointing additional sex-determining genes in this group.

Effects of biological sex and pregnancy in experimental autoimmune encephalomyelitis: It's complicated

Experimental autoimmune encephalomyelitis (EAE) can be induced in many animal strains by inoculation with central nervous system antigens and adjuvant or by the passive transfer of lymphocytes reactive with these antigens and is widely used as an animal model for multiple sclerosis (MS). There are reports that female sex and pregnancy affect EAE. Here we review the effects of biological sex and the effects of pregnancy on the clinical features (including disease susceptibility) and pathophysiology of EAE. We also review reports of the possible mechanisms underlying these differences. These include sex-related differences in the immune system and in the central nervous system, the effects of hormones and the sex chromosomes and molecules unique to pregnancy. We also review sex differences in the response to factors that can modify the course of EAE. Our conclusion is that the effects of biological sex in EAE vary amongst animal models and should not be widely extrapolated. In EAE, it is therefore essential that studies looking at the effects of biological sex or pregnancy give full information about the model that is used (i.e. animal strain, sex, the inducing antigen, timing of EAE induction in relation to pregnancy, etc.). In addition, it would be preferable if more than one EAE model were used, to show if any observed effects are generalizable. This is clearly a field that requires further work. However, understanding of the mechanisms of sex differences could lead to greater understanding of EAE, and suggest possible therapies for MS.

How low can you go? Introducing SeXY: sex identification from low-quantity sequencing data despite lacking assembled sex chromosomes

Accurate sex identification is crucial for elucidating the biology of a species. In the absence of directly observable sexual characteristics, sex identification of wild fauna can be challenging, if not impossible. Molecular sexing offers a powerful alternative to morphological sexing approaches. Here, we present SeXY, a novel sex‐identification pipeline, for very low‐coverage shotgun sequencing data from a single individual. SeXY was designed to utilize low‐effort screening data for sex identification and does not require a conspecific sex‐chromosome assembly as reference. We assess the accuracy of our pipeline to data quantity by downsampling sequencing data from 100,000 to 1000 mapped reads and to reference genome selection by mapping to a variety of reference genomes of various qualities and phylogenetic distance. We show that our method is 100% accurate when mapping to a high‐quality (highly contiguous N50 > 30 Mb) conspecific genome, even down to 1000 mapped reads. For lower‐quality reference assemblies (N50 < 30 Mb), our method is 100% accurate with 50,000 mapped reads, regardless of reference assembly quality or phylogenetic distance. The SeXY pipeline provides several advantages over previously implemented methods; SeXY (i) requires sequencing data from only a single individual, (ii) does not require assembled conspecific sex chromosomes, or even a conspecific reference assembly, (iii) takes into account variation in coverage across the genome, and (iv) is accurate with only 1000 mapped reads in many cases.

A Potential Role of the Translation Elongation Factor eef1a1 in Gonadal High-Temperature Perception in Chinese Tongue Sole (Cynoglossus semilaevis)

Simple Summary

The eukaryotic elongation factor 1 alpha (eef1a) gene is vital for protein translation by delivering aminoacylated tRNAs to the A/P site of the ribosome via the GTP-dependent reaction. Here, the Chinese tongue sole (Cynoglossus semilaevis) eef1a1 gene was identified, and its potential role in gonadal high-temperature perception was assessed. The full-length sequence of eef1a1 cDNA was 1809 base pair (bp) encoding a putative protein of 461 amino acids. The expression levels of eef1a1 in the ovary were significantly higher than that in the testis from 6 mpf to 3 ypf. Under high-temperature induction during sex differentiation, eef1a1 was significantly down-regulated in males, while the difference was not detected in females. Furthermore, the rapid response of eef1a1 to environmental high temperature was assessed in vitro. Our findings suggest that C. semilaevis eef1a1 might be essential for the molecular response regulatory network of external temperature affecting internal sex differentiation.

Abstract

The eukaryotic translation elongation factor 1 alpha (eef1a) gene has a well-defined role in protein synthesis. However, its role in external temperature perception and internal sex differentiation and development is still unclear. In this study, eef1a1 was identified and functionally analyzed in Chinese tongue sole (Cynoglossus semilaevis). The eef1a1 cDNA, 1809 bp in length, had a 1386 bp open reading frame (ORF) that encoded a 461 amino acid polypeptide containing one EF-1_alpha domain. eef1a1 expression levels were investigated across different tissues and during gonadal development. In the gonad, eef1a1 showed a sexually dimorphic expression pattern with a statistically higher expression level in the ovary than in the testis from 6 months postfertilization to 3 years postfertilization. Under high temperature (28 °C) treatment during C. semilaevis sex differentiation (from 30 days postfertilization to 3 months postfertilization), eef1a1 was statistically down-regulated in males, while the difference was not detected in females. In addition, the dual-luciferase assay exhibited that eef1a1 can respond to high temperature rapidly. Based on these results, C. semilaevis eef1a1 might have a dual role in the perception of external temperature changes and sex differentiation regulation.

Identification of ancestral sex chromosomes in the frog Glandirana rugosa bearing XX-XY and ZZ-ZW sex-determining systems

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy.

Sequencing of laser captured Z and W chromosomes of the tocantins paradoxical frog (Pseudis tocantins) provides insights on repeatome and chromosomal homology

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome.

Chromosome-Level Genome Assembly Reveals Dynamic Sex Chromosomes in Neotropical Leaf-Litter Geckos (Sphaerodactylidae: Sphaerodactylus)

Sex determination is a critical element of successful vertebrate development, suggesting that sex chromosome systems might be evolutionarily stable across lineages. For example, mammals and birds have maintained conserved sex chromosome systems over long evolutionary time periods. Other vertebrates, in contrast, have undergone frequent sex chromosome transitions, which is even more amazing considering we still know comparatively little across large swaths of their respective phylogenies. One reptile group in particular, the gecko lizards (infraorder Gekkota), shows an exceptional lability with regard to sex chromosome transitions and may possess the majority of transitions within squamates (lizards and snakes). However, detailed genomic and cytogenetic information about sex chromosomes is lacking for most gecko species, leaving large gaps in our understanding of the evolutionary processes at play. To address this, we assembled a chromosome-level genome for a gecko (Sphaerodactylidae: Sphaerodactylus) and used this assembly to search for sex chromosomes among six closely related species using a variety of genomic data, including whole-genome re-sequencing, RADseq, and RNAseq. Previous work has identified XY systems in two species of Sphaerodactylus geckos. We expand upon that work to identify between two and four sex chromosome cis-transitions (XY to a new XY) within the genus. Interestingly, we confirmed two different linkage groups as XY sex chromosome systems that were previously unknown to act as sex chromosomes in tetrapods (syntenic with Gallus chromosome 3 and Gallus chromosomes 18/30/33), further highlighting a unique and fascinating trend that most linkage groups have the potential to act as sex chromosomes in squamates.

Tenuous Transcriptional Threshold of Human Sex Determination. I. SRY and Swyer Syndrome at the Edge of Ambiguity

Male sex determination in mammals is initiated by SRY, a Y-encoded transcription factor. The protein contains a high-mobility-group (HMG) box mediating sequence-specific DNA bending. Mutations causing XY gonadal dysgenesis (Swyer syndrome) cluster in the box and ordinarily arise de novo. Rare inherited variants lead to male development in one genetic background (the father) but not another (his sterile XY daughter). De novo and inherited mutations occur at an invariant Tyr adjoining the motif's basic tail (box position 72; Y127 in SRY). In SRY-responsive cell lines CH34 and LNCaP, de novo mutations Y127H and Y127C reduced SRY activity (as assessed by transcriptional activation of principal target gene Sox9) by 5- and 8-fold, respectively. Whereas Y127H impaired testis-specific enhancer assembly, Y127C caused accelerated proteasomal proteolysis; activity was in part rescued by proteasome inhibition. Inherited variant Y127F was better tolerated: its expression was unperturbed, and activity was reduced by only twofold, a threshold similar to other inherited variants. Biochemical studies of wild-type (WT) and variant HMG boxes demonstrated similar specific DNA affinities (within a twofold range), with only subtle differences in sharp DNA bending as probed by permutation gel electrophoresis and fluorescence resonance-energy transfer (FRET); thermodynamic stabilities of the free boxes were essentially identical. Such modest perturbations are within the range of species variation. Whereas our cell-based findings rationalize the de novo genotype-phenotype relationships, a molecular understanding of inherited mutation Y127F remains elusive. Our companion study uncovers cryptic biophysical perturbations suggesting that the para-OH group of Y127 anchors a novel water-mediated DNA clamp.

Unusual Mammalian Sex Determination Systems: A Cabinet of Curiosities

Therian mammals have among the oldest and most conserved sex-determining systems known to date. Any deviation from the standard XX/XY mammalian sex chromosome constitution usually leads to sterility or poor fertility, due to the high differentiation and specialization of the X and Y chromosomes. Nevertheless, a handful of rodents harbor so-called unusual sex-determining systems. While in some species, fertile XY females are found, some others have completely lost their Y chromosome. These atypical species have fascinated researchers for over 60 years, and constitute unique natural models for the study of fundamental processes involved in sex determination in mammals and vertebrates. In this article, we review current knowledge of these species, discuss their similarities and differences, and attempt to expose how the study of their exceptional sex-determining systems can further our understanding of general processes involved in sex chromosome and sex determination evolution.

Remnant of Unrelated Amniote Sex Chromosomal Linkage Sharing on the Same Chromosome in House Gecko Lizards, Providing a Better Understanding of the Ancestral Super-Sex Chromosome

Comparative chromosome maps investigating sex chromosomal linkage groups in amniotes and microsatellite repeat motifs of a male house gecko lizard (Hemidactylus frenatus, HFR) and a flat-tailed house gecko lizard (H. platyurus, HPL) of unknown sex were examined using 75 bacterial artificial chromosomes (BACs) from chicken and zebra finch genomes. No massive accumulations of microsatellite repeat motifs were found in either of the gecko lizards, but 10 out of 13 BACs mapped on HPL chromosomes were associated with other amniote sex chromosomes. Hybridization of the same BACs onto multiple different chromosome pairs suggested transitions to sex chromosomes across amniotes. No BAC hybridization signals were found on HFR chromosomes. However, HFR diverged from HPL about 30 million years ago, possibly due to intrachromosomal rearrangements occurring in the HFR lineage. By contrast, heterochromatin likely reshuffled patterns between HPL and HFR, as observed from C-positive heterochromatin distribution. Six out of ten BACs showed partial homology with squamate reptile chromosome 2 (SR2) and snake Z and/or W sex chromosomes. The gecko lizard showed shared unrelated sex chromosomal linkages-the remnants of a super-sex chromosome. A large ancestral super-sex chromosome showed a correlation between SR2 and snake W sex chromosomes.

Against the mainstream: exceptional evolutionary stability of ZW sex chromosomes across the fish families Triportheidae and Gasteropelecidae (Teleostei: Characiformes)

Teleost fishes exhibit a breath-taking diversity of sex determination and differentiation mechanisms. They encompass at least nine sex chromosome systems with often low degree of differentiation, high rate of inter- and intra-specific variability, and frequent turnovers. Nevertheless, several mainly female heterogametic systems at an advanced stage of genetic differentiation and high evolutionary stability have been also found across teleosts, especially among Neotropical characiforms. In this study, we aim to characterize the ZZ/ZW sex chromosome system in representatives of the Triportheidae family (Triportheus auritus, Agoniates halecinus, and the basal-most species Lignobrycon myersi) and its sister clade Gasteropelecidae (Carnegiella strigata, Gasteropelecus levis, and Thoracocharax stellatus). We applied both conventional and molecular cytogenetic approaches including chromosomal mapping of 5S and 18S ribosomal DNA clusters, cross-species chromosome painting (Zoo-FISH) with sex chromosome-derived probes and comparative genomic hybridization (CGH). We identified the ZW sex chromosome system for the first time in A. halecinus and G. levis and also in C. strigata formerly reported to lack sex chromosomes. We also brought evidence for possible mechanisms underlying the sex chromosome differentiation, including inversions, repetitive DNA accumulation, and exchange of genetic material. Our Zoo-FISH experiments further strongly indicated that the ZW sex chromosomes of Triportheidae and Gasteropelecidae are homeologous, suggesting their origin before the split of these lineages (approx. 40-70 million years ago). Such extent of sex chromosome stability is almost exceptional in teleosts, and hence, these lineages afford a special opportunity to scrutinize unique evolutionary forces and pressures shaping sex chromosome evolution in fishes and vertebrates in general.

Sex chromosome evolution among amniotes: is the origin of sex chromosomes non-random?

Sex chromosomes are a great example of a convergent evolution at the genomic level, having evolved dozens of times just within amniotes. An intriguing question is whether this repeated evolution was random, or whether some ancestral syntenic blocks have significantly higher chance to be co-opted for the role of sex chromosomes owing to their gene content related to gonad development. Here, we summarize current knowledge on the evolutionary history of sex determination and sex chromosomes in amniotes and evaluate the hypothesis of non-random emergence of sex chromosomes. The current data on the origin of sex chromosomes in amniotes suggest that their evolution is indeed non-random. However, this non-random pattern is not very strong, and many syntenic blocks representing putatively independently evolved sex chromosomes are unique. Still, repeatedly co-opted chromosomes are an excellent model system, as independent co-option of the same genomic region for the role of sex chromosome offers a great opportunity for testing evolutionary scenarios on the sex chromosome evolution under the explicit control for the genomic background and gene identity. Future studies should use these systems more to explore the convergent/divergent evolution of sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

PAX2 + Mesenchymal Origin of Gonadal Supporting Cells Is Conserved in Birds

During embryonic gonadal development, the supporting cell lineage is the first cell type to differentiate, giving rise to Sertoli cells in the testis and pre-granulosa cells in the ovary. These cells are thought to direct other gonadal cell lineages down the testis or ovarian pathways, including the germline. Recent research has shown that, in contrast to mouse, chicken gonadal supporting cells derive from a PAX2/OSR1/DMRT1/WNT4 positive mesenchymal cell population. These cells colonize the undifferentiated genital ridge during early gonadogenesis, around the time that germ cells migrate into the gonad. During the process of somatic gonadal sex differentiation, PAX2 expression is down-regulated in embryonic chicken gonads just prior to up-regulation of testis- and ovary-specific markers and prior to germ cell differentiation. Most research on avian gonadal development has focused on the chicken model, and related species from the Galloanserae clade. There is a lack of knowledge on gonadal sex differentiation in other avian lineages. Comparative analysis in birds is required to fully understand the mechanisms of avian sex determination and gonadal differentiation. Here we report the first comparative molecular characterization of gonadal supporting cell differentiation in birds from each of the three main clades, Galloanserae (chicken and quail), Neoaves (zebra finch) and Palaeognathe (emu). Our analysis reveals conservation of PAX2+ expression and a mesenchymal origin of supporting cells in each clade. Moreover, down-regulation of PAX2 expression precisely defines the onset of gonadal sex differentiation in each species. Altogether, these results indicate that gonadal morphogenesis is conserved among the major bird clades.

Joint identification of sex and sex-linked scaffolds in non-model organisms using low depth sequencing data

Being able to assign sex to individuals and identify autosomal and sex-linked scaffolds are essential in most population genomic analyses. Non-model organisms often have genome assemblies at scaffold-level and lack characterization of sex-linked scaffolds. Previous methods to identify sex and sex-linked scaffolds have relied on synteny between the non-model organism and a closely related species or prior knowledge about the sex of the samples to identify sex-linked scaffolds. In the latter case, the difference in depth of coverage between the autosomes and the sex chromosomes are used. Here, we present "sex assignment through coverage" (SATC), a method to assign sex to samples and identify sex-linked scaffolds from next generation sequencing (NGS) data. The method works for species with a homogametic/heterogametic sex determination system and only requires a scaffold-level reference assembly and sampling of both sexes with whole genome sequencing (WGS) data. We use the sequencing depth distribution across scaffolds to jointly identify: (i) male and female individuals, and (ii) sex-linked scaffolds. This is achieved through projecting the scaffold depths into a low-dimensional space using principal component analysis (PCA) and subsequent Gaussian mixture clustering. We demonstrate the applicability of our method using data from five mammal species and a bird species complex. The method is freely available at https://github.com/popgenDK/SATC as R code and a graphical user interface (GUI).

Sex chromosome evolution from a heteromorphic to a homomorphic system by inter-population hybridization in a frog

Sex chromosomes generally evolve from a homomorphic to heteromorphic state. Once a heteromorphic system is established, the sex chromosome system may remain stable for an extended period. Here, we show the opposite case of sex chromosome evolution from a heteromorphic to a homomorphic system in the Japanese frog Glandirana rugosa. One geographic group, Neo-ZW, has ZZ-ZW type heteromorphic sex chromosomes. We found that its western edge populations, which are geographically close to another West-Japan group with homomorphic sex chromosomes of XX-XY type, showed homozygous genotypes of sex-linked genes in both sexes. Karyologically, no heteromorphic sex chromosomes were identified. Sex-reversal experiments revealed that the males were heterogametic in sex determination. In addition, we identified another similar population around at the southwestern edge of the Neo-ZW group in the Kii Peninsula: the frogs had homomorphic sex chromosomes under male heterogamety, while shared mitochondrial haplotypes with the XY group, which is located in the east and bears heteromorphic sex chromosomes. In conclusion, our study revealed that the heteromorphic sex chromosome systems independently reversed back to or turned over to a homomorphic system around each of the western and southwestern edges of the Neo-ZW group through hybridization with the West-Japan group bearing homomorphic sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

With or Without W? Molecular and Cytogenetic Markers are Not Sufficient for Identification of Environmentally-Induced Sex Reversal in the Bearded Dragon

Transitions from environmental sex determination (ESD) to genotypic sex determination (GSD) require an intermediate step of sex reversal, i.e., the production of individuals with a mismatch between the ancestral genotypic and the phenotypic sex. Among amniotes, the sole well-documented transition in this direction was shown in the laboratory in the central bearded dragon, Pogona vitticeps, where very high incubation temperatures led to the production of females with the male-typical (ZZ) genotype. These sex-reversed females then produced offspring whose sex depended on the incubation temperature. Sex-reversed animals identified by molecular and cytogenetic markers were also reported in the field, and their increasing incidence was speculated as a climate warming-driven transition in sex determination. We show that the molecular and cytogenetic markers normally sex-linked in P. vitticeps are also sex-linked in P. henrylawsoni and P. minor, which points to quite ancient sex chromosomes in this lineage. Nevertheless, we demonstrate, based on a crossing experiment with a male bearded dragon who possesses a mismatch between phenotypic sex and genotype, that the used cytogenetic and molecular markers might not be reliable for the identification of sex reversal. Sex reversal should not be considered as the only mechanism causing a mismatch between genetic sex-linked markers and phenotypic sex, which can emerge also by other processes, here most likely by a rare recombination between regions of sex chromosomes which are normally sex-linked. We warn that sex-linked, even apparently for a long evolutionary time, and sex-specific molecular and cytogenetic markers are not a reliable tool for the identification of sex-reversed individuals in a population and that sex reversal has to be verified by other approaches, particularly by observation of the sex ratio of the progeny.

The rise and fall of the ancient northern pike master sex-determining gene

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Igf3 is essential for ovary differentiation in zebrafish†

Zebrafish gonadal sexual differentiation is an important but poorly understood subject. Previously, we have identified a novel insulin-like growth factor (Igf) named insulin-like growth factor 3 (Igf3) in teleosts. The importance of Igf3 in oocyte maturation and ovulation has been recently demonstrated by us in zebrafish. In this study, we have further found the essential role of Igf3 in gonadal sexual differentiation of zebrafish. A differential expression pattern of igf3 between ovary and testis during sex differentiation (higher level in ovary than in testis) was found in zebrafish. An igf3 knockout zebrafish line was established using TALENs-mediated gene knockout technique. Intriguingly, all igf3 homozygous mutants were males due to the female-to-male sex reversal occurred during sex differentiation. Further analysis showed that Igf3 did not seem to affect the formation of so-called juvenile ovary and oocyte-like germ cells. Oocyte development was arrested at primary growth stage, and the ovary was gradually sex-reversed to testis before 60 day post fertilization (dpf). Such sex reversal was likely due to decreased germ cell proliferation by suppressing PI3K/Akt pathway in early ovaries of igf3 mutants. Estrogen is considered as a master regulator in fish sex differentiation. Here, we found that igf3 expression could be upregulated by estrogen in early stages of ovarian follicles as evidenced in in vitro treatment assays and cyp19a1a mutant zebrafish, and E2 failed to rescue the defects of igf3 mutants in ovarian development, suggesting that Igf3 may serve as a downstream factor of estrogen signaling in sex differentiation. Taken together, we demonstrated that Igf3 is essential for ovary differentiation in zebrafish.

Sex chromosome evolution in frogs-helpful insights from chromosome painting in the genus Engystomops

The differentiation of sex chromosomes is thought to be interrupted by relatively frequent sex chromosome turnover and/or occasional recombination between sex chromosomes (fountain-of-youth model) in some vertebrate groups as fishes, amphibians, and lizards. As a result, we observe the prevalence of homomorphic sex chromosomes in these groups. Here, we provide evidence for the loss of sex chromosome heteromorphism in the Amazonian frogs of the genus Engystomops, which harbors an intriguing history of sex chromosome evolution. In this species complex composed of two named species, two confirmed unnamed species, and up to three unconfirmed species, highly divergent karyotypes are present, and heteromorphic X and Y chromosomes were previously found in two species. We describe the karyotype of a lineage estimated to be the sister of all remaining Amazonian Engystomops (named Engystomops sp.) and perform chromosome painting techniques using one probe for the Y chromosome and one probe for the non-centromeric heterochromatic bands of the X chromosome of E. freibergi to compare three Engystomops karyotypes. The Y probe detected the Y chromosomes of E. freibergi and E. petersi and one homolog of chromosome pair 11 of Engystomops sp., suggesting their common evolutionary origin. The X probe showed no interspecific hybridization, revealing that X chromosome heterochromatin is strongly divergent among the studied species. In the light of the phylogenetic relationships, our data suggest that sex chromosome heteromorphism may have occurred early in the evolution of the Amazonian Engystomops and have been lost in two unnamed but confirmed candidate species.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

The Y Chromosome: A Complex Locus for Genetic Analyses of Complex Human Traits

The Human Y chromosome (ChrY) has been demonstrated to be a powerful tool for phylogenetics, population genetics, genetic genealogy and forensics. However, the importance of ChrY genetic variation in relation to human complex traits is less clear. In this review, we summarise existing evidence about the inherent complexities of ChrY variation and their use in association studies of human complex traits. We present and discuss the specific particularities of ChrY genetic variation, including Y chromosomal haplogroups, that need to be considered in the design and interpretation of genetic epidemiological studies involving ChrY.

High dynamism for neo-sex chromosomes: satellite DNAs reveal complex evolution in a grasshopper

A common characteristic of sex chromosomes is the accumulation of repetitive DNA, which accounts for their diversification and degeneration. In grasshoppers, the X0 sex-determining system in males is considered ancestral. However, in some species, derived variants like neo-XY in males evolved several times independently by Robertsonian translocation. This is the case of Ronderosia bergii, in which further large pericentromeric inversion in the neo-Y also took place, making this species particularly interesting for investigating sex chromosome evolution. Here, we characterized the satellite DNAs (satDNAs) and transposable elements (TEs) of the species to investigate the quantitative differences in repeat composition between male and female genomes putatively associated with sex chromosomes. We found a total of 53 satDNA families and 56 families of TEs. The satDNAs were 13.5% more abundant in males than in females, while TEs were just 1.02% more abundant in females. These results imply differential amplification of satDNAs on neo-Y chromosome and a minor role of TEs in sex chromosome differentiation. We showed highly differentiated neo-XY sex chromosomes owing to major amplification of satDNAs in neo-Y. Furthermore, chromosomal mapping of satDNAs suggests high turnover of neo-sex chromosomes in R. bergii at the intrapopulation level, caused by multiple paracentric inversions, amplifications, and transpositions. Finally, the species is an example of the action of repetitive DNAs in the generation of variability for sex chromosomes after the suppression of recombination, and helps understand sex chromosome evolution at the intrapopulation level.

Insight into the Possible Formation Mechanism of the Intersex Phenotype of Lanzhou Fat-Tailed Sheep Using Whole-Genome Resequencing

Simple Summary

Individuals with hermaphroditism are a serious hazard to animal husbandry and production due to their abnormal fertility. The molecular mechanism of sheep intersex formation was unclear. This study was the first to locate the homologous sequence of the goat polled intersex syndrome (PIS) region in sheep and found that the intersex traits of Lanzhou fat-tailed sheep were not caused by the lack of this region. By detecting the selective sweep regions, vital genes associated with androgen biosynthesis and the follicle stimulating hormone response entry were found, including steroid 5 alpha-reductase 2 (SRD5A2), and pro-apoptotic WT1 regulator (PAWR). Additionally, the copy number variations of the four regions on chr9, chr1, chr4, and chr16 may affect the expression of the gonadal development genes, zinc finger protein, FOG family member 2 (ZFPM2), LIM homeobox 8 (LHX8), inner mitochondrial membrane peptidase subunit 2 (IMMP2L) and slit guidance ligand 3 (SLIT3), respectively, and further affect the formation of intersex traits.

Abstract

Intersex, also known as hermaphroditism, is a serious hazard to animal husbandry and production. The mechanism of ovine intersex formation is not clear. Therefore, genome-wide resequencing on the only two intersex and two normal Lanzhou fat-tailed (LFT) sheep, an excellent but endangered Chinese indigenous sheep breed, was performed. Herein, the deletion of homologous sequences of the goat polled intersex syndrome (PIS) region (8787 bp, 247747059–247755846) on chromosome 1 of the LFT sheep was not the cause of the ovine intersex trait. By detecting the selective sweep regions, we found that the genes related to androgen biosynthesis and follicle stimulating hormone response items, such as steroid 5 alpha-reductase 2 (SRD5A2), steroid 5 alpha-reductase 3 (SRD5A3), and pro-apoptotic WT1 regulator (PAWR), may be involved in the formation of intersex traits. Furthermore, the copy number variations of the four regions, chr9: 71660801–71662800, chr1: 50776001–50778000, chr4: 58119201–58121600, and chr16: 778801–780800, may affect the expression of the zinc finger protein, FOG family member 2 (ZFPM2), LIM homeobox 8 (LHX8), inner mitochondrial membrane peptidase subunit 2 (IMMP2L) and slit guidance ligand 3 (SLIT3) genes, respectively, which contribute to the appearance of intersex traits. These results may supply a theoretical basis for the timely detection and elimination of intersex individuals in sheep, which could accelerate the healthy development of animal husbandry.

Whole genome-wide chromosome fusion and new gene birth in the Monopterus albus genome

Background

Teleost fishes account for over half of extant vertebrate species. A core question in biology is how genomic changes drive phenotypic diversity that relates to the origin of teleost fishes.

Results

Here, we used comparative genomic analyses with chromosome assemblies of diverse lineages of vertebrates and reconstructed an ancestral vertebrate genome, which revealed phylogenomic trajectories in vertebrates. We found that the whole-genome-wide chromosome fission/fusions took place in the Monopterus albus lineage after the 3-round whole-genome duplication. Four times of genomic fission/fusions events resulted in the whole genome-wide chromosome fusions in the genomic history of the lineage. In addition, abundant recently evolved new genes for reproduction emerged in the Monopterus albus after separated from medaka. Notably, we described evolutionary trajectories of conserved blocks related to sex determination genes in teleosts.

Conclusions

These data pave the way for a better understanding of genomic evolution in extant teleosts.

Sex differences in autophagy-mediated diseases: toward precision medicine

Nearly all diseases in humans, to a certain extent, exhibit sex differences, including differences in the onset, progression, prevention, therapy, and prognosis of diseases. Accumulating evidence shows that macroautophagy/autophagy, as a mechanism for development, differentiation, survival, and homeostasis, is involved in numerous aspects of sex differences in diseases such as cancer, neurodegeneration, and cardiovascular diseases. Advances in our knowledge regarding sex differences in autophagy-mediated diseases have enabled an understanding of their roles in human diseases, although the underlying molecular mechanisms of sex differences in autophagy remain largely unexplored. In this review, we discuss current advances in our insight into the biology of sex differences in autophagy and disease, information that will facilitate precision medicine.Abbreviations: AD: Azheimer disease; AMBRA1: autophagy and beclin 1 regulator 1; APP: amyloid beta precursor protein; AR: androgen receptor; AMPK: AMP-activated protein kinase; ATG: autophagy related; ATP6AP2: ATPase H+ transporting accessory protein 2; BCL2L1: BCL2 like 1; BECN1: beclin 1; CTSD: cathepsin D; CYP19A1: cytochrome P450 family 19 subfamily A member 1; DSD: disorders of sex development; eALDI: enhancer alternate long-distance initiator; ESR1: estrogen receptor 1; ESR2: estrogen receptor 2; FYCO1: FYVE and coiled-coil domain autophagy adaptor 1; GABARAP: GABA type A receptor-associated protein; GLA: galactosidase alpha; GTEx: genotype-tissue expression; HDAC6: histone deacetylase 6; I-R: ischemia-reperfusion; LAMP2: lysosomal associated membrane protein 2; MAP1LC3B/LC3B: microtubule associated protein 1 light chain 3 beta; MTOR: mechanistic target of rapamycin kinase; m6A: N6-methyladenosine; MYBL2: MYB proto-oncogene like 2; PIK3C3: phosphatidylinositol 3-kinase catalytic subunit type 3; PSEN1: presenilin 1; PSEN2: presenilin 2; RAB9A, RAB9A: member RAS oncogene family; RAB9B, RAB9B: member RAS oncogene family; RAB40AL: RAB40A like; SF1: splicing factor 1; SOX9: SRY-box transcription factor 9; SRY: sex determining region Y; TFEB: transcription factor EB; ULK1: unc-51 like autophagy activating kinase 1; UVRAG: UV radiation resistance associated; VDAC2: voltage dependent anion channel 2; WDR45: WD repeat domain 45; XPDS: X-linked parkinsonism and spasticity; YTHDF2: YTH N6-methyladenosine RNA binding protein 2.

Pathways to sex determination in plants: how many roads lead to Rome?

The presence of thousands of independent origins of dioecy in angiosperms provides a unique opportunity to address the parallel evolution of the molecular pathways underlying unisexual flowers. Recent progress towards identifying sex determination genes has identified hormone response pathways, mainly associated with cytokinin and ethylene response pathways, as having been recruited multiple times independently to control unisexuality. Moreover, transcriptomics has begun to identify commonalities among intermediate sections of signal transduction pathways. These recent advances set the stage for development of a comparative evolutionary development research program to identify the shared and unique aspects of the genetic pathways of unisexual flower development in angiosperms.

Heteromorphic Sex Chromosomes and Their DNA Content in Fish: An Insight through Satellite DNA Accumulation in Megaleporinus elongatus

The repetitive DNA content of fish sex chromosomes provides valuable insights into specificities and patterns of their genetic sex determination systems. In this study, we revealed the genomic satellite DNA (satDNA) content of Megaleporinuselongatus, a Neotropical fish species with Z1Z1Z2Z2/Z1W1Z2W2 multiple sex chromosomes, through high-throughput analysis and graph-based clustering, isolating 68 satDNA families. By physically mapping these sequences in female metaphases, we discovered 15 of the most abundant satDNAs clustered in its chromosomes, 9 of which were found exclusively in the highly heterochromatic W1. This heteromorphic sex chromosome showed the highest amount of satDNA accumulations in this species. The second most abundant family, MelSat02-26, shared FISH signals with the NOR-bearing pair in similar patterns and is linked to the multiple sex chromosome system. Our results demonstrate the diverse satDNA content in M. elongatus, especially in its heteromorphic sex chromosome. Additionally, we highlighted the different accumulation patterns and distribution of these sequences across species by physically mapping these satDNAs in other Anostomidae, Megaleporinusmacrocephalus and Leporinusfriderici (a species without differentiated sex chromosomes).

First Report of Sex Chromosomes in Night Lizards (Scincoidea: Xantusiidae)

Squamate reptiles (lizards, snakes, and amphibians) are an outstanding group for studying sex chromosome evolution—they are old, speciose, geographically widespread, and exhibit myriad sex-determining modes. Yet, the vast majority of squamate species lack heteromorphic sex chromosomes. Cataloging the sex chromosome systems of species lacking easily identifiable, heteromorphic sex chromosomes, therefore, is essential before we are to fully understand the evolution of vertebrate sex chromosomes. Here, we use restriction site-associated DNA sequencing (RADseq) to classify the sex chromosome system of the granite night lizard, Xantusia henshawi. RADseq is an effective alternative to traditional cytogenetic methods for determining a species’ sex chromosome system (i.e., XX/XY or ZZ/ZW), particularly in taxa with non-differentiated sex chromosomes. Although many xantusiid lineages have been karyotyped, none possess heteromorphic sex chromosomes. We identified a ZZ/ZW sex chromosome system in X. henshawi—the first such data for this family. Furthermore, we report that the X. henshawi sex chromosome contains fragments of genes found on Gallus gallus chromosomes 7, 12, and 18 (which are homologous to Anolis carolinensis chromosome 2), the first vertebrate sex chromosomes to utilize this linkage group.

Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data

Background

Mammalian X and Y chromosomes share a common evolutionary origin and retain regions of high sequence similarity. Similar sequence content can confound the mapping of short next-generation sequencing reads to a reference genome. It is therefore possible that the presence of both sex chromosomes in a reference genome can cause technical artifacts in genomic data and affect downstream analyses and applications. Understanding this problem is critical for medical genomics and population genomic inference.

Results

Here, we characterize how sequence homology can affect analyses on the sex chromosomes and present XYalign, a new tool that (1) facilitates the inference of sex chromosome complement from next-generation sequencing data; (2) corrects erroneous read mapping on the sex chromosomes; and (3) tabulates and visualizes important metrics for quality control such as mapping quality, sequencing depth, and allele balance. We find that sequence homology affects read mapping on the sex chromosomes and this has downstream effects on variant calling. However, we show that XYalign can correct mismapping, resulting in more accurate variant calling. We also show how metrics output by XYalign can be used to identify XX and XY individuals across diverse sequencing experiments, including low- and high-coverage whole-genome sequencing, and exome sequencing. Finally, we discuss how the flexibility of the XYalign framework can be leveraged for other uses including the identification of aneuploidy on the autosomes. XYalign is available open source under the GNU General Public License (version 3).

Conclusions

Sex chromsome sequence homology causes the mismapping of short reads, which in turn affects downstream analyses. XYalign provides a reproducible framework to correct mismapping and improve variant calling on the sex chromsomes.

Transcriptome of Gonads From High Temperature Induced Sex Reversal During Sex Determination and Differentiation in Chinese Tongue Sole, Cynoglossus semilaevis

The sex of Chinese tongue sole (Cynoglossus semilaevis) is determined by both genetic sex determination (GSD) and environmental sex determination (ESD), making it an ideal model to study the relationship between sex-determination and temperature. In the present study, transcriptomes of undifferentiated gonads from genetic females and males, as well as differentiated gonads from males, females, and pseudomales under high and normal temperature treatments were generated for comparative transcriptomic analysis. A mean of 68.24 M high-quality clean reads was obtained for each library. Differentially expressed genes (DEGs) between different sexes and environmental treatments were identified, revealing that the heat shock protein gene family was involved in the high temperature induced sex reversal. The Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were enriched in pseudomale and genetic female comparison included neuroactive ligand-receptor interaction, cortisol synthesis and secretion, and steroid hormone biosynthesis. Furthermore, weighted gene co-expression network analyses were conducted on all samples, and two modules were positive correlated with pseudomale under high temperature. An illustrated protein-protein interaction map of the module identified a hub gene, hsc70. These findings provide insights into the genetic network that is involved in sex determination and sexual differentiation, and improve our understanding of genes involved in sex reversal under high temperature.

Developmental Systems Drift and the Drivers of Sex Chromosome Evolution

Phenotypic invariance—the outcome of purifying selection—is a hallmark of biological importance. However, invariant phenotypes might be controlled by diverged genetic systems in different species. Here, we explore how an important and invariant phenotype—the development of sexually differentiated individuals—is controlled in over two dozen species in the frog family Pipidae. We uncovered evidence in different species for 1) an ancestral W chromosome that is not found in many females and is found in some males, 2) independent losses and 3) autosomal segregation of this W chromosome, 4) changes in male versus female heterogamy, and 5) substantial variation among species in recombination suppression on sex chromosomes. We further provide evidence of, and evolutionary context for, the origins of at least seven distinct systems for regulating sex determination among three closely related genera. These systems are distinct in their genomic locations, evolutionary origins, and/or male versus female heterogamy. Our findings demonstrate that the developmental control of sexual differentiation changed via loss, sidelining, and empowerment of a mechanistically influential gene, and offer insights into novel factors that impinge on the diverse evolutionary fates of sex chromosomes.

Novel Sex Chromosomes in 3 Cichlid Fishes from Lake Tanganyika

African cichlids are well known for their adaptive radiations, but it is now apparent that they also harbor an extraordinary diversity of sex chromosome systems. In this study, we sequenced pools of males and females from species in 3 different genera of cichlids from Lake Tanganyika. We then searched for regions that were differentiated following the patterns expected for sex chromosomes. We report 2 novel sex chromosomes systems, an XY system on LG19 in Tropheus sp. "black" and a ZW system on LG7 in Hemibates stenosoma. We also identify a ZW system on LG5 in Cyprichromis leptosoma that may be convergent with a system previously described in Lake Malawi cichlids. Our data also identify candidate single nucleotide polymorphisms for the blue/yellow tail color polymorphism observed among male C. leptosoma.