Spatial transcriptomic and single-nucleus analysis reveals heterogeneity in a gigantic single-celled syncytium.

In multicellular organisms, the specification, coordination, and compartmentalization of cell types enable the formation of complex body plans. However, some eukaryotic protists such as slime molds generate diverse and complex structures while remaining in a multinucleate syncytial state. It is unknown if different regions of these giant syncytial cells have distinct transcriptional responses to environmental encounters and if nuclei within the cell diversify into heterogeneous states. Here, we performed spatial transcriptome analysis of the slime mold Physarum polycephalum in the plasmodium state under different environmental conditions and used single-nucleus RNA-sequencing to dissect gene expression heterogeneity among nuclei. Our data identifies transcriptome regionality in the organism that associates with proliferation, syncytial substructures, and localized environmental conditions. Further, we find that nuclei are heterogenous in their transcriptional profile and may process local signals within the plasmodium to coordinate cell growth, metabolism, and reproduction. To understand how nuclei variation within the syncytium compares to heterogeneity in single-nucleus cells, we analyzed states in single Physarum amoebal cells. We observed amoebal cell states at different stages of mitosis and meiosis, and identified cytokinetic features that are specific to nuclei divisions within the syncytium. Notably, we do not find evidence for predefined transcriptomic states in the amoebae that are observed in the syncytium. Our data shows that a single-celled slime mold can control its gene expression in a region-specific manner while lacking cellular compartmentalization and suggests that nuclei are mobile processors facilitating local specialized functions. More broadly, slime molds offer the extraordinary opportunity to explore how organisms can evolve regulatory mechanisms to divide labor, specialize, balance competition with cooperation, and perform other foundational principles that govern the logic of life.

NGN2 induces diverse neuron types from human pluripotency.

Human neurons engineered from induced pluripotent stem cells (iPSCs) through neurogenin 2 (NGN2) overexpression are widely used to study neuronal differentiation mechanisms and to model neurological diseases. However, the differentiation paths and heterogeneity of emerged neurons have not been fully explored. Here, we used single-cell transcriptomics to dissect the cell states that emerge during NGN2 overexpression across a time course from pluripotency to neuron functional maturation. We find a substantial molecular heterogeneity in the neuron types generated, with at least two populations that express genes associated with neurons of the peripheral nervous system. Neuron heterogeneity is observed across multiple iPSC clones and lines from different individuals. We find that neuron fate acquisition is sensitive to NGN2 expression level and the duration of NGN2-forced expression. Our data reveal that NGN2 dosage can regulate neuron fate acquisition, and that NGN2-iN heterogeneity can confound results that are sensitive to neuron type.

Comparison of induced neurons reveals slower structural and functional maturation in humans than in apes.

We generated induced excitatory neurons (iNeurons, iNs) from chimpanzee, bonobo, and human stem cells by expressing the transcription factor neurogenin-2 (NGN2). Single-cell RNA sequencing showed that genes involved in dendrite and synapse development are expressed earlier during iNs maturation in the chimpanzee and bonobo than the human cells. In accordance, during the first 2 weeks of differentiation, chimpanzee and bonobo iNs showed repetitive action potentials and more spontaneous excitatory activity than human iNs, and extended neurites of higher total length. However, the axons of human iNs were slightly longer at 5 weeks of differentiation. The timing of the establishment of neuronal polarity did not differ between the species. Chimpanzee, bonobo, and human neurites eventually reached the same level of structural complexity. Thus, human iNs develop slower than chimpanzee and bonobo iNs, and this difference in timing likely depends on functions downstream of NGN2.

Human Stem Cell Resources Are an Inroad to Neandertal DNA Functions.

Induced pluripotent stem cells (iPSCs) from diverse humans offer the potential to study human functional variation in controlled culture environments. A portion of this variation originates from an ancient admixture between modern humans and Neandertals, which introduced alleles that left a phenotypic legacy on individual humans today. Here, we show that a large iPSC repository harbors extensive Neandertal DNA, including alleles that contribute to human phenotypes and diseases, encode hundreds of amino acid changes, and alter gene expression in specific tissues. We provide a database of the inferred introgressed Neandertal alleles for each individual iPSC line, together with the annotation of the predicted functional variants. We also show that transcriptomic data from organoids generated from iPSCs can be used to track Neandertal-derived RNA over developmental processes. Human iPSC resources provide an opportunity to experimentally explore Neandertal DNA function and its contribution to present-day phenotypes, and potentially study Neandertal traits.

Organoid single-cell genomic atlas uncovers human-specific features of brain development.

The human brain has undergone substantial change since humans diverged from chimpanzees and the other great apes1,2. However, the genetic and developmental programs that underlie this divergence are not fully understood. Here we have analysed stem cell-derived cerebral organoids using single-cell transcriptomics and accessible chromatin profiling to investigate gene-regulatory changes that are specific to humans. We first analysed cell composition and reconstructed differentiation trajectories over the entire course of human cerebral organoid development from pluripotency, through neuroectoderm and neuroepithelial stages, followed by divergence into neuronal fates within the dorsal and ventral forebrain, midbrain and hindbrain regions. Brain-region composition varied in organoids from different iPSC lines, but regional gene-expression patterns remained largely reproducible across individuals. We analysed chimpanzee and macaque cerebral organoids and found that human neuronal development occurs at a slower pace relative to the other two primates. Using pseudotemporal alignment of differentiation paths, we found that human-specific gene expression resolved to distinct cell states along progenitor-to-neuron lineages in the cortex. Chromatin accessibility was dynamic during cortex development, and we identified divergence in accessibility between human and chimpanzee that correlated with human-specific gene expression and genetic change. Finally, we mapped human-specific expression in adult prefrontal cortex using single-nucleus RNA sequencing analysis and identified developmental differences that persist into adulthood, as well as cell-state-specific changes that occur exclusively in the adult brain. Our data provide a temporal cell atlas of great ape forebrain development, and illuminate dynamic gene-regulatory features that are unique to humans.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation.

BACKGROUND: Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes. FINDINGS: Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated. CONCLUSION: The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Convergent evolution of the ladder-like ventral nerve cord in Annelida.

BACKGROUND: A median, segmented, annelid nerve cord has repeatedly been compared to the arthropod and vertebrate nerve cords and became the most used textbook representation of the annelid nervous system. Recent phylogenomic analyses, however, challenge the hypothesis that a subepidermal rope-ladder-like ventral nerve cord (VNC) composed of a paired serial chain of ganglia and somata-free connectives represents either a plesiomorphic or a typical condition in annelids. RESULTS: Using a comparative approach by combining phylogenomic analyses with morphological methods (immunohistochemistry and CLSM, histology and TEM), we compiled a comprehensive dataset to reconstruct the evolution of the annelid VNC. Our phylogenomic analyses generally support previous topologies. However, the so far hard-to-place Apistobranchidae and Psammodrilidae are now incorporated among the basally branching annelids with high support. Based on this topology we reconstruct an intraepidermal VNC as the ancestral state in Annelida. Thus, a subepidermal ladder-like nerve cord clearly represents a derived condition. CONCLUSIONS: Based on the presented data, a ladder-like appearance of the ventral nerve cord evolved repeatedly, and independently of the transition from an intraepidermal to a subepidermal cord during annelid evolution. Our investigations thereby propose an alternative set of neuroanatomical characteristics for the last common ancestor of Annelida or perhaps even Spiralia.

The complete mitochondrial genome of Lacerta bilineata and comparison with its closely related congener L. Viridis.

We sequenced the mitochondrial genome of the Western green lizard (Lacerta bilineata) using Illumina technology and additional Sanger sequencing. The assembled 17 086 bp mitogenome had a GC content of 40.32% and consisted of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region (CR), with a gene order identical to the chordate consensus. In addition, we re-sequenced the mitogenome of the closely related Eastern green lizard L. viridis using the same techniques as for L. bilineata. The mitogenomes of L. bilineata and L. viridis showed a sequence identity of 94.4% and 99.9%, respectively, relative to the previously published L. viridis mitogenome. The phylogenetic reconstruction based on 17 Lacertinae mitogenomes using Anolis carolinensis as the outgroup supported L. bilineata and its sister species L. viridis as distinct lineages.

Ancient horizontal transfers of retrotransposons between birds and ancestors of human pathogenic nematodes.

Parasite host switches may trigger disease emergence, but prehistoric host ranges are often unknowable. Lymphatic filariasis and loiasis are major human diseases caused by the insect-borne filarial nematodes Brugia, Wuchereria and Loa. Here we show that the genomes of these nematodes and seven tropical bird lineages exclusively share a novel retrotransposon, AviRTE, resulting from horizontal transfer (HT). AviRTE subfamilies exhibit 83-99% nucleotide identity between genomes, and their phylogenetic distribution, paleobiogeography and invasion times suggest that HTs involved filarial nematodes. The HTs between bird and nematode genomes took place in two pantropical waves, >25-22 million years ago (Myr ago) involving the Brugia/Wuchereria lineage and >20-17 Myr ago involving the Loa lineage. Contrary to the expectation from the mammal-dominated host range of filarial nematodes, we hypothesize that these major human pathogens may have independently evolved from bird endoparasites that formerly infected the global breadth of avian biodiversity.

Evolution of mitochondrial gene order in Annelida.

Annelida is a highly diverse animal group with over 21,000 described species. As part of Lophotrochozoa, the vast majority of annelids are currently classified into two groups: Errantia and Sedentaria, together forming Pleistoannelida. Besides these taxa, Sipuncula, Amphinomidae, Chaetopteridae, Oweniidae and Magelonidae can be found branching at the base of the tree. Comparisons of mitochondrial genomes have been used to investigate phylogenetic relationship within animal taxa. Complete annelid mitochondrial genomes are available for some Sedentaria and Errantia and in most cases exhibit a highly conserved gene order. Only two complete genomes have been published from the basal branching lineages and these are restricted to Sipuncula. We describe the first complete mitochondrial genome sequences for all other basal branching annelid families: Owenia fusiformis (Oweniidae), Magelona mirabilis (Magelonidae), Eurythoe complanata (Amphinomidae), Chaetopterus variopedatus and Phyllochaetopterus sp. (Chaetopteridae). The mitochondrial gene order of all these taxa is substantially different from the pattern found in Pleistoannelida. Additionally, we report the first mitochondrial genomes in Annelida that encode genes on both strands. Our findings demonstrate that the supposedly highly conserved mitochondrial gene order suggested for Annelida is restricted to Pleistoannelida, representing the ground pattern of this group. All investigated basal branching annelid taxa show a completely different arrangement of genes than observed in Pleistoannelida. The gene order of protein coding and ribosomal genes in Magelona mirabilis differs only in two transposition events from a putative lophotrochozoan ground pattern and might be the closest to an ancestral annelid pattern. The mitochondrial genomes of Myzostomida show the conserved pattern of Pleistoannelida, thereby supporting their inclusion in this taxon.

The evolution of annelids reveals two adaptive routes to the interstitial realm.

Many animals permanently inhabit the marine interstitium, the space between sand grains [1, 2]. Different evolutionary scenarios may explain the existence of interstitial animals [3, 4]. These scenarios include (1) that the interstitial realm is the ancestral habitat of bilaterians [5, 6], (2) that interstitial taxa evolved from larger ancestors by miniaturization, or (3) progenesis [3]. The first view mirrors the former hypothesis that interstitial annelids, called archiannelids, were at the base of the annelid radiation [7]. Based on morphological data, however, progenesis is generally favored for interstitial annelids today [3, 4, 8]. Herein, our phylogenomic approach revealed that interstitial archiannelids are robustly placed into two groups nested within the annelid tree. Evolution of the first group comprising among others Dinophilidae is best explained by progenesis. In contrast, the second group comprising Protodrilida and Polygordiidae appears to have evolved by stepwise miniaturization adapting from coarser to finer sediments. Thus, in addition to progenesis [3, 4], miniaturization, thought to be too slow for an adaptation to the interstitium [3], is an important second route allowing adaptation to interstitial environments. Both progenesis and miniaturization should be considered when investigating evolution of interstitial taxa [1, 3].

The making of a branching annelid: an analysis of complete mitochondrial genome and ribosomal data of Ramisyllis multicaudata.

Ramisyllis multicaudata is a member of Syllidae (Annelida, Errantia, Phyllodocida) with a remarkable branching body plan. Using a next-generation sequencing approach, the complete mitochondrial genomes of R. multicaudata and Trypanobia sp. are sequenced and analysed, representing the first ones from Syllidae. The gene order in these two syllids does not follow the order proposed as the putative ground pattern in Errantia. The phylogenetic relationships of R. multicaudata are discerned using a phylogenetic approach with the nuclear 18S and the mitochondrial 16S and cox1 genes. Ramisyllis multicaudata is the sister group of a clade containing Trypanobia species. Both genera, Ramisyllis and Trypanobia, together with Parahaplosyllis, Trypanosyllis, Eurysyllis, and Xenosyllis are located in a long branched clade. The long branches are explained by an accelerated mutational rate in the 18S rRNA gene. Using a phylogenetic backbone, we propose a scenario in which the postembryonic addition of segments that occurs in most syllids, their huge diversity of reproductive modes, and their ability to regenerate lost parts, in combination, have provided an evolutionary basis to develop a new branching body pattern as realised in Ramisyllis.

Phylogenomic analyses uncover origin and spread of the Wolbachia pandemic.

Of all obligate intracellular bacteria, Wolbachia is probably the most common. In general, Wolbachia are either widespread, opportunistic reproductive parasites of arthropods or essential mutualists in a single group of filarial nematodes, including many species of medical significance. To date, a robust phylogenetic backbone of Wolbachia is lacking and consequently, many Wolbachia-related phenomena cannot be discussed in a broader evolutionary context. Here we present the first comprehensive phylogenomic analysis of Wolbachia supergroup relationships based on new whole-genome-shotgun data. Our results suggest that Wolbachia has switched between its two major host groups at least twice. The ability of some arthropod-infecting Wolbachia to universally infect and to adapt to a broad range of hosts quickly is restricted to a single monophyletic lineage (containing supergroups A and B). Thus, the currently observable pandemic has likely a single evolutionary origin and is unique within the radiation of Wolbachia strains.

Platyzoan paraphyly based on phylogenomic data supports a noncoelomate ancestry of spiralia.

Based on molecular data three major clades have been recognized within Bilateria: Deuterostomia, Ecdysozoa, and Spiralia. Within Spiralia, small-sized and simply organized animals such as flatworms, gastrotrichs, and gnathostomulids have recently been grouped together as Platyzoa. However, the representation of putative platyzoans was low in the respective molecular phylogenetic studies, in terms of both, taxon number and sequence data. Furthermore, increased substitution rates in platyzoan taxa raised the possibility that monophyletic Platyzoa represents an artifact due to long-branch attraction. In order to overcome such problems, we employed a phylogenomic approach, thereby substantially increasing 1) the number of sampled species within Platyzoa and 2) species-specific sequence coverage in data sets of up to 82,162 amino acid positions. Using established and new measures (long-branch score), we disentangled phylogenetic signal from misleading effects such as long-branch attraction. In doing so, our phylogenomic analyses did not recover a monophyletic origin of platyzoan taxa that, instead, appeared paraphyletic with respect to the other spiralians. Platyhelminthes and Gastrotricha formed a monophylum, which we name Rouphozoa. To the exclusion of Gnathifera, Rouphozoa and all other spiralians represent a monophyletic group, which we name Platytrochozoa. Platyzoan paraphyly suggests that the last common ancestor of Spiralia was a simple-bodied organism lacking coelomic cavities, segmentation, and complex brain structures, and that more complex animals such as annelids evolved from such a simply organized ancestor. This conclusion contradicts alternative evolutionary scenarios proposing an annelid-like ancestor of Bilateria and Spiralia and several independent events of secondary reduction.

Illuminating the base of the annelid tree using transcriptomics.

Annelida is one of three animal groups possessing segmentation and is central in considerations about the evolution of different character traits. It has even been proposed that the bilaterian ancestor resembled an annelid. However, a robust phylogeny of Annelida, especially with respect to the basal relationships, has been lacking. Our study based on transcriptomic data comprising 68,750-170,497 amino acid sites from 305 to 622 proteins resolves annelid relationships, including Chaetopteridae, Amphinomidae, Sipuncula, Oweniidae, and Magelonidae in the basal part of the tree. Myzostomida, which have been indicated to belong to the basal radiation as well, are now found deeply nested within Annelida as sister group to Errantia in most analyses. On the basis of our reconstruction of a robust annelid phylogeny, we show that the basal branching taxa include a huge variety of life styles such as tube dwelling and deposit feeding, endobenthic and burrowing, tubicolous and filter feeding, and errant and carnivorous forms. Ancestral character state reconstruction suggests that the ancestral annelid possessed a pair of either sensory or grooved palps, bicellular eyes, biramous parapodia bearing simple chaeta, and lacked nuchal organs. Because the oldest fossil of Annelida is reported for Sipuncula (520 Ma), we infer that the early diversification of annelids took place at least in the Lower Cambrian.

Early divergence, broad distribution, and high diversity of animal chitin synthases.

Even though chitin is one of the most abundant biopolymers in nature, current knowledge on chitin formation is largely based only on data from fungi and insects. This study reveals unanticipated broad taxonomic distribution and extensive diversification of chitin synthases (CSs) in Metazoa, shedding new light on the relevance of chitin in animals and suggesting unforeseen complexity of chitin synthesis in many groups. We uncovered robust orthologs to insect type CSs in several representatives of deuterostomes, which generally are not thought to possess chitin. This suggests a broader distribution and function of chitin in this branch of the animal kingdom. We characterize a new CS type present not only in basal metazoans such as sponges and cnidarians but also in several bilaterian representatives. The most extensive diversification of CSs took place during emergence of lophotrochozoans, the third large group of protostomes next to arthropods and nematodes, resulting in coexistence of up to ten CS paralogs in molluscs. Independent fusion to different kinds of myosin motor domains in fungi and lophotrochozoans points toward high relevance of CS interaction with the cytoskeleton for fine-tuned chitin secretion. Given the fundamental role that chitin plays in the morphology of many animals, the here presented CS diversification reveals many evolutionary complexities. Our findings strongly suggest a very broad and multifarious occurrence of chitin and question an ancestral role as cuticular component. The molecular mechanisms underlying regulation of animal chitin synthesis are most likely far more complex and diverse than existing data from insects suggest.

Expression pattern of Piwi-like genes in adult Myzostoma cirriferum (Annelida).

Piwi-like genes are a subgroup of Argonaute genes which participate as gene regulators by gene silencing. In most bilaterians, such as mouse, human, insects, and zebrafish, their expression is mostly limited to gonadal stem cells. But there are some striking exceptions to this pattern; flatworms and acoels also express Piwi-like genes in somatic stem cells, due to their unique replacement system. Annelid species like Capitella teleta and Platynereis dumerilii express these genes in cells of the posterior growth zone as well as in gonadal stem cells. To investigate the expression pattern of Piwi-like genes in another annelid, we established in situ hybridization for adult Myzostoma cirriferum. Piwi-like gene transcripts recovered in an mRNA-seq library of pooled adult stages of M. cirriferum were expanded using RACE PCR, cloned and sequenced. ML analysis confirmed the identity of both transcripts as part of the Piwi1-like or Piwi2-like subfamily of Argonaute proteins. The results of in situ hybridization studies show that the expression of both Piwi-like genes, Mc-Piwi1 and Mc-Piwi2, is clearly located only in gonadal stem cells, and as such we did not find any evidence for the existence of a posterior growth zone nor expression in somatic stem cells.

Myoanatomy of Myzostoma cirriferum (Annelida, Myzostomida): implications for the evolution of the myzostomid body plan.

Studies of rare genomic marker systems suggest that Myzostomida are a subgroup of Annelida and phylogenomic analyses indicate an early divergence of this taxon within annelids. However, adult myzostomids show a highly specialized body plan, which lacks typical annelid features, such as external body annulation, coelomic cavities with metanephridia, and segmental ganglia of the nervous system. The putative loss of these features might be due to the parasitic/symbiotic lifestyle of myzostomids associated with echinoderms. In contrast, the larval anatomy and adult locomotory system resemble those of annelids. To clarify whether the myoanatomy of myzostomids reflects their relationship to annelids, we analyzed the distribution of f-actin, a common component of muscle fibers, in specimens of Myzostoma cirriferum using phalloidin-rhodamine labeling in conjunction with confocal laser-scanning microscopy. Our data reveal that the musculature of the myzostomid body comprises an outer circular layer, an inner longitudinal layer, numerous dorsoventral muscles, and prominent muscles of the parapodial complex. These features correspond well with the common organization of the muscular system in Annelida. In contrast to other annelids, however, several elements of the muscular system in M. cirriferum, including the musculature of the body wall, and the parapodial flexor muscles, exhibit radial symmetry overlaying a bilateral body plan. These findings are in line with the annelid affinity of myzostomids and suggest that the apparent partial radial symmetry of M. cirriferum arose secondarily in this species. Based on our data, we provide a scenario on the rearrangements of muscle fibers that might have taken place in the lineage leading to this species.