Single-fly 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 (>1Mb 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.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

Urban evolution of thermal physiology in a range-expanding, mycophagous fruit fly, Drosophila tripunctata

In Drosophila spp., their often high number of annual generations, large population sizes and large amounts of standing genetic variation should predispose them to undergo contemporary adaptation to climatic warming. However, a number of laboratory experimental evolution studies in this group of organisms suggest strong limits on the rate and magnitude of contemporary thermal adaptation. Here, we explore this discrepancy by examining the potential for rapid evolutionary divergence between wild populations of Drosophila tripunctata Loew, 1862 from rural and urban sites. We performed a multi-generation common garden study and found evidence for the evolution of higher heat tolerance (critical thermal maximum) in flies from urban populations. We also detected evolutionary divergence in cold resistance (chill coma recovery time), with diminished cold resistance in flies from urban populations, although the effect was weaker than the shift in heat tolerance. Our study provides evidence of contemporary urban thermal adaptation, although the magnitude of phenotypic change lagged the magnitude of environmental temperature change across the urbanization gradient, suggesting potential limits on the evolution of urban thermal physiology.

DNA polymerase theta suppresses mitotic crossing over

Polymerase theta-mediated end joining (TMEJ) is a chromosome break repair pathway that is able to rescue the lethality associated with the loss of proteins involved in early steps in homologous recombination (e.g., BRCA1/2). This is due to the ability of polymerase theta (Pol θ) to use resected, 3' single stranded DNA tails to repair chromosome breaks. These resected DNA tails are also the starting substrate for homologous recombination. However, it remains unknown if TMEJ can compensate for the loss of proteins involved in more downstream steps during homologous recombination. Here we show that the Holliday junction resolvases SLX4 and GEN1 are required for viability in the absence of Pol θ in Drosophila melanogaster, and lack of all three proteins results in high levels of apoptosis. Flies deficient in Pol θ and SLX4 are extremely sensitive to DNA damaging agents, and mammalian cells require either Pol θ or SLX4 to survive. Our results suggest that TMEJ and Holliday junction formation/resolution share a common DNA substrate, likely a homologous recombination intermediate, that when left unrepaired leads to cell death. One major consequence of Holliday junction resolution by SLX4 and GEN1 is cancer-causing loss of heterozygosity due to mitotic crossing over. We measured mitotic crossovers in flies after a Cas9-induced chromosome break, and observed that this mutagenic form of repair is increased in the absence of Pol θ. This demonstrates that TMEJ can function upstream of the Holiday junction resolvases to protect cells from loss of heterozygosity. Our work argues that Pol θ can thus compensate for the loss of the Holliday junction resolvases by using homologous recombination intermediates, suppressing mitotic crossing over and preserving the genomic stability of cells.

Local Inversion Heterozygosity Alters Recombination throughout the Genome

Crossovers (COs) are formed during meiosis by the repair of programmed DNA double-strand breaks (DSBs) and are required for the proper segregation of chromosomes. More DSBs are made than COs, and the remaining DSBs are repaired as noncrossovers (NCOs). The distribution of recombination events along a chromosome occurs in a stereotyped pattern that is shaped by CO-promoting and CO-suppressing forces, collectively referred to as crossover patterning mechanisms. Chromosome inversions are structural aberrations that, when heterozygous, disrupt the recombination landscape by suppressing crossing over. In Drosophila species, the local suppression of COs by heterozygous inversions triggers an increase in crossing over on freely recombining chromosomes termed the interchromosomal (IC) effect [1, 2]. The molecular mechanism(s) by which heterozygous inversions suppress COs, whether noncrossover gene conversions (NCOGCs) are similarly affected, and what mediates the increase in COs in the rest of the genome remain open questions. By sequencing whole genomes of individual offspring from mothers containing heterozygous inversions, we show that, although COs are suppressed by inversions, NCOGCs occur throughout inversions at higher than wild-type frequencies. We confirm that CO frequency increases on the freely recombining chromosomes, yet CO interference remains intact. Intriguingly, NCOGCs do not increase in frequency on the freely recombining chromosomes and the total number of DSBs is approximately the same per genome. Together, our data show that heterozygous inversions change the recombination landscape by altering the relative proportions of COs and NCOGCs and suggest that DSB fate may be plastic until a CO assurance checkpoint has been satisfied.