Highly Contiguous Genome Assembly of Drosophila prolongata - a Model for Evolution of Sexual Dimorphism and Male-specific Innovations.

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2Mb) genome assembly of 220Mb. The repetitive content of the genome is 24.6%, the plurality of which are LTR retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, GWAS, and other genomic analyses.

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.

Highly Contiguous Genome Assembly of Drosophila prolongata - a Model for Evolution of Sexual Dimorphism and Male-specific Innovations.

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2Mb) genome assembly of 220Mb. The repetitive content of the genome is 24.6%, the plurality of which are LTR retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, GWAS, and other genomic analyses.

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.

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.

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.

A single-cell atlas of the sexually dimorphic Drosophila foreleg and its sensory organs during development.

To respond to the world around them, animals rely on the input of a network of sensory organs distributed throughout the body. Distinct classes of sensory organs are specialized for the detection of specific stimuli such as strain, pressure, or taste. The features that underlie this specialization relate both to the neurons that innervate sensory organs and the accessory cells they comprise. To understand the genetic basis of this diversity of cell types, both within and between sensory organs, we performed single-cell RNA sequencing on the first tarsal segment of the male Drosophila melanogaster foreleg during pupal development. This tissue displays a wide variety of functionally and structurally distinct sensory organs, including campaniform sensilla, mechanosensory bristles, and chemosensory taste bristles, as well as the sex comb, a recently evolved male-specific structure. In this study, we characterize the cellular landscape in which the sensory organs reside, identify a novel cell type that contributes to the construction of the neural lamella, and resolve the transcriptomic differences among support cells within and between sensory organs. We identify the genes that distinguish between mechanosensory and chemosensory neurons, resolve a combinatorial transcription factor code that defines 4 distinct classes of gustatory neurons and several types of mechanosensory neurons, and match the expression of sensory receptor genes to specific neuron classes. Collectively, our work identifies core genetic features of a variety of sensory organs and provides a rich, annotated resource for studying their development and function.

Evolution and development of male-specific leg brushes in Drosophilidae.

The origin, diversification, and secondary loss of sexually dimorphic characters are common in animal evolution. In some cases, structurally and functionally similar traits have evolved independently in multiple lineages. Prominent examples of such traits include the male-specific grasping structures that develop on the front legs of many dipteran insects. In this report, we describe the evolution and development of one of these structures, the male-specific "sex brush." The sex brush is composed of densely packed, irregularly arranged modified bristles and is found in several distantly related lineages in the family Drosophilidae. Phylogenetic analysis using 250 genes from over 200 species provides modest support for a single origin of the sex brush followed by many secondary losses; however, independent origins of the sex brush cannot be ruled out completely. We show that sex brushes develop in very similar ways in all brush-bearing lineages. The dense packing of brush hairs is explained by the specification of bristle precursor cells at a near-maximum density permitted by the lateral inhibition mechanism, as well as by the reduced size of the surrounding epithelial cells. In contrast to the female and the ancestral male condition, where bristles are arranged in stereotypical, precisely spaced rows, cell migration does not contribute appreciably to the formation of the sex brush. The complex phylogenetic history of the sex brush can make it a valuable model for investigating coevolution of sex-specific morphology and mating behavior.

Secondary reversion to sexual monomorphism associated with tissue-specific loss of doublesex expression.

Animal evolution is characterized by frequent turnover of sexually dimorphic traits-new sex-specific characters are gained, and some ancestral sex-specific characters are lost, in many lineages. In insects, sexual differentiation is predominantly cell autonomous and depends on the expression of the doublesex (dsx) transcription factor. In most cases, cells that transcribe dsx have the potential to undergo sex-specific differentiation, while those that lack dsx expression do not. Consistent with this mode of development, comparative research has shown that the origin of new sex-specific traits can be associated with the origin of new spatial domains of dsx expression. In this report, we examine the opposite situation-a secondary loss of the sex comb, a male-specific grasping structure that develops on the front legs of some drosophilid species. We show that while the origin of the sex comb is linked to an evolutionary gain of dsx expression in the leg, sex comb loss in a newly identified species of Lordiphosa (Drosophilidae) is associated with a secondary loss of dsx expression. We discuss how the developmental control of sexual dimorphism affects the mechanisms by which sex-specific traits can evolve.

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.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes.

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

DrosoPhyla: Resources for Drosophilid Phylogeny and Systematics.

The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4,400 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the genera, subgenera, and species groups in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred a species tree with highly supported groups for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.

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.

Evolution of sexual development and sexual dimorphism in insects.

Most animal species consist of two distinct sexes. At the morphological, physiological, and behavioral levels the differences between males and females are numerous and dramatic, yet at the genomic level they are often slight or absent. This disconnect is overcome because simple genetic differences or environmental signals are able to direct the sex-specific expression of a shared genome. A canonical picture of how this process works in insects emerged from decades of work on Drosophila. But recent years have seen an explosion of molecular-genetic and developmental work on a broad range of insects. Drawing these studies together, we describe the evolution of sexual dimorphism from a comparative perspective and argue that insect sex determination and differentiation systems are composites of rapidly evolving and highly conserved elements.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses.

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

A hierarchical Bayesian mixture model for inferring the expression state of genes in transcriptomes.

Transcriptomes are key to understanding the relationship between genotype and phenotype. The ability to infer the expression state (active or inactive) of genes in the transcriptome offers unique benefits for addressing this issue. For example, qualitative changes in gene expression may underly the origin of novel phenotypes, and expression states are readily comparable between tissues and species. However, inferring the expression state of genes is a surprisingly difficult problem, owing to the complex biological and technical processes that give rise to observed transcriptomic datasets. Here, we develop a hierarchical Bayesian mixture model that describes this complex process and allows us to infer expression state of genes from replicate transcriptomic libraries. We explore the statistical behavior of this method with analyses of simulated datasets-where we demonstrate its ability to correctly infer true (known) expression states-and empirical-benchmark datasets, where we demonstrate that the expression states inferred from RNA-sequencing (RNA-seq) datasets using our method are consistent with those based on independent evidence. The power of our method to correctly infer expression states is generally high and remarkably, approaches the maximum possible power for this inference problem. We present an empirical analysis of primate-brain transcriptomes, which identifies genes that have a unique expression state in humans. Our method is implemented in the freely available R package zigzag.

Evolution of sexually dimorphic pheromone profiles coincides with increased number of male-specific chemosensory organs in Drosophila prolongata.

Binary communication systems that involve sex-specific signaling and sex-specific signal perception play a key role in sexual selection and in the evolution of sexually dimorphic traits. The driving forces and genetic changes underlying such traits can be investigated in systems where sex-specific signaling and perception have emerged recently and show evidence of potential coevolution. A promising model is found in Drosophila prolongata, which exhibits a species-specific increase in the number of male chemosensory bristles. We show that this transition coincides with recent evolutionary changes in cuticular hydrocarbon (CHC) profiles. Long-chain CHCs that are sexually monomorphic in the closest relatives of D. prolongata (D. rhopaloa, D. carrolli, D. kurseongensis, and D. fuyamai) are strongly male-biased in this species. We also identify an intraspecific female-limited polymorphism, where some females have male-like CHC profiles. Both the origin of sexually dimorphic CHC profiles and the female-limited polymorphism in D. prolongata involve changes in the relative amounts of three mono-alkene homologs, 9-tricosene, 9-pentacosene, and 9-heptacosene, all of which share a common biosynthetic origin and point to a potentially simple genetic change underlying these traits. Our results suggest that pheromone synthesis may have coevolved with chemosensory perception and open the way for reconstructing the origin of sexual dimorphism in this communication system.

Sex-specific evolution of relative leg size in Drosophila prolongata results from changes in the intersegmental coordination of tissue growth.

Evolution of relative organ size is the most prolific source of morphological diversity, yet the underlying molecular mechanisms that modify growth control are largely unknown. Models where organ proportions have undergone recent evolutionary changes hold the greatest promise for understanding this process. Uniquely among Drosophila species, Drosophila prolongata displays a dramatic, male-specific increase in the size of its forelegs relative to other legs. By comparing leg development between males and females of D. prolongata and its closest relative Drosophila carrolli, we show that the exaggerated male forelegs are produced by a sex- and segment-specific increase in mitosis during the final larval instar. Intersegmental compensatory control, where smaller leg primordia grow at a faster rate, is observed in both species and sexes. However, the equilibrium growth rates that determine the final relative proportion between the first and second legs have shifted in male D. prolongata compared both to conspecific females and to D. carrolli. We suggest that the observed developmental changes that produce new adult proportions reflect an interplay between conserved growth coordination mechanisms and evolving organ-specific growth targets.

Hemimetabolous insects elucidate the origin of sexual development via alternative splicing.

Insects are the only known animals in which sexual differentiation is controlled by sex-specific splicing. The doublesex transcription factor produces distinct male and female isoforms, which are both essential for sex-specific development. dsx splicing depends on transformer, which is also alternatively spliced such that functional Tra is only present in females. This pathway has evolved from an ancestral mechanism where dsx was independent of tra and expressed and required only in males. To reconstruct this transition, we examined three basal, hemimetabolous insect orders: Hemiptera, Phthiraptera, and Blattodea. We show that tra and dsx have distinct functions in these insects, reflecting different stages in the changeover from a transcription-based to a splicing-based mode of sexual differentiation. We propose that the canonical insect tra-dsx pathway evolved via merger between expanding dsx function (from males to both sexes) and narrowing tra function (from a general splicing factor to dedicated regulator of dsx).