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

Sex chromosome evolution in haploid plants: Microchromosomes, disappearing chromosomes, and giant chromosomes

As in many diploid organisms with genetic sex determination, haploid-dominant organisms have also evolved sex chromosomes or extensive genomic regions that lack genetic recombination. An understanding of sex chromosome evolution should explain the causes and consequences of such regions in both diploids and haploids. However, haploids have been little studied, even though differences from sex chromosomes in diploids carry implications concerning the evolution of suppressed recombination in diploid organisms, and make predictions about genome evolution in the sex-linked regions of haploids that can now be tested by approaches using genome sequences. I review these ideas, and the current empirical evidence concerning them, in more detail than recent reviews focusing on progress in understanding the mechanisms involved in sex determination. I also discuss evidence that one specific prediction, that genetic degeneration should be minor in haploids, is not upheld. I suggest that this prediction does not take account of all processes leading to gene loss from sex-linked regions and that profound degeneration may evolve if sex-linked genes become duplicated to autosomes, a process that also appears to occur in diploids. I emphasize types of data that are needed to make progress in testing several of the ideas described.

Small RNA-mediated suppression of sex chromosome meiotic conflicts during Drosophila male gametogenesis

Meiosis is an evolutionarily conserved process in eukaryotes that ensures equal segregation of alleles and chromosomes during reproduction. Although parity in allelic transmission is the norm, selfish genes such as meiotic drivers can violate Mendel's first law of segregation. Sex chromosome drive is a form of meiotic drive that leads to unequal segregation of sex chromosomes, resulting in sex-ratio distortion and/or sterility in the offspring. Adverse fitness effects due to sex chromosome drive trigger the evolution of suppressors to restore Mendelian segregation. However, the molecular mechanisms by which suppressors emerge and counteract meiotic drive genes remain unclear. Recent studies from Drosophila have shed light on the critical roles of small RNA-mediated post-transcriptional silencing in mitigating sex chromosome meiotic conflicts. This review highlights the recruitment of two distinct small RNA pathways to combat intragenomic conflicts during male gametogenesis and seeks to reveal the impact of molecular arms races between meiotic drivers and their suppressors in shaping genome and sex chromosome evolution.

Transcription of a centromere-enriched retroelement and local retention of its RNA are significant features of the CENP-A chromatin landscape

BACKGROUND

Centromeres depend on chromatin containing the conserved histone H3 variant CENP-A for function and inheritance, while the role of centromeric DNA repeats remains unclear. Retroelements are prevalent at centromeres across taxa and represent a potential mechanism for promoting transcription to aid in CENP-A incorporation or for generating RNA transcripts to maintain centromere integrity.

RESULTS

In this study, we probe into the transcription and RNA localization of the centromere-enriched retroelement G2/Jockey-3 (hereafter referred to as Jockey-3) in Drosophila melanogaster, currently the only in vivo model with assembled centromeres. We find that Jockey-3 is a major component of the centromeric transcriptome and produces RNAs that localize to centromeres in metaphase. Leveraging the polymorphism of Jockey-3 and a de novo centromere system, we show that these RNAs remain associated with their cognate DNA sequences in cis, suggesting they are unlikely to perform a sequence-specific function at all centromeres. We show that Jockey-3 transcription is positively correlated with the presence of CENP-A and that recent Jockey-3 transposition events have occurred preferentially at CENP-A-containing chromatin.

CONCLUSIONS

We propose that Jockey-3 preferentially inserts at the centromere to ensure its own selfish propagation, while contributing to transcription across these regions. Given the conservation of retroelements as centromere components through evolution, our findings may offer a basis for understanding similar associations in other species.

Turnover of retroelements and satellite DNA drives centromere reorganization over short evolutionary timescales in Drosophila

Centromeres reside in rapidly evolving, repeat-rich genomic regions, despite their essential function in chromosome segregation. Across organisms, centromeres are rich in selfish genetic elements such as transposable elements and satellite DNAs that can bias their transmission through meiosis. However, these elements still need to cooperate at some level and contribute to, or avoid interfering with, centromere function. To gain insight into the balance between conflict and cooperation at centromeric DNA, we take advantage of the close evolutionary relationships within the Drosophila simulans clade-D. simulans, D. sechellia, and D. mauritiana-and their relative, D. melanogaster. Using chromatin profiling combined with high-resolution fluorescence in situ hybridization on stretched chromatin fibers, we characterize all centromeres across these species. We discovered dramatic centromere reorganization involving recurrent shifts between retroelements and satellite DNAs over short evolutionary timescales. We also reveal the recent origin (<240 Kya) of telocentric chromosomes in D. sechellia, where the X and fourth centromeres now sit on telomere-specific retroelements. Finally, the Y chromosome centromeres, which are the only chromosomes that do not experience female meiosis, do not show dynamic cycling between satDNA and TEs. The patterns of rapid centromere turnover in these species are consistent with genetic conflicts in the female germline and have implications for centromeric DNA function and karyotype evolution. Regardless of the evolutionary forces driving this turnover, the rapid reorganization of centromeric sequences over short evolutionary timescales highlights their potential as hotspots for evolutionary innovation.

A germline PAF1 paralog complex ensures cell type-specific gene expression

Animal germline development and fertility rely on paralogs of general transcription factors that recruit RNA polymerase II to ensure cell type-specific gene expression. It remains unclear whether gene expression processes downstream from such paralog-based transcription is distinct from that of canonical RNA polymerase II genes. In Drosophila, the testis-specific TBP-associated factors (tTAFs) activate over a thousand spermatocyte-specific gene promoters to enable meiosis and germ cell differentiation. Here, we show that efficient termination of tTAF-activated transcription relies on testis-specific paralogs of canonical polymerase-associated factor 1 complex (PAF1C) proteins, which form a testis-specific PAF1C (tPAF). Consequently, tPAF mutants show aberrant expression of hundreds of downstream genes due to read-in transcription. Furthermore, tPAF facilitates expression of Y-linked male fertility factor genes and thus serves to maintain spermatocyte-specific gene expression. Consistently, tPAF is required for the segregation of meiotic chromosomes and male fertility. Supported by comparative in vivo protein interaction assays, we provide a mechanistic model for the functional divergence of tPAF and the PAF1C and identify transcription termination as a developmentally regulated process required for germline-specific gene expression.

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.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

Single-fly genome assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Genetic variation in recalcitrant repetitive regions of the Drosophila melanogaster genome

Many essential functions of organisms are encoded in highly repetitive genomic regions, including histones involved in DNA packaging, centromeres that are core components of chromosome segregation, ribosomal RNA comprising the protein translation machinery, telomeres that ensure chromosome integrity, piRNA clusters encoding host defenses against selfish elements, and virtually the entire Y chromosome. These regions, formed by highly similar tandem arrays, pose significant challenges for experimental and informatic study, impeding sequence-level descriptions essential for understanding genetic variation. Here, we report the assembly and variation analysis of such repetitive regions in Drosophila melanogaster, offering significant improvements to the existing community reference assembly. Our work successfully recovers previously elusive segments, including complete reconstructions of the histone locus and the pericentric heterochromatin of the X chromosome, spanning the Stellate locus to the distal flank of the rDNA cluster. To infer structural changes in these regions where alignments are often not practicable, we introduce landmark anchors based on unique variants that are putatively orthologous. These regions display considerable structural variation between different D. melanogaster strains, exhibiting differences in copy number and organization of homologous repeat units between haplotypes. In the histone cluster, although we observe minimal genetic exchange indicative of crossing over, the variation patterns suggest mechanisms such as unequal sister chromatid exchange. We also examine the prevalence and scale of concerted evolution in the histone and Stellate clusters and discuss the mechanisms underlying these observed patterns.

Spatially revealed roles for lncRNAs in Drosophila spermatogenesis, Y chromosome function and evolution

Unlike coding genes, the number of lncRNA genes in organism genomes is relatively proportional to organism complexity. From plants to humans, the tissues with highest numbers and levels of lncRNA gene expression are the male reproductive organs. To learn why, we initiated a genome-wide analysis of Drosophila lncRNA spatial expression patterns in these tissues. The numbers of genes and levels of expression observed greatly exceed those previously reported, due largely to a preponderance of non-polyadenylated transcripts. In stark contrast to coding genes, the highest numbers of lncRNAs expressed are in post-meiotic spermatids. Correlations between expression levels, localization and previously performed genetic analyses indicate high levels of function and requirement. More focused analyses indicate that lncRNAs play major roles in evolution by controlling transposable element activities, Y chromosome gene expression and sperm construction. A new type of lncRNA-based particle found in seminal fluid may also contribute to reproductive outcomes.

Transcriptionally active chromatin loops contain both 'active' and 'inactive' histone modifications that exhibit exclusivity at the level of nucleosome clusters

Chromatin state is thought to impart regulatory function to the underlying DNA sequence. This can be established through histone modifications and chromatin organisation, but exactly how these factors relate to one another to regulate gene expression is unclear. In this study, we have used super-resolution microscopy to image the Y loops of Drosophila melanogaster primary spermatocytes, which are enormous transcriptionally active chromatin fibres, each representing single transcription units that are individually resolvable in the nuclear interior. We previously found that the Y loops consist of regular clusters of nucleosomes, with an estimated median of 54 nucleosomes per cluster with wide variation.In this study, we report that the histone modifications H3K4me3, H3K27me3, and H3K36me3 are also clustered along the Y loops, with H3K4me3 more associated with diffuse chromatin compared to H3K27me3. These histone modifications form domains that can be stretches of Y loop chromatin micrometres long, or can be in short alternating domains. The different histone modifications are associated with different sizes of chromatin clusters and unique morphologies. Strikingly, a single chromatin cluster almost always only contains only one type of the histone modifications that were labelled, suggesting exclusivity, and therefore regulation at the level of individual chromatin clusters. The active mark H3K36me3 is more associated with actively elongating RNA polymerase II than H3K27me3, with polymerase often appearing on what are assumed to be looping regions on the periphery of chromatin clusters.These results provide a foundation for understanding the relationship between chromatin state, chromatin organisation, and transcription regulation - with potential implications for pause-release dynamics, splicing complex organisation and chromatin dynamics during polymerase progression along a gene.

The implications of satellite DNA instability on cellular function and evolution

Abundant tandemly repeated satellite DNA is present in most eukaryotic genomes. Previous limitations including a pervasive view that it was uninteresting junk DNA, combined with challenges in studying it, are starting to dissolve - and recent studies have found important functions for satellite DNAs. The observed rapid evolution and implied instability of satellite DNA now has important significance for their functions and maintenance within the genome. In this review, we discuss the processes that lead to satellite DNA copy number instability, and the importance of mechanisms to manage the potential negative effects of instability. Satellite DNA is vulnerable to challenges during replication and repair, since it forms difficult-to-process secondary structures and its homology within tandem arrays can result in various types of recombination. Satellite DNA instability may be managed by DNA or chromatin-binding proteins ensuring proper nuclear localization and repair, or by proteins that process aberrant structures that satellite DNAs tend to form. We also discuss the pattern of satellite DNA mutations from recent mutation accumulation (MA) studies that have tracked changes in satellite DNA for up to 1000 generations with minimal selection. Finally, we highlight examples of satellite evolution from studies that have characterized satellites across millions of years of Drosophila fruit fly evolution, and discuss possible ways that selection might act on the satellite DNA composition.

The Drosophila melanogaster Y-linked gene, WDY, is required for sperm to swim in the female reproductive tract

Unique patterns of inheritance and selection on Y chromosomes have led to the evolution of specialized gene functions. We report CRISPR mutants in Drosophila of the Y-linked gene, WDY, which is required for male fertility. We demonstrate that the sperm tails of WDY mutants beat approximately half as fast as those of wild-type and that mutant sperm do not propel themselves within the male ejaculatory duct or female reproductive tract. Therefore, although mature sperm are produced by WDY mutant males, and are transferred to females, those sperm fail to enter the female sperm storage organs. We report genotype-dependent and regional differences in sperm motility that appear to break the correlation between sperm tail beating and propulsion. Furthermore, we identify a significant change in hydrophobicity at a residue at a putative calcium-binding site in WDY orthologs at the split between the melanogaster and obscura species groups, when WDY first became Y-linked. This suggests that a major functional change in WDY coincided with its appearance on the Y chromosome. Finally, we show that mutants for another Y-linked gene, PRY, also show a sperm storage defect that may explain their subfertility. Overall, we provide direct evidence for the long-held presumption that protein-coding genes on the Drosophila Y regulate sperm motility.

High levels of intra-strain structural variation in Drosophila simulans X pericentric heterochromatin

Large genome structural variations can impact genome regulation and integrity. Repeat-rich regions like pericentric heterochromatin are vulnerable to structural rearrangements although we know little about how often these rearrangements occur over evolutionary time. Repetitive genome regions are particularly difficult to study with genomic approaches, as they are missing from most genome assemblies. However, cytogenetic approaches offer a direct way to detect large rearrangements involving pericentric heterochromatin. Here, we use a cytogenetic approach to reveal large structural rearrangements associated with the X pericentromeric region of Drosophila simulans. These rearrangements involve large blocks of satellite DNA-the 500-bp and Rsp-like satellites-which colocalize in the X pericentromeric heterochromatin. We find that this region is polymorphic not only among different strains, but between isolates of the same strain from different labs, and even within individual isolates. On the one hand, our observations raise questions regarding the potential impact of such variation at the phenotypic level and our ability to control for such genetic variability. On the other hand, this highlights the very rapid turnover of the pericentric heterochromatin most likely associated with genomic instability of the X pericentromere. It represents a unique opportunity to study the dynamics of pericentric heterochromatin, the evolution of associated satellites on a very short time scale, and to better understand how structural variation arises.

Chromosome-level genome assembly of an agricultural pest Zeugodacus tau (Diptera: Tephritidae)

The fruit fly Zeugodacus tau (Diptera: Tephritidae) is a major pest of melons and other cucurbits in Southeast Asia. In this study, we used Illumina, Nanopore, and Hi-C sequencing technologies to assemble a reference genome of Z. tau at the chromosomal level. The assembled genome was 421.79 Mb and consisted of six chromosomes (one X-chromosome + five autosomes). The contig N50 was 4.23 Mb. We identified 20,922 protein-coding genes, of which 17,251 (82.45%) were functionally annotated. Additionally, we found 247 rRNAs, 435 tRNAs, 67 small nuclear RNAs, and 829 small RNAs in the genome. Repetitive elements accounted for 55.30 Mb (13.15%) of the genome. This high-quality genome assembly is valuable for evolutionary and genetic studies of Z. tau and its relative species.

Chromosome-specific maturation of the epigenome in the Drosophila male germline

Spermatogenesis in the Drosophila male germline proceeds through a unique transcriptional program controlled both by germline-specific transcription factors and by testis-specific versions of core transcriptional machinery. This program includes the activation of genes on the heterochromatic Y chromosome, and reduced transcription from the X chromosome, but how expression from these sex chromosomes is regulated has not been defined. To resolve this, we profiled active chromatin features in the testes from wildtype and meiotic arrest mutants and integrate this with single-cell gene expression data from the Fly Cell Atlas. These data assign the timing of promoter activation for genes with germline-enriched expression throughout spermatogenesis, and general alterations of promoter regulation in germline cells. By profiling both active RNA polymerase II and histone modifications in isolated spermatocytes, we detail widespread patterns associated with regulation of the sex chromosomes. Our results demonstrate that the X chromosome is not enriched for silencing histone modifications, implying that sex chromosome inactivation does not occur in the Drosophila male germline. Instead, a lack of dosage compensation in spermatocytes accounts for the reduced expression from this chromosome. Finally, profiling uncovers dramatic ubiquitinylation of histone H2A and lysine-16 acetylation of histone H4 across the Y chromosome in spermatocytes that may contribute to the activation of this heterochromatic chromosome.

Genetic conflict and the origin of multigene families: implications for sex chromosome evolution

Sex chromosomes are havens for intragenomic conflicts. The absence of recombination between sex chromosomes creates the opportunity for the evolution of segregation distorters: selfish genetic elements that hijack different aspects of an individual's reproduction to increase their own transmission. Biased (non-Mendelian) segregation, however, often occurs at a detriment to their host's fitness, and therefore can trigger evolutionary arms races that can have major consequences for genome structure and regulation, gametogenesis, reproductive strategies and even speciation. Here, we review an emerging feature from comparative genomic and sex chromosome evolution studies suggesting that meiotic drive is pervasive: the recurrent evolution of paralogous sex-linked gene families. Sex chromosomes of several species independently acquire and co-amplify rapidly evolving gene families with spermatogenesis-related functions, consistent with a history of intragenomic conflict over transmission. We discuss Y chromosome features that might contribute to the tempo and mode of evolution of X/Y co-amplified gene families, as well as their implications for the evolution of complexity in the genome. Finally, we propose a framework that explores the conditions that might allow for recurrent bouts of fixation of drivers and suppressors, in a dosage-sensitive fashion, and therefore the co-amplification of multigene families on sex chromosomes.

The chromatin source-sink hypothesis: a shared mode of chromatin-mediated regulations

We propose that several chromatin-mediated regulatory processes are dominated by source-sink relationships in which factors operate as 'sources' to produce or provide a resource and compete with each other to occupy separate 'sinks'. In this model, large portions of genomic DNA operate as 'sinks', which are filled by 'sources', such as available histone variants, covalent modifications to histones, the readers of these modifications and non-coding RNAs. Competing occupation for the sinks by different sources leads to distinct states of genomic equilibrium in differentiated cells. During dynamic developmental events, such as sexual reproduction, we propose that dramatic and rapid reconfiguration of source-sink relationships modifies chromatin states. We envision that re-routing of sources could occur by altering the dimensions of the sink, by reconfiguration of existing sink occupation or by varying the size of the source, providing a central mechanism to explain a plethora of epigenetic phenomena, which contribute to phenotypic variegation, zygotic genome activation and nucleolar dominance.

Genomes of the cosmopolitan fruit pest Bactrocera dorsalis (Diptera: Tephritidae) reveal its global invasion history and thermal adaptation

INTRODUCTION

The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural pests worldwide, with highly debated species delimitation, origin, and global spread routes.

OBJECTIVES

Our study intended to (i) resolve the taxonomic uncertainties between B. dorsalis and B. carambolae, (ii) reveal the population structure and global invasion routes of B. dorsalis across Asia, Africa, and Oceania, and (iii) identify genomic regions that are responsible for the thermal adaptation of B. dorsalis.

METHODS

Based on a high-quality chromosome-level reference genome assembly, we explored the population relationship using a genome-scale single nucleotide polymorphism dataset generated from the resequencing data of 487 B. dorsalis genomes and 25 B. carambolae genomes. Genome-wide association studies and silencing using RNA interference were used to identify and verify the candidate genes associated with extreme thermal stress.

RESULTS

We showed that B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide: (i) from Northern India to Northern Southeast Asia, then to Southern Southeast Asia; (ii) from Northern India to Northern Southeast Asian, then to China and Hawaii; and (iii) from Southern India toward the African mainland, then to Madagascar, which is mainly facilitated by human activities including trade and immigration. Twenty-seven genes were identified by a genome-wide association study to be associated with 11 temperature bioclimatic variables. The Cyp6a9 gene may enhance the thermal adaptation of B. dorsalis and thus boost its invasion, which tended to be upregulated at a hardening temperature of 38 °C. Functional verification using RNA interference silencing against Cyp6a9, led to the specific decrease in Cyp6a9 expression, reducing the survival rate of dsRNA-feeding larvae exposed to extreme thermal stress of 45 °C after heat hardening treatments in B. dorsalis.

CONCLUSION

This study provides insights into the evolutionary history and genetic basis of temperature adaptation in B. dorsalis.

Chromosome-specific maturation of the epigenome in the Drosophila male germline

Spermatogenesis in the Drosophila male germline proceeds through a unique transcriptional program controlled both by germline-specific transcription factors and by testis-specific versions of core transcriptional machinery. This program includes the activation of genes on the heterochromatic Y chromosome, and reduced transcription from the X chromosome, but how expression from these sex chromosomes is regulated has not been defined. To resolve this, we profiled active chromatin features in the testes from wildtype and meiotic arrest mutants and integrate this with single-cell gene expression data from the Fly Cell Atlas. These data assign the timing of promoter activation for genes with germline-enriched expression throughout spermatogenesis, and general alterations of promoter regulation in germline cells. By profiling both active RNA polymerase II and histone modifications in isolated spermatocytes, we detail widespread patterns associated with regulation of the sex chromosomes. Our results demonstrate that the X chromosome is not enriched for silencing histone modifications, implying that sex chromosome inactivation does not occur in the Drosophila male germline. Instead, a lack of dosage compensation in spermatocytes accounts for the reduced expression from this chromosome. Finally, profiling uncovers dramatic ubiquitinylation of histone H2A and lysine-16 acetylation of histone H4 across the Y chromosome in spermatocytes that may contribute to the activation of this heterochromatic chromosome.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

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

Y chromosome toxicity does not contribute to sex-specific differences in longevity

While sex chromosomes carry sex-determining genes, they also often differ from autosomes in size and composition, consisting mainly of silenced heterochromatic repetitive DNA. Even though Y chromosomes show structural heteromorphism, the functional significance of such differences remains elusive. Correlative studies suggest that the amount of Y chromosome heterochromatin might be responsible for several male-specific traits, including sex-specific differences in longevity observed across a wide spectrum of species, including humans. However, experimental models to test this hypothesis have been lacking. Here we use the Drosophila melanogaster Y chromosome to investigate the relevance of sex chromosome heterochromatin in somatic organs in vivo. Using CRISPR-Cas9, we generated a library of Y chromosomes with variable levels of heterochromatin. We show that these different Y chromosomes can disrupt gene silencing in trans, on other chromosomes, by sequestering core components of the heterochromatin machinery. This effect is positively correlated to the level of Y heterochromatin. However, we also find that the ability of the Y chromosome to affect genome-wide heterochromatin does not generate physiological sex differences, including sexual dimorphism in longevity. Instead, we discovered that it is the phenotypic sex, female or male, that controls sex-specific differences in lifespan, rather than the presence of a Y chromosome. Altogether, our findings dismiss the 'toxic Y' hypothesis that postulates that the Y chromosome leads to reduced lifespan in XY individuals.

Escaping but not the inactive X-linked protein complex coding genes may achieve X-chromosome dosage compensation and underlie X chromosome inactivation-related diseases

X chromosome dosage compensation (XDC) refers to the process by which X-linked genes acquire expression equivalence between two sexes. Ohno proposed that XDC is achieved by two-fold upregulations of X-linked genes in both sexes and by silencing one X chromosome (X chromosome inactivation, XCI) in females. However, genes subject to two-fold upregulations as well as the underlying mechanism remain unclear. It's reported that gene dosage changes may only affect X-linked dosage-sensitive genes, such as protein complex coding genes (PCGs). Our results showed that in human PCGs are more likely to escape XCI and escaping PCGs (EsP) show two-fold higher expression than inactivated PCGs (InP) or other X-linked genes at RNA and protein levels in both sexes, which suggest that EsP may achieve upregulations and XDC. The higher expressions of EsP possibly result from the upregulations of the single active X chromosome (Xa), rather than escaping expressions from the inactive X chromosome (Xi). EsP genes have relatively high expression levels in humans and lower dN/dS ratios, suggesting that they are likely under stronger selection pressure over evolutionary time. Our study also suggests that SP1 transcription factor is significantly enriched in EsP and may be involved in the up-regulations of EsP on the active X. Finally, human EsP genes in this study are enriched in the toll-like receptor pathway, NF-kB pathway, apoptotic pathway, and abnormal mental, developmental and reproductive phenotypes. These findings suggest misregulations of EsP may be involved in autoimmune, reproductive, and neurological diseases, providing insight for the diagnosis and treatment of these diseases.

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.

Y chromosome-linked variation affects locomotor activity in male Drosophila melanogaster and is robust to differences in thermal environment

Although containing genes important for sex determination, genetic variation within the Y chromosome was traditionally predicted to contribute little to the expression of sexually dimorphic traits. This prediction was shaped by the assumption that the chromosome harbours few protein-coding genes, and that capacity for Y-linked variation to shape adaptation would be hindered by the chromosome's lack of recombination and holandric inheritance. Consequently, most studies exploring the genotypic contributions to sexually dimorphic traits have focused on the autosomes and X chromosome. Yet, several studies have now demonstrated that the Y chromosome harbours variation affecting male fitness, moderating the expression of hundreds of genes across the nuclear genome. Furthermore, emerging results have shown that expression of this Y-linked variation may be sensitive to environmental heterogeneity, leading to the prediction that Y-mediated gene-by-environment interactions will shape the expression of sexually dimorphic phenotypes. We tested this prediction, investigating whether genetic variation across six distinct Y chromosome haplotypes affects the expression of locomotor activity, at each of two temperatures (20 and 28 °C) in male fruit flies (Drosophila melanogaster). Locomotor activity is a sexually dimorphic trait in this species, previously demonstrated to be under intralocus sexual conflict. We demonstrate Y haplotype effects on male locomotor activity, but the rank order and magnitude of these effects were unaltered by differences in temperature. Our study contributes to a growing number of studies demonstrating Y-linked effects moderating expression of traits evolving under sexually antagonistic selection, suggesting a role for the Y chromosome in shaping outcomes of sexual conflict.

Long-read-based Genome Assembly of Drosophila gunungcola Reveals Fewer Chemosensory Genes in Flower-breeding Species

Drosophila gunungcola exhibits reproductive activities on the fresh flowers of several plant species and is an emerging model to study the co-option of morphological and behavioral traits in male courtship display. Here, we report a near-chromosome-level genome assembly that was constructed based on long-read PacBio sequencing data (with ∼66× coverage) and annotated with the assistant from RNA-seq transcriptome data of whole organisms at various developmental stages. A nuclear genome of 189 Mb with 13,950 protein-coding genes and a mitogenome of 17.5 kb were acquired. Few interchromosomal rearrangements were found in the comparisons of synteny with Drosophila elegans, its sister species, and Drosophila melanogaster, suggesting that the gene compositions on each Muller element are evolutionarily conserved. Loss events of several OR and IR genes in D. gunungcola and D. elegans were revealed when orthologous genomic regions were compared across species in the D. melanogaster species group. This high-quality reference genome will facilitate further comparative studies on traits related to the evolution of sexual behavior and diet specialization.

The Drosophila melanogaster Y-linked gene, WDY, is required for sperm to swim in the female reproductive tract

Unique patterns of inheritance and selection on Y chromosomes lead to the evolution of specialized gene functions. Yet characterizing the function of genes on Y chromosomes is notoriously difficult. We report CRISPR mutants in Drosophila of the Y-linked gene, WDY, which is required for male fertility. WDY mutants produce mature sperm with beating tails that can be transferred to females but fail to enter the female sperm storage organs. We demonstrate that the sperm tails of WDY mutants beat approximately half as fast as wild-type sperm's and that the mutant sperm do not propel themselves within the male ejaculatory duct or female reproductive tract (RT). These specific motility defects likely cause the sperm storage defect and sterility of the mutants. Regional and genotype-dependent differences in sperm motility suggest that sperm tail beating and propulsion do not always correlate. Furthermore, we find significant differences in the hydrophobicity of key residues of a putative calcium-binding domain between orthologs of WDY that are Y-linked and those that are autosomal. Given that WDY appears to be evolving under positive selection, our results suggest that WDY's functional evolution coincides with its transition from autosomal to Y-linked in Drosophila melanogaster and its most closely related species. Finally, we show that mutants for another Y-linked gene, PRY, also show a sperm storage defect that may explain their subfertility. In contrast to WDY, PRY mutants do swim in the female RT, suggesting they are defective in yet another mode of motility, navigation, or a necessary interaction with the female RT. Overall, we provide direct evidence for the long-held presumption that protein-coding genes on the Drosophila Y regulate sperm motility.

Expansion and loss of sperm nuclear basic protein genes in Drosophila correspond with genetic conflicts between sex chromosomes

Many animal species employ sperm nuclear basic proteins (SNBPs) or protamines to package sperm genomes tightly. SNBPs vary across animal lineages and evolve rapidly in mammals. We used a phylogenomic approach to investigate SNBP diversification in Drosophila species. We found that most SNBP genes in Drosophila melanogaster evolve under positive selection except for genes essential for male fertility. Unexpectedly, evolutionarily young SNBP genes are more likely to be critical for fertility than ancient, conserved SNBP genes. For example, CG30056 is dispensable for male fertility despite being one of three SNBP genes universally retained in Drosophila species. We found 19 independent SNBP gene amplification events that occurred preferentially on sex chromosomes. Conversely, the montium group of Drosophila species lost otherwise-conserved SNBP genes, coincident with an X-Y chromosomal fusion. Furthermore, SNBP genes that became linked to sex chromosomes via chromosomal fusions were more likely to degenerate or relocate back to autosomes. We hypothesize that autosomal SNBP genes suppress meiotic drive, whereas sex-chromosomal SNBP expansions lead to meiotic drive. X-Y fusions in the montium group render autosomal SNBPs dispensable by making X-versus-Y meiotic drive obsolete or costly. Thus, genetic conflicts between sex chromosomes may drive SNBP rapid evolution during spermatogenesis in Drosophila species.

The Drosophila simulans Genome Lacks the crystal-Stellate System

The cry-Ste system is a genetic interaction system between heterochromatin and euchromatin in Drosophila melanogaster, regulated via the piRNA pathway. Deregulation of this system leads to meiotic defects and male sterility. Although the cry-Ste system is peculiar to D. melanogaster, ancestors of Ste and Su(Ste) elements are present in the three closely related species, D. simulans, D. sechellia, and D. mauritiana. The birth, evolution, and maintenance of this genetic system in Drosophila melanogaster are of interest. We investigate the presence of sequences homologous to cry and Ste elements in the simulans complex and describe their chromosomal distribution. The organization and expression of cry- and Ste-like sequences were further characterized in the D. simulans genome. Our results allow us to conclude that the cry-Ste genetic interaction system is absent in the D. simulans genome.

The Rattlesnake W Chromosome: A GC-Rich Retroelement Refugium with Retained Gene Function Across Ancient Evolutionary Strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a "refugium" for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

The nanoCUT&RUN technique visualizes telomeric chromatin in Drosophila

Advances in genomic technology led to a more focused pattern for the distribution of chromosomal proteins and a better understanding of their functions. The recent development of the CUT&RUN technique marks one of the important such advances. Here we develop a modified CUT&RUN technique that we termed nanoCUT&RUN, in which a high affinity nanobody to GFP is used to bring micrococcal nuclease to the binding sites of GFP-tagged chromatin proteins. Subsequent activation of the nuclease cleaves the chromatin, and sequencing of released DNA identifies binding sites. We show that nanoCUT&RUN efficiently produces high quality data for the TRL transcription factor in Drosophila embryos, and distinguishes binding sites specific between two TRL isoforms. We further show that nanoCUT&RUN dissects the distributions of the HipHop and HOAP telomere capping proteins, and uncovers unexpected binding of telomeric proteins at centromeres. nanoCUT&RUN can be readily applied to any system in which a chromatin protein of interest, or its isoforms, carries the GFP tag.

The method of chromatin immunoprecipitation followed by genomic sequencing (ChIP-seq) has been employed to study the distribution of chromatin binding proteins genome-wide. Such studies have greatly enhanced our understanding of the function of the target proteins. However, the uses of chemical crosslinking combined with the procedure of antibody-medicated precipitation of the protein-DNA complex have limited the efficiency of ChIP-seq. The recently developed CUT&RUN method has greatly improved that efficiency. We here developed the “nanoCUT&RUN” extension of CUT&RUN, which can be readily applied to any target protein with a GFP tag. Using nanoCUT&RUN, we profiled the HipHop and HOAP proteins that protect telomeric chromatin in Drosophila. We uncovered sequence-independent binding of both proteins predominantly to telomeres. Interestingly, HipHop binding can also be detected at centromeric chromatin suggestive of a novel function of a telomere capping protein.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

Validation of reference-assisted assembly using existing and novel Heliothine genomes

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

Enrichment of Non-B-Form DNA at D. melanogaster Centromeres

Centromeres are essential chromosomal regions that mediate the accurate inheritance of genetic information during eukaryotic cell division. Despite their conserved function, centromeres do not contain conserved DNA sequences and are instead epigenetically marked by the presence of the centromere-specific histone H3 variant centromeric protein A. The functional contribution of centromeric DNA sequences to centromere identity remains elusive. Previous work found that dyad symmetries with a propensity to adopt noncanonical secondary DNA structures are enriched at the centromeres of several species. These findings lead to the proposal that noncanonical DNA structures may contribute to centromere specification. Here, we analyze the predicted secondary structures of the recently identified centromere DNA sequences of Drosophila melanogaster. Although dyad symmetries are only enriched on the Y centromere, we find that other types of noncanonical DNA structures, including melted DNA and G-quadruplexes, are common features of all D. melanogaster centromeres. Our work is consistent with previous models suggesting that noncanonical DNA secondary structures may be conserved features of centromeres with possible implications for centromere specification.

Transposable element accumulation drives size differences among polymorphic Y Chromosomes in Drosophila

Y Chromosomes of many species are gene poor and show low levels of nucleotide variation, yet often display high amounts of structural diversity. Dobzhansky cataloged several morphologically distinct Y Chromosomes in Drosophila pseudoobscura that differ in size and shape, but the molecular causes of their dramatic size differences are unclear. Here we use cytogenetics and long-read sequencing to study the sequence content of polymorphic Y Chromosomes in D. pseudoobscura We show that Y Chromosomes differ almost 2-fold in size, ranging from 30 to 60 Mb. Most of this size difference is caused by a handful of active transposable elements (TEs) that have recently expanded on the largest Y Chromosome, with different elements being responsible for Y expansion on differently sized D. pseudoobscura Y's. We show that Y Chromosomes differ in their heterochromatin enrichment, expression of Y-enriched TEs, and also influence expression of dozens of autosomal and X-linked genes. The same helitron element that showed the most drastic amplification on the largest Y in D. pseudoobscura independently amplified on a polymorphic large Y Chromosome in D. affinis, suggesting that some TEs are inherently more prone to become deregulated on Y Chromosomes.

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

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.

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.

Blm helicase facilitates rapid replication of repetitive DNA sequences in early Drosophila development

The absence of functional BLM DNA helicase, a member of the RecQ family of helicases, is responsible for the rare human disorder Bloom Syndrome, which results in developmental abnormalities, DNA repair defects, genomic instability, and a predisposition to cancer. In Drosophila melanogaster, the orthologous Blm protein is essential during early development when the embryo is under the control of maternal gene products. We show that lack of functional maternal Blm during the syncytial cell cycles of Drosophila embryonic development results in severe nuclear defects and lethality. Amongst the small fraction of embryos from Blm mutant mothers that survive to adulthood, a prominent sex-bias favors the class that inherits less repetitive DNA content, which serves as an endogenous source of replication stress. This selection against repetitive DNA content reflects a role for Blm in facilitating replication through repetitive sequences during the rapid S-phases of syncytial cell cycles. During these syncytial cycles, Blm is not required for complex DNA double-strand break repair; however, the progeny sex-bias resulting from the absence of maternal Blm is exacerbated by repetitive DNA sequences and by the slowing of replication fork progression, suggesting that the essential role for Blm during this stage is to manage replication fork stress brought about by impediments to fork progression. Additionally, our data suggest that Blm is only required to manage this replication stress during embryonic development, and likely only during the early, rapid syncytial cell cycles, and not at later developmental stages. These results provide novel insights into Blm function throughout development.

A transposon expression burst accompanies the activation of Y-chromosome fertility genes during Drosophila spermatogenesis

Transposable elements (TEs) must replicate in germline cells to pass novel insertions to offspring. In Drosophila melanogaster ovaries, TEs can exploit specific developmental windows of opportunity to evade host silencing and increase their copy numbers. However, TE activity and host silencing in the distinct cell types of Drosophila testis are not well understood. Here, we reanalyze publicly available single-cell RNA-seq datasets to quantify TE expression in the distinct cell types of the Drosophila testis. We develop a method for identification of TE and host gene expression modules and find that a distinct population of early spermatocytes expresses a large number of TEs at much higher levels than other germline and somatic components of the testes. This burst of TE expression coincides with the activation of Y chromosome fertility factors and spermatocyte-specific transcriptional regulators, as well as downregulation of many components of the piRNA pathway. The TEs expressed by this cell population are specifically enriched on the Y chromosome and depleted on the X chromosome, relative to other active TEs. These data suggest that some TEs may achieve high insertional activity in males by exploiting a window of opportunity for mobilization created by the activation of spermatocyte-specific and Y chromosome-specific transcriptional programs.

The Organization of Pericentromeric Heterochromatin in Polytene Chromosome 3 of the Drosophilamelanogaster Line with the Rif11; SuURES Su(var)3-906 Mutations Suppressing Underreplication

Although heterochromatin makes up 40% of the Drosophila melanogaster genome, its organization remains little explored, especially in polytene chromosomes, as it is virtually not represented in them due to underreplication. Two all-new approaches were used in this work: (i) with the use of a newly synthesized Drosophila line that carries three mutations, Rif11, SuURES and Su(var)3-906, suppressing the underreplication of heterochromatic regions, we obtained their fullest representation in polytene chromosomes and described their structure; (ii) 20 DNA fragments with known positions on the physical map as well as molecular genetic features of the genome (gene density, histone marks, heterochromatin proteins, origin recognition complex proteins, replication timing sites and satellite DNAs) were mapped in the newly polytenized heterochromatin using FISH and bioinformatics data. The borders of the heterochromatic regions and variations in their positions on arm 3L have been determined for the first time. The newly polytenized heterochromatic material exhibits two main types of morphology: a banding pattern (locations of genes and short satellites) and reticular chromatin (locations of large blocks of satellite DNA). The locations of the banding and reticular polytene heterochromatin was determined on the physical map.

Heterochromatin-dependent transcription of satellite DNAs in the Drosophila melanogaster female germline

Large blocks of tandemly repeated DNAs-satellite DNAs (satDNAs)-play important roles in heterochromatin formation and chromosome segregation. We know little about how satDNAs are regulated; however, their misregulation is associated with genomic instability and human diseases. We use the Drosophila melanogaster germline as a model to study the regulation of satDNA transcription and chromatin. Here we show that complex satDNAs (>100-bp repeat units) are transcribed into long noncoding RNAs and processed into piRNAs (PIWI interacting RNAs). This satDNA piRNA production depends on the Rhino-Deadlock-Cutoff complex and the transcription factor Moonshiner-a previously described non-canonical pathway that licenses heterochromatin-dependent transcription of dual-strand piRNA clusters. We show that this pathway is important for establishing heterochromatin at satDNAs. Therefore, satDNAs are regulated by piRNAs originating from their own genomic loci. This novel mechanism of satDNA regulation provides insight into the role of piRNA pathways in heterochromatin formation and genome stability.

Improved contiguity of the threespine stickleback genome using long-read sequencing

While the cost and time for assembling a genome has drastically decreased, it still remains a challenge to assemble a highly contiguous genome. These challenges are rapidly being overcome by the integration of long-read sequencing technologies. Here, we use long-read sequencing to improve the contiguity of the threespine stickleback fish (Gasterosteus aculeatus) genome, a prominent genetic model species. Using Pacific Biosciences sequencing, we assembled a highly contiguous genome of a freshwater fish from Paxton Lake. Using contigs from this genome, we were able to fill over 76.7% of the gaps in the existing reference genome assembly, improving contiguity over fivefold. Our gap filling approach was highly accurate, validated by 10X Genomics long-distance linked-reads. In addition to closing a majority of gaps, we were able to assemble segments of telomeres and centromeres throughout the genome. This highlights the power of using long sequencing reads to assemble highly repetitive and difficult to assemble regions of genomes. This latest genome build has been released through a newly designed community genome browser that aims to consolidate the growing number of genomics datasets available for the threespine stickleback fish.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

Establishment of H3K9me3-dependent heterochromatin during embryogenesis in Drosophila miranda

Heterochromatin is a key architectural feature of eukaryotic genomes crucial for silencing of repetitive elements. During Drosophila embryonic cellularization, heterochromatin rapidly appears over repetitive sequences, but the molecular details of how heterochromatin is established are poorly understood. Here, we map the genome-wide distribution of H3K9me3-dependent heterochromatin in individual embryos of Drosophila miranda at precisely staged developmental time points. We find that canonical H3K9me3 enrichment is established prior to cellularization and matures into stable and broad heterochromatin domains through development. Intriguingly, initial nucleation sites of H3K9me3 enrichment appear as early as embryonic stage 3 over transposable elements (TEs) and progressively broaden, consistent with spreading to neighboring nucleosomes. The earliest nucleation sites are limited to specific regions of a small number of recently active retrotransposon families and often appear over promoter and 5' regions of LTR retrotransposons, while late nucleation sites develop broadly across the entirety of most TEs. Interestingly, early nucleating TEs are strongly associated with abundant maternal piRNAs and show early zygotic transcription. These results support a model of piRNA-associated co-transcriptional silencing while also suggesting additional mechanisms for site-restricted H3K9me3 nucleation at TEs in pre-cellular Drosophila embryos.

Toxic Y chromosome: Increased repeat expression and age-associated heterochromatin loss in male Drosophila with a young Y chromosome

Sex-specific differences in lifespan are prevalent across the tree of life and influenced by heteromorphic sex chromosomes. In species with XY sex chromosomes, females often outlive males. Males and females can differ in their overall repeat content due to the repetitive Y chromosome, and repeats on the Y might lower survival of the heterogametic sex (toxic Y effect). Here, we take advantage of the well-assembled young Y chromosome of Drosophila miranda to study the sex-specific dynamics of chromatin structure and repeat expression during aging in male and female flies. Male D. miranda have about twice as much repetitive DNA compared to females, and live shorter than females. Heterochromatin is crucial for silencing of repetitive elements, yet old D. miranda flies lose H3K9me3 modifications in their pericentromere, with heterochromatin loss being more severe during aging in males than females. Satellite DNA becomes de-repressed more rapidly in old vs. young male flies relative to females. In contrast to what is observed in D. melanogaster, we find that transposable elements (TEs) are expressed at higher levels in male D. miranda throughout their life. We show that epigenetic silencing via heterochromatin formation is ineffective on the TE-rich neo-Y chromosome, presumably due to active transcription of a large number of neo-Y linked genes, resulting in up-regulation of Y-linked TEs already in young males. This is consistent with an interaction between the evolutionary age of the Y chromosome and the genomic effects of aging. Our data support growing evidence that "toxic Y chromosomes" can diminish male fitness and a reduction in heterochromatin can contribute to sex-specific aging.

Evolutionary Dynamics of Abundant 7-bp Satellites in the Genome of Drosophila virilis

The factors that drive the rapid changes in abundance of tandem arrays of highly repetitive sequences, known as satellite DNA, are not well understood. Drosophila virilis has one of the highest relative amounts of simple satellites of any organism that has been studied, with an estimated >40% of its genome composed of a few related 7-bp satellites. Here, we use D. virilis as a model to understand technical biases affecting satellite sequencing and the evolutionary processes that drive satellite composition. By analyzing sequencing data from Illumina, PacBio, and Nanopore platforms, we identify platform-specific biases and suggest best practices for accurate characterization of satellites by sequencing. We use comparative genomics and cytogenetics to demonstrate that the highly abundant AAACTAC satellite family arose from a related satellite in the branch leading to the virilis phylad 4.5-11 Ma before exploding in abundance in some species of the clade. The most abundant satellite is conserved in sequence and location in the pericentromeric region but has diverged widely in abundance among species, whereas the satellites nearest the centromere are rapidly turning over in sequence composition. By analyzing multiple strains of D. virilis, we saw that the abundances of two centromere-proximal satellites are anticorrelated along a geographical gradient, which we suggest could be caused by ongoing conflicts at the centromere. In conclusion, we illuminate several key attributes of satellite evolutionary dynamics that we hypothesize to be driven by processes including selection, meiotic drive, and constraints on satellite sequence and abundance.

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.