Transcriptomic and proteomic elucidation of Z chromosome dosage compensation in Helicoverpa armigera.

Transcriptomic data have been used to study sex chromosome dosage compensation (SCDC) in approximately 10 Lepidoptera ZW species, yielding a consensus compensation pattern of Z ≈ ZZ < AA . $$ \approx \mathrm{ZZ}<\mathrm{AA}. $$ It remains unclear whether this compensation pattern holds when examining more Lepidoptera ZW species and/or using proteomic data to analyse SCDC. Here we combined transcriptomic and proteomic data as well as transcriptional level of six individual Z genes to reveal the SCDC pattern in Helicoverpa armigera, a polyphagous lepidopteran pest of economic importance. Transcriptomic analysis showed that the Z chromosome expression of H. armigera was balanced between male and female but substantially reduced relative to autosome expression, exhibiting an SCDC pattern of Z ≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ . When using H. amigera midgut proteomic data, the SCDC pattern of this species changed from Z ≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ at transcriptomic level to Z = ZZ = AA at the proteomic level. RT-qPCR analysis of transcript abundance of six Z genes found that compensation for each Z gene could vary from no compensation to overcompensation, depending on the individual genes and tissues tested. These results demonstrate for the first time the existence of a translational compensation mechanism, which is operating in addition to a translational mechanism, such as has been reported in other lepidopteran species. And the transcriptional compensation mechanism functions to accomplish Z chromosome dosage balance between the sexes (M = F on the Z chromosome), whereas the translation compensation mechanism operates to achieve dosage compensation between Z chromosome and autosome (Z = AA).

Genome-wide association mapping uncovers sex-associated copy number variation markers and female hemizygous regions on the W chromosome in Salix viminalis.

BACKGROUND: Sex chromosomes are in some species largely undifferentiated (homomorphic) with restricted sex determination regions. Homomorphic but different sex chromosomes are found in the closely related genera Populus and Salix indicating flexible sex determination systems, ideal for studies of processes involved in sex chromosome evolution. We have performed genome-wide association studies of sex and analysed sex chromosomes in a population of 265 wild collected Salix viminalis accessions and studied the sex determining locus. RESULTS: A total of 19,592 markers were used in association analyses using both Fisher's exact tests and a single-marker mixed linear model, which resulted in 48 and 41 sex-associated (SA) markers respectively. Across all 48 SA markers, females were much more often heterozygous than males, which is expected if females were the heterogametic sex. The majority of the SA markers were, based on positions in the S. purpurea genome, located on chromosome 15, previously demonstrated to be the sex chromosome. Interestingly, when mapping the genotyping-by-sequencing sequence tag harbouring the two SA markers with the highest significance to the S. viminalis genomic scaffolds, five regions of very high similarity were found: three on a scaffold that represents a part of chromosome 15, one on a scaffold that represents a part of chromosome 9 and one on a scaffold not anchored to the genome. Based on segregation differences of the alleles at the two marker positions and on differences in PCR amplification between females and males we conclude that females had multiple copies of this DNA fragment (chromosome 9 and 15), whereas males only had one (chromosome 9). We therefore postulate that the female specific sequences have been copied from chromosome 9 and inserted on chromosome 15, subsequently developing into a hemizygous W chromosome linked region. CONCLUSIONS: Our results support that sex determination in S. viminalis is controlled by one locus on chromosome 15. The segregation patterns observed at the SA markers furthermore confirm that S. viminalis females are the heterogametic sex. We also identified a translocation from chromosome 9 to the W chromosome.

Epistatic interactions between sex chromosomes and autosomes can affect the stability of sex determination systems.

Sex determination (SD) is an essential and ancient developmental process, but the genetic systems that regulate this process are surprisingly variable. Why SD mechanisms vary so much is a longstanding question in evolutionary biology. SD genes are generally located on sex chromosomes which also carry genes that interact epistatically with autosomes to affect fitness. How this affects the evolutionary stability of SD mechanisms is still unknown. Here, we explore how epistatic interactions between a sexually antagonistic (SA) non-SD gene, located on either an ancestral or novel sex chromosome, and an autosomal gene affect the conditions under which an evolutionary transition to a new SD system occurs. We find that when the SD gene is linked to an ancestral sex-chromosomal gene which engages in epistatic interactions, epistasis enhances the stability of the sex chromosomes so that they are retained under conditions where transitions would otherwise occur. This occurs both when weaker fitness effects are associated with the ancestral sex chromosome pair or stronger fitness effects associated with a newly evolved SD gene. However, the probability that novel SD genes spread is unaffected if they arise near genes involved in epistasis. This discrepancy occurs because, on autosomes, SA allele frequencies are typically lower than on sex chromosomes. In our model, increased frequencies of these alleles contribute to a higher frequency of epistasis which may therefore more readily occur on sex chromosomes. Because sex chromosome-autosome interactions are abundant and can take several forms, they may play a large role in maintaining sex chromosomes.

The Genome of Armadillidium vulgare (Crustacea, Isopoda) Provides Insights into Sex Chromosome Evolution in the Context of Cytoplasmic Sex Determination.

The terrestrial isopod Armadillidium vulgare is an original model to study the evolution of sex determination and symbiosis in animals. Its sex can be determined by ZW sex chromosomes, or by feminizing Wolbachia bacterial endosymbionts. Here, we report the sequence and analysis of the ZW female genome of A. vulgare. A distinguishing feature of the 1.72 gigabase assembly is the abundance of repeats (68% of the genome). We show that the Z and W sex chromosomes are essentially undifferentiated at the molecular level and the W-specific region is extremely small (at most several hundreds of kilobases). Our results suggest that recombination suppression has not spread very far from the sex-determining locus, if at all. This is consistent with A. vulgare possessing evolutionarily young sex chromosomes. We characterized multiple Wolbachia nuclear inserts in the A. vulgare genome, none of which is associated with the W-specific region. We also identified several candidate genes that may be involved in the sex determination or sexual differentiation pathways. The A. vulgare genome serves as a resource for studying the biology and evolution of crustaceans, one of the most speciose and emblematic metazoan groups.

Global Dosage Compensation Is Ubiquitous in Lepidoptera, but Counteracted by the Masculinization of the Z Chromosome.

While chromosome-wide dosage compensation of the X chromosome has been found in many species, studies in ZW clades have indicated that compensation of the Z is more localized and/or incomplete. In the ZW Lepidoptera, some species show complete compensation of the Z chromosome, while others lack full equalization, but what drives these inconsistencies is unclear. Here, we compare patterns of male and female gene expression on the Z chromosome of two closely related butterfly species, Papilio xuthus and Papilio machaon, and in multiple tissues of two moths species, Plodia interpunctella and Bombyx mori, which were previously found to differ in the extent to which they equalize Z-linked gene expression between the sexes. We find that, while some species and tissues seem to have incomplete dosage compensation, this is in fact due to the accumulation of male-biased genes and the depletion of female-biased genes on the Z chromosome. Once this is accounted for, the Z chromosome is fully compensated in all four species, through the up-regulation of Z expression in females and in some cases additional down-regulation in males. We further find that both sex-biased genes and Z-linked genes have increased rates of expression divergence in this clade, and that this can lead to fast shifts in patterns of gene expression even between closely related species. Taken together, these results show that the uneven distribution of sex-biased genes on sex chromosomes can confound conclusions about dosage compensation and that Z chromosome-wide dosage compensation is not only possible but ubiquitous among Lepidoptera.

The Z chromosome is enriched for sperm proteins in two divergent species of Lepidoptera.

Genes that promote sexual conflict, such as those with a sex-limited fitness benefit, are expected to accumulate differentially on sex chromosomes relative to autosomes. Few tests of this hypothesis exist for male homogametic (ZZ) taxa, however, and most use RNA expression data to identify such genes. Here, we employ a different identification method by using proteomic analysis of sperm cells to identify genes with a sex-limited benefit. We tested for a bias in genomic location of sperm protein genes in two species of Lepidoptera. An excess of sperm protein genes was identified on the Z chromosomes of both the Carolina sphinx moth (Manduca sexta) and the monarch butterfly (Danaus plexippus). Taking into consideration a Z-autosome fusion in monarchs, we discover that the ancestrally sex-linked portion of the genome is the source of this enrichment, while the newly sex-linked portion still appears similar to autosomes in relative abundance of sperm protein genes. Together, these results point to an enrichment of male-beneficial genes on the Z chromosome and demonstrate the usefulness of proteomic datasets in sexual conflict research.

Positive Selection Underlies Faster-Z Evolution of Gene Expression in Birds.

The elevated rate of evolution for genes on sex chromosomes compared with autosomes (Fast-X or Fast-Z evolution) can result either from positive selection in the heterogametic sex or from nonadaptive consequences of reduced relative effective population size. Recent work in birds suggests that Fast-Z of coding sequence is primarily due to relaxed purifying selection resulting from reduced relative effective population size. However, gene sequence and gene expression are often subject to distinct evolutionary pressures; therefore, we tested for Fast-Z in gene expression using next-generation RNA-sequencing data from multiple avian species. Similar to studies of Fast-Z in coding sequence, we recover clear signatures of Fast-Z in gene expression; however, in contrast to coding sequence, our data indicate that Fast-Z in expression is due to positive selection acting primarily in females. In the soma, where gene expression is highly correlated between the sexes, we detected Fast-Z in both sexes, although at a higher rate in females, suggesting that many positively selected expression changes in females are also expressed in males. In the gonad, where intersexual correlations in expression are much lower, we detected Fast-Z for female gene expression, but crucially, not males. This suggests that a large amount of expression variation is sex-specific in its effects within the gonad. Taken together, our results indicate that Fast-Z evolution of gene expression is the product of positive selection acting on recessive beneficial alleles in the heterogametic sex. More broadly, our analysis suggests that the adaptive potential of Z chromosome gene expression may be much greater than that of gene sequence, results which have important implications for the role of sex chromosomes in speciation and sexual selection.

ZW sex-chromosome evolution and contagious parthenogenesis in Artemia brine shrimp.

Eurasian brine shrimp (genus Artemia) have closely related sexual and asexual lineages of parthenogenetic females, which produce rare males at low frequencies. Although they are known to have ZW chromosomes, these are not well characterized, and it is unclear whether they are shared across the clade. Furthermore, the underlying genetic architecture of the transmission of asexuality, which can occur when rare males mate with closely related sexual females, is not well understood. We produced a chromosome-level assembly for the sexual Eurasian species Artemia sinica and characterized in detail the pair of sex chromosomes of this species. We combined this new assembly with short-read genomic data for the sexual species Artemia sp. Kazakhstan and several asexual lineages of Artemia parthenogenetica, allowing us to perform an in-depth characterization of sex-chromosome evolution across the genus. We identified a small differentiated region of the ZW pair that is shared by all sexual and asexual lineages, supporting the shared ancestry of the sex chromosomes. We also inferred that recombination suppression has spread to larger sections of the chromosome independently in the American and Eurasian lineages. Finally, we took advantage of a rare male, which we backcrossed to sexual females, to explore the genetic basis of asexuality. Our results suggest that parthenogenesis is likely partly controlled by a locus on the Z chromosome, highlighting the interplay between sex determination and asexuality.

Dosage Compensation throughout the Schistosoma mansoni Lifecycle: Specific Chromatin Landscape of the Z Chromosome.

Differentiated sex chromosomes are accompanied by a difference in gene dose between X/Z-specific and autosomal genes. At the transcriptomic level, these sex-linked genes can lead to expression imbalance, or gene dosage can be compensated by epigenetic mechanisms and results into expression level equalization. Schistosoma mansoni has been previously described as a ZW species (i.e., female heterogamety, in opposition to XY male heterogametic species) with a partial dosage compensation, but underlying mechanisms are still unexplored. Here, we combine transcriptomic (RNA-Seq) and epigenetic data (ChIP-Seq against H3K4me3, H3K27me3, and H4K20me1 histone marks) in free larval cercariae and intravertebrate parasitic stages. For the first time, we describe differences in dosage compensation status in ZW females, depending on the parasitic status: free cercariae display global dosage compensation, whereas intravertebrate stages show a partial dosage compensation. We also highlight regional differences of gene expression along the Z chromosome in cercariae, but not in the intravertebrate stages. Finally, we feature a consistent permissive chromatin landscape of the Z chromosome in both sexes and stages. We argue that dosage compensation in schistosomes is characterized by chromatin remodeling mechanisms in the Z-specific region.

Multilayered Tuning of Dosage Compensation and Z-Chromosome Masculinization in the Wood White (Leptidea sinapis) Butterfly.

In species with genetic sex determination, dosage compensation can evolve to equal expression levels of sex-linked and autosomal genes. Current knowledge about dosage compensation has mainly been derived from male-heterogametic (XX/XY) model organisms, whereas less is understood about the process in female-heterogametic systems (ZZ/ZW). In moths and butterflies, downregulation of Z-linked expression in males (ZZ) to match the expression level in females (ZW) is often observed. However, little is known about the underlying regulatory mechanisms, or if dosage compensation patterns vary across ontogenetic stages. In this study, we assessed dynamics of Z-linked and autosomal expression levels across developmental stages in the wood white (Leptidea sinapis). We found that although expression of Z-linked genes in general was reduced compared with autosomal genes, dosage compensation was actually complete for some categories of genes, in particular sex-biased genes, but equalization in females was constrained to a narrower gene set. We also observed a noticeable convergence in Z-linked expression between males and females after correcting for sex-biased genes. Sex-biased expression increased successively across developmental stages, and male-biased genes were enriched on the Z-chromosome. Finally, all five core genes associated with the ribonucleoprotein dosage compensation complex male-specific lethal were detected in adult females, in correspondence with a reduction in the expression difference between autosomes and the single Z-chromosome. We show that tuning of gene dosage is multilayered in Lepidoptera and argue that expression balance across chromosomal classes may predominantly be driven by enrichment of male-biased genes on the Z-chromosome and cooption of available dosage regulators.

Triploid Colubrid Snake Provides Insight into the Mechanism of Sex Determination in Advanced Snakes.

The advanced snakes (Caenophidia), the important amniote lineage encompassing more than 3,000 living species, possess highly conserved female heterogamety across all families. However, we still lack any knowledge on the gene(s) and the molecular mechanism controlling sex determination. Triploid individuals spontaneously appear in populations of diploid species and can provide an important insight into the evolution of sex determination. Here, we report a case of spontaneous triploidy in a male of the twin-spotted ratsnake (Elaphe bimaculata) with ZZW sex chromosomes. We speculate that as both ZZ and ZZW individuals develop male gonads, the ratio between the number of Z chromosomes and autosomes, and not the presence of the W chromosome in the genome, drives sex determination in the advanced snakes.

Sex Chromosomes of the Iconic Moth Abraxas grossulariata (Lepidoptera, Geometridae) and Its Congener A. sylvata.

The magpie moth, Abraxas grossulariata, is an iconic species in which female heterogamety was discovered at the beginning of the 20th century. However, the sex chromosomes of this species have not yet been cytologically identified. We describe the sex chromosomes of A. grossulariata and its congener, A. sylvata. Although these species split only around 9.5 million years ago, and both species have the expected WZ/ZZ chromosomal system of sex determination and their sex chromosomes share the major ribosomal DNA (rDNA) representing the nucleolar organizer region (NOR), we found major differences between their karyotypes, including between their sex chromosomes. The species differ in chromosome number, which is 2n = 56 in A. grossularita and 2n = 58 in A. sylvata. In addition, A. grossularita autosomes exhibit massive autosomal blocks of heterochromatin, which is a very rare phenomenon in Lepidoptera, whereas the autosomes of A. sylvata are completely devoid of distinct heterochromatin. Their W chromosomes differ greatly. Although they are largely composed of female-specific DNA sequences, as shown by comparative genomic hybridization, cross-species W-chromosome painting revealed considerable sequence differences between them. The results suggest a relatively rapid molecular divergence of Abraxas W chromosomes by the independent spreading of female-specific repetitive sequences.

RNA sequencing reveals a complete but an unconventional type of dosage compensation in the domestic silkworm Bombyx mori.

Sex chromosomal dose difference between sexes is often normalized by a gene regulatory mechanism called dosage compensation (DC). Studies indicate that DC mechanisms are generally effective in XY rather than ZW systems. However, DC studies in lepidopterans (ZW system) gave bewildering results. In Manduca sexta, DC was complete and in Plodia interpunctella, it was incomplete. In Heliconius species, dosage was found to be partly incomplete. In domesticated silkmoth Bombyx mori, DC studies have yielded contradictory results thus far, showing incomplete DC based on microarray data and a possible existence of DC based on recent reanalysis of same data. In this study, analysis of B. mori sexed embryos (78, 96 and 120 h) and larval heads using RNA sequencing suggest an onset of DC at 120 h. The average Z-linked expression is substantially less than autosomes, and the male-biased Z-linked expression observed at initial stages (78 and 96 h) gets almost compensated at 120 h embryonic stage and perfectly compensated in heads. Based on these findings, we suggest a complete but an unconventional type of DC, which may be achieved by reduced Z-linked expression in males (ZZ). To our knowledge, this is the first next-generation sequencing report showing DC in B. mori, clarifying the previous contradictions.

ZZ/ZW sex chromosome system in the endangered fish Lignobrycon myersi Miranda-Ribeiro, 1956 (Teleostei, Characiformes, Triportheidae).

Lignobrycon myersi is an endemic fish species from a few coastal rivers in northeastern Brazil. Based on molecular evidence, Lignobrycon myersi and genera Triportheus Cope, 1872, Agoniates Müller & Troschel, 1845, Clupeacharax Pearson, 1924 and Engraulisoma Castro, 1981 were placed in the family Triportheidae. In the present work, we report the first cytogenetic data for Lignobrycon myersi to test the hypothesis that Lignobrycon and Triportheus are closely related. Studied specimens presented 2n=52 with 28 metacentric (m), 18 submetacentric (sm) and six subtelocentric (st) chromosomes for males and 27 m, 19 sm and 6 st for females, characterizing a ZZ/ZW sex chromosome system. The Z chromosome corresponds to the largest chromosome in karyotype while the W is about 50% smaller than the Z and largely heterochromatic. Terminal nucleolus organizer regions, GC-rich sites and 18S rDNA signals were detected on pair 14. However, additional 18S rDNA sites were observed in the W chromosome. The 5S rDNA was mainly detected on long arms of pair 7. The apparent synapomorphic chromosomal traits of Triportheus and Lignobrycon myersi reinforce their close phylogenetic relationship, suggesting that the ZZ/ZW chromosome system in both genera has arisen before cladogenic events.

Independent stratum formation on the avian sex chromosomes reveals inter-chromosomal gene conversion and predominance of purifying selection on the W chromosome.

We used a comparative approach spanning three species and 90 million years to study the evolutionary history of the avian sex chromosomes. Using whole transcriptomes, we assembled the largest cross-species dataset of W-linked coding content to date. Our results show that recombination suppression in large portions of the avian sex chromosomes has evolved independently, and that long-term sex chromosome divergence is consistent with repeated and independent inversions spreading progressively to restrict recombination. In contrast, over short-term periods we observe heterogeneous and locus-specific divergence. We also uncover four instances of gene conversion between both highly diverged and recently evolved gametologs, suggesting a complex mosaic of recombination suppression across the sex chromosomes. Lastly, evidence from 16 gametologs reveal that the W chromosome is evolving with a significant contribution of purifying selection, consistent with previous findings that W-linked genes play an important role in encoding sex-specific fitness.