DrosoPhyla: Resources for Drosophilid Phylogeny and Systematics

Gene regulatory network co-option is sufficient to induce a morphological novelty in Drosophila

Identifying the molecular origins by which new morphological structures evolve is one of the long standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR/Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, mis-expression of Shavenbaby in the sheath of D. melanogaster , a naïve species that lacks these extensions, is sufficient to induce small trichomes. These induced extensions rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionarily novelty, our work shows that the trichome-forming network is flexible enough that it can be co-opted in a new context, and subsequently refined to produce unique apical projections that are barely recognizable compared to their simpler ancestral beginnings.

Mating dynamics of a sperm-limited drosophilid, Zaprionus indianus

When males have large sperm, they may become sperm limited and mating dynamics may be affected. One such species is Zaprionus indianus, a drosophilid that is an introduced pest species in the Americas. We examined aspects of mating behavior in Z. indianus to determine the senses necessary for mating and measure female and male remating habits. We found that vision is necessary for successful copulation, but wings, which produce courtship song, are not needed. Males need their foretarsi to successfully copulate and although the foretarsi may be needed for chemoreception, their role in hanging on to the female during copulation may be more important for successful mating. Females that mate once run out of sperm in approximately five days, although mating a second time greatly increases offspring production. Females do not seem to exert pre-mating choice among males with respect to mating with a familiar versus a novel male. Males are not capable of mating continuously and fail to produce offspring in many copulations. Overall, females of this species benefit from polyandry, providing an opportunity to study sexual selection in females. In addition, the dynamics of male competition for fertilizing eggs needs to be studied.

Skeleton phylogeny reconstructed with transcriptomes for the tribe Drosophilini (Diptera: Drosophilidae)

The family Drosophilidae is one of the most important model systems in evolutionary biology. Thanks to advances in high-throughput sequencing technology, a number of molecular phylogenetic analyses have been undertaken by using large data sets of many genes and many species sampled across this family. Especially, recent analyses using genome sequences have depicted the family-wide skeleton phylogeny with high confidence. However, the taxon sampling is still insufficient for minor lineages and non-Drosophila genera. In this study, we carried out phylogenetic analyses using a large number of transcriptome-based nucleotide sequences, focusing on the largest, core tribe Drosophilini in the Drosophilidae. In our analyses, some noise factors against phylogenetic reconstruction were taken into account by removing putative paralogy from the datasets and examining the effects of missing data, i.e. gene occupancy and site coverage, and incomplete lineage sorting. The inferred phylogeny has newly resolved the following phylogenetic positions/relationships at the genomic scale: (i) the monophyly of the subgenus Siphlodora including Zaprionus flavofasciatus to be transferred therein; (ii) the paraphyly of the robusta and melanica species groups within a clade comprised of the robusta, melanica and quadrisetata groups and Z. flavofasciatus; (iii) Drosophila curviceps (representing the curviceps group), D. annulipes (the quadrilineata subgroup of the immigrans group) and D. maculinotata clustered into a clade sister to the Idiomyia + Scaptomyza clade, forming together the expanded Hawaiian drosophilid lineage; (iv) Dichaetophora tenuicauda (representing the lineage comprised of the Zygothrica genus group and Dichaetophora) placed as the sister to the clade of the expanded Hawaiian drosophilid lineage and Siphlodora; and (v) relationships of the subgenus Drosophila and the genus Zaprionus as follows: (Zaprionus, (the quadrilineata subgroup, ((D. sternopleuralis, the immigrans group proper), (the quinaria radiation, the tripunctata radiation)))). These results are to be incorporated into the so-far published phylogenomic tree as a backbone (constraint) tree for grafting much more species based on sequences of a limited number of genes. Such a comprehensive, highly confident phylogenetic tree with extensive and dense taxon sampling will provide an essential framework for comparative studies of the Drosophilidae.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has since followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appears to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, Sex Peptide Receptor, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appear to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, SPR, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Unistrand piRNA clusters are an evolutionarily conserved mechanism to suppress endogenous retroviruses across the Drosophila genus

The PIWI-interacting RNA (piRNA) pathway prevents endogenous genomic parasites, i.e. transposable elements, from damaging the genetic material of animal gonadal cells. Specific regions in the genome, called piRNA clusters, are thought to define each species' piRNA repertoire and therefore its capacity to recognize and silence specific transposon families. The unistrand cluster flamenco (flam) is essential in the somatic compartment of the Drosophila ovary to restrict Gypsy-family transposons from infecting the neighbouring germ cells. Disruption of flam results in transposon de-repression and sterility, yet it remains unknown whether this silencing mechanism is present more widely. Here, we systematically characterise 119 Drosophila species and identify five additional flam-like clusters separated by up to 45 million years of evolution. Small RNA-sequencing validated these as bona-fide unistrand piRNA clusters expressed in somatic cells of the ovary, where they selectively target transposons of the Gypsy family. Together, our study provides compelling evidence of a widely conserved transposon silencing mechanism that co-evolved with virus-like Gypsy-family transposons.

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.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Evolution of a fatty acyl-CoA elongase underlies desert adaptation in Drosophila

Traits that allow species to survive in extreme environments such as hot-arid deserts have independently evolved in multiple taxa. However, the genetic and evolutionary mechanisms underlying these traits have thus far not been elucidated. Here, we show that Drosophila mojavensis, a desert-adapted fruit fly species, has evolved high desiccation resistance by producing long-chain methyl-branched cuticular hydrocarbons (mbCHCs) that contribute to a cuticular lipid layer reducing water loss. We show that the ability to synthesize these longer mbCHCs is due to evolutionary changes in a fatty acyl-CoA elongase (mElo). mElo knockout in D. mojavensis led to loss of longer mbCHCs and reduction of desiccation resistance at high temperatures but did not affect mortality at either high temperatures or desiccating conditions individually. Phylogenetic analysis showed that mElo is a Drosophila-specific gene, suggesting that while the physiological mechanisms underlying desert adaptation may be similar between species, the genes involved in these mechanisms may be species or lineage specific.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila , many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila . We first inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of these species could develop on diet with toxin. Third, we examined fly development on a diet containing a natural mix of toxins extracted from the Death Cap Amanita phalloides mushroom. Both tolerant and susceptible species developed on diet with this mix, though tolerant species survived at significantly higher concentrations. Finally, we asked how cyclopeptide tolerance might have evolved across the immigrans-tripunctata radiation and inferred that toxin tolerance was ancestral and subsequently lost multiple times. Our results suggest the evolutionary history of cyclopeptide tolerance is complex, and simply describing this trait as present or absent does not fully capture the occurrence or impact on this adaptive radiation. More broadly, the evolution of novelty can be more complex than previously thought, and that accurate descriptions of such novelties are critical in studies examining their evolution.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

Desiccation resistance differences in Drosophila species can be largely explained by variations in cuticular hydrocarbons

Maintaining water balance is a universal challenge for organisms living in terrestrial environments, especially for insects, which have essential roles in our ecosystem. Although the high surface area to volume ratio in insects makes them vulnerable to water loss, insects have evolved different levels of desiccation resistance to adapt to diverse environments. To withstand desiccation, insects use a lipid layer called cuticular hydrocarbons (CHCs) to reduce water evaporation from the body surface. It has long been hypothesized that the water-proofing capability of this CHC layer, which can confer different levels of desiccation resistance, depends on its chemical composition. However, it is unknown which CHC components are important contributors to desiccation resistance and how these components can determine differences in desiccation resistance. In this study, we used machine-learning algorithms, correlation analyses, and synthetic CHCs to investigate how different CHC components affect desiccation resistance in 50 Drosophila and related species. We showed that desiccation resistance differences across these species can be largely explained by variation in CHC composition. In particular, length variation in a subset of CHCs, the methyl-branched CHCs (mbCHCs), is a key determinant of desiccation resistance. There is also a significant correlation between the evolution of longer mbCHCs and higher desiccation resistance in these species. Given that CHCs are almost ubiquitous in insects, we suggest that evolutionary changes in insect CHC components can be a general mechanism for the evolution of desiccation resistance and adaptation to diverse and changing environments.

Trichomes on female reproductive tract: rapid diversification and underlying gene regulatory network in Drosophila suzukii and its related species

Background

The ovipositors of some insects are external female genitalia, which have their primary function to deliver eggs. Drosophila suzukii and its sibling species D. subpulchrella are known to have acquired highly sclerotized and enlarged ovipositors upon their shifts in oviposition sites from rotting to ripening fruits. Inside the ovipositor plates, there are scale-like polarized protrusions termed “oviprovector scales” that are likely to aid the mechanical movement of the eggs. The size and spatial distribution of the scales need to be rearranged following the divergence of the ovipositors. In this study, we examined the features of the oviprovector scales in D. suzukii and its closely related species. We also investigated whether the scales are single-cell protrusions comprised of F-actin under the same conserved gene regulatory network as the well-characterized trichomes on the larval cuticular surface.

Results

The oviprovector scales of D. suzukii and D. subpulchrella were distinct in size and spatial arrangement compared to those of D. biarmipes and other closely related species. The scale numbers also varied greatly among these species. The comparisons of the size of the scales suggested a possibility that the apical cell area of the oviprovector has expanded upon the elongation of the ovipositor plates in these species. Our transcriptome analysis revealed that 43 out of the 46 genes known to be involved in the trichome gene regulatory network are expressed in the developing female genitalia of D. suzukii and D. subpulchrella. The presence of Shavenbaby (Svb) or svb was detected in the inner cavity of the developing ovipositors of D. melanogaster, D. suzukii, and D. subpulchrella. Also, shavenoid (sha) was expressed in the corresponding patterns in the developing ovipositors and showed differential expression levels between D. suzukii and D. subpulchrella at 48 h APF.

Conclusions

The oviprovector scales have divergent size and spatial arrangements among species. Therefore, these scales may represent a rapidly diversifying morphological trait of the female reproductive tract reflecting ecological contexts. Furthermore, our results showed that the gene regulatory network underlying trichome formation is also utilized to develop the rapidly evolving trichomes on the oviprovectors of these flies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02046-1.

Evolution and genomic basis of the plant-penetrating ovipositor: a key morphological trait in herbivorous Drosophilidae

Herbivorous insects are extraordinarily diverse, yet are found in only one-third of insect orders. This skew may result from barriers to plant colonization, coupled with phylogenetic constraint on plant-colonizing adaptations. The plant-penetrating ovipositor, however, is one trait that surmounts host plant physical defences and may be evolutionarily labile. Ovipositors densely lined with hard bristles have evolved repeatedly in herbivorous lineages, including within the Drosophilidae. However, the evolution and genetic basis of this innovation has not been well studied. Here, we focused on the evolution of this trait in Scaptomyza, a genus sister to Hawaiian Drosophila, that contains a herbivorous clade. Our phylogenetic approach revealed that ovipositor bristle number increased as herbivory evolved in the Scaptomyza lineage. Through a genome-wide association study, we then dissected the genomic architecture of variation in ovipositor bristle number within S. flava. Top-associated variants were enriched for transcriptional repressors, and the strongest associations included genes contributing to peripheral nervous system development. Individual genotyping supported the association at a variant upstream of Gαi, a neural development gene, contributing to a gain of 0.58 bristles/major allele. These results suggest that regulatory variation involving conserved developmental genes contributes to this key morphological trait involved in plant colonization.

Phylogenetic resolution of the fly superfamily Ephydroidea-Molecular systematics of the enigmatic and diverse relatives of Drosophilidae

The schizophoran superfamily Ephydroidea (Diptera: Cyclorrhapha) includes eight families, ranging from the well-known vinegar flies (Drosophilidae) and shore flies (Ephydridae), to several small, relatively unusual groups, the phylogenetic placement of which has been particularly challenging for systematists. An extraordinary diversity in life histories, feeding habits and morphology are a hallmark of fly biology, and the Ephydroidea are no exception. Extreme specialization can lead to "orphaned" taxa with no clear evidence for their phylogenetic position. To resolve relationships among a diverse sample of Ephydroidea, including the highly modified flies in the families Braulidae and Mormotomyiidae, we conducted phylogenomic sampling. Using exon capture from Anchored Hybrid Enrichment and transcriptomics to obtain 320 orthologous nuclear genes sampled for 32 species of Ephydroidea and 11 outgroups, we evaluate a new phylogenetic hypothesis for representatives of the superfamily. These data strongly support monophyly of Ephydroidea with Ephydridae as an early branching radiation and the placement of Mormotomyiidae as a family-level lineage sister to all remaining families. We confirm placement of Cryptochetidae as sister taxon to a large clade containing both Drosophilidae and Braulidae-the latter a family of honeybee ectoparasites. Our results reaffirm that sampling of both taxa and characters is critical in hyperdiverse clades and that these factors have a major influence on phylogenomic reconstruction of the history of the schizophoran fly radiation.

Interspecific variation in sex-specific gustatory organs in Drosophila

Drosophila males use leg gustatory bristles to discriminate between male and female cuticular pheromones as an important part of courtship behavior. In Drosophila melanogaster, several male-specific gustatory bristles are present on the anterior surface of the first tarsal segment of the prothoracic leg, in addition to a larger set of gustatory bristles found in both sexes. These bristles are thought to be specialized for pheromone detection. Here, we report the number and location of sex-specific gustatory bristles in 27 other Drosophila species. Although some species have a pattern similar to D. melanogaster, others lack anterior male-specific bristles but have many dorsal male-specific gustatory bristles instead. Some species have both anterior and dorsal male-specific bristles, while others lack sexual dimorphism entirely. In several distantly related species, the number of gustatory bristles is much greater in males than in females due to a male-specific transformation of ancestrally mechanosensory bristles to a chemosensory identity. This variation in the extent and pattern of sexual dimorphism may affect the formation and function of neuronal circuits that control Drosophila courtship and contribute to the evolution of mating behavior.

Hematopoietic plasticity mapped in Drosophila and other insects

Hemocytes, similar to vertebrate blood cells, play important roles in insect development and immunity, but it is not well understood how they perform their tasks. New technology, in particular single-cell transcriptomic analysis in combination with Drosophila genetics, may now change this picture. This review aims to make sense of recently published data, focusing on Drosophila melanogaster and comparing to data from other drosophilids, the malaria mosquito, Anopheles gambiae, and the silkworm, Bombyx mori. Basically, the new data support the presence of a few major classes of hemocytes: (1) a highly heterogenous and plastic class of professional phagocytes with many functions, called plasmatocytes in Drosophila and granular cells in other insects. (2) A conserved class of cells that control melanin deposition around parasites and wounds, called crystal cells in D. melanogaster, and oenocytoids in other insects. (3) A new class of cells, the primocytes, so far only identified in D. melanogaster. They are related to cells of the so-called posterior signaling center of the larval hematopoietic organ, which controls the hematopoiesis of other hemocytes. (4) Different kinds of specialized cells, like the lamellocytes in D. melanogaster, for the encapsulation of parasites. These cells undergo rapid evolution, and the homology relationships between such cells in different insects are uncertain. Lists of genes expressed in the different hemocyte classes now provide a solid ground for further investigation of function.

Evolution of the Colocasiomyia gigantea Species Group (Diptera: Drosophilidae): Phylogeny, Biogeography and Shift of Host Use

Simple Summary

All the species in the Colocasiomyia gigantea group breed on monsteroid host plants (aroids in the subfamily Monsteroideae). So far, we have not resolved the phylogenetic relationship among these fly species, making it difficult to trace the origin and history of the species diversification, biogeography and host plant selection. In this study, we reconstructed the evolutionary relationships between these species using multilocus DNA sequence data, and we inferred their ancestral areas and host plants. According to the results, this group diverged from its sister taxon through a split between the northeastern Oriental region and Sundaland + Wallacea, with the subsequent diversification occurring largely in the first region. We inferred the most likely ancestral host genus of this group to be Rhaphidophora Hassk, with possible subsequent shifts to Scindapsus Schott and/or Epipremnum Schott. We discuss the potential of the group as a model system for studies in evolutionary ecology and developmental genetics.

Abstract

The gigantea species group of the genus Colocasiomyia de Meijere (Diptera: Drosophilidae) is among the four aroid-breeding species groups in this genus; however, it differs from the remaining three groups in the host use: all the flies in this group use plants from the subfamily Monsteroideae instead of from the subfamily Aroideae. So far, we have not resolved the phylogenetic relationship within this group, making it difficult to trace its geographical origin, pattern of species diversification and history of host plant use. In this study, we reconstructed the phylogenetic relationships within the C. gigantea group using DNA sequences of eight (two mitochondrial and six nuclear) gene markers, and we inferred the ancestral areas and host plants of the group based on the resulting phylogeny. According to the results, the C. gigantea group may have diverged from its sister group (i.e., the C. cristata group) through vicariance between the northeastern Oriental region and Sundaland + Wallacea, and the subsequent diversification of the C. gigantea group occurred mostly in the northeastern Oriental region, although an Oriental-to-Sundaland dispersal was followed by vicariance between these two areas, which finally gave rise to the C. gigantea-C. scindapsae lineage in the latter area. We inferred the most likely ancestral host plant of the C. gigantea group to be of the genus Rhaphidophora Hassk, with possible subsequent shifts to Scindapsus Schott and/or Epipremnum Schott plants. We discuss the potential for the egg filaments in the C. gigantea group to be used as a model system for comparative studies in pollination mutualism and developmental genetics concerning tubulogenesis.

Dietary Utilization Drives the Differentiation of Gut Bacterial Communities between Specialist and Generalist Drosophilid Flies

Gut bacteria play vital roles in the dietary detoxification, digestion, and nutrient supplementation of hosts during dietary specialization. The roles of gut bacteria in the host can be unveiled by comparing communities of specialist and generalist bacterial species. However, these species usually have a long evolutionary history, making it difficult to determine whether bacterial community differentiation is due to host dietary adaptation or phylogenetic divergence. In this regard, we investigated the bacterial communities from two Araceae-feeding Colocasiomyia species and further performed a meta-analysis by incorporating the published data from Drosophila bacterial community studies. The compositional and functional differentiation of bacterial communities was uncovered by comparing three (Araceae-feeding, mycophagous, and cactophilic) specialists with generalist flies. The compositional differentiation showed that Bacteroidetes and Firmicutes inhabited specialists, while more Proteobacteria lived in generalists. The functional prediction based on the bacterial community compositions suggested that amino acid metabolism and energy metabolism are overrepresented pathways in specialists and generalists, respectively. The differences were mainly associated with the higher utilization of structural complex carbohydrates, protein utilization, vitamin B₁₂ acquisition, and demand for detoxification in specialists than in generalists. The complementary roles of bacteria reveal a connection between gut bacterial communities and fly dietary specialization.

IMPORTANCE Gut bacteria may play roles in the dietary utilization of hosts, especially in specialist animals, during long-term host-microbe interaction. By comparing the gut bacterial communities between specialist and generalist drosophilid flies, we found that specialists harbor more bacteria linked to complex carbohydrate degradation, amino acid metabolism, vitamin B₁₂ formation, and detoxification than do generalists. This study reveals the roles of gut bacteria in drosophilid species in dietary utilization.

Biology and ecology of the Oriental flower-breeding Drosophila elegans and related species

Animals adapt to their environments in the course of evolution. One effective approach to elucidate mechanisms of adaptive evolution is to compare closely related species with model organisms in which knowledge of the molecular and physiological bases of various traits has been accumulated. Drosophila elegans and its close relatives, belonging to the same species group as the model organism D. melanogaster, exhibit various unique characteristics such as flower-breeding habit, courtship display, territoriality, sexual dimorphism, and colour polymorphism. Their ease of culturing and availability of genomic information makes them a useful model for understanding mechanisms of adaptive evolution. Here, we review the morphology, distribution, and phylogenetic relationships of D. elegans and related species, as well as their characteristic flower-dependent biology, food habits, and life-history traits. We also describe their unique mating and territorial behaviours and note their distinctive karyotype and the genetic mechanisms of morphological diversity that have recently been revealed.

The Physiological and Evolutionary Ecology of Sperm Thermal Performance

Ongoing anthropogenic climate change has increased attention on the ecological and evolutionary consequences of thermal variation. Most research in this field has focused on the physiology and behavior of diploid whole organisms. The thermal performance of haploid gamete stages directly tied to reproductive success has received comparatively little attention, especially in the context of the evolutionary ecology of wild (i.e., not domesticated) organisms. Here, we review evidence for the effects of temperature on sperm phenotypes, emphasizing data from wild organisms whenever possible. We find that temperature effects on sperm are pervasive, and that above normal temperatures in particular are detrimental. That said, there is evidence that sperm traits can evolve adaptively in response to temperature change, and that adaptive phenotypic plasticity in sperm traits is also possible. We place results in the context of thermal performance curves, and encourage this framework to be used as a guide for experimental design to maximize ecological relevance as well as the comparability of results across studies. We also highlight gaps in our understanding of sperm thermal performance that require attention to more fully understand thermal adaptation and the consequences of global change.

Interpopulation variation of transposable elements of the hAT superfamily in Drosophila willistoni (Diptera: Drosophilidae): in-situ approach

Transposable elements are abundant and dynamic part of the genome, influencing organisms in different ways through their presence or mobilization, or by acting directly on pre- and post-transcriptional regulatory regions. We compared and evaluated the presence, structure, and copy number of three hAT superfamily transposons (hobo, BuT2, and mar) in five strains of Drosophila willistoni species. These D. willistoni strains are of different geographical origins, sampled across the north-south occurrence of this species. We used sequenced clones of the hAT elements in fluorescence in-situ hybridizations in the polytene chromosomes of three strains of D. willistoni. We also analyzed the structural characteristics and number of copies of these hAT elements in the 10 currently available sequenced genomes of the willistoni group. We found that hobo, BuT2, and mar were widely distributed in D. willistoni polytene chromosomes and sequenced genomes of the willistoni group, except for mar, which is restricted to the subgroup willistoni. Furthermore, the elements hobo, BuT2, and mar have different evolutionary histories. The transposon differences among D. willistoni strains, such as variation in the number, structure, and chromosomal distribution of hAT transposons, could reflect the genomic and chromosomal plasticity of D. willistoni species in adapting to highly variable environments.

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high‐quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome‐wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert‐adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate‐continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics‐restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large‐scale and multifaceted genomic changes.

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

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

Repression precedes independent evolutionary gains of a highly specific gene expression pattern

Highly specific expression patterns can be caused by the overlapping activities of activator and repressor sequences in enhancers. However, few studies illuminate how these sequences evolve in the origin of new enhancers. Here, we show that expression of the bond gene in the semicircular wall epithelium (swe) of the Drosophila melanogaster male ejaculatory bulb (EB) is controlled by an enhancer consisting of an activator region that requires Abdominal-B driving expression in the entire EB and a repressor region that restricts this expression to the EB swe. Although this expression pattern is independently gained in the distantly related Scaptodrosophila lebanonensis and does not require Abdominal-B, we show that functionally similar repressor sequences are present in Scaptodrosophila and also in species that do not express bond in the EB. We suggest that during enhancer evolution, repressor sequences can precede the evolution of activator sequences and may lead to similar but independently evolved expression patterns.

Pu et al. show that the independent gain of a highly specific expression pattern across distantly related species may be because of the preexistence of repressor sequences that precedes the diversification of these species. This may reflect a general mechanism underlying the evolution of highly specific enhancers.

Highly contiguous assemblies of 101 drosophilid genomes

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.