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.

Integrating phylogenies into single-cell RNA sequencing analysis allows comparisons across species, genes, and cells.

Comparisons of single-cell RNA sequencing (scRNA-seq) data across species can reveal links between cellular gene expression and the evolution of cell functions, features, and phenotypes. These comparisons evoke evolutionary histories, as depicted by phylogenetic trees, that define relationships between species, genes, and cells. This Essay considers each of these in turn, laying out challenges and solutions derived from a phylogenetic comparative approach and relating these solutions to previously proposed methods for the pairwise alignment of cellular dimensional maps. This Essay contends that species trees, gene trees, cell phylogenies, and cell lineages can all be reconciled as descriptions of the same concept-the tree of cellular life. By integrating phylogenetic approaches into scRNA-seq analyses, challenges for building informed comparisons across species can be overcome, and hypotheses about gene and cell evolution can be robustly tested.

Cell cycle variants during Drosophila male accessory gland development.

The Drosophila melanogaster male accessory gland (AG) is a functional analog of the mammalian prostate and seminal vesicles containing two secretory epithelial cell types, termed main and secondary cells. This tissue is responsible for making and secreting seminal fluid proteins and other molecules that contribute to successful reproduction. The cells of this tissue are binucleate and polyploid, due to variant cell cycles that include endomitosis and endocycling during metamorphosis. Here, we provide evidence of additional cell cycle variants in this tissue. We show that main cells of the gland are connected by ring canals that form after the penultimate mitosis, and we describe an additional post-eclosion endocycle required for gland maturation that is dependent on juvenile hormone signaling. We present evidence that the main cells of the D. melanogaster AG undergo a unique cell cycle reprogramming throughout organ development that results in step-wise cell cycle truncations culminating in cells containing two octoploid nuclei with under-replicated heterochromatin in the mature gland. We propose this tissue as a model to study developmental and hormonal temporal control of cell cycle variants in terminally differentiating tissues.

Normalizing need not be the norm: count-based math for analyzing single-cell data.

Counting transcripts of mRNA are a key method of observation in modern biology. With advances in counting transcripts in single cells (single-cell RNA sequencing or scRNA-seq), these data are routinely used to identify cells by their transcriptional profile, and to identify genes with differential cellular expression. Because the total number of transcripts counted per cell can vary for technical reasons, the first step of many commonly used scRNA-seq workflows is to normalize by sequencing depth, transforming counts into proportional abundances. The primary objective of this step is to reshape the data such that cells with similar biological proportions of transcripts end up with similar transformed measurements. But there is growing concern that normalization and other transformations result in unintended distortions that hinder both analyses and the interpretation of results. This has led to an intense focus on optimizing methods for normalization and transformation of scRNA-seq data. Here, we take an alternative approach, by avoiding normalization and transformation altogether. We abandon the use of distances to compare cells, and instead use a restricted algebra, motivated by measurement theory and abstract algebra, that preserves the count nature of the data. We demonstrate that this restricted algebra is sufficient to draw meaningful and practical comparisons of gene expression through the use of the dot product and other elementary operations. This approach sidesteps many of the problems with common transformations, and has the added benefit of being simpler and more intuitive. We implement our approach in the package countland, available in python and R.

Single-fly assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life.

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

The evolution of ovary-biased gene expression in Hawaiian Drosophila.

With detailed data on gene expression accessible from an increasingly broad array of species, we can test the extent to which our developmental genetic knowledge from model organisms predicts expression patterns and variation across species. But to know when differences in gene expression across species are significant, we first need to know how much evolutionary variation in gene expression we expect to observe. Here we provide an answer by analyzing RNAseq data across twelve species of Hawaiian Drosophilidae flies, focusing on gene expression differences between the ovary and other tissues. We show that over evolutionary time, there exists a cohort of ovary specific genes that is stable and that largely corresponds to described expression patterns from laboratory model Drosophila species. Our results also provide a demonstration of the prediction that, as phylogenetic distance increases, variation between species overwhelms variation between tissue types. Using ancestral state reconstruction of expression, we describe the distribution of evolutionary changes in tissue-biased expression, and use this to identify gains and losses of ovary-biased expression across these twelve species. We then use this distribution to calculate the evolutionary correlation in expression changes between genes, and demonstrate that genes with known interactions in D. melanogaster are significantly more correlated in their evolution than genes with no or unknown interactions. Finally, we use this correlation matrix to infer new networks of genes that share evolutionary trajectories, and we present these results as a dataset of new testable hypotheses about genetic roles and interactions in the function and evolution of the Drosophila ovary.

Phylotranscriptomics Reveals Discordance in the Phylogeny of Hawaiian Drosophila and Scaptomyza (Diptera: Drosophilidae).

Island radiations present natural laboratories for studying the evolutionary process. The Hawaiian Drosophilidae are one such radiation, with nearly 600 described species and substantial morphological and ecological diversification. These species are largely divided into a few major clades, but the relationship between clades remains uncertain. Here, we present new assembled transcriptomes from 12 species across these clades, and use these transcriptomes to resolve the base of the evolutionary radiation. We recover a new hypothesis for the relationship between clades, and demonstrate its support over previously published hypotheses. We then use the evolutionary radiation to explore dynamics of concordance in phylogenetic support, by analyzing the gene and site concordance factors for every possible topological combination of major groups. We show that high bootstrap values mask low evolutionary concordance, and we demonstrate that the most likely topology is distinct from the topology with the highest support across gene trees and from the topology with highest support across sites. We then combine all previously published genetic data for the group to estimate a time-calibrated tree for over 300 species of drosophilids. Finally, we digitize dozens of published Hawaiian Drosophilidae descriptions, and use this to pinpoint probable evolutionary shifts in reproductive ecology as well as body, wing, and egg size. We show that by examining the entire landscape of tree and trait space, we can gain a more complete understanding of how evolutionary dynamics play out across an island radiation.

Repeated loss of variation in insect ovary morphology highlights the role of development in life-history evolution.

The number of offspring an organism can produce is a key component of its evolutionary fitness and life history. Here we perform a test of the hypothesized trade-off between the number and size of offspring using thousands of descriptions of the number of egg-producing compartments in the insect ovary (ovarioles), a common proxy for potential offspring number in insects. We find evidence of a negative relationship between egg size and ovariole number when accounting for adult body size. However, in contrast to prior claims, we note that this relationship is not generalizable across all insect clades, and we highlight several factors that may have contributed to this size-number trade-off being stated as a general rule in previous studies. We reconstruct the evolution of the arrangement of cells that contribute nutrients and patterning information during oogenesis (nurse cells), and show that the diversification of ovariole number and egg size have both been largely independent of their presence or position within the ovariole. Instead, we show that ovariole number evolution has been shaped by a series of transitions between variable and invariant states, with multiple independent lineages evolving to have almost no variation in ovariole number. We highlight the implications of these invariant lineages on our understanding of the specification of ovariole number during development, as well as the importance of considering developmental processes in theories of life-history evolution.

Antisense oligonucleotides and nucleic acids generate hypersensitive platelets.

INTRODUCTION: Despite the great promise for therapies using antisense oligonucleotides (ASOs), their adverse effects, which include pro-inflammatory effects and thrombocytopenia, have limited their use. Previously, these effects have been linked to the phosphorothioate (PS) backbone necessary to prevent rapid ASO degradation in plasma. The main aim of this study was to assess the impact of the nucleic acid portion of an ASO-type drug on platelets and determine if it may contribute to thrombosis or thrombocytopenia. METHODS: Platelets were isolated from healthy donors and men with advanced prostate cancer. Effects of antisense oligonucleotides (ASO), oligonucleotides, gDNA, and microRNA on platelet activation and aggregation were evaluated. A mouse model of lung thrombosis was used to confirm the effects of PS-modified oligonucleotides in vivo. RESULTS: Platelet exposure to gDNA, miRNA, and oligonucleotides longer than 16-mer at a concentration above 8 mM resulted in the formation of hypersensitive platelets, characterized by an increased sensitivity to low-dose thrombin (0.1 nM) and increase in p-Selectin expression (6-8 fold greater than control; p < 0.001). The observed nucleic acid (NA) effects on platelets were toll-like receptor (TLR) -7 subfamily dependent. Injection of a p-Selectin inhibitor significantly (p = 0.02) reduced the formation of oligonucleotide-associated pulmonary microthrombosis in vivo. CONCLUSION: Our results suggest that platelet exposure to nucleic acids independent of the presence of a PS modification leads to a generation of hypersensitive platelets and requires TLR-7 subfamily receptors. ASO studies conducted in cancer patients may benefit from testing the ASO effects on platelets ex vivo before initiation of patient treatment.

Null hypotheses for developmental evolution.

How much evolutionary change in development do we expect? In this Spotlight, we argue that, as developmental biologists, we are in a prime position to contribute to the definition of a null hypothesis for developmental evolution: in other words, a hypothesis for how much developmental evolution we expect to observe over time. Today, we have access to an unprecedented array of developmental data from across the tree of life. Using these data, we can now consider development in the light of evolution, and vice versa, more deeply than ever before. As we do this, we may need to re-examine previous assumptions that appeared to serve us well when data points were fewer. Specifically, we think it is important to challenge assumptions that change is very rare for all developmental traits, especially if this assumption is used to sustain an erroneous view that evolution always optimizes adaptive traits toward increasing complexity.

A dataset of egg size and shape from more than 6,700 insect species.

Offspring size is a fundamental trait in disparate biological fields of study. This trait can be measured as the size of plant seeds, animal eggs, or live young, and it influences ecological interactions, organism fitness, maternal investment, and embryonic development. Although multiple evolutionary processes have been predicted to drive the evolution of offspring size, the phylogenetic distribution of this trait remains poorly understood, due to the difficulty of reliably collecting and comparing offspring size data from many species. Here we present a dataset of 10,449 morphological descriptions of insect eggs, with records for 6,706 unique insect species and representatives from every extant hexapod order. The dataset includes eggs whose volumes span more than eight orders of magnitude. We created this dataset by partially automating the extraction of egg traits from the primary literature. In the process, we overcame challenges associated with large-scale phenotyping by designing and employing custom bioinformatic solutions to common problems. We matched the taxa in this dataset to the currently accepted scientific names in taxonomic and genetic databases, which will facilitate the use of these data for testing pressing evolutionary hypotheses in offspring size evolution.

Insect egg size and shape evolve with ecology but not developmental rate.

Over the course of evolution, organism size has diversified markedly. Changes in size are thought to have occurred because of developmental, morphological and/or ecological pressures. To perform phylogenetic tests of the potential effects of these pressures, here we generated a dataset of more than ten thousand descriptions of insect eggs, and combined these with genetic and life-history datasets. We show that, across eight orders of magnitude of variation in egg volume, the relationship between size and shape itself evolves, such that previously predicted global patterns of scaling do not adequately explain the diversity in egg shapes. We show that egg size is not correlated with developmental rate and that, for many insects, egg size is not correlated with adult body size. Instead, we find that the evolution of parasitoidism and aquatic oviposition help to explain the diversification in the size and shape of insect eggs. Our study suggests that where eggs are laid, rather than universal allometric constants, underlies the evolution of insect egg size and shape.

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

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

Improved phylogenetic resolution within Siphonophora (Cnidaria) with implications for trait evolution.

Siphonophores are a diverse group of hydrozoans (Cnidaria) that are found at most depths of the ocean - from the surface, like the familiar Portuguese man of war, to the deep sea. They play important roles in ocean ecosystems, and are among the most abundant gelatinous predators. A previous phylogenetic study based on two ribosomal RNA genes provided insight into the internal relationships between major siphonophore groups. There was, however, little support for many deep relationships within the clade Codonophora. Here, we present a new siphonophore phylogeny based on new transcriptome data from 29 siphonophore species analyzed in combination with 14 publicly available genomic and transcriptomic datasets. We use this new phylogeny to reconstruct several traits that are central to siphonophore biology, including sexual system (monoecy vs. dioecy), gain and loss of zooid types, life history traits, and habitat. The phylogenetic relationships in this study are largely consistent with the previous phylogeny, but we find strong support for new clades within Codonophora that were previously unresolved. These results have important implications for trait evolution within Siphonophora, including favoring the hypothesis that monoecy arose at least twice.