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.

The Transmission Patterns of the Endosymbiont Wolbachia within the Hawaiian Drosophilidae Adaptive Radiation.

The evolution of endosymbionts and their hosts can lead to highly dynamic interactions with varying fitness effects for both the endosymbiont and host species. Wolbachia, a ubiquitous endosymbiont of arthropods and nematodes, can have both beneficial and detrimental effects on host fitness. We documented the occurrence and patterns of transmission of Wolbachia within the Hawaiian Drosophilidae and examined the potential contributions of Wolbachia to the rapid diversification of their hosts. Screens for Wolbachia infections across a minimum of 140 species of Hawaiian Drosophila and Scaptomyza revealed species-level infections of 20.0%, and across all 399 samples, a general infection rate of 10.3%. Among the 44 Wolbachia strains we identified using a modified Wolbachia multi-locus strain typing scheme, 30 (68.18%) belonged to supergroup B, five (11.36%) belonged to supergroup A, and nine (20.45%) had alleles with conflicting supergroup assignments. Co-phylogenetic reconciliation analysis indicated that Wolbachia strain diversity within their endemic Hawaiian Drosophilidae hosts can be explained by vertical (e.g., co-speciation) and horizontal (e.g., host switch) modes of transmission. Results from stochastic character trait mapping suggest that horizontal transmission is associated with the preferred oviposition substrate of the host, but not the host's plant family or island of occurrence. For Hawaiian Drosophilid species of conservation concern, with 13 species listed as endangered and 1 listed as threatened, knowledge of Wolbachia strain types, infection status, and potential for superinfection could assist with conservation breeding programs designed to bolster population sizes, especially when wild populations are supplemented with laboratory-reared, translocated individuals. Future research aimed at improving the understanding of the mechanisms of Wolbachia transmission in nature, their impact on the host, and their role in host species formation may shed light on the influence of Wolbachia as an evolutionary driver, especially in Hawaiian ecosystems.

Reproductive Capacity Evolves in Response to Ecology through Common Changes in Cell Number in Hawaiian Drosophila.

Lifetime reproductive capacity is a critical fitness component. In insects, female reproductive capacity is largely determined by the number of ovarioles, the egg-producing subunits of the ovary [e.g., 1]. Recent work has provided insights into ovariole number regulation in Drosophila melanogaster. However, whether mechanisms discovered under laboratory conditions explain evolutionary variation in natural populations is an outstanding question. We investigated potential effects of ecology on the developmental processes underlying ovariole number evolution among Hawaiian Drosophila, a large adaptive radiation wherein the highest and lowest ovariole numbers of the family have evolved within 25 million years. Previous studies proposed that ovariole number correlated with oviposition substrate [2-4] but sampled largely one clade of these flies and were limited by a provisional phylogeny and the available comparative methods. We test this hypothesis by applying phylogenetic modeling to an expanded sampling of ovariole numbers and substrate types and show support for these predictions across all major groups of Hawaiian Drosophila, wherein ovariole number variation is best explained by adaptation to specific substrates. Furthermore, we show that oviposition substrate evolution is linked to changes in the allometric relationship between body size and ovariole number. Finally, we provide evidence that the major changes in ovarian cell number that regulate D. melanogaster ovariole number also regulate ovariole number in Hawaiian drosophilids. Thus, we provide evidence that this remarkable adaptive radiation is linked to evolutionary changes in a key reproductive trait regulated at least partly by variation in the same developmental parameters that operate in the model species D. melanogaster.

The Native Hawaiian Insect Microbiome Initiative: A Critical Perspective for Hawaiian Insect Evolution.

Insects associate with a diversity of microbes that can shape host ecology and diversity by providing essential biological and adaptive services. For most insect groups, the evolutionary implications of host-microbe interactions remain poorly understood. Geographically discrete areas with high biodiversity offer powerful, simplified model systems to better understand insect-microbe interactions. Hawaii boasts a diverse endemic insect fauna (~6000 species) characterized by spectacular adaptive radiations. Despite this, little is known about the role of bacteria in shaping this diversity. To address this knowledge gap, we inaugurate the Native Hawaiian Insect Microbiome Initiative (NHIMI). The NHIMI is an effort intended to develop a framework for informing evolutionary and biological studies in Hawaii. To initiate this effort, we have sequenced the bacterial microbiomes of thirteen species representing iconic, endemic Hawaiian insect groups. Our results show that native Hawaiian insects associate with a diversity of bacteria that exhibit a wide phylogenetic breadth. Several groups show predictable associations with obligate microbes that permit diet specialization. Others exhibit unique ecological transitions that are correlated with shifts in their microbiomes (e.g., transition to carrion feeding from plant-feeding in Nysius wekiuicola). Finally, some groups, such as the Hawaiian Drosophila, have relatively diverse microbiomes with a conserved core of bacterial taxa across multiple species and islands.

Rapid adaptive radiation and host plant conservation in the Hawaiian picture wing Drosophila (Diptera: Drosophilidae).

The Hawaiian picture wing Drosophila are a striking example of adaptive radiation in specialist saprophages on an island system. We use DNA sequences from five nuclear genes with a total of 4260 nucleotides to provide a comprehensive phylogeny and biogeographic analysis of 90 species in the Hawaiian Drosophila picture wing clade. The current analysis indicates that the evolution of the picture wing clade took place more recently than previously suggested. The relationships of several morphologically anomalous taxa are resolved with strong support. Biogeography and host plant analyses show two periods of rapid divergence occurred when Kauai and Oahu were the main high islands, indicating that a combination of complex topographical features of islands and development of novel host plant associations was key to the rapid diversification of these lineages. For the past 2 million years, host associations within lineages have been largely stable, and speciation has occurred primarily due to the establishment of populations on newer islands as they arose followed by divergence by isolation. The existence of several apparently relictual taxa suggests that extinction has also played a major role in assembly of the present Hawaiian Drosophila fauna.

Phylogenetics of the antopocerus-modified tarsus clade of Hawaiian Drosophila: diversification across the Hawaiian Islands.

The Hawaiian Drosophilidae radiation is an ecologically and morphologically diverse clade of almost 700 described species. A phylogenetic approach is key to understanding the evolutionary forces that have given rise to this diverse lineage. Here we infer the phylogeny for the antopocerus, modified tarsus and ciliated tarsus (AMC) clade, a lineage comprising 16% (91 of 687 species) of the described Hawaiian Drosophilidae. To improve on previous analyses we constructed the largest dataset to date for the AMC, including a matrix of 15 genes for 68 species. Results strongly support most of the morphologically defined species groups as monophyletic. We explore the correlation of increased diversity in biogeography, sexual selection and ecology on the present day diversity seen in this lineage using a combination of dating methods, rearing records, and distributional data. Molecular dating analyses indicate that AMC lineage started diversifying about 4.4 million years ago, culminating in the present day AMC diversity. We do not find evidence that ecological speciation or sexual selection played a part in generating this diversity, but given the limited number of described larval substrates and secondary sexual characters analyzed we can not rule these factors out entirely. An increased rate of diversification in the AMC is found to overlap with the emergence of multiple islands in the current chain of high islands, specifically Oahu and Kauai.

DNA barcoding a regional fauna: Irish solitary bees.

As the globally dominant group of pollinators, bees provide a key ecosystem service for natural and agricultural landscapes. Their corresponding global decline thus poses an important threat to plant populations and the ecosystems they support. Bee conservation requires rapid and effective tools to identify and delineate species. Here, we apply DNA barcoding to Irish solitary bees as the first step towards a DNA barcode library for European solitary bees. Using the standard barcoding sequence, we were able to identify 51 of 55 species. Potential problems included a suite of species in the genus Andrena, which were recalcitrant to sequencing, mitochondrial heteroplasmy and parasitic flies, which led to the production of erroneous sequences from DNA extracts. DNA barcoding enabled the assignment of morphologically unidentifiable females of the parasitic genus Sphecodes to their nominal taxa. It also enabled correction of the Irish bee list for morphologically inaccurately identified specimens. However, the standard COI barcode was unable to differentiate the recently diverged taxa Sphecodes ferruginatus and S. hyalinatus. Overall, our results show that DNA barcoding provides an excellent identification tool for Irish solitary bees and should be rolled out to provide a database for solitary bees globally.

Mitochondrial heteroplasmy and DNA barcoding in Hawaiian Hylaeus (Nesoprosopis) bees (Hymenoptera: Colletidae).

BACKGROUND: The past several years have seen a flurry of papers seeking to clarify the utility and limits of DNA barcoding, particularly in areas such as species discovery and paralogy due to nuclear pseudogenes. Heteroplasmy, the coexistence of multiple mitochondrial haplotypes in a single organism, has been cited as a potentially serious problem for DNA barcoding but its effect on identification accuracy has not been tested. In addition, few studies of barcoding have tested a large group of closely-related species with a well-established morphological taxonomy. In this study we examine both of these issues, by densely sampling the Hawaiian Hylaeus bee radiation. RESULTS: Individuals from 21 of the 49 a priori morphologically-defined species exhibited coding sequence heteroplasmy at levels of 1-6% or more. All homoplasmic species were successfully identified by COI using standard methods of analysis, but only 71% of heteroplasmic species. The success rate in identifying heteroplasmic species was increased to 86% by treating polymorphisms as character states rather than ambiguities. Nuclear pseudogenes (numts) were also present in four species, and were distinguishable from heteroplasmic sequences by patterns of nucleotide and amino acid change. CONCLUSIONS: Heteroplasmy significantly decreased the reliability of species identification. In addition, the practical issue of dealing with large numbers of polymorphisms- and resulting increased time and labor required - makes the development of DNA barcode databases considerably more complex than has previously been suggested. The impact of heteroplasmy on the utility of DNA barcoding as a bulk specimen identification tool will depend upon its frequency across populations, which remains unknown. However, DNA barcoding is still likely to remain an important identification tool for those species that are difficult or impossible to identify through morphology, as is the case for the ecologically important solitary bee fauna.

Tissue segregation of mitochondrial haplotypes in heteroplasmic Hawaiian bees: implications for DNA barcoding.

The issue of mitochondrial heteroplasmy has been cited as a theoretical problem for DNA barcoding but is only beginning to be examined in natural systems. We sequenced multiple DNA extractions from 20 individuals of four Hawaiian Hylaeus bee species known to be heteroplasmic. All species showed strong differences at polymorphic sites between abdominal and muscle tissue in most individuals, and only two individuals had no obvious segregation. Two specimens produced completely clean sequences from abdominal DNA. The fact that these differences are clearly visible by direct sequencing indicates that substantial intra-individual mtDNA diversity may be overlooked when DNA is taken from small tissue fragments. At the same time, differences in haplotype distribution among individuals may result in incorrect recognition of cryptic species. Because DNA barcoding studies typically use only a small fragment of an organism, they are particularly vulnerable to sequencing bias where heteroplasmy and haplotype segregation are present. It is important to anticipate this possibility prior to undertaking large-scale barcoding projects to reduce the likelihood of haplotype segregation confounding the results.

Low nuclear DNA variation supports a recent origin of Hawaiian Hylaeus bees (Hymenoptera: Colletidae).

Previous phylogenetic work on the Hawaiian bees of the genus Hylaeus, based on mitochondrial DNA and morphology, appeared to support a recent origin for the group, but support for the resulting tree was weak. Four nuclear genes with varying evolutionary rates -- arginine kinase, EF-1alpha, opsin, and wingless -- were sequenced for a reduced taxon set in an attempt to find one or more data set that would provide better support. All showed very low variation (<2%) in the ingroup. Comparison among genes revealed a much higher than expected rate of evolution in mtDNA, especially at first and second positions. While the data from the nuclear genes showed insufficient variation for phylogenetic analysis, the strong sequence similarity among the Hawaiian species supports the previous hypothesis of a recent origin for the group.