High rate of translocation-based gene birth on the Drosophila Y chromosome

Hiding in plain sight: the Y chromosome and its reinvigorated role in evolutionary processes

Recent methodological approaches have expanded our understanding of Y chromosome sequence, revealed unexpected Y diversity, and sparked a growing realization of its importance in evolutionary processes. To fully understand the diversity and importance of the Y chromosome, we suggest the need to move from a holotype Y chromosome sequence, based on a single individual and meant to represent the species, to a thorough understanding of Y chromosome haplotype diversity, its phenotypic implications, and its phylogeographic distribution. Additionally, the Y chromosome may play an important role in two key rules of speciation that have otherwise been attributed to the X, namely Haldane's Rule and the Large-X Effect. Emerging genomic tools and analytical approaches are just now giving us the means to ask how important this small, often forgotten region of the genome is in evolutionary processes.

Bdyof is a Y-chromosome-specific gene required for male development in Bactrocera dorsalis

BACKGROUND

In many organisms, the Y chromosome contains important genes associated with sex determination and male reproductive development. However, there have been few studies of Y-chromosome-specific genes in non-model species due to the incomplete information of Y chromosome genome and difficulty in sequencing. Here, we screened 90 candidate Y-specific sequences in a constructed transcriptome assembly library by using the chromosome quotient method, among which 11 were unreported sequences associated with male reproductive development, including Bactrocera dorsalis Y-specific Oligozoospermia factor (Bdyof) with the highest expression in the testis.

RESULTS

CRISPR/Cas9-mediated knockout of Bdyof resulted in abnormal testis development, significantly reduced sperm count, and obviously lower egg hatching rate in homozygous mutant flies. In addition, Bdyof knockout decreased the expression of dsx-M.

CONCLUSION

This results provides new insights into the biological processes related to male reproductive development controlled by the Y-chromosome-specific gene Bdyof, thus providing a promising molecular target for the study of agricultural pests. © 2024 Society of Chemical Industry.

Subcellular Enrichment Patterns of New Genes in Drosophila Evolution

The evolutionary patterns of proteins within subcellular compartments underlie the innovation and diversification foundation of the living eukaryotic organism. The location of proteins in subcellular compartments promotes the formation of network interaction modules, which in turn reshape the architecture of higher-level protein-protein interaction networks. Here, we conducted the most up-to-date gene age dating of Drosophila melanogaster by employing recently available long-read sequencing genomes as references. We found that an elevated gene fixation in the most recent common ancestor of Drosophila genus predated the divergence between two Drosophila subgenera, and a significant tendency of these genes in D. melanogaster encode proteins that localize to the extracellular matrix, accompanying the adaptive radiation of Drosophila species. Proteins encoded by genes located in the extracellular space exhibit higher sequence divergence, suggesting a rapid evolutionary process. We also observed that proteins encoded by genes originating from the same evolutionary branches tend to co-localize in the same subcellular compartments, and proteins in the same subcellular compartment tend to interact with each other. The proteins encoded by genes that have persisted through deeper branches exhibit broader localization across multiple subcellular compartments, enhancing the likelihood of their integration into various protein or gene regulatory networks, thereby increasing functional diversity. These evolutionary patterns not only contribute to understanding the evolution of subcellular localization in proteins encoded by genes originating from different branches, but also provide insights into the evolution of protein-protein networks driven by the emergence of new genes.

Identification and functional verification of Y-chromosome-specific gene typo-gyf in Bactrocera dorsalis

Genes on the Y chromosome play important roles in male sex determination and development. The identification of Y-chromosome-specific genes not only provides a theoretical basis for the study of male reproductive development, but also offers genetic control targets for agricultural pests. However, Y-chromosome genes are rarely characterized due to their high repeatability and high heterochromatinization, especially in the oriental fruit fly. In this study, 1 011 Y-chromosome-specific candidate sequences were screened from 2 to 4 h Bactrocera dorsalis embryo datasets with the chromosome quotient method, 6 of which were identified as Y-chromosome-specific sequences by polymerase chain reaction, including typo-gyf, a 19 126-bp DNA sequence containing a 575-amino acid open reading frame. Testicular deformation and a significant reduction in sperm number were observed after typo-gyf knockdown with RNA interference in embryos. After typo-gyf knockout with clustered regularly interspaced palindromic repeats (CRISPR) / CRISPR-associated protein 9 in the embryonic stage, the sex ratio of the emergent adults was unbalanced, with far more females than males. A genotype analysis of these females with the Y-chromosome gene MoY revealed no sex reversal. Typo-gyf knockout led to the death of XY individuals in the embryonic stage. We conclude that typo-gyf is an essential gene for male survival, and is also involved in testicular development and spermatogenesis. The identification of typo-gyf and its functional verification provide insight into the roles of Y-chromosome genes in male development.

On the resolution of sexual conflict over shared traits

Anisogamy, different-sized male and female gametes, sits at the heart of sexual selection and conflict between the sexes. Sperm producers (males) and egg producers (females) of the same species generally share most, if not all, of the same genome, but selection frequently favours different trait values in each sex for traits common to both. The extent to which this conflict might be resolved, and the potential mechanisms by which this can occur, have been widely debated. Here, we summarize recent findings and emphasize that once the sexes evolve, sexual selection is ongoing, and therefore new conflict is always possible. In addition, sexual conflict is largely a multivariate problem, involving trait combinations underpinned by networks of interconnected genes. Although these complexities can hinder conflict resolution, they also provide multiple possible routes to decouple male and female phenotypes and permit sex-specific evolution. Finally, we highlight difficulty in the study of sexual conflict over shared traits and promising directions for future research.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

The genomic distribution of transposable elements is driven by spatially variable purifying selection

It is widely accepted that the genomic distribution of transposable elements (TEs) mainly reflects the outcome of purifying selection and insertion bias (1). Nevertheless, the relative importance of these two evolutionary forces could not be tested thoroughly. Here, we introduce an experimental system, which allows separating purifying selection from TE insertion bias. We used experimental evolution to study the TE insertion patterns in Drosophila simulans founder populations harboring 1040 insertions of an active P-element. After 10 generations at a large population size, we detected strong selection against P-element insertions. The exception were P-element insertions in genomic regions for which a strong insertion bias has been proposed (2-4). Because recurrent P-element insertions cannot explain this pattern, we conclude that purifying selection, with variable strength along the chromosomes, is the major determinant of the genomic distribution of P-elements. Genomic regions with relaxed purifying selection against P-element insertions exhibit normal levels of purifying selection against base substitutions. This suggests that different types of purifying selection operate on base substitutions and P-element insertions. Our results highlight the power of experimental evolution to understand basic evolutionary processes, which are difficult to infer from patterns of natural variation alone.

Sex-linked gene traffic underlies the acquisition of sexually dimorphic UV color vision in Heliconius butterflies

The acquisition of novel sexually dimorphic traits poses an evolutionary puzzle: How do new traits arise and become sex-limited? Recently acquired color vision, sexually dimorphic in animals like primates and butterflies, presents a compelling model for understanding how traits become sex-biased. For example, some Heliconius butterflies uniquely possess UV (ultraviolet) color vision, which correlates with the expression of two differentially tuned UV-sensitive rhodopsins, UVRh1 and UVRh2. To discover how such traits become sexually dimorphic, we studied Heliconius charithonia, which exhibits female-specific UVRh1 expression. We demonstrate that females, but not males, discriminate different UV wavelengths. Through whole-genome shotgun sequencing and assembly of the H. charithonia genome, we discovered that UVRh1 is present on the W chromosome, making it obligately female-specific. By knocking out UVRh1, we show that UVRh1 protein expression is absent in mutant female eye tissue, as in wild-type male eyes. A PCR survey of UVRh1 sex-linkage across the genus shows that species with female-specific UVRh1 expression lack UVRh1 gDNA in males. Thus, acquisition of sex linkage is sufficient to achieve female-specific expression of UVRh1, though this does not preclude other mechanisms, like cis-regulatory evolution from also contributing. Moreover, both this event, and mutations leading to differential UV opsin sensitivity, occurred early in the history of Heliconius. These results suggest a path for acquiring sexual dimorphism distinct from existing mechanistic models. We propose a model where gene traffic to heterosomes (the W or the Y) genetically partitions a trait by sex before a phenotype shifts (spectral tuning of UV sensitivity).

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Searching for signatures of sexually antagonistic selection on stickleback sex chromosomes

Intralocus sexually antagonistic selection occurs when an allele is beneficial to one sex but detrimental to the other. This form of selection is thought to be key to the evolution of sex chromosomes but is hard to detect. Here we perform an analysis of phased young sex chromosomes to look for signals of sexually antagonistic selection in the Japan Sea stickleback (Gasterosteus nipponicus). Phasing allows us to date the suppression of recombination on the sex chromosome and provides unprecedented resolution to identify sexually antagonistic selection in the recombining region of the chromosome. We identify four windows with elevated divergence between the X and Y in the recombining region, all in or very near genes associated with phenotypes potentially under sexually antagonistic selection in humans. We are unable, however, to rule out the alternative hypothesis that the peaks of divergence result from demographic effects. Thus, although sexually antagonistic selection is a key hypothesis for the formation of supergenes on sex chromosomes, it remains challenging to detect.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

Evolution of the sex-determining region in Ginkgo biloba

Sex chromosomes or sex-determining regions (SDR) have been discovered in many dioecious plant species, including the iconic 'living fossil' Ginkgo biloba, though the location and size of the SDR in G. biloba remain contradictory. Here we resolve these controversies and analyse the evolution of the SDR in this species. Based on transcriptome sequencing data from four genetic crosses we reconstruct male- and female-specific genetic maps and locate the SDR to the middle of chromosome 2. Integration of the genetic maps with the genome sequence reveals that recombination in and around the SDR is suppressed in a region of about 50 Mb in both males and females. However, occasional recombination does occur except a small, less than 5 Mb long region that does not recombine in males. Based on synonymous divergence between homologous X- and Y-linked genes in this region, we infer that the Ginkgo SDR is fairly old-at least of Cretaceous origin. The analysis of substitution rates and gene expression reveals only slight Y-degeneration. These results are consistent with findings in other dioecious plants with homomorphic sex chromosomes, where the SDR is typically small and evolves in a region with pre-existing reduced recombination, surrounded by long actively recombining pseudoautosomal regions. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Gene duplication to the Y chromosome in Trinidadian Guppies

Differences in allele frequencies at autosomal genes between males and females in a population can result from two scenarios. First, unresolved sexual conflict over survival can produce allelic differentiation between the sexes. However, given the substantial mortality costs required to produce allelic differences between males and females at each generation, it remains unclear how many loci within the genome experience significant sexual conflict over survival. Alternatively, recent studies have shown that similarity between autosomal and Y sequences can create perceived allelic differences between the sexes. However, Y duplications are most likely in species with large non-recombining regions, in part because they simply represent larger targets for duplications. We assessed the genomes of 120 wild-caught guppies, which experience extensive predation- and pathogen-induced mortality and have a relatively small ancestral Y chromosome. We identified seven autosomal genes that show allelic differences between male and female adults. Five of these genes show clear evidence of whole or partial gene duplication between the Y chromosome and the autosomes. The remaining two genes show evidence of partial homology to the Y. Overall, our findings suggest that the guppy genome experiences a very low level of unresolved sexual conflict over survival, and instead the Y chromosome, despite its small ancestral size and recent origin, may nonetheless accumulate genes with male-specific functions.

Unique structure and positive selection promote the rapid divergence of Drosophila Y chromosomes

Y chromosomes across diverse species convergently evolve a gene-poor, heterochromatic organization enriched for duplicated genes, LTR retrotransposons, and satellite DNA. Sexual antagonism and a loss of recombination play major roles in the degeneration of young Y chromosomes. However, the processes shaping the evolution of mature, already degenerated Y chromosomes are less well-understood. Because Y chromosomes evolve rapidly, comparisons between closely related species are particularly useful. We generated de novo long-read assemblies complemented with cytological validation to reveal Y chromosome organization in three closely related species of the Drosophila simulans complex, which diverged only 250,000 years ago and share >98% sequence identity. We find these Y chromosomes are divergent in their organization and repetitive DNA composition and discover new Y-linked gene families whose evolution is driven by both positive selection and gene conversion. These Y chromosomes are also enriched for large deletions, suggesting that the repair of double-strand breaks on Y chromosomes may be biased toward microhomology-mediated end joining over canonical non-homologous end-joining. We propose that this repair mechanism contributes to the convergent evolution of Y chromosome organization across organisms.

New Genes in the Drosophila Y Chromosome: Lessons from D. willistoni

Y chromosomes play important roles in sex determination and male fertility. In several groups (e.g., mammals) there is strong evidence that they evolved through gene loss from a common X-Y ancestor, but in Drosophila the acquisition of new genes plays a major role. This conclusion came mostly from studies in two species. Here we report the identification of the 22 Y-linked genes in D. willistoni. They all fit the previously observed pattern of autosomal or X-linked testis-specific genes that duplicated to the Y. The ratio of gene gains to gene losses is ~25 in D. willistoni, confirming the prominent role of gene gains in the evolution of Drosophila Y chromosomes. We also found four large segmental duplications (ranging from 62 kb to 303 kb) from autosomal regions to the Y, containing ~58 genes. All but four of these duplicated genes became pseudogenes in the Y or disappeared. In the GK20609 gene the Y-linked copy remained functional, whereas its original autosomal copy degenerated, demonstrating how autosomal genes are transferred to the Y chromosome. Since the segmental duplication that carried GK20609 contained six other testis-specific genes, it seems that chance plays a significant role in the acquisition of new genes by the Drosophila Y chromosome.

The genetic architecture of temperature adaptation is shaped by population ancestry and not by selection regime

BACKGROUND

Understanding the genetic architecture of temperature adaptation is key for characterizing and predicting the effect of climate change on natural populations. One particularly promising approach is Evolve and Resequence, which combines advantages of experimental evolution such as time series, replicate populations, and controlled environmental conditions, with whole genome sequencing. Recent analysis of replicate populations from two different Drosophila simulans founder populations, which were adapting to the same novel hot environment, uncovered very different architectures-either many selection targets with large heterogeneity among replicates or fewer selection targets with a consistent response among replicates.

RESULTS

Here, we expose the founder population from Portugal to a cold temperature regime. Although almost no selection targets are shared between the hot and cold selection regime, the adaptive architecture was similar. We identify a moderate number of targets under strong selection (19 selection targets, mean selection coefficient = 0.072) and parallel responses in the cold evolved replicates. This similarity across different environments indicates that the adaptive architecture depends more on the ancestry of the founder population than the specific selection regime.

CONCLUSIONS

These observations will have broad implications for the correct interpretation of the genomic responses to a changing climate in natural populations.

Extreme Y chromosome polymorphism corresponds to five male reproductive morphs of a freshwater fish

Loss of recombination between sex chromosomes often depletes Y chromosomes of functional content and genetic variation, which might limit their potential to generate adaptive diversity. Males of the freshwater fish Poecilia parae occur as one of five discrete morphs, all of which shoal together in natural populations where morph frequency has been stable for over 50 years. Each morph uses a different complex reproductive strategy and morphs differ dramatically in colour, body size and mating behaviour. Morph phenotype is passed perfectly from father to son, indicating there are five Y haplotypes segregating in the species, which encode the complex male morph characteristics. Here, we examine Y diversity in natural populations of P. parae. Using linked-read sequencing on multiple P. parae females and males of all five morphs, we find that the genetic architecture of the male morphs evolved on the Y chromosome after recombination suppression had occurred with the X. Comparing Y chromosomes between each of the morphs, we show that, although the Ys of the three minor morphs that differ in colour are highly similar, there are substantial amounts of unique genetic material and divergence between the Ys of the three major morphs that differ in reproductive strategy, body size and mating behaviour. Altogether, our results suggest that the Y chromosome is able to overcome the constraints of recombination loss to generate extreme diversity, resulting in five discrete Y chromosomes that control complex reproductive strategies.

Evolution of genome structure in the Drosophila simulans species complex

The rapid evolution of repetitive DNA sequences, including satellite DNA, tandem duplications, and transposable elements, underlies phenotypic evolution and contributes to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented and misassembled in most contemporary genome assemblies. We generated highly contiguous de novo reference genomes for the Drosophila simulans species complex (D. simulans, D. mauritiana, and D. sechellia), which speciated ∼250,000 yr ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive sequences between these four species. We find that at least 15% of the D. simulans complex species genomes fail to align uniquely to D. melanogaster owing to structural divergence-twice the number of single-nucleotide substitutions. We also find rapid turnover of satellite DNA and extensive structural divergence in heterochromatic regions, whereas the euchromatic gene content is mostly conserved. Despite the overall preservation of gene synteny, euchromatin in each species has been shaped by clade- and species-specific inversions, transposable elements, expansions and contractions of satellite and tRNA tandem arrays, and gene duplications. We also find rapid divergence among Y-linked genes, including copy number variation and recent gene duplications from autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.

Fitness effects for Ace insecticide resistance mutations are determined by ambient temperature

BACKGROUND

Insect pest control programs often use periods of insecticide treatment with intermittent breaks, to prevent fixing of mutations conferring insecticide resistance. Such mutations are typically costly in an insecticide-free environment, and their frequency is determined by the balance between insecticide treatment and cost of resistance. Ace, a key gene in neuronal signaling, is a prominent target of many insecticides and across several species, three amino acid replacements (I161V, G265A, and F330Y) provide resistance against several insecticides. Because temperature disturbs neuronal signaling homeostasis, we reasoned that the cost of insecticide resistance could be modulated by ambient temperature.

RESULTS

Experimental evolution of a natural Drosophila simulans population at hot and cold temperature regimes uncovered a surprisingly strong effect of ambient temperature. In the cold temperature regime, the resistance mutations were strongly counter selected (s = - 0.055), but in a hot environment, the fitness costs of resistance mutations were reduced by almost 50% (s = - 0.031). We attribute this unexpected observation to the advantage of the reduced enzymatic activity of resistance mutations in hot environments.

CONCLUSION

We show that fitness costs of insecticide resistance genes are temperature-dependent and suggest that the duration of insecticide-free periods need to be adjusted for different climatic regions to reflect these costs. We suggest that such environment-dependent fitness effects may be more common than previously assumed and pose a major challenge for modeling climate change.

Sequence Transpositions Restore Genes on the Highly Degenerated W Chromosomes of Songbirds

The female-specific W chromosomes of most Neognathae birds are highly degenerated and gene-poor. Previous studies have demonstrated that the gene repertoires of the Neognathae bird W chromosomes, despite being in small numbers, are conserved across bird species, likely due to purifying selection maintaining the regulatory and dosage-sensitive genes. Here we report the discovery of DNA-based sequence duplications from the Z to the W chromosome in birds-of-paradise (Paradisaeidae, Passeriformes), through sequence transposition. The original transposition involved nine genes, but only two of them (ANXA1 and ALDH1A1) survived on the W chromosomes. Both ANXA1 and ALDH1A1 are predicted to be dosage-sensitive, and the expression of ANXA1 is restricted to ovaries in all the investigated birds. These analyses suggest the newly transposed gene onto the W chromosomes can be favored for their role in restoring dosage imbalance or through female-specific selection. After examining seven additional songbird genomes, we further identified five other transposed genes on the W chromosomes of Darwin's finches and one in the great tit, expanding the observation of the Z-to-W transpositions to a larger range of bird species, but not all transposed genes exhibit dosage-sensitivity or ovary-biased expression We demonstrate a new mechanism by which the highly degenerated W chromosomes of songbirds can acquire genes from the homologous Z chromosomes, but further functional investigations are needed to validate the evolutionary forces underlying the transpositions.

Patterns of Sex Chromosome Differentiation in Spiders: Insights from Comparative Genomic Hybridisation

Spiders are an intriguing model to analyse sex chromosome evolution because of their peculiar multiple X chromosome systems. Y chromosomes were considered rare in this group, arising after neo-sex chromosome formation by X chromosome-autosome rearrangements. However, recent findings suggest that Y chromosomes are more common in spiders than previously thought. Besides neo-sex chromosomes, they are also involved in the ancient X1X2Y system of haplogyne spiders, whose origin is unknown. Furthermore, spiders seem to exhibit obligatorily one or two pairs of cryptic homomorphic XY chromosomes (further cryptic sex chromosome pairs, CSCPs), which could represent the ancestral spider sex chromosomes. Here, we analyse the molecular differentiation of particular types of spider Y chromosomes in a representative set of ten species by comparative genomic hybridisation (CGH). We found a high Y chromosome differentiation in haplogyne species with X1X2Y system except for Loxosceles spp. CSCP chromosomes exhibited generally low differentiation. Possible mechanisms and factors behind the observed patterns are discussed. The presence of autosomal regions marked predominantly or exclusively with the male or female probe was also recorded. We attribute this pattern to intraspecific variability in the copy number and distribution of certain repetitive DNAs in spider genomes, pointing thus to the limits of CGH in this arachnid group. In addition, we confirmed nonrandom association of chromosomes belonging to particular CSCPs at spermatogonial mitosis and spermatocyte meiosis and their association with multiple Xs throughout meiosis. Taken together, our data suggest diverse evolutionary pathways of molecular differentiation in different types of spider Y chromosomes.

Assembly of the threespine stickleback Y chromosome reveals convergent signatures of sex chromosome evolution

BACKGROUND

Heteromorphic sex chromosomes have evolved repeatedly across diverse species. Suppression of recombination between X and Y chromosomes leads to degeneration of the Y chromosome. The progression of degeneration is not well understood, as complete sequence assemblies of heteromorphic Y chromosomes have only been generated across a handful of taxa with highly degenerate sex chromosomes. Here, we describe the assembly of the threespine stickleback (Gasterosteus aculeatus) Y chromosome, which is less than 26 million years old and at an intermediate stage of degeneration. Our previous work identified that the non-recombining region between the X and the Y spans approximately 17.5 Mb on the X chromosome.

RESULTS

We combine long-read sequencing with a Hi-C-based proximity guided assembly to generate a 15.87 Mb assembly of the Y chromosome. Our assembly is concordant with cytogenetic maps and Sanger sequences of over 90 Y chromosome BAC clones. We find three evolutionary strata on the Y chromosome, consistent with the three inversions identified by our previous cytogenetic analyses. The threespine stickleback Y shows convergence with more degenerate sex chromosomes in the retention of haploinsufficient genes and the accumulation of genes with testis-biased expression, many of which are recent duplicates. However, we find no evidence for large amplicons identified in other sex chromosome systems. We also report an excellent candidate for the master sex-determination gene: a translocated copy of Amh (Amhy).

CONCLUSIONS

Together, our work shows that the evolutionary forces shaping sex chromosomes can cause relatively rapid changes in the overall genetic architecture of Y chromosomes.

Evolution of novel genes in three-spined stickleback populations

Eukaryotic genomes frequently acquire new protein-coding genes which may significantly impact an organism’s fitness. Novel genes can be created, for example, by duplication of large genomic regions or de novo, from previously non-coding DNA. Either way, creation of a novel transcript is an essential early step during novel gene emergence. Most studies on the gain-and-loss dynamics of novel genes so far have compared genomes between species, constraining analyses to genes that have remained fixed over long time scales. However, the importance of novel genes for rapid adaptation among populations has recently been shown. Therefore, since little is known about the evolutionary dynamics of transcripts across natural populations, we here study transcriptomes from several tissues and nine geographically distinct populations of an ecological model species, the three-spined stickleback. Our findings suggest that novel genes typically start out as transcripts with low expression and high tissue specificity. Early expression regulation appears to be mediated by gene-body methylation. Although most new and narrowly expressed genes are rapidly lost, those that survive and subsequently spread through populations tend to gain broader and higher expression levels. The properties of the encoded proteins, such as disorder and aggregation propensity, hardly change. Correspondingly, young novel genes are not preferentially under positive selection but older novel genes more often overlap with F_ST outlier regions. Taken together, expression of the surviving novel genes is rapidly regulated, probably via epigenetic mechanisms, while structural properties of encoded proteins are non-debilitating and might only change much later.

Dissecting Fertility Functions of Drosophila Y Chromosome Genes with CRISPR

Gene-poor, repeat-rich regions of the genome are poorly understood and have been understudied due to technical challenges and the misconception that they are degenerating "junk." Yet multiple lines of evidence indicate these regions may be an important source of variation that could drive adaptation and species divergence, particularly through regulation of fertility. The ∼40 Mb Y chromosome of Drosophila melanoga st er contains only 16 known protein-coding genes, and is highly repetitive and entirely heterochromatic. Most of the genes originated from duplication of autosomal genes and have reduced nonsynonymous substitution rates, suggesting functional constraint. We devised a genetic strategy for recovering and retaining stocks with sterile Y-linked mutations and combined it with CRISPR to create mutants with deletions that disrupt three Y-linked genes. Two genes, PRY and FDY, had no previously identified functions. We found that PRY mutant males are subfertile, but FDY mutant males had no detectable fertility defects. FDY, the newest known gene on the Y chromosome, may have fertility effects that are conditional or too subtle to detect. The third gene, CCY, had been predicted but never formally shown to be required for male fertility. CRISPR targeting and RNA interference of CCY caused male sterility. Surprisingly, however, our CCY mutants were sterile even in the presence of an extra wild-type Y chromosome, suggesting that perturbation of the Y chromosome can lead to dominant sterility. Our approach provides an important step toward understanding the complex functions of the Y chromosome and parsing which functions are accomplished by genes vs. repeat elements.

Sex-Ratio Meiotic Drive Shapes the Evolution of the Y Chromosome in Drosophila simulans

The recent emergence and spread of X-linked segregation distorters-called "Paris" system-in the worldwide species Drosophila simulans has elicited the selection of drive-resistant Y chromosomes. Here, we investigate the evolutionary history of 386 Y chromosomes originating from 29 population samples collected over a period of 20 years, showing a wide continuum of phenotypes when tested against the Paris distorters, from high sensitivity to complete resistance (males sire ∼95% to ∼40% female progeny). Analyzing around 13 kb of Y-linked gene sequences in a representative subset of nine Y chromosomes, we identified only three polymorphic sites resulting in three haplotypes. Remarkably, one of the haplotypes is associated with resistance. This haplotype is fixed in all samples from Sub-Saharan Africa, the region of origin of the drivers. Exceptionally, with the spread of the drivers in Egypt and Morocco, we were able to record the replacement of the sensitive lineage by the resistant haplotype in real time, within only a few years. In addition, we performed in situ hybridization, using satellite DNA probes, on a subset of 21 Y chromosomes from six locations. In contrast to the low molecular polymorphism, this revealed extensive structural variation suggestive of rapid evolution, either neutral or adaptive. Moreover, our results show that intragenomic conflicts can drive astonishingly rapid replacement of Y chromosomes and suggest that the emergence of Paris segregation distorters in East Africa occurred less than half a century ago.

Signature of sexual conflict is actually conflict resolved

There has been substantial interest of late in using population genetic methods to study sexual conflict, where an allele increases the fitness of one sex at some cost to the other (Mank, 2017). Population genomic scans for sexual conflict offer an important advance given the difficulties of identifying antagonistic alleles from more traditional methods, and could greatly increase our understanding of the extent and loci of sexual conflict. This is particularly true for studies in natural populations, for which obtaining accurate fitness measurements for each sex can be challenging. In this issue of Molecular Ecology, Bissegger, Laurentino, Roesti, and Berner (2019) present a cautionary tale about how to interpret these population genomic data.

Identification of Y-chromosome scaffolds of the Queensland fruit fly reveals a duplicated gyf gene paralogue common to many Bactrocera pest species

Bactrocera tryoni (Queensland fruit fly) are polyphagous horticultural pests of eastern Australia. Heterogametic males contain a sex-determining Y-chromosome thought to be gene poor and repetitive. Here, we report 39 Y-chromosome scaffolds (~700 kb) from B. tryoni identified using genotype-by-sequencing data and whole-genome resequencing. Male diagnostic PCR assays validated eight Y-scaffolds, and one (Btry4096) contained a novel gene with five exons that encode a predicted 575 amino acid protein. The Y-gene, referred to as typo-gyf, is a truncated Y-chromosome paralogue of X-chromosome gene gyf (1773 aa). The Y-chromosome contained ~41 copies of typo-gyf, and expression occurred in male flies and embryos. Analysis of 13 tephritid transcriptomes confirmed typo-gyf expression in six additional Bactrocera species, including Bactrocera latifrons, Bactrocera dorsalis and Bactrocera zonata. Molecular dating estimated typo-gyf evolved within the past 8.02 million years (95% highest posterior density 10.56-5.52 million years), after the split with Bactrocera oleae. Phylogenetic analysis also highlighted complex evolutionary histories among several Bactrocera species, as discordant nuclear (116 genes) and mitochondrial (13 genes) topologies were observed. B. tryoni Y-sequences may provide useful sites for future transgene insertions, and typo-gyf could act as a Y-chromosome diagnostic marker for many Bactrocera species, although its function is unknown.

Identification and expression profiling of MSY genes of yak for bull fertility

Yak (Bos grunniens) is a unique bovine species and considered as lifeline of highlanders. The male subfertility in yak is a matter of concern that causes huge economic loses. The spermatogenesis and male reproduction machinery are critically governed by Y-linked genes which tend to acquire necessary information in the course of evolution. The Y-linked fertility genes are present in multiple copies with testis-limited expression. To understand this novel complexity, 12 male-specific region of Y chromosome (MSY) genes have been studied in the yak. Targeted genes are amplified in male and female genomic DNA and confirmed the male derived specificity. Moreover, testis and sperm-specific expressions of MSY genes are distinct among different tissues. The quantitative polymerase chain reaction results validate the expression pattern of these genes in various tissues with predominant expression intestis and sperm. The sequencing of resultant yak MSY genes gives significant result and shows similarity with cattle (Bos indicus), but few nucleotide mismatches define the proposition of infertile male in the F₁ hybrid of cattle and yak. The identified MSY genes can be used to establish male-specific characteristics and to differentiate male and female yak genotypically. Further, these genes may act as valuable resources to understand the capacity of spermatogenesis, embryogenesis, cellular growth, azoospermia and malesubfertility in the yak.

Genetic redundancy fuels polygenic adaptation in Drosophila

The genetic architecture of adaptive traits is of key importance to predict evolutionary responses. Most adaptive traits are polygenic-i.e., result from selection on a large number of genetic loci-but most molecularly characterized traits have a simple genetic basis. This discrepancy is best explained by the difficulty in detecting small allele frequency changes (AFCs) across many contributing loci. To resolve this, we use laboratory natural selection to detect signatures for selective sweeps and polygenic adaptation. We exposed 10 replicates of a Drosophila simulans population to a new temperature regime and uncovered a polygenic architecture of an adaptive trait with high genetic redundancy among beneficial alleles. We observed convergent responses for several phenotypes-e.g., fitness, metabolic rate, and fat content-and a strong polygenic response (99 selected alleles; mean s = 0.059). However, each of these selected alleles increased in frequency only in a subset of the evolving replicates. We discerned different evolutionary paradigms based on the heterogeneous genomic patterns among replicates. Redundancy and quantitative trait (QT) paradigms fitted the experimental data better than simulations assuming independent selective sweeps. Our results show that natural D. simulans populations harbor a vast reservoir of adaptive variation facilitating rapid evolutionary responses using multiple alternative genetic pathways converging at a new phenotypic optimum. This key property of beneficial alleles requires the modification of testing strategies in natural populations beyond the search for convergence on the molecular level.

Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome

Heterochromatic regions of the genome are repeat-rich and poor in protein coding genes, and are therefore underrepresented in even the best genome assemblies. One of the most difficult regions of the genome to assemble are sex-limited chromosomes. The Drosophila melanogaster Y chromosome is entirely heterochromatic, yet has wide-ranging effects on male fertility, fitness, and genome-wide gene expression. The genetic basis of this phenotypic variation is difficult to study, in part because we do not know the detailed organization of the Y chromosome. To study Y chromosome organization in D. melanogaster, we develop an assembly strategy involving the in silico enrichment of heterochromatic long single-molecule reads and use these reads to create targeted de novo assemblies of heterochromatic sequences. We assigned contigs to the Y chromosome using Illumina reads to identify male-specific sequences. Our pipeline extends the D. melanogaster reference genome by 11.9 Mb, closes 43.8% of the gaps, and improves overall contiguity. The addition of 10.6 MB of Y-linked sequence permitted us to study the organization of repeats and genes along the Y chromosome. We detected a high rate of duplication to the pericentric regions of the Y chromosome from other regions in the genome. Most of these duplicated genes exist in multiple copies. We detail the evolutionary history of one sex-linked gene family, crystal-Stellate While the Y chromosome does not undergo crossing over, we observed high gene conversion rates within and between members of the crystal-Stellate gene family, Su(Ste), and PCKR, compared to genome-wide estimates. Our results suggest that gene conversion and gene duplication play an important role in the evolution of Y-linked genes.