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.

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.

A phylogenetic and proteomic reconstruction of eukaryotic chromatin evolution.

Histones and associated chromatin proteins have essential functions in eukaryotic genome organization and regulation. Despite this fundamental role in eukaryotic cell biology, we lack a phylogenetically comprehensive understanding of chromatin evolution. Here, we combine comparative proteomics and genomics analysis of chromatin in eukaryotes and archaea. Proteomics uncovers the existence of histone post-translational modifications in archaea. However, archaeal histone modifications are scarce, in contrast with the highly conserved and abundant marks we identify across eukaryotes. Phylogenetic analysis reveals that chromatin-associated catalytic functions (for example, methyltransferases) have pre-eukaryotic origins, whereas histone mark readers and chaperones are eukaryotic innovations. We show that further chromatin evolution is characterized by expansion of readers, including capture by transposable elements and viruses. Overall, our study infers detailed evolutionary history of eukaryotic chromatin: from its archaeal roots, through the emergence of nucleosome-based regulation in the eukaryotic ancestor, to the diversification of chromatin regulators and their hijacking by genomic parasites.

Functional and genomic comparative study of the bitter taste receptor family TAS2R: Insight into the role of human TAS2R5.

Bitterness is perceived in humans by 25 subtypes of bitter taste receptors (hTAS2R) that range from broadly tuned to more narrowly tuned receptors. hTAS2R5 is one of the most narrowly tuned bitter taste receptors in humans. In this study, we review the literature on this receptor and show there is no consensus about its role. We then compare the possible role of hTAS2R5 with that of the proteins of the TAS2R family in rat, mouse, and pig. A phylogenetic tree of all mammalian TAS2R domain-containing proteins showed that human hTAS2R5 has no ortholog in pig, mouse, or rat genomes. By comparing the agonists that are common to hTAS2R5 and other members of the family, we observed that hTAS2R39 is the receptor that shares most agonists with hTAS2R5. In mouse, some of these agonists activate mTas2r105 and mTas2r144, which are distant paralogs of hTAS2R5. mTas2r144 seems to be the receptor that is most similar to hTAS2R5 because they are both activated by the same agonists and have affinities in the same range of values. Then, we can conclude that hTAS2R5 has a unique functional specificity in humans as it is activated by selective agonists and that its closest functional homolog in mouse is the phylogenetically distant mTas2r144.

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.

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.

Resistance to pirimiphos-methyl in West African Anopheles is spreading via duplication and introgression of the Ace1 locus.

Vector population control using insecticides is a key element of current strategies to prevent malaria transmission in Africa. The introduction of effective insecticides, such as the organophosphate pirimiphos-methyl, is essential to overcome the recurrent emergence of resistance driven by the highly diverse Anopheles genomes. Here, we use a population genomic approach to investigate the basis of pirimiphos-methyl resistance in the major malaria vectors Anopheles gambiae and A. coluzzii. A combination of copy number variation and a single non-synonymous substitution in the acetylcholinesterase gene, Ace1, provides the key resistance diagnostic in an A. coluzzii population from Côte d'Ivoire that we used for sequence-based association mapping, with replication in other West African populations. The Ace1 substitution and duplications occur on a unique resistance haplotype that evolved in A. gambiae and introgressed into A. coluzzii, and is now common in West Africa primarily due to selection imposed by other organophosphate or carbamate insecticides. Our findings highlight the predictive value of this complex resistance haplotype for phenotypic resistance and clarify its evolutionary history, providing tools to for molecular surveillance of the current and future effectiveness of pirimiphos-methyl based interventions.

Conservative route to genome compaction in a miniature annelid.

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

RNA editing: an overlooked source of fine-scale adaptation in insect vectors?

RNA editing is a source of molecular diversity that regulates the functional repertoire of animal transcriptomes. Multiple studies in Drosophila have revealed that conserved editing events can be a source of evolutionary adaptations, and there is a solid body of evidence linking editing and the fine-tuning of neural genes, which are often targeted by insecticides used in vector control. Yet, despite these suggestive connections, genome-wide analyses of editing in insect vectors are conspicuously lacking. Future advances will require complementing the growing wealth of vector genomes with targeted transcriptome analyses. Here, we review recent investigations of the genetic footprints of adaptive RNA editing in insects and provide an overview of new methodologies applicable to studies of RNA editing in insect vectors.

Isolation and transcriptomic analysis of Anopheles gambiae oenocytes enables the delineation of hydrocarbon biosynthesis.

The surface of insects is coated in cuticular hydrocarbons (CHCs); variations in the composition of this layer affect a range of traits including adaptation to arid environments and defence against pathogens and toxins. In the African malaria vector, Anopheles gambiae quantitative and qualitative variance in CHC composition have been associated with speciation, ecological habitat and insecticide resistance. Understanding how these modifications arise will inform us of how mosquitoes are responding to climate change and vector control interventions. CHCs are synthesised in sub-epidermal cells called oenocytes that are very difficult to isolate from surrounding tissues. Here we utilise a transgenic line with fluorescent oenocytes to purify these cells for the first time. Comparative transcriptomics revealed the enrichment of biological processes related to long chain fatty acyl-CoA biosynthesis and elongation of mono-, poly-unsaturated and saturated fatty acids and enabled us to delineate, and partially validate, the hydrocarbon biosynthetic pathway in An. gambiae.

The bodies of insects are encased in an exoskeleton or cuticle that is key for their survival. The cuticle helps protect insects against damage, prevents water loss and can defend against pesticides. A better understanding of the role of the cuticle for survival in mosquitoes and other insects could lead to new ways to prevent the spread of diseases such as malaria. The cuticle is coated with various molecules from a group of chemicals called hydrocarbons. This coating is made by specialized cells called oenocytes and helps to protect insects. Hydrocarbons can also influence communications between certain insects by acting as recognition signals. In mosquitoes, oenocytes make several hydrocarbons using a set of processes that are not well understood, and the types of hydrocarbons they make can vary between individuals of the same species. It is unclear how this mixture of hydrocarbons is generated and how differences in the mixture can determine how mosquitoes adapt to their surroundings. Grigoraki et al. studied the genes that were active in isolated oenocytes from the mosquito Anopheles gambiae, which carries the parasite that causes malaria. The study revealed a set of genes which are highly active in oenocytes and control the production of fatty acids, a group of molecules used to make hydrocarbons. Other genes involved in creating hydrocarbons were also found. Grigoraki et al. further investigated a specific gene called FAS1899 and showed that loss of this gene reduces overall hydrocarbon production by 25%. Additionally, genes for transporting and recycling molecules and for producing fats were also shown to be active, which may indicate that oenocytes have a variety of unexplored roles besides making hydrocarbons. Grigoraki et al. identify the genes involved in producing the hydrocarbon coating of mosquitoes and demonstrate their significance. Further work is needed to understand the precise roles of each of these genes and how they are regulated to adapt the hydrocarbon coating to different situations. This can help explain how the hydrocarbon coating changes in mosquitoes, for example in response to the use of insecticides or climate change. This information is important to adapt and develop new tools to improve mosquito control.

Evolution of the Insecticide Target Rdl in African Anopheles Is Driven by Interspecific and Interkaryotypic Introgression.

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across nonconcordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors.

A unicellular relative of animals generates a layer of polarized cells by actomyosin-dependent cellularization.

In animals, cellularization of a coenocyte is a specialized form of cytokinesis that results in the formation of a polarized epithelium during early embryonic development. It is characterized by coordinated assembly of an actomyosin network, which drives inward membrane invaginations. However, whether coordinated cellularization driven by membrane invagination exists outside animals is not known. To that end, we investigate cellularization in the ichthyosporean Sphaeroforma arctica, a close unicellular relative of animals. We show that the process of cellularization involves coordinated inward plasma membrane invaginations dependent on an actomyosin network and reveal the temporal order of its assembly. This leads to the formation of a polarized layer of cells resembling an epithelium. We show that this stage is associated with tightly regulated transcriptional activation of genes involved in cell adhesion. Hereby we demonstrate the presence of a self-organized, clonally-generated, polarized layer of cells in a unicellular relative of animals.

Reconstruction of protein domain evolution using single-cell amplified genomes of uncultured choanoflagellates sheds light on the origin of animals.

Understanding the origins of animal multicellularity is a fundamental biological question. Recent genome data have unravelled the role that co-option of pre-existing genes played in the origin of animals. However, there were also some important genetic novelties at the onset of Metazoa. To have a clear understanding of the specific genetic innovations and how they appeared, we need the broadest taxon sampling possible, especially among early-branching animals and their unicellular relatives. Here, we take advantage of single-cell genomics to expand our understanding of the genomic diversity of choanoflagellates, the sister-group to animals. With these genomes, we have performed an updated and taxon-rich reconstruction of protein evolution from the Last Eukaryotic Common Ancestor (LECA) to animals. Our novel data re-defines the origin of some genes previously thought to be metazoan-specific, like the POU transcription factor, which we show appeared earlier in evolution. Moreover, our data indicate that the acquisition of new genes at the stem of Metazoa was mainly driven by duplications and protein domain rearrangement processes at the stem of Metazoa. Furthermore, our analysis allowed us to reveal protein domains that are essential to the maintenance of animal multicellularity. Our analyses also demonstrate the utility of single-cell genomics from uncultured taxa to address evolutionary questions. This article is part of a discussion meeting issue 'Single cell ecology'.

The Expansion of Inosine at the Wobble Position of tRNAs, and Its Role in the Evolution of Proteomes.

The modification of adenosine to inosine at the first position of transfer RNA (tRNA) anticodons (I34) is widespread among bacteria and eukaryotes. In bacteria, the modification is found in tRNAArg and is catalyzed by tRNA adenosine deaminase A, a homodimeric enzyme. In eukaryotes, I34 is introduced in up to eight different tRNAs by the heterodimeric adenosine deaminase acting on tRNA. This substrate expansion significantly influenced the evolution of eukaryotic genomes in terms of codon usage and tRNA gene composition. However, the selective advantages driving this process remain unclear. Here, we have studied the evolution of I34, tRNA adenosine deaminase A, adenosine deaminase acting on tRNA, and their relevant codons in a large set of bacterial and eukaryotic species. We show that a functional expansion of I34 to tRNAs other than tRNAArg also occurred within bacteria, in a process likely initiated by the emergence of unmodified A34-containing tRNAs. In eukaryotes, we report on a large variability in the use of I34 in protists, in contrast to a more uniform presence in fungi, plans, and animals. Our data support that the eukaryotic expansion of I34-tRNAs was driven by the improvement brought by these tRNAs to the synthesis of proteins highly enriched in certain amino acids.

Global transcriptome analysis of the aphelid Paraphelidium tribonemae supports the phagotrophic origin of fungi.

Aphelids are little-known phagotrophic parasites of algae whose life cycle and morphology resemble those of the parasitic rozellids (Cryptomycota, Rozellomycota). In previous phylogenetic analyses of RNA polymerase and rRNA genes, aphelids, rozellids and Microsporidia (parasites of animals) formed a clade, named Opisthosporidia, which appeared as the sister group to Fungi. However, the statistical support for the Opisthosporidia was always moderate. Here, we generated full life-cycle transcriptome data for the aphelid species Paraphelidium tribonemae. In-depth multi-gene phylogenomic analyses using several protein datasets place this aphelid as the closest relative of fungi to the exclusion of rozellids and Microsporidia. In contrast with the comparatively reduced Rozella allomycis genome, we infer a rich, free-living-like aphelid proteome, with a metabolism similar to fungi, including cellulases likely involved in algal cell-wall penetration and enzymes involved in chitin biosynthesis. Our results suggest that fungi evolved from complex aphelid-like ancestors that lost phagotrophy and became osmotrophic.

Origin of exon skipping-rich transcriptomes in animals driven by evolution of gene architecture.

BACKGROUND: Alternative splicing, particularly through intron retention and exon skipping, is a major layer of pre-translational regulation in eukaryotes. While intron retention is believed to be the most prevalent mode across non-animal eukaryotes, animals have unusually high rates of exon skipping. However, when and how this high prevalence of exon skipping evolved is unknown. Since exon skipping can greatly expand proteomes, answering these questions sheds light on the evolution of higher organismal complexity in metazoans. RESULTS: We used RNA-seq data to quantify exon skipping and intron retention frequencies across 65 eukaryotic species, with particular focus on early branching animals and unicellular holozoans. We found that only bilaterians have significantly increased their exon skipping frequencies compared to all other eukaryotic groups. Unlike in other eukaryotes, however, exon skipping in nearly all animals, including non-bilaterians, is strongly enriched for frame-preserving sequences, suggesting that exon skipping involvement in proteome expansion predated the increase in frequency. We also identified architectural features consistently associated with higher exon skipping rates within all studied eukaryotic genomes. Remarkably, these architectures became more prevalent during animal evolution, indicating co-evolution between genome architectures and exon skipping frequencies. CONCLUSION: We suggest that the increase of exon skipping rates in animals followed a two-step process. First, exon skipping in early animals became enriched for frame-preserving events. Second, bilaterian ancestors dramatically increased their exon skipping frequencies, likely driven by the interplay between a shift in their genome architectures towards more exon definition and recruitment of frame-preserving exon skipping events to functionally diversify their cell-specific proteomes.

Comparative genomic analysis of the 'pseudofungus' Hyphochytrium catenoides.

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete 'pseudofungus'; Hyphochytrium catenoides Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.