Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage

SMC-mediated dosage compensation in C. elegans evolved in the presence of an ancestral nematode mechanism

Mechanisms of X chromosome dosage compensation have been studied extensively in a few model species representing clades of shared sex chromosome ancestry. However, the diversity within each clade as a function of sex chromosome evolution is largely unknown. Here, we anchor ourselves to the nematode Caenorhabditis elegans, for which a well-studied mechanism of dosage compensation occurs through a specialized structural maintenance of chromosomes (SMC) complex, and explore the diversity of dosage compensation in the surrounding phylogeny of nematodes. Through phylogenetic analysis of the C. elegans dosage compensation complex and a survey of its epigenetic signatures, including X-specific topologically associating domains (TADs) and X-enrichment of H4K20me1, we found that the condensin-mediated mechanism evolved recently in the lineage leading to Caenorhabditis through an SMC-4 duplication. Intriguingly, an independent duplication of SMC-4 and the presence of X-specific TADs in Pristionchus pacificus suggest that condensin-mediated dosage compensation arose more than once. mRNA-seq analyses of gene expression in several nematode species indicate that dosage compensation itself is ancestral, as expected from the ancient XO sex determination system. Indicative of the ancestral mechanism, H4K20me1 is enriched on the X chromosomes in Oscheius tipulae, which does not contain X-specific TADs or SMC-4 paralogs. Together, our results indicate that the dosage compensation system in C. elegans is surprisingly new, and condensin may have been co-opted repeatedly in nematodes, suggesting that the process of evolving a chromosome-wide gene regulatory mechanism for dosage compensation is constrained.

Significance statement

X chromosome dosage compensation mechanisms evolved in response to Y chromosome degeneration during sex chromosome evolution. However, establishment of dosage compensation is not an endpoint. As sex chromosomes change, dosage compensation strategies may have also changed. In this study, we performed phylogenetic and epigenomic analyses surrounding Caenorhabditis elegans and found that the condensin-mediated dosage compensation mechanism in C. elegans is surprisingly new, and has evolved in the presence of an ancestral mechanism. Intriguingly, condensin-based dosage compensation may have evolved more than once in the nematode lineage, the other time in Pristionchus. Together, our work highlights a previously unappreciated diversity of dosage compensation mechanisms within a clade, and suggests constraints in evolving new mechanisms in the presence of an existing one.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Chromosome-Level Genome Assembly of the Blacktail Brush Lizard, Urosaurus nigricaudus, Reveals Dosage Compensation in an Endemic Lizard

Urosaurus nigricaudus is a phrynosomatid lizard endemic to the Baja California Peninsula in Mexico. This work presents a chromosome-level genome assembly and annotation from a male individual. We used PacBio long reads and HiRise scaffolding to generate a high-quality genomic assembly of 1.87 Gb distributed in 327 scaffolds, with an N50 of 279 Mb and an L50 of 3. Approximately 98.4% of the genome is contained in 14 scaffolds, with 6 large scaffolds (334-127 Mb) representing macrochromosomes and 8 small scaffolds (63-22 Mb) representing microchromosomes. Using standard gene modeling and transcriptomic data, we predicted 17,902 protein-coding genes on the genome. The repeat content is characterized by a large proportion of long interspersed nuclear elements that are relatively old. Synteny analysis revealed some microchromosomes with high repeat content are more prone to rearrangements but that both macro- and microchromosomes are well conserved across reptiles. We identified scaffold 14 as the X chromosome. This microchromosome presents perfect dosage compensation where the single X of males has the same expression levels as two X chromosomes in females. Finally, we estimated the effective population size for U. nigricaudus was extremely low, which may reflect a reduction in polymorphism related to it becoming a peninsular endemic.

Transposon wave remodeled the epigenomic landscape in the rapid evolution of X-Chromosome dosage compensation

Sex chromosome dosage compensation is a model to understand the coordinated evolution of transcription; however, the advanced age of the sex chromosomes in model systems makes it difficult to study how the complex regulatory mechanisms underlying chromosome-wide dosage compensation can evolve. The sex chromosomes of Poecilia picta have undergone recent and rapid divergence, resulting in widespread gene loss on the male Y, coupled with complete X Chromosome dosage compensation, the first case reported in a fish. The recent de novo origin of dosage compensation presents a unique opportunity to understand the genetic and evolutionary basis of coordinated chromosomal gene regulation. By combining a new chromosome-level assembly of P. picta with whole-genome bisulfite sequencing and RNA-seq data, we determine that the YY1 transcription factor (YY1) DNA binding motif is associated with male-specific hypomethylated regions on the X, but not the autosomes. These YY1 motifs are the result of a recent and rapid repetitive element expansion on the P. picta X Chromosome, which is absent in closely related species that lack dosage compensation. Taken together, our results present compelling support that a disruptive wave of repetitive element insertions carrying YY1 motifs resulted in the remodeling of the X Chromosome epigenomic landscape and the rapid de novo origin of a dosage compensation system.

The molecular evolution of mammalian spermatogenesis

The testis is a key male reproductive organ that produces gametes through the process of spermatogenesis. Testis morphologies, sperm phenotypes, and the process of spermatogenesis evolve rapidly in mammals, presumably due to the evolutionary pressure on males to give rise to their own offspring. Here, we review studies illuminating the molecular evolution of the testis, in particular large-scale transcriptomic studies, which were based on bulk tissue samples and, more recently, individual cells. Together with various genomic and epigenomic data, these studies have unveiled the cellular source, molecular mechanisms, and evolutionary forces that underlie the rapid phenotypic evolution of the testis. They also revealed shared (ancestral) and species-specific spermatogenic gene expression programs. The insights and available data that have accumulated also provide a valuable resource for the investigation and treatment of male fertility disorders - a dramatically increasing problem in modern industrial societies.

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.

Genome-wide analysis of RNA-chromatin interactions in lizards as a mean for functional lncRNA identification

BACKGROUND

Long non-coding RNAs (lncRNAs) are defined as transcribed molecules longer than 200 nucleotides with little to no protein-coding potential. LncRNAs can regulate gene expression of nearby genes (cis-acting) or genes located on other chromosomes (trans-acting). Several methodologies have been developed to capture lncRNAs associated with chromatin at a genome-wide level. Analysis of RNA-DNA contacts can be combined with epigenetic and RNA-seq data to define potential lncRNAs involved in the regulation of gene expression.

RESULTS

We performed Chromatin Associated RNA sequencing (ChAR-seq) in Anolis carolinensis to obtain the genome-wide map of the associations that RNA molecules have with chromatin. We analyzed the frequency of DNA contacts for different classes of RNAs and were able to define cis- and trans-acting lncRNAs. We integrated the ChAR-seq map of RNA-DNA contacts with epigenetic data for the acetylation of lysine 16 on histone H4 (H4K16ac), a mark connected to actively transcribed chromatin in lizards. We successfully identified three trans-acting lncRNAs significantly associated with the H4K16ac signal, which are likely involved in the regulation of gene expression in A. carolinensis.

CONCLUSIONS

We show that the ChAR-seq method is a powerful tool to explore the RNA-DNA map of interactions. Moreover, in combination with epigenetic data, ChAR-seq can be applied in non-model species to establish potential roles for predicted lncRNAs that lack functional annotations.

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

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals^1-6, probably owing to the evolutionary pressure on males to be reproductively successful⁷. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals.

Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei – a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

A highly-complete chromosome-scale genome assembly of the brown anole, Anolis sagrei, provides insight into the evolution of sex chromosomes and is a crucial resource for this model lizard species.

The HASTER lncRNA promoter is a cis-acting transcriptional stabilizer of HNF1A

The biological purpose of long non-coding RNAs (lncRNAs) is poorly understood. Haploinsufficient mutations in HNF1A homeobox A (HNF1A), encoding a homeodomain transcription factor, cause diabetes mellitus. Here, we examine HASTER, the promoter of an lncRNA antisense to HNF1A. Using mouse and human models, we show that HASTER maintains cell-specific physiological HNF1A concentrations through positive and negative feedback loops. Pancreatic β cells from Haster mutant mice consequently showed variegated HNF1A silencing or overexpression, resulting in hyperglycaemia. HASTER-dependent negative feedback was essential to prevent HNF1A binding to inappropriate genomic regions. We demonstrate that the HASTER promoter DNA, rather than the lncRNA, modulates HNF1A promoter-enhancer interactions in cis and thereby regulates HNF1A transcription. Our studies expose a cis-regulatory element that is unlike classic enhancers or silencers, it stabilizes the transcription of its target gene and ensures the fidelity of a cell-specific transcription factor program. They also show that disruption of a mammalian lncRNA promoter can cause diabetes mellitus.

Monotreme-Specific Conserved Putative Proteins Derived from Retroviral Reverse Transcriptase

Endogenous retroviruses (ERVs) have played an essential role in the evolution of mammals. ERV-derived genes are reported in the therians, many of which are involved in placental development; however, the contribution of the ERV-derived genes in monotremes, which are oviparous mammals, remains to be uncovered. Here, we conducted a comprehensive search for possible ERV-derived genes in platypus and echidna genomes and identified three reverse transcriptase-like genes named RTOM1, RTOM2, and RTOM3 clustered in the GRIP2 intron. Comparative genomic analyses revealed that RTOM1, RTOM2, and RTOM3 are strongly conserved and are under purifying selection between these species. These could be generated by tandem duplications before the divergence of platypus and echidna. All RTOM transcripts were specifically expressed in the testis, possibly suggesting their physiological importance. This is the first study reporting monotreme-specific de novo gene candidates derived from ERVs, which provides new insights into the unique evolution of monotremes.

Genetic determination and JARID2 over-expression in a thermal incubation experiment in Casque-Headed Lizard

Non-avian reptiles, unlike mammals and birds, have undergone numerous sex determination changes. Casque-Headed Lizards have replaced the ancestral XY system shared across pleurodonts with a new pair of XY chromosomes. However, the evolutionary forces that triggered this transition have remained unclear. An interesting hypothesis suggests that species with intermediate states, with sex chromosomes but also thermal-induced sex reversal at specific incubation temperatures, could be more susceptible to sex determination turnovers. We contrasted genotypic data (presence/absence of the Y chromosome) against the histology of gonads of embryos from stages 35–37 incubated at various temperatures, including typical male-producing (26°C) and female-producing (32°C) temperatures. Our work apparently reports for the first time the histology of gonads, including morphological changes, from stages 35–37 of development in the family Corytophanidae. We also observed that all embryos developed hemipenes, suggesting sex-linked developmental heterochrony. We observed perfect concordance between genotype and phenotype at all temperatures. However, analysis of transcriptomic data from embryos incubated at 26°C and 32°C identified transcript variants of the chromatin modifiers JARID2 and KDM6B that have been linked to temperature-dependent sex determination in other reptiles. Our work tested the validity of a mixed sex determination system in the family Corytophanidae. We found that XY chromosomes are dominant; however, our work supports the hypothesis of a conserved transcriptional response to incubation temperatures across non-avian reptiles that could be a reminiscence of an ancestral sex determination system.

Evolution of the short form of DNMT3A, DNMT3A2, occurred in the common ancestor of mammals

Genomic imprinting is found in marsupial and eutherian mammals, but not in monotremes. While the primary regulator of genomic imprinting in eutherians is differential DNA methylation between parental alleles, conserved imprinted genes in marsupials tend to lack DNA methylation at their promoters. DNA methylation at eutherian imprinted genes is mainly catalysed by a DNA methyltransferase (DNMT) enzyme, DNMT3A. There are two isoforms of eutherian DNMT3A: DNMT3A and DNMT3A2. DNMT3A2 is the primary isoform for establishing DNA methylation at eutherian imprinted genes and is essential for eutherian genomic imprinting. In this study, we investigated whether DNMT3A2 is also present in the two other mammalian lineages, marsupials and monotremes. We identified DNMT3A2 in both marsupials and monotremes, although imprinting has not been identified in monotremes. By analysing genomic sequences and transcriptome data across vertebrates, we concluded that the evolution of DNMT3A2 occurred in the common ancestor of mammals. In addition, DNMT3A/3A2 gene and protein expression during gametogenesis showed distinct sexual dimorphisms in a marsupial, the tammar wallaby, and this pattern coincided with the sex-specific DNA methylation reprogramming in this species as it does in mice. Our results show that DNMT3A2 is present in all mammalian groups and suggests that the basic DNMT3A/3A2-based DNA methylation mechanism is conserved at least in therian mammals.

Dosage compensation evolution in plants: theories, controversies and mechanisms

In a minority of flowering plants, separate sexes are genetically determined by sex chromosomes. The Y chromosome has a non-recombining region that degenerates, causing a reduced expression of Y genes. In some species, the lower Y expression is accompanied by dosage compensation (DC), a mechanism that re-equalizes male and female expression and/or brings XY male expression back to its ancestral level. Here, we review work on DC in plants, which started as early as the late 1960s with cytological approaches. The use of transcriptomics fired a controversy as to whether DC existed in plants. Further work revealed that various plants exhibit partial DC, including a few species with young and homomorphic sex chromosomes. We are starting to understand the mechanisms responsible for DC in some plants, but in most species, we lack the data to differentiate between global and gene-by-gene DC. Also, it is unknown why some species evolve many dosage compensated genes while others do not. Finally, the forces that drive DC evolution remain mysterious, both in plants and animals. We review the multiple evolutionary theories that have been proposed to explain DC patterns in eukaryotes with XY or ZW sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Lack of age-related mosaic loss of W chromosome in long-lived birds

Females and males often exhibit different survival in nature, and it has been hypothesized that sex chromosomes may play a role in driving differential survival rates. For instance, the Y chromosome in mammals and the W chromosome in birds are often degenerated, with reduced numbers of genes, and loss of the Y chromosome in old men is associated with shorter life expectancy. However, mosaic loss of sex chromosomes has not been investigated in any non-human species. Here, we tested whether mosaic loss of the W chromosome (LOW) occurs with ageing in wild birds as a natural consequence of cellular senescence. Using loci-specific PCR and a target sequencing approach we estimated LOW in both young and adult individuals of two long-lived bird species and showed that the copy number of W chromosomes remains constant across age groups. Our results suggest that LOW is not a consequence of cellular ageing in birds. We concluded that the inheritance of the W chromosome in birds, unlike the Y chromosome in mammals, is more stable.

Evolutionary Characterization of the Short Protein SPAAR

Microproteins (<100 amino acids) are receiving increasing recognition as important participants in numerous biological processes, but their evolutionary dynamics are poorly understood. SPAAR is a recently discovered microprotein that regulates muscle regeneration and angiogenesis through interactions with conserved signaling pathways. Interestingly, SPAAR does not belong to any known protein family and has known homologs exclusively among placental mammals. This lack of distant homology could be caused by challenges in homology detection of short sequences, or it could indicate a recent de novo emergence from a noncoding sequence. By integrating syntenic alignments and homology searches, we identify SPAAR orthologs in marsupials and monotremes, establishing that SPAAR has existed at least since the emergence of mammals. SPAAR shows substantial primary sequence divergence but retains a conserved protein structure. In primates, we infer two independent evolutionary events leading to the de novo origination of 5' elongated isoforms of SPAAR from a noncoding sequence and find evidence of adaptive evolution in this extended region. Thus, SPAAR may be of ancient origin, but it appears to be experiencing continual evolutionary innovation in mammals.

Sex Chromosomes and Master Sex-Determining Genes in Turtles and Other Reptiles

Among tetrapods, the well differentiated heteromorphic sex chromosomes of birds and mammals have been highly investigated and their master sex-determining (MSD) gene, Dmrt1 and SRY, respectively, have been identified. The homomorphic sex chromosomes of reptiles have been the least studied, but the gap with birds and mammals has begun to fill. This review describes our current knowledge of reptilian sex chromosomes at the cytogenetic and molecular level. Most of it arose recently from various studies comparing male to female gene content. This includes restriction site-associated DNA sequencing (RAD-Seq) experiments in several male and female samples, RNA sequencing and identification of Z- or X-linked genes by male/female comparative transcriptome coverage, and male/female transcriptomic or transcriptome/genome substraction approaches allowing the identification of Y- or W-linked transcripts. A few putative master sex-determining (MSD) genes have been proposed, but none has been demonstrated yet. Lastly, future directions in the field of reptilian sex chromosomes and their MSD gene studies are considered.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. 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 I)'.

Thermosensitive sex chromosome dosage compensation in ZZ/ZW softshell turtles, Apalone spinifera

Sex chromosome dosage compensation (SCDC) overcomes gene-dose imbalances that disturb transcriptional networks, as when ZW females or XY males are hemizygous for Z/X genes. Mounting data from non-model organisms reveal diverse SCDC mechanisms, yet their evolution remains obscure, because most informative lineages with variable sex chromosomes are unstudied. Here, we discovered SCDC in turtles and an unprecedented thermosensitive SCDC in eukaryotes. We contrasted RNA-seq expression of Z-genes, their autosomal orthologues, and control autosomal genes in Apalone spinifera (ZZ/ZW) and Chrysemys picta turtles with temperature-dependent sex determination (TSD) (proxy for ancestral expression). This approach disentangled chromosomal context effects on Z-linked and autosomal expression, from lineage effects owing to selection or drift. Embryonic Apalone SCDC is tissue- and age-dependent, regulated gene-by-gene, complete in females via Z-upregulation in both sexes (Type IV) but partial and environmentally plastic via Z-downregulation in males (accentuated at colder temperature), present in female hatchlings and a weakly suggestive in adult liver (Type I). Results indicate that embryonic SCDC evolved with/after sex chromosomes in Apalone's family Tryonichidae, while co-opting Z-gene upregulation present in the TSD ancestor. Notably, Apalone's SCDC resembles pygmy snake's, and differs from the full-SCDC of Anolis lizards who share homologous sex chromosomes (XY), advancing our understanding of how XX/XY and ZZ/ZW systems compensate gene-dose imbalance. 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)'.

Whole-chromosome fusions in the karyotype evolution of Sceloporus (Iguania, Reptilia) are more frequent in sex chromosomes than autosomes

Whole-chromosome fusions play a major role in the karyotypic evolution of reptiles. It has been suggested that certain chromosomes tend to fuse with sex chromosomes more frequently than others. However, the comparative genomic synteny data are too scarce to draw strong conclusions. We obtained and sequenced chromosome-specific DNA pools of Sceloporus malachiticus, an iguanian species which has experienced many chromosome fusions. We found that four of seven lineage-specific fusions involved sex chromosomes, and that certain syntenic blocks which constitute the sex chromosomes, such as the homologues of the Anolis carolinensis chromosomes 11 and 16, are repeatedly involved in sex chromosome formation in different squamate species. To test the hypothesis that the karyotypic shift could be associated with changes in recombination patterns, we performed a synaptonemal complex analysis in this species and in Sceloporus variabilis (2n = 34). It revealed that the sex chromosomes in S. malachiticus had two distal pseudoautosomal regions and a medial differentiated region. We found that multiple fusions little affected the recombination rate in S. malachiticus. Our data confirm more frequent involvement of certain chromosomes in sex chromosome formation, but do not reveal a connection between the gonosome-autosome fusions and the evolution of recombination rate. 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)'.

Distinct mechanisms mediate X chromosome dosage compensation in Anopheles and Drosophila

CRISPR knockout of msl-2 and epigenome analyses in Anopheles reveal that X chromosome dosage compensation in mosquitos and Drosophila is achieved by two different molecular mechanisms.

Sex chromosomes induce potentially deleterious gene expression imbalances that are frequently corrected by dosage compensation (DC). Three distinct molecular strategies to achieve DC have been previously described in nematodes, fruit flies, and mammals. Is this a consequence of distinct genomes, functional or ecological constraints, or random initial commitment to an evolutionary trajectory? Here, we study DC in the malaria mosquito Anopheles gambiae. The Anopheles and Drosophila X chromosomes evolved independently but share a high degree of homology. We find that Anopheles achieves DC by a mechanism distinct from the Drosophila MSL complex–histone H4 lysine 16 acetylation pathway. CRISPR knockout of Anopheles msl-2 leads to embryonic lethality in both sexes. Transcriptome analyses indicate that this phenotype is not a consequence of defective X chromosome DC. By immunofluorescence and ChIP, H4K16ac does not preferentially enrich on the male X. Instead, the mosquito MSL pathway regulates conserved developmental genes. We conclude that a novel mechanism confers X chromosome up-regulation in Anopheles. Our findings highlight the pluralism of gene-dosage buffering mechanisms even under similar genomic and functional constraints.

Viviparous Reptile Regarded to Have Temperature-Dependent Sex Determination Has Old XY Chromosomes

The water skinks Eulamprus tympanum and Eulamprus heatwolei show thermally induced sex determination where elevated temperatures give rise to male offspring. Paradoxically, Eulamprus species reproduce in temperatures of 12–15 °C making them outliers when compared with reptiles that use temperature as a cue for sex determination. Moreover, these two species are among the very few viviparous reptiles reported to have thermally induced sex determination. Thus, we tested whether these skinks possess undetected sex chromosomes with thermal override. We produced transcriptome and genome data for E. heatwolei. We found that E. heatwolei presents XY chromosomes that include 14 gametologs with regulatory functions. The Y chromosomal region is 79–116 Myr old and shared between water and spotted skinks. Our work provides clear evidence that climate could be useful to predict the type of sex determination systems in reptiles and it also indicates that viviparity is strictly associated with sex chromosomes.

Rapid Evolution of Complete Dosage Compensation in Poecilia

Dosage compensation balances gene expression between the sexes in systems with diverged heterogametic sex chromosomes. Theory predicts that dosage compensation should rapidly evolve in tandem with the divergence of sex chromosomes to prevent the deleterious effects of dosage imbalances that occur as a result of sex chromosome divergence. Examples of complete dosage compensation, where gene expression of the entire sex chromosome is compensated, are rare, and have only been found in relatively ancient sex chromosome systems. Consequently, very little is known about the evolutionary dynamics of complete dosage compensation systems. Within the family Poeciliidae the subgenus Lebistes share the same sex chromosome system which originated 18.48–26.08 Ma. In Poecilia reticulata and P. wingei, the Y chromosome has been largely maintained, whereas the Y in the closely related species P. picta and P. parae has rapidly degraded. We recently found P. picta to be the first example of complete dosage compensation in a fish. Here, we show that P. parae also has complete dosage compensation, thus complete dosage compensation likely evolved in the short (∼3.7 Myr) interval after the split of the ancestor of these two species from P. reticulata, but before they diverged from each other. These data suggest that novel dosage compensation mechanisms can evolve rapidly, thus supporting the longstanding theoretical prediction that such mechanisms arise in tandem with rapidly diverging sex chromosomes.

Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae)

Heteromorphic W and Y sex chromosomes often experience gene loss and heterochromatinization, which is frequently viewed as their "degeneration". However, the evolutionary trajectories of the heterochromosomes are in fact more complex since they may not only lose but also acquire new sequences. Previously, we found that the heterochromatic W chromosome of a lizard Eremias velox (Lacertidae) is decondensed and thus transcriptionally active during the lampbrush stage. To determine possible sources of this transcription, we sequenced DNA from a microdissected W chromosome sample and a total female DNA sample and analyzed the results of reference-based and de novo assembly. We found a new repetitive sequence, consisting of fragments of an autosomal protein-coding gene ATF7IP2, several SINE elements, and sequences of unknown origin. This repetitive element is distributed across the whole length of the W chromosome, except the centromeric region. Since it retained only 3 out of 10 original ATF7IP2 exons, it remains unclear whether it is able to produce a protein product. Subsequent studies are required to test the presence of this element in other species of Lacertidae and possible functionality. Our results provide further evidence for the view of W and Y chromosomes as not just "degraded" copies of Z and X chromosomes but independent genomic segments in which novel genetic elements may arise.

Epigenetics drive the evolution of sex chromosomes in animals and plants

We review how epigenetics affect sex chromosome evolution in animals and plants. In a few species, sex is determined epigenetically through the action of Y-encoded small RNAs. Epigenetics is also responsible for changing the sex of individuals through time, even in species that carry sex chromosomes, and could favour species adaptation through breeding system plasticity. The Y chromosome accumulates repeats that become epigenetically silenced which leads to an epigenetic conflict with the expression of Y genes and could accelerate Y degeneration. Y heterochromatin can be lost through ageing, which activates transposable elements and lowers male longevity. Y chromosome degeneration has led to the evolution of meiotic sex chromosome inactivation in eutherians (placentals) and marsupials, and dosage compensation mechanisms in animals and plants. X-inactivation convergently evolved in eutherians and marsupials via two independently evolved non-coding RNAs. In Drosophila, male X upregulation by the male specific lethal (MSL) complex can spread to neo-X chromosomes through the transposition of transposable elements that carry an MSL-binding motif. We discuss similarities and possible differences between plants and animals and suggest future directions for this dynamic field of research. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Evolution of dosage compensation does not depend on genomic background

Organisms have evolved various mechanisms to cope with the differences in the gene copy numbers between sexes caused by degeneration of Y and W sex chromosomes. Complete dosage compensation or at least expression balance between sexes has been reported predominantly in XX/XY systems, but rarely in ZZ/ZW systems. However, this often-reported pattern is based on comparisons of lineages where sex chromosomes evolved from nonhomologous genomic regions, potentially differing in sensitivity to differences in gene copy numbers. Here we document that two reptilian lineages (XX/XY iguanas and ZZ/ZW softshell turtles), which independently co-opted the same ancestral genomic region for the function of sex chromosomes, evolved different gene dose regulatory mechanisms. The independent co-option of the same genomic region for the role of sex chromosomes as in the iguanas and the softshell turtles offers great opportunity for testing evolutionary scenarios on sex chromosome evolution under the explicit control of the genomic background and gene identity. We show that the parallel loss of functional genes from the Y chromosome of the green anole and the W chromosome of the Florida softshell turtle led to different dosage compensation mechanisms. Our approach controlling for genetic background thus does not support that the variability in the regulation of gene dose differences is a consequence of ancestral autosomal gene content.

Poorly differentiated XX/XY sex chromosomes are widely shared across skink radiation

Differentiated sex chromosomes are believed to be evolutionarily stable, while poorly differentiated sex chromosomes are considered to be prone to turnovers. With around 1700 currently known species forming ca 15% of reptile species diversity, skinks (family Scincidae) are a very diverse group of squamates known for their large ecological and morphological variability. Skinks generally have poorly differentiated and cytogenetically indistinguishable sex chromosomes, and their sex determination was suggested to be highly variable. Here, we determined X-linked genes in the common sandfish (Scincus scincus) and demonstrate that skinks have shared the same homologous XX/XY sex chromosomes across their wide phylogenetic spectrum for at least 85 million years, approaching the age of the highly differentiated ZZ/ZW sex chromosomes of birds and advanced snakes. Skinks thus demonstrate that even poorly differentiated sex chromosomes can be evolutionarily stable. The conservation of sex chromosomes across skinks allows us to introduce the first molecular sexing method widely applicable in this group.

MeDAS: a Metazoan Developmental Alternative Splicing database

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Ion channel profiling of the Lymnaea stagnalis ganglia via transcriptome analysis

BACKGROUND

The pond snail Lymnaea stagnalis (L. stagnalis) has been widely used as a model organism in neurobiology, ecotoxicology, and parasitology due to the relative simplicity of its central nervous system (CNS). However, its usefulness is restricted by a limited availability of transcriptome data. While sequence information for the L. stagnalis CNS transcripts has been obtained from EST libraries and a de novo RNA-seq assembly, the quality of these assemblies is limited by a combination of low coverage of EST libraries, the fragmented nature of de novo assemblies, and lack of reference genome.

RESULTS

In this study, taking advantage of the recent availability of a preliminary L. stagnalis genome, we generated an RNA-seq library from the adult L. stagnalis CNS, using a combination of genome-guided and de novo assembly programs to identify 17,832 protein-coding L. stagnalis transcripts. We combined our library with existing resources to produce a transcript set with greater sequence length, completeness, and diversity than previously available ones. Using our assembly and functional domain analysis, we profiled L. stagnalis CNS transcripts encoding ion channels and ionotropic receptors, which are key proteins for CNS function, and compared their sequences to other vertebrate and invertebrate model organisms. Interestingly, L. stagnalis transcripts encoding numerous putative Ca2+ channels showed the most sequence similarity to those of Mus musculus, Danio rerio, Xenopus tropicalis, Drosophila melanogaster, and Caenorhabditis elegans, suggesting that many calcium channel-related signaling pathways may be evolutionarily conserved.

CONCLUSIONS

Our study provides the most thorough characterization to date of the L. stagnalis transcriptome and provides insights into differences between vertebrates and invertebrates in CNS transcript diversity, according to function and protein class. Furthermore, this study provides a complete characterization of the ion channels of Lymnaea stagnalis, opening new avenues for future research on fundamental neurobiological processes in this model system.

Meiotic Cells Counteract Programmed Retrotransposon Activation via RNA-Binding Translational Repressor Assemblies

Retrotransposon proliferation poses a threat to germline integrity. While retrotransposons must be activated in developing germ cells in order to survive and propagate, how they are selectively activated in the context of meiosis is unclear. We demonstrate that the transcriptional activation of Ty3/Gypsy retrotransposons and host defense are controlled by master meiotic regulators. We show that budding yeast Ty3/Gypsy co-opts binding sites of the essential meiotic transcription factor Ndt80 upstream of the integration site, thereby tightly linking its transcriptional activation to meiotic progression. We also elucidate how yeast cells thwart Ty3/Gypsy proliferation by blocking translation of the retrotransposon mRNA using amyloid-like assemblies of the RNA-binding protein Rim4. In mammals, several inactive Ty3/Gypsy elements are undergoing domestication. We show that mammals utilize equivalent master meiotic regulators (Stra8, Mybl1, Dazl) to regulate Ty3/Gypsy-derived genes in developing gametes. Our findings inform how genes that are evolving from retrotransposons can build upon existing regulatory networks during domestication.

The evolution of sex chromosome dosage compensation in animals

The evolution of heteromorphic sex chromosomes shall lead to gene expression dosage problems, as in at least one of the sexes, the sex-linked gene dose has been reduced by half. It has been proposed that the transcriptional output of the whole X or Z chromosome should be doubled for complete dosage compensation in heterogametic sex. However, owing to the variability of the existing methods to determine the transcriptional differences between sex chromosomes and autosomes (S:A ratios) in different studies, we collected more than 500 public RNA-Seq data set from multiple tissues and species in major clades and proposed a unified computational framework for unbiased and comparable measurement of the S:A ratios of multiple species. We also tested the evolution of dosage compensation more directly by assessing changes in the expression levels of the current sex-linked genes relative to those of the ancestral sex-linked genes. We found that in mammals and birds, the S:A ratio is approximately 0.5, whereas in insects, fishes, and flatworms, the S:A ratio is approximately 1.0. Further analysis showed that the fraction of dosage-sensitive housekeeping genes on the X/Z chromosome is significantly correlated with the S:A ratio. In addition, the degree of degeneration of the Y chromosome may be responsible for the change in the S:A ratio in mammals without a dosage compensation mechanism. Our observations offer unequivocal support for the sex chromosome insensitivity hypothesis in animals and suggest that dosage sensitivity states of sex chromosomes are a major factor underlying different evolutionary strategies of dosage compensation.

Transcriptome and translatome co-evolution in mammals

Gene expression programs define shared and species-specific phenotypes, but their evolution remains largely uncharacterized beyond the transcriptome layer¹. Here we report an analysis of the co-evolution of translatomes and transcriptomes using ribosome-profling and matched RNA-sequencing data for three organs (brain, liver and testis) in fve mammals (human, macaque, mouse, opossum and platypus) and a bird (chicken). Our within-species analyses reveal that translational regulation is widespread in the diferent organs, in particular across the spermatogenic cell types of the testis. The between-species divergence in gene expression is around 20% lower at the translatome layer than at the transcriptome layer owing to extensive buffering between the expression layers, which especially preserved old, essential and housekeeping genes. Translational upregulation specifcally counterbalanced global dosage reductions during the evolution of sex chromosomes and the efects of meiotic sex-chromosome inactivation during spermatogenesis. Despite the overall prevalence of bufering, some genes evolved faster at the translatome layer—potentially indicating adaptive changes in expression; testis tissue shows the highest fraction of such genes. Further analyses incorporating mass spectrometry proteomics data establish that the co-evolution of transcriptomes and translatomes is refected at the proteome layer. Together, our work uncovers co-evolutionary patterns and associated selective forces across the expression layers, and provides a resource for understanding their interplay in mammalian organs.

Expression Evolution of Ancestral XY Gametologs across All Major Groups of Placental Mammals

Placental mammals present 180 million-year-old Y chromosomes that have retained a handful of dosage-sensitive genes. However, the expression evolution of Y-linked genes across placental groups has remained largely unexplored. Here, we expanded the number of Y gametolog sequences by analyzing ten additional species from previously unexplored groups. We detected seven remarkably conserved genes across 25 placental species with known Y repertoires. We then used RNA-seq data from 17 placental mammals to unveil the expression evolution of XY gametologs. We found that Y gametologs followed, on average, a 3-fold expression loss and that X gametologs also experienced some expression reduction, particularly in primates. Y gametologs gained testis specificity through an accelerated expression decay in somatic tissues. Moreover, despite the substantial expression decay of Y genes, the combined expression of XY gametologs in males is higher than that of both X gametologs in females. Finally, our work describes several features of the Y chromosome in the last common mammalian ancestor.

Eutherian-Specific Gene TRIML2 Attenuates Inflammation in the Evolution of Placentation

Evolution of highly invasive placentation in the stem lineage of eutherians and subsequent extension of pregnancy set eutherians apart from other mammals, that is, marsupials with short-lived placentas, and oviparous monotremes. Recent studies suggest that eutherian implantation evolved from marsupial attachment reaction, an inflammatory process induced by the direct contact of fetal placenta with maternal endometrium after the breakdown of the shell coat, and shortly before the onset of parturition. Unique to eutherians, a dramatic downregulation of inflammation after implantation prevents the onset of premature parturition, and is critical for the maintenance of gestation. This downregulation likely involved evolutionary changes on maternal as well as fetal/placental side. Tripartite-motif family-like2 (TRIML2) only exists in eutherian genomes and shows preferential expression in preimplantation embryos, and trophoblast-derived structures, such as chorion and placental disc. Comparative genomic evidence supports that TRIML2 originated from a gene duplication event in the stem lineage of Eutheria that also gave rise to eutherian TRIML1. Compared with TRIML1, TRIML2 lost the catalytic RING domain of E3 ligase. However, only TRIML2 is induced in human choriocarcinoma cell line JEG3 with poly(I:C) treatment to simulate inflammation during viral infection. Its knockdown increases the production of proinflammatory cytokines and reduces trophoblast survival during poly(I:C) stimulation, while its overexpression reduces proinflammatory cytokine production, supporting TRIML2’s role as a regulatory inhibitor of the inflammatory pathways in trophoblasts. TRIML2’s potential virus-interacting PRY/SPRY domain shows significant signature of selection, suggesting its contribution to the evolution of eutherian-specific inflammation regulation during placentation.

Conservation, acquisition, and functional impact of sex-biased gene expression in mammals

Sex differences abound in human health and disease, as they do in other mammals used as models. The extent to which sex differences are conserved at the molecular level across species and tissues is unknown. We surveyed sex differences in gene expression in human, macaque, mouse, rat, and dog, across 12 tissues. In each tissue, we identified hundreds of genes with conserved sex-biased expression-findings that, combined with genomic analyses of human height, explain ~12% of the difference in height between females and males. We surmise that conserved sex biases in expression of genes otherwise operating equivalently in females and males contribute to sex differences in traits. However, most sex-biased expression arose during the mammalian radiation, which suggests that careful attention to interspecies divergence is needed when modeling human sex differences.

Deciphering Ancestral Sex Chromosome Turnovers Based on Analysis of Male Mutation Bias

The age of sex chromosomes is commonly obtained by comparing the substitution rates of XY gametologs. Coupled with phylogenetic reconstructions, one can refine the origin of a sex chromosome system relative to specific speciation events. However, these approaches are insufficient to determine the presence and duration of ancestral sex chromosome systems that were lost in some species. In this study, we worked with genomic and transcriptomic data from mammals and squamates and analyzed the effect of male mutation bias on X-linked sequences in these groups. We searched for signatures indicating whether monotremes shared the same sex chromosomes with placental mammals or whether pleurodonts and acrodonts had a common ancestral sex chromosome system. Our analyses indicate that platypus did not share the XY chromosomes with placental mammals, in agreement with previous work. In contrast, analyses of agamids showed that this lineage maintained the pleurodont XY chromosomes for several million years. We performed multiple simulations using different strengths of male mutation bias to confirm the results. Overall, our work shows that variations in substitution rates due to male mutation bias could be applied to uncover signatures of ancestral sex chromosome systems.

Corytophanids Replaced the Pleurodont XY System with a New Pair of XY Chromosomes

Almost all lizard families in the pleurodont clade share the same XY system. This system was meticulously studied in Anolis carolinensis, where it shows a highly degenerated Y chromosome and a male-specific X chromosome dosage compensation mechanism. Corytophanids (casque-headed lizards) have been proposed as the only family in the pleurodont clade to lack the XY system. In this study, we worked with extensive genomic and transcriptomic data from Basiliscus vittatus, a member of the Corytophanidae family that inhabits the tropical rainforests of Mexico. We confirmed that B. vittatus underwent a sex chromosome system turnover, which consisted in the loss of the pleurodont XY system and the gain of a new pair of XY chromosomes that are orthologous to chicken chromosome 17. We estimated the origin of the sex chromosome system to have occurred ∼63 Ma in the ancestor of corytophanids. Moreover, we identified 12 XY gametologues with particular attributes, such as functions related to the membrane and intracellular trafficking, very low expression levels, blood specificity, and incomplete dosage compensation in males.

Programmed DNA elimination of germline development genes in songbirds

In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline–soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits.

Songbirds have extensive germline–soma genome differences due to developmental elimination of a germline-specific chromosome (GRC). Here, the authors show that the GRC contains dozens of expressed developmental genes, some of which have been on the GRC since the ancestor of all songbirds.

Ancestral male recombination in Drosophila albomicans produced geographically restricted neo-Y chromosome haplotypes varying in age and onset of decay

Male Drosophila typically have achiasmatic meiosis, and fusions between autosomes and the Y chromosome have repeatedly created non-recombining neo-Y chromosomes that degenerate. Intriguingly, Drosophila nasuta males recombine, but their close relative D. albomicans reverted back to achiasmy after evolving neo-sex chromosomes. Here we use genome-wide polymorphism data to reconstruct the complex evolutionary history of neo-sex chromosomes in D. albomicans and examine the effect of recombination and its cessation on the initiation of neo-Y decay. Population and phylogenomic analyses reveal three distinct neo-Y types that are geographically restricted. Due to ancestral recombination with the neo-X, overall nucleotide diversity on the neo-Y is similar to the neo-X but severely reduced within neo-Y types. Consistently, the neo-Y chromosomes fail to form a monophyletic clade in sliding window trees outside of the region proximal to the fusion. Based on tree topology changes, we inferred the recombination breakpoints that produced haplotypes specific to each neo-Y type. We show that recombination became suppressed at different time points for the different neo-Y haplotypes. Haplotype age correlates with onset of neo-Y decay, and older neo-Y haplotypes show more fixed gene disruption via frameshift indels and down-regulation of neo-Y alleles. Genes are downregulated independently on the different neo-Ys, but are depleted of testes-expressed genes across all haplotypes. This indicates that genes important for male function are initially shielded from degeneration. Our results offer a time course of the early progression of Y chromosome evolution, showing how the suppression of recombination, through the reversal to achiasmy in D. albomicans males, initiates the process of degeneration.

Dichotomy of Dosage Compensation along the Neo Z Chromosome of the Monarch Butterfly

Mechanisms of sex chromosome dosage compensation (SCDC) differ strikingly among animals. In Drosophila flies, chromosome-wide transcription is doubled from the single X chromosome in hemizygous (XY) males, whereas in Caenorhabditis nematodes, expression is halved for both X copies in homozygous (XX) females [1, 2]. Unlike other female-heterogametic (WZ female and ZZ male) animals, moths and butterflies exhibit sex chromosome dosage compensation patterns typically seen only in male-heterogametic species [3]. The monarch butterfly carries a newly derived Z chromosome segment that arose from an autosomal fusion with the ancestral Z [4]. Using a highly contiguous genome assembly, we show that gene expression is balanced between sexes along the entire Z chromosome but with distinct modes of compensation on the two segments. On the ancestral Z segment, depletion of H4K16ac corresponds to nearly halving of biallelic transcription in males, a pattern convergent to nematodes. Conversely, the newly derived Z segment shows a Drosophila-like mode of compensation, with enriched H4K16ac levels corresponding to doubled monoallelic transcription in females. Our work reveals that, contrary to the expectation of co-opting regulatory mechanisms readily in place, the evolution of plural modes of dosage compensation is also possible along a single sex chromosome within a species.

Escaping the evolutionary trap? Sex chromosome turnover in basilisks and related lizards (Corytophanidae: Squamata)

Most pleurodont lizard families (anoles, iguanas and their relatives), with the exception of the basilisks and casquehead lizards (family Corytophanidae), share homologous XX/XY sex chromosomes, syntenic with chicken chromosome 15. Here, we used a suite of methods (i.e. RADseq, RNAseq and qPCR) to identify corytophanid sex chromosomes for the first time. We reveal that all examined corytophanid species have partially degenerated XX/XY sex chromosomes, syntenic with chicken chromosome 17. Transcriptomic analyses showed that the expression of X-linked genes in the corytophanid, Basiliscus vittatus, is not balanced between the sexes, which is rather exceptional under male heterogamety, and unlike the dosage-balanced sex chromosomes in other well-studied XX/XY systems, including the green anole, Anolis carolinensis. Corytophanid sex chromosomes may represent a rare example of a turnover away from stable, differentiated sex chromosomes. However, because of poor phylogenetic resolution among pleurodont families, we cannot reject the alternative hypothesis that corytophanid sex chromosomes evolved independently from an unknown ancestral system.

Genomic Data Reveal Conserved Female Heterogamety in Giant Salamanders with Gigantic Nuclear Genomes

Systems of genetic sex determination and the homology of sex chromosomes in different taxa vary greatly across vertebrates. Much progress remains to be made in understanding systems of genetic sex determination in non-model organisms, especially those with homomorphic sex chromosomes and/or large genomes. We used reduced representation genome sequencing to investigate genetic sex determination systems in the salamander family Cryptobranchidae (genera Cryptobranchus and Andrias), which typifies both of these inherent difficulties. We tested hypotheses of male- or female-heterogamety by sequencing hundreds of thousands of anonymous genomic regions in a panel of known-sex cryptobranchids and characterized patterns of presence/absence, inferred zygosity, and depth of coverage to identify sex-linked regions of these 56 gigabase genomes. Our results strongly support the hypothesis that all cryptobranchid species possess homologous systems of female heterogamety, despite maintenance of homomorphic sex chromosomes over nearly 60 million years. Additionally, we report a robust, non-invasive genetic assay for sex diagnosis in Cryptobranchus and Andrias which may have great utility for conservation efforts with these endangered salamanders. Co-amplification of these W-linked markers in both cryptobranchid genera provides evidence for long-term sex chromosome stasis in one of the most divergent salamander lineages. These findings inform hypotheses about the ancestral mode of sex determination in salamanders, but suggest that comparative data from other salamander families are needed. Our results further demonstrate that massive genomes are not necessarily a barrier to effective genome-wide sequencing and that the resulting data can be highly informative about sex determination systems in taxa with homomorphic sex chromosomes.

Regional epigenetic differentiation of the Z Chromosome between sexes in a female heterogametic system

In male heterogametic systems, the X Chromosome is epigenetically differentiated between males and females, to facilitate dosage compensation. For example, the X Chromosome in female mammals is largely inactivated. Relative to well-studied male heterogametic systems, the extent of epigenetic differentiation between male and female Z Chromosomes in female heterogametic species, which often lack complete dosage compensation, is poorly understood. Here, we examined the chromosomal DNA methylation landscapes of male and female Z Chromosomes in two distantly related avian species, namely chicken and white-throated sparrow. We show that, in contrast to the pattern in mammals, male and female Z Chromosomes in these species exhibit highly similar patterns of DNA methylation, which is consistent with weak or absent dosage compensation. We further demonstrate that the epigenetic differences between male and female chicken Z Chromosomes are localized to a few regions, including a previously identified male hypermethylated region 1 (MHM1; CGNC: 80601). We discovered a novel region with elevated male-to-female methylation ratios on the chicken Z Chromosome (male hypermethylated region 2 [MHM2]; CGNC: 80602). The MHM1 and MHM2, despite little sequence similarity between them, bear similar molecular features that are likely associated with their functions. We present evidence consistent with female hypomethylation of MHMs and up-regulation of nearby genes. Therefore, despite little methylation differentiation between sexes, extremely localized DNA methylation differences between male and female chicken Z Chromosomes have evolved and affect expression of nearby regions. Our findings offer new insights into epigenetic regulation of gene expression between sexes in female heterogametic systems.

Gene expression across mammalian organ development

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified developmental stage correspondences across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

The eXceptional nature of the X chromosome

The X chromosome is unique in the genome. In this review we discuss recent advances in our understanding of the genetics and epigenetics of the X chromosome. The X chromosome shares limited conservation with its ancestral homologue the Y chromosome and the resulting difference in X-chromosome dosage between males and females is largely compensated for by X-chromosome inactivation. The process of inactivation is initiated by the long non-coding RNA X-inactive specific transcript (XIST) and achieved through interaction with multiple synergistic silencing pathways. Identification of Xist-interacting proteins has given insight into these processes yet the cascade of events from initiation to maintenance have still to be resolved. In particular, the initiation of inactivation in humans has been challenging to study as: it occurs very early in development; most human embryonic stem cell lines already have an inactive X; and the process seems to differ from mouse. Another difference between human and mouse X inactivation is the larger number of human genes that escape silencing. In humans over 20% of X-linked genes continue to be expressed from the otherwise inactive X chromosome. We are only beginning to understand how such escape occurs but there is growing recognition that escapees contribute to sexually dimorphic traits. The unique biology and epigenetics of the X chromosome have often led to its exclusion from disease studies, yet the X constitutes 5% of the genome and is an important contributor to disease, often in a sex-specific manner.

The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes

Here we use a chromosome-level genome assembly of a prairie rattlesnake (Crotalus viridis), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies.

Guy1, a Y-linked embryonic signal, regulates dosage compensation in Anopheles stephensi by increasing X gene expression

We previously showed that Guy1, a primary signal expressed from the Y chromosome, is a strong candidate for a male-determining factor that confers female-specific lethality in Anopheles stephensi (Criscione et al., 2016). Here, we present evidence that Guy1 increases X gene expression in Guy1-transgenic females from two independent lines, providing a mechanism underlying the Guy1-conferred female lethality. The median level gene expression (MGE) of X-linked genes is significantly higher than autosomal genes in Guy1-transgenic females while there is no significant difference in MGE between X and autosomal genes in wild-type females. Furthermore, Guy1 significantly upregulates at least 40% of the 996 genes across the X chromosome in transgenic females. Guy1-conferred female-specific lethality is remarkably stable and completely penetrant. These findings indicate that Guy1 regulates dosage compensation in An. stephensi and components of dosage compensation may be explored to develop novel strategies to control mosquito-borne diseases.