Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony.

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

Stepwise emergence of the neuronal gene expression program in early animal evolution.

The assembly of the neuronal and other major cell type programs occurred early in animal evolution. We can reconstruct this process by studying non-bilaterians like placozoans. These small disc-shaped animals not only have nine morphologically described cell types and no neurons but also show coordinated behaviors triggered by peptide-secreting cells. We investigated possible neuronal affinities of these peptidergic cells using phylogenetics, chromatin profiling, and comparative single-cell genomics in four placozoans. We found conserved cell type expression programs across placozoans, including populations of transdifferentiating and cycling cells, suggestive of active cell type homeostasis. We also uncovered fourteen peptidergic cell types expressing neuronal-associated components like the pre-synaptic scaffold that derive from progenitor cells with neurogenesis signatures. In contrast, earlier-branching animals like sponges and ctenophores lacked this conserved expression. Our findings indicate that key neuronal developmental and effector gene modules evolved before the advent of cnidarian/bilaterian neurons in the context of paracrine cell signaling.

cGLRs are a diverse family of pattern recognition receptors in innate immunity.

Cyclic GMP-AMP synthase (cGAS) is an enzyme in human cells that controls an immune response to cytosolic DNA. Upon binding DNA, cGAS synthesizes a nucleotide signal 2'3'-cGAMP that activates STING-dependent downstream immunity. Here, we discover that cGAS-like receptors (cGLRs) constitute a major family of pattern recognition receptors in innate immunity. Building on recent analysis in Drosophila, we identify >3,000 cGLRs present in nearly all metazoan phyla. A forward biochemical screening of 150 animal cGLRs reveals a conserved mechanism of signaling including response to dsDNA and dsRNA ligands and synthesis of isomers of the nucleotide signals cGAMP, c-UMP-AMP, and c-di-AMP. Combining structural biology and in vivo analysis in coral and oyster animals, we explain how synthesis of distinct nucleotide signals enables cells to control discrete cGLR-STING signaling pathways. Our results reveal cGLRs as a widespread family of pattern recognition receptors and establish molecular rules that govern nucleotide signaling in animal immunity.

A stony coral cell atlas illuminates the molecular and cellular basis of coral symbiosis, calcification, and immunity.

Stony corals are colonial cnidarians that sustain the most biodiverse marine ecosystems on Earth: coral reefs. Despite their ecological importance, little is known about the cell types and molecular pathways that underpin the biology of reef-building corals. Using single-cell RNA sequencing, we define over 40 cell types across the life cycle of Stylophora pistillata. We discover specialized immune cells, and we uncover the developmental gene expression dynamics of calcium-carbonate skeleton formation. By simultaneously measuring the transcriptomes of coral cells and the algae within them, we characterize the metabolic programs involved in symbiosis in both partners. We also trace the evolution of these coral cell specializations by phylogenetic integration of multiple cnidarian cell type atlases. Overall, this study reveals the molecular and cellular basis of stony coral biology.

Cnidarian Cell Type Diversity and Regulation Revealed by Whole-Organism Single-Cell RNA-Seq.

The emergence and diversification of cell types is a leading factor in animal evolution. So far, systematic characterization of the gene regulatory programs associated with cell type specificity was limited to few cell types and few species. Here, we perform whole-organism single-cell transcriptomics to map adult and larval cell types in the cnidarian Nematostella vectensis, a non-bilaterian animal with complex tissue-level body-plan organization. We uncover eight broad cell classes in Nematostella, including neurons, cnidocytes, and digestive cells. Each class comprises different subtypes defined by the expression of multiple specific markers. In particular, we characterize a surprisingly diverse repertoire of neurons, which comparative analysis suggests are the result of lineage-specific diversification. By integrating transcription factor expression, chromatin profiling, and sequence motif analysis, we identify the regulatory codes that underlie Nematostella cell-specific expression. Our study reveals cnidarian cell type complexity and provides insights into the evolution of animal cell-specific genomic regulation.

The Dynamic Regulatory Genome of Capsaspora and the Origin of Animal Multicellularity.

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

Evolutionary cell type mapping with single-cell genomics.

A fundamental characteristic of animal multicellularity is the spatial coexistence of functionally specialized cell types that are all encoded by a single genome sequence. Cell type transcriptional programs are deployed and maintained by regulatory mechanisms that control the asymmetric, differential access to genomic information in each cell. This genome regulation ultimately results in specific cellular phenotypes. However, the emergence, diversity, and evolutionary dynamics of animal cell types remain almost completely unexplored beyond a few species. Single-cell genomics is emerging as a powerful tool to build comprehensive catalogs of cell types and their associated gene regulatory programs in non-traditional model species. We review the current state of sampling efforts across the animal tree of life and challenges ahead for the comparative study of cell type programs. We also discuss how the phylogenetic integration of cell atlases can lead to the development of models of cell type evolution and a phylogenetic taxonomy of cells.

The origin of Metazoa: a unicellular perspective.

The first animals evolved from an unknown single-celled ancestor in the Precambrian period. Recently, the identification and characterization of the genomic and cellular traits of the protists most closely related to animals have shed light on the origin of animals. Comparisons of animals with these unicellular relatives allow us to reconstruct the first evolutionary steps towards animal multicellularity. Here, we review the results of these investigations and discuss their implications for understanding the earliest stages of animal evolution, including the origin of metazoan genes and genome function.

Orthology Clusters from Gene Trees with Possvm.

Possvm (Phylogenetic Ortholog Sorting with Species oVerlap and MCL [Markov clustering algorithm]) is a tool that automates the process of identifying clusters of orthologous genes from precomputed phylogenetic trees and classifying gene families. It identifies orthology relationships between genes using the species overlap algorithm to infer taxonomic information from the gene tree topology, and then uses the MCL to identify orthology clusters and provide annotated gene families. Our benchmarking shows that this approach, when provided with accurate phylogenies, is able to identify manually curated orthogroups with very high precision and recall. Overall, Possvm automates the routine process of gene tree inspection and annotation in a highly interpretable manner, and provides reusable outputs and phylogeny-aware gene annotations that can be used to inform comparative genomics and gene family evolution analyses.

The eukaryotic ancestor had a complex ubiquitin signaling system of archaeal origin.

The origin of the eukaryotic cell is one of the most important transitions in the history of life. However, the emergence and early evolution of eukaryotes remains poorly understood. Recent data have shown that the last eukaryotic common ancestor (LECA) was much more complex than previously thought. The LECA already had the genetic machinery encoding the endomembrane apparatus, spliceosome, nuclear pore, and myosin and kinesin cytoskeletal motors. It is unclear, however, when the functional regulation of these cellular components evolved. Here, we address this question by analyzing the origin and evolution of the ubiquitin (Ub) signaling system, one of the most important regulatory layers in eukaryotes. We delineated the evolution of the whole Ub, Small-Ub-related MOdifier (SUMO), and Ub-fold modifier 1 (Ufm1) signaling networks by analyzing representatives from all major eukaryotic, bacterial, and archaeal lineages. We found that the Ub toolkit had a pre-eukaryotic origin and is present in three extant archaeal groups. The pre-eukaryotic Ub toolkit greatly expanded during eukaryogenesis, through massive gene innovation and diversification of protein domain architectures. This resulted in a LECA with essentially all of the Ub-related genes, including the SUMO and Ufm1 Ub-like systems. Ub and SUMO signaling further expanded during eukaryotic evolution, especially labeling and delabeling enzymes responsible for substrate selection. Additionally, we analyzed protein domain architecture evolution and found that multicellular lineages have the most complex Ub systems in terms of domain architectures. Together, we demonstrate that the Ub system predates the origin of eukaryotes and that a burst of innovation during eukaryogenesis led to a LECA with complex posttranslational regulation.

Tracing the Evolutionary History of Inositol, 1, 4, 5-Trisphosphate Receptor: Insights from Analyses of Capsaspora owczarzaki Ca2+ Release Channel Orthologs.

Cellular Ca(2+) homeostasis is tightly regulated and is pivotal to life. Inositol 1,4,5-trisphosphate receptors (IP3Rs) and ryanodine receptors (RyRs) are the major ion channels that regulate Ca(2+) release from intracellular stores. Although these channels have been extensively investigated in multicellular organisms, an appreciation of their evolution and the biology of orthologs in unicellular organisms is largely lacking. Extensive phylogenetic analyses reveal that the IP3R gene superfamily is ancient and diverged into two subfamilies, IP3R-A and IP3R-B/RyR, at the dawn of Opisthokonta. IP3R-B/RyR further diversified into IP3R-B and RyR at the stem of Filozoa. Subsequent evolution and speciation of Holozoa is associated with duplication of IP3R-A and RyR genes, and loss of IP3R-B in the vertebrate lineages. To gain insight into the properties of IP3R important for the challenges of multicellularity, the IP3R-A and IP3R-B family orthologs were cloned from Capsaspora owczarzaki, a close unicellular relative to Metazoa (designated as CO.IP3R-A and CO.IP3R-B). Both proteins were targeted to the endoplasmic reticulum. However, CO.IP3R-A, but strikingly not CO.IP3R-B, bound IP3, exhibited robust Ca(2+) release activity and associated with mammalian IP3Rs. These data indicate strongly that CO.IP3R-A as an exemplar of ancestral IP3R-A orthologs forms bona fide IP3-gated channels. Notably, however, CO.IP3R-A appears not to be regulated by Ca(2+), ATP or Protein kinase A-phosphorylation. Collectively, our findings explore the origin, conservation, and diversification of IP3R gene families and provide insight into the functionality of ancestral IP3Rs and the added specialization of these proteins in Metazoa.

Transcription factor evolution in eukaryotes and the assembly of the regulatory toolkit in multicellular lineages.

Transcription factors (TFs) are the main players in transcriptional regulation in eukaryotes. However, it remains unclear what role TFs played in the origin of all of the different eukaryotic multicellular lineages. In this paper, we explore how the origin of TF repertoires shaped eukaryotic evolution and, in particular, their role into the emergence of multicellular lineages. We traced the origin and expansion of all known TFs through the eukaryotic tree of life, using the broadest possible taxon sampling and an updated phylogenetic background. Our results show that the most complex multicellular lineages (i.e., those with embryonic development, Metazoa and Embryophyta) have the most complex TF repertoires, and that these repertoires were assembled in a stepwise manner. We also show that a significant part of the metazoan and embryophyte TF toolkits evolved earlier, in their respective unicellular ancestors. To gain insights into the role of TFs in the development of both embryophytes and metazoans, we analyzed TF expression patterns throughout their ontogeny. The expression patterns observed in both groups recapitulate those of the whole transcriptome, but reveal some important differences. Our comparative genomics and expression data reshape our view on how TFs contributed to eukaryotic evolution and reveal the importance of TFs to the origins of multicellularity and embryonic development.

Early evolution of the T-box transcription factor family.

Developmental transcription factors are key players in animal multicellularity, being members of the T-box family that are among the most important. Until recently, T-box transcription factors were thought to be exclusively present in metazoans. Here, we report the presence of T-box genes in several nonmetazoan lineages, including ichthyosporeans, filastereans, and fungi. Our data confirm that Brachyury is the most ancient member of the T-box family and establish that the T-box family diversified at the onset of Metazoa. Moreover, we demonstrate functional conservation of a homolog of Brachyury of the protist Capsaspora owczarzaki in Xenopus laevis. By comparing the molecular phenotype of C. owczarzaki Brachyury with that of homologs of early branching metazoans, we define a clear difference between unicellular holozoan and metazoan Brachyury homologs, suggesting that the specificity of Brachyury emerged at the origin of Metazoa. Experimental determination of the binding preferences of the C. owczarzaki Brachyury results in a similar motif to that of metazoan Brachyury and other T-box classes. This finding suggests that functional specificity between different T-box classes is likely achieved by interaction with alternative cofactors, as opposed to differences in binding specificity.

Insights into the origin of metazoan filopodia and microvilli.

Filopodia are fine actin-based cellular projections used for both environmental sensing and cell motility, and they are essential organelles for metazoan cells. In this study, we reconstruct the origin of metazoan filopodia and microvilli. We first report on the evolutionary assembly of the filopodial molecular toolkit and show that homologs of many metazoan filopodial components, including fascin and myosin X, were already present in the unicellular or colonial progenitors of metazoans. Furthermore, we find that the actin crosslinking protein fascin localizes to filopodia-like structures and microvilli in the choanoflagellate Salpingoeca rosetta. In addition, homologs of filopodial genes in the holozoan Capsaspora owczarzaki are upregulated in filopodia-bearing cells relative to those that lack them. Therefore, our findings suggest that proteins essential for metazoan filopodia and microvilli are functionally conserved in unicellular and colonial holozoans and that the last common ancestor of metazoans bore a complex and specific filopodial machinery.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative.

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

An amphioxus neurula stage cell atlas supports a complex scenario for the emergence of vertebrate head mesoderm.

The emergence of new structures can often be linked to the evolution of novel cell types that follows the rewiring of developmental gene regulatory subnetworks. Vertebrates are characterized by a complex body plan compared to the other chordate clades and the question remains of whether and how the emergence of vertebrate morphological innovations can be related to the appearance of new embryonic cell populations. We previously proposed, by studying mesoderm development in the cephalochordate amphioxus, a scenario for the evolution of the vertebrate head mesoderm. To further test this scenario at the cell population level, we used scRNA-seq to construct a cell atlas of the amphioxus neurula, stage at which the main mesodermal compartments are specified. Our data allowed us to validate the presence of a prechordal-plate like territory in amphioxus. Additionally, the transcriptomic profile of somite cell populations supports the homology between specific territories of amphioxus somites and vertebrate cranial/pharyngeal and lateral plate mesoderm. Finally, our work provides evidence that the appearance of the specific mesodermal structures of the vertebrate head was associated to both segregation of pre-existing cell populations, and co-option of new genes for the control of myogenesis.

Single cell atlas of Xenoturbella bocki highlights limited cell-type complexity.

Phylogenetic analyses over the last two decades have united a few small, and previously orphan clades, the nematodermatids, acoels and xenoturbelids, into the phylum Xenacoelomorpha. Some phylogenetic analyses support a sister relationship between Xenacoelomorpha and Ambulacraria (Xenambulacraria), while others suggest that Xenacoelomorpha may be sister to the rest of the Bilateria (Nephrozoa). An understanding of the cell type complements of Xenacoelomorphs is essential to assessing these alternatives as well as to our broader understanding of bilaterian cell type evolution. Employing whole organism single-cell RNA-seq in the marine xenacoelomorph worm Xenoturbella bocki, we show that Xenambulacrarian nerve nets share regulatory features and a peptidergic identity with those found in cnidarians and protostomes and more broadly share muscle and gland cell similarities with other metazoans. Taken together, these data are consistent with broad homologies of animal gland, muscle, and neurons as well as more specific affinities between Xenoturbella and acoel gut and epidermal tissues, consistent with the monophyly of Xenacoelomorpha.

Ancient origin of the integrin-mediated adhesion and signaling machinery.

The evolution of animals (metazoans) from their unicellular ancestors required the emergence of novel mechanisms for cell adhesion and cell-cell communication. One of the most important cell adhesion mechanisms for metazoan development is integrin-mediated adhesion and signaling. The integrin adhesion complex mediates critical interactions between cells and the extracellular matrix, modulating several aspects of cell physiology. To date this machinery has been considered strictly metazoan specific. Here we report the results of a comparative genomic analysis of the integrin adhesion machinery, using genomic data from several unicellular relatives of Metazoa and Fungi. Unexpectedly, we found that core components of the integrin adhesion complex are encoded in the genome of the apusozoan protist Amastigomonas sp., and therefore their origins predate the divergence of Opisthokonta, the clade that includes metazoans and fungi. Furthermore, our analyses suggest that key components of this apparatus have been lost independently in fungi and choanoflagellates. Our data highlight the fact that many of the key genes that had formerly been cited as crucial for metazoan origins have a much earlier origin. This underscores the importance of gene cooption in the unicellular-to-multicellular transition that led to the emergence of the Metazoa.

Systematic benchmarking of single-cell ATAC-sequencing protocols.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Unexpected repertoire of metazoan transcription factors in the unicellular holozoan Capsaspora owczarzaki.

How animals (metazoans) originated from their single-celled ancestors remains a major question in biology. As transcriptional regulation is crucial to animal development, deciphering the early evolution of associated transcription factors (TFs) is critical to understanding metazoan origins. In this study, we uncovered the repertoire of 17 metazoan TFs in the amoeboid holozoan Capsaspora owczarzaki, a representative of a unicellular lineage that is closely related to choanoflagellates and metazoans. Phylogenetic and comparative genomic analyses with the broadest possible taxonomic sampling allowed us to formulate new hypotheses regarding the origin and evolution of developmental metazoan TFs. We show that the complexity of the TF repertoire in C. owczarzaki is strikingly high, pushing back further the origin of some TFs formerly thought to be metazoan specific, such as T-box or Runx. Nonetheless, TF families whose beginnings antedate the origin of the animal kingdom, such as homeodomain or basic helix-loop-helix, underwent significant expansion and diversification along metazoan and eumetazoan stems.