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.

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.

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.

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.

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.

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.

Phylogenomic Analyses Support Traditional Relationships within Cnidaria.

Cnidaria, the sister group to Bilateria, is a highly diverse group of animals in terms of morphology, lifecycles, ecology, and development. How this diversity originated and evolved is not well understood because phylogenetic relationships among major cnidarian lineages are unclear, and recent studies present contrasting phylogenetic hypotheses. Here, we use transcriptome data from 15 newly-sequenced species in combination with 26 publicly available genomes and transcriptomes to assess phylogenetic relationships among major cnidarian lineages. Phylogenetic analyses using different partition schemes and models of molecular evolution, as well as topology tests for alternative phylogenetic relationships, support the monophyly of Medusozoa, Anthozoa, Octocorallia, Hydrozoa, and a clade consisting of Staurozoa, Cubozoa, and Scyphozoa. Support for the monophyly of Hexacorallia is weak due to the equivocal position of Ceriantharia. Taken together, these results further resolve deep cnidarian relationships, largely support traditional phylogenetic views on relationships, and provide a historical framework for studying the evolutionary processes involved in one of the most ancient animal radiations.

Automation and Evaluation of the SOWH Test with SOWHAT.

The Swofford-Olsen-Waddell-Hillis (SOWH) test evaluates statistical support for incongruent phylogenetic topologies. It is commonly applied to determine if the maximum likelihood tree in a phylogenetic analysis is significantly different than an alternative hypothesis. The SOWH test compares the observed difference in log-likelihood between two topologies to a null distribution of differences in log-likelihood generated by parametric resampling. The test is a well-established phylogenetic method for topology testing, but it is sensitive to model misspecification, it is computationally burdensome to perform, and its implementation requires the investigator to make several decisions that each have the potential to affect the outcome of the test. We analyzed the effects of multiple factors using seven data sets to which the SOWH test was previously applied. These factors include a number of sample replicates, likelihood software, the introduction of gaps to simulated data, the use of distinct models of evolution for data simulation and likelihood inference, and a suggested test correction wherein an unresolved "zero-constrained" tree is used to simulate sequence data. To facilitate these analyses and future applications of the SOWH test, we wrote SOWHAT, a program that automates the SOWH test. We find that inadequate bootstrap sampling can change the outcome of the SOWH test. The results also show that using a zero-constrained tree for data simulation can result in a wider null distribution and higher p-values, but does not change the outcome of the SOWH test for most of the data sets tested here. These results will help others implement and evaluate the SOWH test and allow us to provide recommendations for future applications of the SOWH test. SOWHAT is available for download from https://github.com/josephryan/SOWHAT.

Stem cells in Nanomia bijuga (Siphonophora), a colonial animal with localized growth zones.

BACKGROUND: Siphonophores (Hydrozoa) have unparalleled colony-level complexity, precision of colony organization, and functional specialization between zooids (i.e., the units that make up colonies). Previous work has shown that, unlike other colonial animals, most growth in siphonophores is restricted to one or two well-defined growth zones that are the sites of both elongation and zooid budding. It remained unknown, however, how this unique colony growth and development is realized at the cellular level. RESULTS: To understand the colony-level growth and development of siphonophores at the cellular level, we characterize the distribution of proliferating cells and interstitial stem cells (i-cells) in the siphonophore Nanomia bijuga. Within the colony, we find evidence that i-cells are present at the tip of the horn, the structure within the growth zone that gives rise to new zooids. Co-localized gene expression of vasa-1, pl10, piwi, nanos-1, and nanos-2 suggests that i-cells persist in the youngest zooid buds and that i-cells become progressively restricted to specific regions within the zooids until they are mostly absent from the oldest zooids. The examined genes remain expressed in gametogenic regions. No evidence for i-cells is found in the stem between maturing zooids. Domains of high cell proliferation include regions where the examined genes are expressed, but also include some areas in which the examined genes were not expressed such as the stem within the growth zones. Cell proliferation in regions devoid of vasa-1, pl10, piwi, nanos-1, and nanos-2 expression indicates the presence of mitotically active epithelial cell lineages and, potentially, progenitor cell populations. CONCLUSIONS: We provide the first evidence for i-cells in a siphonophore. Our findings suggest maintenance of i-cell populations at the sites of growth zones and that these sites are the main source of i-cells. This restriction of stem cells to particular regions in the colony, in combination with localized budding and spatial patterning during pro-bud subdivision, may play a major role in facilitating the precision of siphonophore growth. Spatially restricted maintenance of i-cells in mature zooids and absence of i-cells along the stem may explain the reduced developmental plasticity in older parts of the colony.

The histology of Nanomia bijuga (Hydrozoa: Siphonophora).

The siphonophore Nanomia bijuga is a pelagic hydrozoan (Cnidaria) with complex morphological organization. Each siphonophore is made up of many asexually produced, genetically identical zooids that are functionally specialized and morphologically distinct. These zooids predominantly arise by budding in two growth zones, and are arranged in precise patterns. This study describes the cellular anatomy of several zooid types, the stem, and the gas-filled float, called the pneumatophore. The distribution of cellular morphologies across zooid types enhances our understanding of zooid function. The unique absorptive cells in the palpon, for example, indicate specialized intracellular digestive processing in this zooid type. Though cnidarians are usually thought of as mono-epithelial, we characterize at least two cellular populations in this species which are not connected to a basement membrane. This work provides a greater understanding of epithelial diversity within the cnidarians, and will be a foundation for future studies on N. bijuga, including functional assays and gene expression analyses.