Birth of a new gene on the Y chromosome of Drosophila melanogaster

Co-transcriptional splicing facilitates transcription of gigantic genes

Although introns are typically tens to thousands of nucleotides, there are notable exceptions. In flies as well as humans, a small number of genes contain introns that are more than 1000 times larger than typical introns, exceeding hundreds of kilobases (kb) to megabases (Mb). It remains unknown why gigantic introns exist and how cells overcome the challenges associated with their transcription and RNA processing. The Drosophila Y chromosome contains some of the largest genes identified to date: multiple genes exceed 4Mb, with introns accounting for over 99% of the gene span. Here we demonstrate that co-transcriptional splicing of these gigantic Y-linked genes is important to ensure successful transcription: perturbation of splicing led to the attenuation of transcription, leading to a failure to produce mature mRNA. Cytologically, defective splicing of the Y-linked gigantic genes resulted in disorganization of transcripts within the nucleus suggestive of entanglement of transcripts, likely resulting from unspliced long RNAs. We propose that co-transcriptional splicing maintains the length of nascent transcripts of gigantic genes under a critical threshold, preventing their entanglement and ensuring proper gene expression. Our study reveals a novel biological significance of co-transcriptional splicing.

Cis-regulatory polymorphism at fiz ecdysone oxidase contributes to polygenic evolutionary response to malnutrition in Drosophila

We investigate the contribution of a candidate gene, fiz (fezzik), to complex polygenic adaptation to juvenile malnutrition in Drosophila melanogaster. Experimental populations maintained for >250 generations of experimental evolution to a nutritionally poor larval diet (Selected populations) evolved several-fold lower fiz expression compared to unselected Control populations. Here we show that this divergence in fiz expression is mediated by a cis-regulatory polymorphism. This polymorphism, originally sampled from a natural population in Switzerland, is distinct from a second cis-regulatory SNP previously identified in non-African D. melanogaster populations, implying that two independent cis-regulatory variants promoting high fiz expression segregate in non-African populations. Enzymatic analyses of Fiz protein expressed in E. coli demonstrate that it has ecdysone oxidase activity acting on both ecdysone and 20-hydroxyecdysone. Four of five fiz paralogs annotated to ecdysteroid metabolism also show reduced expression in Selected larvae, implying that malnutrition-driven selection favored general downregulation of ecdysone oxidases. Finally, as an independent test of the role of fiz in poor diet adaptation, we show that fiz knockdown by RNAi results in faster larval growth on the poor diet, but at the cost of greatly reduced survival. These results imply that downregulation of fiz in Selected populations was favored by selection on the nutritionally poor diet because of its role in suppressing growth in response to nutrient shortage. However, they suggest that fiz downregulation is only adaptive in combination with other changes evolved by Selected populations, which ensure that the organism can sustain the faster growth promoted by fiz downregulation.

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

Gene-poor Y-chromosomes substantially impact male trait heritabilities and may help shape sexually dimorphic evolution

How natural selection facilitates sexually dimorphic evolution despite a shared genome is unclear. The patrilineal inheritance of Y-chromosomes makes them an appealing solution. However, they have largely been dismissed due to their gene-poor, heterochromatic nature and because the additive genetic variation necessary for adaptive evolution is theoretically difficult to maintain. Further, previous empirical work has revealed mostly Y-linked sign epistatic variance segregating within populations, which can often impede adaptive evolution. To assess the evolutionary impact of Y-linked variation, we established replicate populations in Drosophila simulans containing multiple Y-chromosomes (YN populations) or a single Y-chromosome variant (Y1 populations) drawn from a single population. We estimated male and female heritabilities for several traits known to be influenced by Y-chromosomes, including the number of sternopleural bristles, abdominal bristles, sex comb teeth, and tibia length. A decrease in YN heritabilities compared with Y1 would be consistent with Y-chromosome variation being sign epistatic. A decrease in Y1 heritabilities would be consistent with Y-chromosome variation being additive, though additive-by-additive epistatic variation cannot be entirely dismissed. Female heritability estimates served as controls and were not expected to differ. We found male Y1 populations exhibited lower heritabilities for all traits except tibia length; consistent with Y-linked additivity (on average YN trait heritabilities were 25% greater than Y1). Female estimates showed no difference. These data suggest Y-chromosomes should play an important role in male trait evolution and may even influence sexually dimorphic evolution by shaping traits shared by both sexes.

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.

The regulation and potential functions of intronic satellite DNA

Satellite DNAs are arrays of tandem repeats found in the eukaryotic genome. They are mainly found in pericentromeric heterochromatin and have been believed to be mostly inert, leading satellite DNAs to be erroneously regarded as junk. Recent studies have started to elucidate the function of satellite DNA, yet little is known about the peculiar case where satellite DNA is found within the introns of protein coding genes, resulting in incredibly large introns, a phenomenon termed intron gigantism. Studies in Drosophila demonstrated that satellite DNA-containing introns are transcribed with the gene and require specialized mechanisms to overcome the burdens imposed by the extremely long stretches of repetitive DNA. Whether intron gigantism confers any benefit or serves any functional purpose for cells and/or organisms remains elusive. Here we review our current understanding of intron gigantism: where it is found, the challenges it imposes, how it is regulated and what purpose it may serve.

Drosophila as a Model System for Studying of the Evolution and Functional Specialization of the Y Chromosome

The Y chromosome is one of the sex chromosomes found in males of animals of different taxa, including insects and mammals. Among all chromosomes, the Y chromosome is characterized by a unique chromatin landscape undergoing dynamic evolutionary change. Being entirely heterochromatic, the Y chromosome as a rule preserves few functional genes, but is enriched in tandem repeats and transposons. Due to difficulties in the assembly of the highly repetitive Y chromosome sequence, deep analyses of Y chromosome evolution, structure, and functions are limited to a few species, one of them being Drosophila melanogaster. Despite Y chromosomes exhibiting high structural divergence between even closely related species, Y-linked genes have evolved convergently and are mainly associated with spermatogenesis-related activities. This indicates that male-specific selection is a dominant force shaping evolution of Y chromosomes across species. This review presents our analysis of current knowledge concerning Y chromosome functions, focusing on recent findings in Drosophila. Here we dissect the experimental and bioinformatics data about the Y chromosome accumulated to date in Drosophila species, providing comparative analysis with mammals, and discussing the relevance of our analysis to a wide range of eukaryotic organisms, including humans.

Novel patterns of expression and recruitment of new genes on the t-haplotype, a mouse selfish chromosome

The t-haplotype of mice is a classical model for autosomal transmission distortion. A largely non-recombining variant of the proximal region of chromosome 17, it is transmitted to more than 90% of the progeny of heterozygous males through the disabling of sperm carrying a standard chromosome. While extensive genetic and functional work has shed light on individual genes involved in drive, much less is known about the evolution and function of the rest of its hundreds of genes. Here, we characterize the sequence and expression of dozens of t-specific transcripts and of their chromosome 17 homologues. Many genes showed reduced expression of the t-allele, but an equal number of genes showed increased expression of their t-copy, consistent with increased activity or a newly evolved function. Genes on the t-haplotype had a significantly higher non-synonymous substitution rate than their homologues on the standard chromosome, with several genes harbouring dN/dS ratios above 1. Finally, the t-haplotype has acquired at least two genes from other chromosomes, which show high and tissue-specific expression. These results provide a first overview of the gene content of this selfish element, and support a more dynamic evolutionary scenario than expected of a large genomic region with almost no recombination.

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.

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.

piRNA-mediated gene regulation and adaptation to sex-specific transposon expression in D. melanogaster male germline

Small noncoding piRNAs act as sequence-specific guides to repress complementary targets in Metazoa. Prior studies in Drosophila ovaries have demonstrated the function of the piRNA pathway in transposon silencing and therefore genome defense. However, the ability of the piRNA program to respond to different transposon landscapes and the role of piRNAs in regulating host gene expression remain poorly understood. Here, we comprehensively analyzed piRNA expression and defined the repertoire of their targets in Drosophila melanogaster testes. Comparison of piRNA programs between sexes revealed sexual dimorphism in piRNA programs that parallel sex-specific transposon expression. Using a novel bioinformatic pipeline, we identified new piRNA clusters and established complex satellites as dual-strand piRNA clusters. While sharing most piRNA clusters, the two sexes employ them differentially to combat the sex-specific transposon landscape. We found two piRNA clusters that produce piRNAs antisense to four host genes in testis, including CG12717/pirate, a SUMO protease gene. piRNAs encoded on the Y chromosome silence pirate, but not its paralog, to exert sex- and paralog-specific gene regulation. Interestingly, pirate is targeted by endogenous siRNAs in a sibling species, Drosophila mauritiana, suggesting distinct but related silencing strategies invented in recent evolution to regulate a conserved protein-coding gene.

Evolution of a New Function by Fusion between Phage DNA and a Bacterial Gene

Mobile genetic elements, such as plasmids, phages, and transposons, are important sources for evolution of novel functions. In this study, we performed a large-scale screening of metagenomic phage libraries for their ability to suppress temperature-sensitivity in Salmonella enterica serovar Typhimurium strain LT2 mutants to examine how phage DNA could confer evolutionary novelty to bacteria. We identified an insert encoding 23 amino acids from a phage that when fused with a bacterial DNA-binding repressor protein (LacI) resulted in the formation of a chimeric protein that localized to the outer membrane. This relocalization of the chimeric protein resulted in increased membrane vesicle formation and an associated suppression of the temperature sensitivity of the bacterium. Both the host LacI protein and the extracellular 23-amino acid stretch are necessary for the generation of the novel phenotype. Furthermore, mutational analysis of the chimeric protein showed that although the native repressor function of the LacI protein is maintained in this chimeric structure, it is not necessary for the new function. Thus, our study demonstrates how a gene fusion between foreign DNA and bacterial DNA can generate novelty without compromising the native function of a given gene.

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.

Drosophila melanogaster Y Chromosome Genes Affect Male Sensitivity to Microbial Infections

The genders of Drosophila melanogaster vary in their sensitivities to microbial pathogens. While many of the immunity-related genes are located on the X chromosome, the polymorphisms within the Y chromosome were also shown to affect the immunity of flies. In this study, we investigated the necessity of individual genes on the Y chromosome (Y-genes) for male sensitivity to microbes. We identified several Y-genes whose genetic inactivation either increases or decreases the sensitivity of males to gastrointestinal infections with fungal Saccharomyces cerevisiae and bacterial Serratia liquefaciens. Specifically, the loss of function mutations in fly kl-5 and Ppr-Y Y-genes lead to increased and decreased sensitivity of males to fungal challenge, respectively, compared to female sensitivity. In contrast, mutations in Drosophila Pp1-Y1, kl-5, kl-3, Ppr-Y, CCY, and FDY Y-genes lead to increased sensitivity of males to bacterial infection, compared to females. Moreover, while these Y-genes are necessary, the Y chromosome is not sufficient for the sensitivity of males to microbes, since the sensitivity of XXY females to fungal and bacterial challenges was not different from the sensitivity of wild-type female flies, compared to males. This study assigns a new immunity-related function to numerous Y-genes in D. melanogaster.

Transcribed sex-specific markers on the Y chromosome of the oriental fruit fly, Bactrocera dorsalis

BACKGROUND

The Oriental fruit fly, Bactrocera dorsalis, is a highly polyphagous invasive species with a high reproductive potential. In many tropical and subtropical parts of the world it ranks as one of the major pests of fruits and vegetables. Due to its economic importance, genetic, cytogenetic, genomic and biotechnological approaches have been applied to understand its biology and to implement the Sterile Insect Technique, currently a part of area-wide control programmes against this fly. Its chromosome complement includes five pairs of autosomes and the sex chromosomes. The X and Y sex chromosomes are heteromorphic and the highly heterochromatic and degenerate Y harbours the male factor BdMoY. The characterization of the Y chromosome in this fly apart from elucidating its role as primary sex determination system, it is also of crucial importance to understand its role in male biology. The repetitive nature of the Y chromosome makes it challenging to sequence and characterise.

RESULTS

Using Representational Difference Analysis, fluorescent in situ hybridisation on mitotic chromosomes and in silico genome resources, we show that the B. dorsalis Y chromosome harbours transcribed sequences of gyf, (typo-gyf) a homologue of the Drosophila melanogaster Gigyf gene, and of a non-LTR retrotransposon R1. Similar sequences are also transcribed on the X chromosome. Paralogues of the Gigyf gene are also present on the Y and X chromosomes of the related species B. tryoni. Another identified Y-specific repetitive sequence linked to BdMoY appears to be specific to B. dorsalis.

CONCLUSIONS

Our random scan of the Y chromosome provides a broad picture of its general composition and represents a starting point for further applicative and evolutionary studies. The identified repetitive sequences can provide a useful Y-marking system for molecular karyotyping of single embryos. Having a robust diagnostic marker associated with BdMoY will facilitate studies on how BdMoY regulates the male sex determination cascade during the embryonic sex-determination window. The Y chromosome, despite its high degeneracy and heterochromatic nature, harbours transcribed sequences of typo-gyf that may maintain their important function in post-transcriptional mRNA regulation. That transcribed paralogous copies of Gigyf are present also on the X and that this genomic distribution is maintained also in B. tryoni raises questions on the evolution of sex chromosomes in Bactrocera and other tephritids.

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.

The Y Chromosome as a Battleground for Intragenomic Conflict

Y chromosomes are typically viewed as genetic wastelands with few intact genes. Recent genomic analyses in Drosophila, however, show that gene gain is prominent on young Y chromosomes. Meiosis- and RNAi-related genes often coamplify on recently formed X and Y chromosomes, are testis-expressed, and produce antisense transcripts and short RNAs. RNAi pathways are also involved in suppressing sex ratio drive in Drosophila. These observations paint a dynamic picture of sex chromosome differentiation, suggesting that rapidly evolving genomic battles over segregation are rampant on young sex chromosomes and utilize RNAi to defend the genome against selfish elements that manipulate fair meiosis. Recurrent sex chromosome drive can have profound ecological, evolutionary, and cellular impacts and account for unique features of sex chromosomes.

The Role of Y Chromosome Genes in Male Fertility in Drosophila melanogaster

The Y chromosome of Drosophila melanogaster is pivotal for male fertility. Yet, only 16 protein-coding genes reside on this chromosome. The Y chromosome is comprised primarily of heterochromatic sequences, including DNA repeats and satellite DNA, and most of the Y chromosome is still missing from the genome sequence. Furthermore, the functions of the majority of genes on the Y chromosome remain elusive. Through multiple genetic strategies, six distinct segments on the Y chromosome have been identified as "male fertility factors," and candidate gene sequences corresponding to each of these loci have been ascribed. In one case, kl-3, a specific protein coding sequence for a fertility factor has been confirmed molecularly. Here, we employed CRISPR/Cas9 to generate mutations, and RNAi, to interrogate the requirements of protein coding sequences on the Y chromosome for male fertility. We show that CRISPR/Cas9-mediated editing of kl-2 and kl-5 causes male sterility, supporting the model that these gene sequences correspond to the cognate fertility factors. We show that another gene, CCY, also functions in male fertility and may be the ks-2 fertility factor. We demonstrate that editing of kl-2, kl-3, and kl-5, and RNAi knockdown of CCY, disrupts nuclear elongation, and leads to defects in sperm individualization, including impairments in the individualization complex (IC) and synchronization. However, CRISPR/Cas9 mediated knockout of some genes on the Y chromosome, such as FDY, Ppr-Y, and Pp1-Y2 do not cause sterility, indicating that not all Y chromosome genes are essential for male fertility.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly

BACKGROUND

The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use.

RESULTS

Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR.

CONCLUSIONS

The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

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.

Unpredictable Effects of the Genetic Background of Transgenic Lines in Physiological Quantitative Traits

Physiology, fitness and disease phenotypes are complex traits exhibiting continuous variation in natural populations. To understand complex trait gene functions transgenic lines of undefined genetic background are often combined to assess quantitative phenotypes ignoring the impact of genetic polymorphisms. Here, we used inbred wild-type strains of the Drosophila Genetics Reference Panel to assess the phenotypic variation of six physiological and fitness traits, namely, female fecundity, survival and intestinal mitosis upon oral infection, defecation rate and fecal pH upon oral infection, and terminal tracheal cell branching in hypoxia. We found continuous variation in the approximately 150 strains tested for each trait, with extreme values differing by more than four standard deviations for all traits. In addition, we assessed the effects of commonly used Drosophila UAS-RNAi transgenic strains and their backcrossed isogenized counterparts, in the same traits plus baseline intestinal mitosis and tracheal branching in normoxia, in heterozygous conditions, when only half of the genetic background was different among strains. We tested 20 non-isogenic strains (10 KK and 10 GD) from the Vienna Drosophila Resource Center and their isogenized counterparts without Gal4 induction. Survival upon infection and female fecundity exhibited differences in 50% and 40% of the tested isogenic vs. non-isogenic pairs, respectively, whereas all other traits were affected in only 10–25% of the cases. When 11 isogenic and their corresponding non-isogenic UAS-RNAi lines were expressed ubiquitously with Gal4, 4 isogenic vs. non-isogenic pairs exhibited differences in survival to infection. Furthermore, when a single UAS-RNAi line was crossed with the same Gal4 transgene inserted in different genetic backgrounds, the quantitative variations observed were unpredictable on the basis of pure line performance. Thus, irrespective of the trait of interest, the genetic background of commonly used transgenic strains needs to be considered carefully during experimentation.

Sex and the TEs: transposable elements in sexual development and function in animals

Transposable elements are endogenous DNA sequences able to integrate into and multiply within genomes. They constitute a major source of genetic innovations, as they can not only rearrange genomes but also spread ready-to-use regulatory sequences able to modify host gene expression, and even can give birth to new host genes. As their evolutionary success depends on their vertical transmission, transposable elements are intrinsically linked to reproduction. In organisms with sexual reproduction, this implies that transposable elements have to manifest their transpositional activity in germ cells or their progenitors. The control of sexual development and function can be very versatile, and several studies have demonstrated the implication of transposable elements in the evolution of sex. In this review, we report the functional and evolutionary relationships between transposable elements and sexual reproduction in animals. In particular, we highlight how transposable elements can influence expression of sexual development genes, and how, reciprocally, they are tightly controlled in gonads. We also review how transposable elements contribute to the organization, expression and evolution of sexual development genes and sex chromosomes. This underscores the intricate co-evolution between host functions and transposable elements, which regularly shift from a parasitic to a domesticated status useful to the host.

A New Portrait of Constitutive Heterochromatin: Lessons from Drosophila melanogaster

Constitutive heterochromatin represents a significant portion of eukaryotic genomes, but its functions still need to be elucidated. Even in the most updated genetics and molecular biology textbooks, constitutive heterochromatin is portrayed mainly as the 'silent' component of eukaryotic genomes. However, there may be more complexity to the relationship between heterochromatin and gene expression. In the fruit fly Drosophila melanogaster, a model for heterochromatin studies, about one-third of the genome is heterochromatic and is concentrated in the centric, pericentric, and telomeric regions of the chromosomes. Recent findings indicate that hundreds of D. melanogaster genes can 'live and work' properly within constitutive heterochromatin. The genomic size of these genes is generally larger than that of euchromatic genes and together they account for a significant fraction of the entire constitutive heterochromatin. Thus, this peculiar genome component in spite its ability to induce silencing, has in fact the means for being quite dynamic. A major scope of this review is to revisit the 'dogma of silent heterochromatin'.

Satellite DNA-containing gigantic introns in a unique gene expression program during Drosophila spermatogenesis

Intron gigantism, where genes contain megabase-sized introns, is observed across species, yet little is known about its purpose or regulation. Here we identify a unique gene expression program utilized for the proper expression of genes with intron gigantism. We find that two Drosophila genes with intron gigantism, kl-3 and kl-5, are transcribed in a spatiotemporal manner over the course of spermatocyte differentiation, which spans ~90 hours. The introns of these genes contain megabases of simple satellite DNA repeats that comprise over 99% of the gene loci, and these satellite-DNA containing introns are transcribed. We identify two RNA-binding proteins that specifically localize to kl-3 and kl-5 transcripts and are needed for the successful transcription or processing of these genes. We propose that genes with intron gigantism require a unique gene expression program, which may serve as a platform to regulate gene expression during cellular differentiation.

Introns are non-coding elements of eukaryotic genes, often containing important regulatory sequences. Curiously, some genes contain introns so large that more than 99% of the gene locus is non-coding. One of the best-studied large genes, Dystrophin, a causative gene for Duchenne Muscular Dystrophy, spans 2.2Mb, only 11kb of which is coding. This phenomenon, ‘intron gigantism’, is observed across species, yet little is known about its purpose or regulation. Here we identify a unique gene expression program utilized for the proper expression of genes with intron gigantism using Drosophila spermatogenic genes as a model system. We show that the gigantic introns of these genes are transcribed in line with the exons, likely as a single transcript. We identify two RNA-binding proteins that specifically localize to the site of transcription and are needed for the successful transcription or processing of these genes. We propose that genes with intron gigantism require a unique gene expression program, which may serve as a platform to regulate gene expression during cellular differentiation.

Genetic and Molecular Analysis of Essential Genes in Centromeric Heterochromatin of the Left Arm of Chromosome 3 in Drosophila melanogaster

A large portion of the Drosophila melanogaster genome is contained within heterochromatic regions of chromosomes, predominantly at centromeres and telomeres. The remaining euchromatic portions of the genome have been extensively characterized with respect to gene organization, function and regulation. However, it has been difficult to derive similar data for sequences within centromeric (centric) heterochromatin because these regions have not been as amenable to analysis by standard genetic and molecular tools. Here we present an updated genetic and molecular analysis of chromosome 3L centric heterochromatin (3L Het). We have generated and characterized a number of new, overlapping deficiencies (Dfs) which remove regions of 3L Het. These Dfs were critically important reagents in our subsequent genetic analysis for the isolation and characterization of lethal point mutations in the region. The assignment of these mutations to genetically-defined essential loci was followed by matching them to gene models derived from genome sequence data: this was done by using molecular mapping plus sequence analysis of mutant alleles, thereby aligning genetic and physical maps of the region. We also identified putative essential gene sequences in 3L Het by using RNA interference to target candidate gene sequences. We report that at least 25, or just under 2/3 of loci in 3L Het, are essential for viability and/or fertility. This work contributes to the functional annotation of centric heterochromatin in Drosophila, and the genetic and molecular tools generated should help to provide important insights into the organization and functions of gene sequences in 3L Het.

Novel Y Chromosome Retrocopies in Canids Revealed through a Genome-Wide Association Study for Sex

The lack of an annotated reference sequence for the canine Y chromosome has limited evolutionary studies, as well as our understanding of the role of Y-linked sequences in phenotypes with a sex bias. In genome-wide association studies (GWASs), we observed spurious associations with autosomal SNPs when sex was unbalanced in case-control cohorts and hypothesized that a subset of SNPs mapped to autosomes are in fact sex-linked. Using the Illumina 230K CanineHD array in a GWAS for sex, we identified SNPs that amplify in both sexes but possess significant allele frequency differences between males and females. We found 48 SNPs mapping to 14 regions of eight autosomes and the X chromosome that are Y-linked, appearing heterozygous in males and monomorphic in females. Within these 14 regions are eight genes: three autosomal and five X-linked. We investigated the autosomal genes (MITF, PPP2CB, and WNK1) and determined that the SNPs are diverged nucleotides in retrocopies that have transposed to the Y chromosome. MITFY and WNK1Y are expressed and appeared recently in the Canidae lineage, whereas PPP2CBY represents a much older insertion with no evidence of expression in the dog. This work reveals novel canid Y chromosome sequences and provides evidence for gene transposition to the Y from autosomes and the X.

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.

Sex-specific adaptation and genomic responses to Y chromosome presence in female reproductive and neural tissues

Y chromosomes typically harbour a small number of genes and an abundance of repetitive sequences. In Drosophila, the Y chromosome comprises multimegabase long segments of repetitive DNA and a handful of protein-coding genes. In mammals, the Y chromosome also harbours a disproportionally high abundance of repeats. Here, we built on a Drosophila melanogaster model in which the Y chromosome is decoupled from sexual determination. Genotypes were genetically identical for the autosomes, X chromosome, and mitochondria, but differ by the presence or dose of the Y chromosome. Addition of an extra Y chromosome had limited impact in males. However, the presence of a Y chromosome in females induced a disproportionate response in genes expressed in the ovaries as well as genes encoded by the mitochondrial genome. Furthermore, the data revealed significant consequences of Y chromosome presence in larvae neuronal tissue. This included the repression of genes implicated in reproductive behaviour, courtship, mating and synaptic function. Our findings exhibit the Y chromosome as a hotspot for sex-specific adaptation. They suggest roles for natural selection on Y-linked genetic elements exerting impact on sex-specific tissues as well as somatic tissues shared by males and females.

Regulation of the X Chromosome in the Germline and Soma of Drosophila melanogaster Males

During the evolution of heteromorphic sex chromosomes, the sex-specific Y chromosome degenerates, while the X chromosome evolves new mechanisms of regulation. Using bioinformatic and experimental approaches, we investigate the expression of the X chromosome in Drosophila melanogaster. We observe nearly complete X chromosome dosage compensation in male somatic tissues, but not in testis. The X chromosome contains disproportionately fewer genes with high expression in testis than the autosomes, even after accounting for the lack of dosage compensation, which suggests that another mechanism suppresses their expression in the male germline. This is consistent with studies of reporter genes and transposed genes, which find that the same gene has higher expression when autosomal than when X-linked. Using a new reporter gene that is expressed in both testis and somatic tissues, we find that the suppression of X-linked gene expression is limited to genes with high expression in testis and that the extent of the suppression is positively correlated with expression level.

Shared and Species-Specific Patterns of Nascent Y Chromosome Evolution in Two Guppy Species

Sex chromosomes form once recombination is halted around the sex-determining locus between a homologous pair of chromosomes, resulting in a male-limited Y chromosome. We recently characterized the nascent sex chromosome system in the Trinidadian guppy (Poeciliareticulata). The guppy Y is one of the youngest animal sex chromosomes yet identified, and therefore offers a unique window into the early evolutionary forces shaping sex chromosome formation, particularly the rate of accumulation of repetitive elements and Y-specific sequence. We used comparisons between male and female genomes in P. reticulata and its sister species, Endler’s guppy (P. wingei), which share an ancestral sex chromosome, to identify male-specific sequences and to characterize the degree of differentiation between the X and Y chromosomes. We identified male-specific sequence shared between P. reticulata and P. wingei consistent with a small ancestral non-recombining region. Our assembly of this Y-specific sequence shows substantial homology to the X chromosome, and appears to be significantly enriched for genes implicated in pigmentation. We also found two plausible candidates that may be involved in sex determination. Furthermore, we found that the P. wingei Y chromosome exhibits a greater signature of repetitive element accumulation than the P. reticulata Y chromosome. This suggests that Y chromosome divergence does not necessarily correlate with the time since recombination suppression. Overall, our results reveal the early stages of Y chromosome divergence in the guppy.

dCas9-targeted locus-specific protein isolation method identifies histone gene regulators

Eukaryotic gene regulation is a complex process, often coordinated by the action of tens to hundreds of proteins. Although previous biochemical studies have identified many components of the basal machinery and various ancillary factors involved in gene regulation, numerous gene-specific regulators remain undiscovered. To comprehensively survey the proteome directing gene expression at a specific genomic locus of interest, we developed an in vitro nuclease-deficient Cas9 (dCas9)-targeted chromatin-based purification strategy, called "CLASP" (Cas9 locus-associated proteome), to identify and functionally test associated gene-regulatory factors. Our CLASP method, coupled to mass spectrometry and functional screens, can be efficiently adapted for isolating associated regulatory factors in an unbiased manner targeting multiple genomic loci across different cell types. Here, we applied our method to isolate the Drosophila melanogaster histone cluster in S2 cells to identify several factors including Vig and Vig2, two proteins that bind and regulate core histone H2A and H3 mRNA via interaction with their 3' UTRs.

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

Using a powerful method that uses inexpensive short reads to detect Y-linked transfers, we show that gene traffic onto the Drosophila Y chromosome is 10 times more frequent than previously thought and includes the first Y-linked retrocopies discovered in these taxa. All 25 identified Y-linked gene transfers were relatively young (<1 million years old), although most appear to be pseudogenes because only three of these transfers show signs of purifying selection. Our method provides compelling evidence that the Drosophila Y chromosome is a highly challenging and dynamic genetic environment that is capable of rapidly diverging between species and promises to reveal fundamental insights into Y chromosome evolution across many taxa.

The Y chromosome is a unique genetic environment defined by a lack of recombination and male-limited inheritance. The Drosophila Y chromosome has been gradually acquiring genes from the rest of the genome, with only seven Y-linked genes being gained over the past 63 million years (0.12 gene gains per million years). Using a next-generation sequencing (NGS)-powered genomic scan, we show that gene transfers to the Y chromosome are much more common than previously suspected: at least 25 have arisen across three Drosophila species over the past 5.4 million years (1.67 per million years for each lineage). The gene transfer rate is significantly lower in Drosophila melanogaster than in the Drosophila simulans clade, primarily due to Y-linked retrotranspositions being significantly more common in the latter. Despite all Y-linked gene transfers being evolutionarily recent (<1 million years old), only three showed evidence for purifying selection (ω ≤ 0.14). Thus, although the resulting Y-linked functional gene acquisition rate (0.25 new genes per million years) is double the longer-term estimate, the fate of most new Y-linked genes is defined by rapid degeneration and pseudogenization. Our results show that Y-linked gene traffic, and the molecular mechanisms governing these transfers, can diverge rapidly between species, revealing the Drosophila Y chromosome to be more dynamic than previously appreciated. Our analytical method provides a powerful means to identify Y-linked gene transfers and will help illuminate the evolutionary dynamics of the Y chromosome in Drosophila and other species.

Neo-sex Chromosomes in the Monarch Butterfly, Danaus plexippus

We report the discovery of a neo-sex chromosome in the monarch butterfly, Danaus plexippus, and several of its close relatives. Z-linked scaffolds in the D. plexippus genome assembly were identified via sex-specific differences in Illumina sequencing coverage. Additionally, a majority of the D. plexippus genome assembly was assigned to chromosomes based on counts of one-to-one orthologs relative to the butterfly Melitaea cinxia (with replication using two other lepidopteran species), in which genome scaffolds have been mapped to linkage groups. Sequencing coverage-based assessments of Z linkage combined with homology-based chromosomal assignments provided strong evidence for a Z-autosome fusion in the Danaus lineage, involving the autosome homologous to chromosome 21 in M. cinxia. Coverage analysis also identified three notable assembly errors resulting in chimeric Z-autosome scaffolds. Cytogenetic analysis further revealed a large W chromosome that is partially euchromatic, consistent with being a neo-W chromosome. The discovery of a neo-Z and the provisional assignment of chromosome linkage for >90% of D. plexippus genes lays the foundation for novel insights concerning sex chromosome evolution in this female-heterogametic model species for functional and evolutionary genomics.

Convergent evolution of Y chromosome gene content in flies

Sex-chromosomes have formed repeatedly across Diptera from ordinary autosomes, and X-chromosomes mostly conserve their ancestral genes. Y-chromosomes are characterized by abundant gene-loss and an accumulation of repetitive DNA, yet the nature of the gene repertoire of fly Y-chromosomes is largely unknown. Here we trace gene-content evolution of Y-chromosomes across 22 Diptera species, using a subtraction pipeline that infers Y genes from male and female genome, and transcriptome data. Few genes remain on old Y-chromosomes, but the number of inferred Y-genes varies substantially between species. Young Y-chromosomes still show clear evidence of their autosomal origins, but most genes on old Y-chromosomes are not simply remnants of genes originally present on the proto-sex-chromosome that escaped degeneration, but instead were recruited secondarily from autosomes. Despite almost no overlap in Y-linked gene content in different species with independently formed sex-chromosomes, we find that Y-linked genes have evolved convergent gene functions associated with testis expression. Thus, male-specific selection appears as a dominant force shaping gene-content evolution of Y-chromosomes across fly species.

While X-chromosome gene content tends to be conserved, Y-chromosome evolution is dynamic and difficult to reconstruct. Here, Mahajan and Bachtrog use a subtraction pipeline to identify Y-linked genes in 22 Diptera species, revealing patterns of Y-chromosome gene-content evolution.

Cross-Species Y Chromosome Function Between Malaria Vectors of the Anopheles gambiae Species Complex

Y chromosome function, structure and evolution is poorly understood in many species, including the Anopheles genus of mosquitoes-an emerging model system for studying speciation that also represents the major vectors of malaria. While the Anopheline Y had previously been implicated in male mating behavior, recent data from the Anopheles gambiae complex suggests that, apart from the putative primary sex-determiner, no other genes are conserved on the Y. Studying the functional basis of the evolutionary divergence of the Y chromosome in the gambiae complex is complicated by complete F1 male hybrid sterility. Here, we used an F1 × F0 crossing scheme to overcome a severe bottleneck of male hybrid incompatibilities that enabled us to experimentally purify a genetically labeled A. gambiae Y chromosome in an A. arabiensis background. Whole genome sequencing (WGS) confirmed that the A. gambiae Y retained its original sequence content in the A. arabiensis genomic background. In contrast to comparable experiments in Drosophila, we find that the presence of a heterospecific Y chromosome has no significant effect on the expression of A. arabiensis genes, and transcriptional differences can be explained almost exclusively as a direct consequence of transcripts arising from sequence elements present on the A. gambiae Y chromosome itself. We find that Y hybrids show no obvious fertility defects, and no substantial reduction in male competitiveness. Our results demonstrate that, despite their radically different structure, Y chromosomes of these two species of the gambiae complex that diverged an estimated 1.85 MYA function interchangeably, thus indicating that the Y chromosome does not harbor loci contributing to hybrid incompatibility. Therefore, Y chromosome gene flow between members of the gambiae complex is possible even at their current level of divergence. Importantly, this also suggests that malaria control interventions based on sex-distorting Y drive would be transferable, whether intentionally or contingent, between the major malaria vector species.

The house fly Y Chromosome is young and minimally differentiated from its ancient X Chromosome partner

Canonical ancient sex chromosome pairs consist of a gene rich X (or Z) Chromosome and a male-limited (or female-limited) Y (or W) Chromosome that is gene poor. In contrast to highly differentiated sex chromosomes, nascent sex chromosome pairs are homomorphic or very similar in sequence content. Nascent sex chromosomes can arise if an existing sex chromosome fuses to an autosome or an autosome acquires a new sex-determining locus/allele. Sex chromosomes often differ between closely related species and can even be polymorphic within species, suggesting that nascent sex chromosomes arise frequently over the course of evolution. Previously documented sex chromosome transitions involve changes to both members of the sex chromosome pair (X and Y, or Z and W). The house fly has sex chromosomes that resemble the ancestral fly karyotype that originated ∼100 million yr ago; therefore, the house fly is expected to have X and Y Chromosomes with different gene content. We tested this hypothesis using whole-genome sequencing and transcriptomic data, and we discovered little evidence for genetic differentiation between the X and Y in house fly. We propose that the house fly has retained the ancient X Chromosome, but the ancestral Y was replaced by an X Chromosome carrying a new male determining gene. Our proposed hypothesis provides a mechanism for how one member of a sex chromosome pair can experience evolutionary turnover while the other member remains unaffected.

Multiple gene movements into and out of haploid sex chromosomes

BACKGROUND

Long-term evolution of sex chromosomes is a dynamic process shaped by gene gain and gene loss. Sex chromosome gene traffic has been studied in XY and ZW systems but no detailed analyses have been carried out for haploid phase UV sex chromosomes. Here, we explore sex-specific sequences of seven brown algal species to understand the dynamics of the sex-determining region (SDR) gene content across 100 million years of evolution.

RESULTS

A core set of sex-linked genes is conserved across all the species investigated, but we also identify modifications of both the U and the V SDRs that occurred in a lineage-specific fashion. These modifications involve gene loss, gene gain and relocation of genes from the SDR to autosomes. Evolutionary analyses suggest that the SDR genes are evolving rapidly and that this is due to relaxed purifying selection. Expression analysis indicates that genes that were acquired from the autosomes have been retained in the SDR because they confer a sex-specific role in reproduction. By examining retroposed genes in Saccharina japonica, we demonstrate that UV sex chromosomes have generated a disproportionate number of functional orphan retrogenes compared with autosomes. Movement of genes out of the UV sex chromosome could be a means to compensate for gene loss from the non-recombining region, as has been suggested for Y-derived retrogenes in XY sexual systems.

CONCLUSION

This study provides the first analysis of gene traffic in a haploid UV system and identifies several features of general relevance to the evolution of sex chromosomes.

Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster

Highly repetitive satellite DNA (satDNA) repeats are found in most eukaryotic genomes. SatDNAs are rapidly evolving and have roles in genome stability and chromosome segregation. Their repetitive nature poses a challenge for genome assembly and makes progress on the detailed study of satDNA structure difficult. Here, we use single-molecule sequencing long reads from Pacific Biosciences (PacBio) to determine the detailed structure of all major autosomal complex satDNA loci in Drosophila melanogaster, with a particular focus on the 260-bp and Responder satellites. We determine the optimal de novo assembly methods and parameter combinations required to produce a high-quality assembly of these previously unassembled satDNA loci and validate this assembly using molecular and computational approaches. We determined that the computationally intensive PBcR-BLASR assembly pipeline yielded better assemblies than the faster and more efficient pipelines based on the MHAP hashing algorithm, and it is essential to validate assemblies of repetitive loci. The assemblies reveal that satDNA repeats are organized into large arrays interrupted by transposable elements. The repeats in the center of the array tend to be homogenized in sequence, suggesting that gene conversion and unequal crossovers lead to repeat homogenization through concerted evolution, although the degree of unequal crossing over may differ among complex satellite loci. We find evidence for higher-order structure within satDNA arrays that suggest recent structural rearrangements. These assemblies provide a platform for the evolutionary and functional genomics of satDNAs in pericentric heterochromatin.

Improved assembly of noisy long reads by k-mer validation

Genome assembly depends critically on read length. Two recent technologies, from Pacific Biosciences (PacBio) and Oxford Nanopore, produce read lengths >20 kb, which yield de novo genome assemblies with vastly greater contiguity than those based on Sanger, Illumina, or other technologies. However, the very high error rates of these two new technologies (∼15% per base) makes assembly imprecise at repeats longer than the read length and computationally expensive. Here we show that the contiguity and quality of the assembly of these noisy long reads can be significantly improved at a minimal cost, by leveraging on the low error rate and low cost of Illumina short reads. Namely, k-mers from the PacBio raw reads that are not present in Illumina reads (which account for ∼95% of the distinct k-mers) are deemed sequencing errors and ignored at the seed alignment step. By focusing on the ∼5% of k-mers that are error free, read overlap sensitivity is dramatically increased. Of equal importance, the validation procedure can be extended to exclude repetitive k-mers, which prevents read miscorrection at repeats and further improves the resulting assemblies. We tested the k-mer validation procedure using one long-read technology (PacBio) and one assembler (MHAP/Celera Assembler), but it is very likely to yield analogous improvements with alternative long-read technologies and assemblers, such as Oxford Nanopore and BLASR/DALIGNER/Falcon, respectively.

Radical remodeling of the Y chromosome in a recent radiation of malaria mosquitoes

Y chromosomes control essential male functions in many species, including sex determination and fertility. However, because of obstacles posed by repeat-rich heterochromatin, knowledge of Y chromosome sequences is limited to a handful of model organisms, constraining our understanding of Y biology across the tree of life. Here, we leverage long single-molecule sequencing to determine the content and structure of the nonrecombining Y chromosome of the primary African malaria mosquito, Anopheles gambiae We find that the An. gambiae Y consists almost entirely of a few massively amplified, tandemly arrayed repeats, some of which can recombine with similar repeats on the X chromosome. Sex-specific genome resequencing in a recent species radiation, the An. gambiae complex, revealed rapid sequence turnover within An. gambiae and among species. Exploiting 52 sex-specific An. gambiae RNA-Seq datasets representing all developmental stages, we identified a small repertoire of Y-linked genes that lack X gametologs and are not Y-linked in any other species except An. gambiae, with the notable exception of YG2, a candidate male-determining gene. YG2 is the only gene conserved and exclusive to the Y in all species examined, yet sequence similarity to YG2 is not detectable in the genome of a more distant mosquito relative, suggesting rapid evolution of Y chromosome genes in this highly dynamic genus of malaria vectors. The extensive characterization of the An. gambiae Y provides a long-awaited foundation for studying male mosquito biology, and will inform novel mosquito control strategies based on the manipulation of Y chromosomes.