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.

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.

The evolutionary trajectory of drosophilid walking.

Neural circuits must both execute the behavioral repertoire of individuals and account for behavioral variation across species. Understanding how this variation emerges over evolutionary time requires large-scale phylogenetic comparisons of behavioral repertoires. Here, we describe the evolution of walking in fruit flies by capturing high-resolution, unconstrained movement from 13 species and 15 strains of drosophilids. We find that walking can be captured in a universal behavior space, the structure of which is evolutionarily conserved. However, the occurrence of and transitions between specific movements have evolved rapidly, resulting in repeated convergent evolution in the temporal structure of locomotion. Moreover, a meta-analysis demonstrates that many behaviors evolve more rapidly than other traits. Thus, the architecture and physiology of locomotor circuits can execute precise individual movements in one species and simultaneously support rapid evolutionary changes in the temporal ordering of these modular elements across clades.

Deep Ancestral Introgression Shapes Evolutionary History of Dragonflies and Damselflies.

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here, we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, nonmodel insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression. [Gene flow; Odonata; phylogenomics; reticulate evolution.].

Widespread introgression across a phylogeny of 155 Drosophila genomes.

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

Molecular Evolution of Phototransduction Pathway Genes in Nocturnal and Diurnal Fireflies (Coleoptera: Lampyridae).

Most organisms are dependent on sensory cues from their environment for survival and reproduction. Fireflies (Coleoptera: Lampyridae) represent an ideal system for studying sensory niche adaptation due to many species relying on bioluminescent communication; as well as a diversity of ecologies. Here; using transcriptomics; we examine the phototransduction pathway in this non-model organism; and provide some of the first evidence for positive selection in the phototransduction pathway beyond opsins in beetles. Evidence for gene duplications within Lampyridae are found in inactivation no afterpotential C and inactivation no afterpotential D. We also find strong support for positive selection in arrestin-2; inactivation no afterpotential D; and transient receptor potential-like; with weak support for positive selection in guanine nucleotide-binding protein G(q) subunit alpha and neither inactivation nor afterpotential C. Taken with other recent work in flies; butterflies; and moths; this represents an exciting new avenue of study as we seek to further understand diversification and constraint on the phototransduction pathway in light of organism ecology.

Phylogeny and classification of Odonata using targeted genomics.

Dragonflies and damselflies are a charismatic, medium-sized insect order (~6300 species) with a unique potential to approach comparative research questions. Their taxonomy and many ecological traits for a large fraction of extant species are relatively well understood. However, until now, the lack of a large-scale phylogeny based on high throughput data with the potential to connect both perspectives has precluded comparative evolutionary questions for these insects. Here, we provide an ordinal hypothesis of classification based on anchored hybrid enrichment using a total of 136 species representing 46 of the 48 families or incertae sedis, and a total of 478 target loci. Our analyses recovered the monophyly for all three suborders: Anisoptera, Anisozygoptera and Zygoptera. Although the backbone of the topology was reinforced and showed the highest support values to date, our genomic data was unable to stronglyresolve portions of the topology. In addition, a quartet sampling approach highlights the potential evolutionary scenarios that may have shaped evolutionary phylogeny (e.g., incomplete lineage sorting and introgression) of this taxon. Finally, in light of our phylogenomic reconstruction and previous morphological and molecular information we proposed an updated odonate classification and define five new families (Amanipodagrionidae fam. nov., Mesagrionidae fam. nov., Mesopodagrionidae fam. nov., Priscagrionidae fam. nov., Protolestidae fam. nov.) and reinstate another two (Rhipidolestidae stat. res., Tatocnemididae stat. res.). Additionally, we feature the problematic taxonomic groupings for examination in future studies to improve our current phylogenetic hypothesis.

Genome analysis of Plectus murrayi, a nematode from continental Antarctica.

Plectus murrayi is one of the most common and locally abundant invertebrates of continental Antarctic ecosystems. Because it is readily cultured on artificial medium in the laboratory and highly tolerant to an extremely harsh environment, P. murrayi is emerging as a model organism for understanding the evolutionary origin and maintenance of adaptive responses to multiple environmental stressors, including freezing and desiccation. The de novo assembled genome of P. murrayi contains 225.741 million base pairs and a total of 14,689 predicted genes. Compared to Caenorhabditis elegans, the architectural components of P. murrayi are characterized by a lower number of protein-coding genes, fewer transposable elements, but more exons, than closely related taxa from less harsh environments. We compared the transcriptomes of lab-reared P. murrayi with wild-caught P. murrayi and found genes involved in growth and cellular processing were up-regulated in lab-cultured P. murrayi, while a few genes associated with cellular metabolism and freeze tolerance were expressed at relatively lower levels. Preliminary comparative genomic and transcriptomic analyses suggest that the observed constraints on P. murrayi genome architecture and functional gene expression, including genome decay and intron retention, may be an adaptive response to persisting in a biotically simplified, yet consistently physically harsh environment.

Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome.

Gut microbial diversity is thought to reflect the co-evolution of microbes and their hosts as well as current host-specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial-nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome-wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (ß-diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among-population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).

Accurate Inference of Tree Topologies from Multiple Sequence Alignments Using Deep Learning.

Reconstructing the phylogenetic relationships between species is one of the most formidable tasks in evolutionary biology. Multiple methods exist to reconstruct phylogenetic trees, each with their own strengths and weaknesses. Both simulation and empirical studies have identified several "zones" of parameter space where accuracy of some methods can plummet, even for four-taxon trees. Further, some methods can have undesirable statistical properties such as statistical inconsistency and/or the tendency to be positively misleading (i.e. assert strong support for the incorrect tree topology). Recently, deep learning techniques have made inroads on a number of both new and longstanding problems in biological research. In this study, we designed a deep convolutional neural network (CNN) to infer quartet topologies from multiple sequence alignments. This CNN can readily be trained to make inferences using both gapped and ungapped data. We show that our approach is highly accurate on simulated data, often outperforming traditional methods, and is remarkably robust to bias-inducing regions of parameter space such as the Felsenstein zone and the Farris zone. We also demonstrate that the confidence scores produced by our CNN can more accurately assess support for the chosen topology than bootstrap and posterior probability scores from traditional methods. Although numerous practical challenges remain, these findings suggest that the deep learning approaches such as ours have the potential to produce more accurate phylogenetic inferences.

Identification and evolutionary characterization of salt-responsive transcription factors in the succulent halophyte Suaeda fruticosa.

Transcription factors are key regulatory elements that affect gene expression in response to specific signals, including environmental stresses such as salinity. Halophytes are specialized plants that have the ability to complete their life cycle in saline environments. In this study we have identified and characterized the evolutionary relationships of putative transcription factors (TF) in an obligate succulent halophyte, Suaeda fruticosa, that are involved in conferring salt tolerance. Using RNA-seq data we have analyzed the expression patterns of certain TF families, predicted protein-protein interactions, and analyzed evolutionary trajectories to elucidate their possible roles in salt tolerance. We have detected the top differentially expressed (DE) transcription factor families (MYB, CAMTA, MADS-box and bZIP) that show the most pronounced response to salinity. The majority of DE genes in the four aforementioned TF families cluster together on TF phylogenetic trees, which suggests common evolutionary origins and trajectories. This research represents the first comprehensive TF study of a leaf succulent halophyte including their evolutionary relationships with TFs in other halophyte and salt-senstive plants. These findings provide a foundation for understanding the function of salt-responsive transcription factors in salt tolerance and associated gene regulation in plants.

Firefly genomes illuminate parallel origins of bioluminescence in beetles.

Fireflies and their luminous courtships have inspired centuries of scientific study. Today firefly luciferase is widely used in biotechnology, but the evolutionary origin of bioluminescence within beetles remains unclear. To shed light on this long-standing question, we sequenced the genomes of two firefly species that diverged over 100 million-years-ago: the North American Photinus pyralis and Japanese Aquatica lateralis. To compare bioluminescent origins, we also sequenced the genome of a related click beetle, the Caribbean Ignelater luminosus, with bioluminescent biochemistry near-identical to fireflies, but anatomically unique light organs, suggesting the intriguing hypothesis of parallel gains of bioluminescence. Our analyses support independent gains of bioluminescence in fireflies and click beetles, and provide new insights into the genes, chemical defenses, and symbionts that evolved alongside their luminous lifestyle.

Overcoming the loss of blue sensitivity through opsin duplication in the largest animal group, beetles.

Opsin proteins are fundamental components of animal vision whose structure largely determines the sensitivity of visual pigments to different wavelengths of light. Surprisingly little is known about opsin evolution in beetles, even though they are the most species rich animal group on Earth and exhibit considerable variation in visual system sensitivities. We reveal the patterns of opsin evolution across 62 beetle species and relatives. Our results show that the major insect opsin class (SW) that typically confers sensitivity to "blue" wavelengths was lost ~300 million years ago, before the origin of modern beetles. We propose that UV and LW opsin gene duplications have restored the potential for trichromacy (three separate channels for colour vision) in beetles up to 12 times and more specifically, duplications within the UV opsin class have likely led to the restoration of "blue" sensitivity up to 10 times. This finding reveals unexpected plasticity within the insect visual system and highlights its remarkable ability to evolve and adapt to the available light and visual cues present in the environment.

The OGCleaner: filtering false-positive homology clusters.

Detecting homologous sequences in organisms is an essential step in protein structure and function prediction, gene annotation and phylogenetic tree construction. Heuristic methods are often employed for quality control of putative homology clusters. These heuristics, however, usually only apply to pairwise sequence comparison and do not examine clusters as a whole. We present the Orthology Group Cleaner (the OGCleaner), a tool designed for filtering putative orthology groups as homology or non-homology clusters by considering all sequences in a cluster. The OGCleaner relies on high-quality orthologous groups identified in OrthoDB to train machine learning algorithms that are able to distinguish between true-positive and false-positive homology groups. This package aims to improve the quality of phylogenetic tree construction especially in instances of lower-quality transcriptome assemblies. AVAILABILITY AND IMPLEMENTATION: https://github.com/byucsl/ogcleaner CONTACT: sfujimoto@gmail.comSupplementary information: Supplementary data are available at Bioinformatics online.

Opsins have evolved under the permanent heterozygote model: insights from phylotranscriptomics of Odonata.

Gene duplication plays a central role in adaptation to novel environments by providing new genetic material for functional divergence and evolution of biological complexity. Several evolutionary models have been proposed for gene duplication to explain how new gene copies are preserved by natural selection, but these models have rarely been tested using empirical data. Opsin proteins, when combined with a chromophore, form a photopigment that is responsible for the absorption of light, the first step in the phototransduction cascade. Adaptive gene duplications have occurred many times within the animal opsins' gene family, leading to novel wavelength sensitivities. Consequently, opsins are an attractive choice for the study of gene duplication evolutionary models. Odonata (dragonflies and damselflies) have the largest opsin repertoire of any insect currently known. Additionally, there is tremendous variation in opsin copy number between species, particularly in the long-wavelength-sensitive (LWS) class. Using comprehensive phylotranscriptomic and statistical approaches, we tested various evolutionary models of gene duplication. Our results suggest that both the blue-sensitive (BS) and LWS opsin classes were subjected to strong positive selection that greatly weakens after multiple duplication events, a pattern that is consistent with the permanent heterozygote model. Due to the immense interspecific variation and duplicability potential of opsin genes among odonates, they represent a unique model system to test hypotheses regarding opsin gene duplication and diversification at the molecular level.

Odonata (dragonflies and damselflies) as a bridge between ecology and evolutionary genomics.

Odonata (dragonflies and damselflies) present an unparalleled insect model to integrate evolutionary genomics with ecology for the study of insect evolution. Key features of Odonata include their ancient phylogenetic position, extensive phenotypic and ecological diversity, several unique evolutionary innovations, ease of study in the wild and usefulness as bioindicators for freshwater ecosystems worldwide. In this review, we synthesize studies on the evolution, ecology and physiology of odonates, highlighting those areas where the integration of ecology with genomics would yield significant insights into the evolutionary processes that would not be gained easily by working on other animal groups. We argue that the unique features of this group combined with their complex life cycle, flight behaviour, diversity in ecological niches and their sensitivity to anthropogenic change make odonates a promising and fruitful taxon for genomics focused research. Future areas of research that deserve increased attention are also briefly outlined.

A cure for the blues: opsin duplication and subfunctionalization for short-wavelength sensitivity in jewel beetles (Coleoptera: Buprestidae).

BACKGROUND: Arthropods have received much attention as a model for studying opsin evolution in invertebrates. Yet, relatively few studies have investigated the diversity of opsin proteins that underlie spectral sensitivity of the visual pigments within the diverse beetles (Insecta: Coleoptera). Previous work has demonstrated that beetles appear to lack the short-wavelength-sensitive (SWS) opsin class that typically confers sensitivity to the "blue" region of the light spectrum. However, this is contrary to established physiological data in a number of Coleoptera. To explore potential adaptations at the molecular level that may compensate for the loss of the SWS opsin, we carried out an exploration of the opsin proteins within a group of beetles (Buprestidae) where short-wave sensitivity has been demonstrated. RNA-seq data were generated to identify opsin proteins from nine taxa comprising six buprestid species (including three male/female pairs) across four subfamilies. Structural analyses of recovered opsins were conducted and compared to opsin sequences in other insects across the main opsin classes-ultraviolet, short-wavelength, and long-wavelength. RESULTS: All nine buprestids were found to express two opsin copies in each of the ultraviolet and long-wavelength classes, contrary to the single copies recovered in all other molecular studies of adult beetle opsin expression. No SWS opsin class was recovered. Furthermore, the male Agrilus planipennis (emerald ash borer-EAB) expressed a third LWS opsin at low levels that is presumed to be a larval copy. Subsequent homology and structural analyses identified multiple amino acid substitutions in the UVS and LWS copies that could confer short-wavelength sensitivity. CONCLUSIONS: This work is the first to compare expressed opsin genes against known electrophysiological data that demonstrate multiple peak sensitivities in Coleoptera. We report the first instance of opsin duplication in adult beetles, which occurs in both the UVS and LWS opsin classes. Through structural comparisons of known insect opsins, we suggest that opsin duplication and amino acid variation within the chromophore binding pocket explains sensitivity in the short-wavelength portion of the visible light spectrum in these species. These findings are the first to reveal molecular complexity of the color vision system within beetles.

Detecting false positive sequence homology: a machine learning approach.

BACKGROUND: Accurate detection of homologous relationships of biological sequences (DNA or amino acid) amongst organisms is an important and often difficult task that is essential to various evolutionary studies, ranging from building phylogenies to predicting functional gene annotations. There are many existing heuristic tools, most commonly based on bidirectional BLAST searches that are used to identify homologous genes and combine them into two fundamentally distinct classes: orthologs and paralogs. Due to only using heuristic filtering based on significance score cutoffs and having no cluster post-processing tools available, these methods can often produce multiple clusters constituting unrelated (non-homologous) sequences. Therefore sequencing data extracted from incomplete genome/transcriptome assemblies originated from low coverage sequencing or produced by de novo processes without a reference genome are susceptible to high false positive rates of homology detection. RESULTS: In this paper we develop biologically informative features that can be extracted from multiple sequence alignments of putative homologous genes (orthologs and paralogs) and further utilized in context of guided experimentation to verify false positive outcomes. We demonstrate that our machine learning method trained on both known homology clusters obtained from OrthoDB and randomly generated sequence alignments (non-homologs), successfully determines apparent false positives inferred by heuristic algorithms especially among proteomes recovered from low-coverage RNA-seq data. Almost ~42 % and ~25 % of predicted putative homologies by InParanoid and HaMStR respectively were classified as false positives on experimental data set. CONCLUSIONS: Our process increases the quality of output from other clustering algorithms by providing a novel post-processing method that is both fast and efficient at removing low quality clusters of putative homologous genes recovered by heuristic-based approaches.

Intra-specific regulatory variation in Drosophila pseudoobscura.

It is generally accepted that gene regulation serves an important role in determining the phenotype. To shed light on the evolutionary forces operating on gene regulation, previous studies mainly focused on the expression differences between species and their inter-specific hybrids. Here, we use RNA-Seq to study the intra-specific distribution of cis- and trans-regulatory variation in Drosophila pseudoobscura. Consistent with previous results, we find almost twice as many genes (26%) with significant trans-effects than genes with significant cis-effects (18%). While this result supports the previous suggestion of a larger mutational target of trans-effects, we also show that trans-effects may be subjected to purifying selection. Our results underline the importance of intra-specific analyses for the understanding of the evolution of gene expression.