The genomic and cellular foundations of animal origins.

The first animals arose more than six hundred million years ago, yet they left little impression in the fossil record. Nonetheless, the cell biology and genome composition of the first animal, the Urmetazoan, can be reconstructed through the study of phylogenetically relevant living organisms. Comparisons among animals and their unicellular and colonial relatives reveal that the Urmetazoan likely possessed a layer of epithelium-like collar cells, preyed on bacteria, reproduced by sperm and egg, and developed through cell division, cell differentiation, and invagination. Although many genes involved in development, body patterning, immunity, and cell-type specification evolved in the animal stem lineage or after animal origins, several gene families critical for cell adhesion, signaling, and gene regulation predate the origin of animals. The ancestral functions of these and other genes may eventually be revealed through studies of gene and genome function in early-branching animals and their closest non-animal relatives.

Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes.

This revision of the classification of eukaryotes follows that of Adl et al., 2012 [J. Euk. Microbiol. 59(5)] and retains an emphasis on protists. Changes since have improved the resolution of many nodes in phylogenetic analyses. For some clades even families are being clearly resolved. As we had predicted, environmental sampling in the intervening years has massively increased the genetic information at hand. Consequently, we have discovered novel clades, exciting new genera and uncovered a massive species level diversity beyond the morphological species descriptions. Several clades known from environmental samples only have now found their home. Sampling soils, deeper marine waters and the deep sea will continue to fill us with surprises. The main changes in this revision are the confirmation that eukaryotes form at least two domains, the loss of monophyly in the Excavata, robust support for the Haptista and Cryptista. We provide suggested primer sets for DNA sequences from environmental samples that are effective for each clade. We have provided a guide to trophic functional guilds in an appendix, to facilitate the interpretation of environmental samples, and a standardized taxonomic guide for East Asian users.

The Evolution of Silicon Transport in Eukaryotes.

Biosilicification (the formation of biological structures from silica) occurs in diverse eukaryotic lineages, plays a major role in global biogeochemical cycles, and has significant biotechnological applications. Silicon (Si) uptake is crucial for biosilicification, yet the evolutionary history of the transporters involved remains poorly known. Recent evidence suggests that the SIT family of Si transporters, initially identified in diatoms, may be widely distributed, with an extended family of related transporters (SIT-Ls) present in some nonsilicified organisms. Here, we identify SITs and SIT-Ls in a range of eukaryotes, including major silicified lineages (radiolarians and chrysophytes) and also bacterial SIT-Ls. Our evidence suggests that the symmetrical 10-transmembrane-domain SIT structure has independently evolved multiple times via duplication and fusion of 5-transmembrane-domain SIT-Ls. We also identify a second gene family, similar to the active Si transporter Lsi2, that is broadly distributed amongst siliceous and nonsiliceous eukaryotes. Our analyses resolve a distinct group of Lsi2-like genes, including plant and diatom Si-responsive genes, and sequences unique to siliceous sponges and choanoflagellates. The SIT/SIT-L and Lsi2 transporter families likely contribute to biosilicification in diverse lineages, indicating an ancient role for Si transport in eukaryotes. We propose that these Si transporters may have arisen initially to prevent Si toxicity in the high Si Precambrian oceans, with subsequent biologically induced reductions in Si concentrations of Phanerozoic seas leading to widespread losses of SIT, SIT-L, and Lsi2-like genes in diverse lineages. Thus, the origin and diversification of two independent Si transporter families both drove and were driven by ancient ocean Si levels.

A six-gene phylogeny provides new insights into choanoflagellate evolution.

Recent studies have shown that molecular phylogenies of the choanoflagellates (Class Choanoflagellatea) are in disagreement with their traditional taxonomy, based on morphology, and that Choanoflagellatea requires considerable taxonomic revision. Furthermore, phylogenies suggest that the morphological and ecological evolution of the group is more complex than has previously been recognized. Here we address the taxonomy of the major choanoflagellate order Craspedida, by erecting four new genera. The new genera are shown to be morphologically, ecologically and phylogenetically distinct from other choanoflagellate taxa. Furthermore, we name five novel craspedid species, as well as formally describe ten species that have been shown to be either misidentified or require taxonomic revision. Our revised phylogeny, including 18 new species and sequence data for two additional genes, provides insights into the morphological and ecological evolution of the choanoflagellates. We examine the distribution within choanoflagellates of these two additional genes, EF-1A and EFL, closely related translation GTPases which are required for protein synthesis. Mapping the presence and absence of these genes onto the phylogeny highlights multiple events of gene loss within the choanoflagellates.

Two levels of protection for the B cell genome during somatic hypermutation.

Somatic hypermutation introduces point mutations into immunoglobulin genes in germinal centre B cells during an immune response. The reaction is initiated by cytosine deamination by the activation-induced deaminase (AID) and completed by error-prone processing of the resulting uracils by mismatch and base excision repair factors. Somatic hypermutation represents a threat to genome integrity and it is not known how the B cell genome is protected from the mutagenic effects of somatic hypermutation nor how often these protective mechanisms fail. Here we show, by extensive sequencing of murine B cell genes, that the genome is protected by two distinct mechanisms: selective targeting of AID and gene-specific, high-fidelity repair of AID-generated uracils. Numerous genes linked to B cell tumorigenesis, including Myc, Pim1, Pax5, Ocab (also called Pou2af1), H2afx, Rhoh and Ebf1, are deaminated by AID but escape acquisition of most mutations through the combined action of mismatch and base excision repair. However, approximately 25% of expressed genes analysed were not fully protected by either mechanism and accumulated mutations in germinal centre B cells. Our results demonstrate that AID acts broadly on the genome, with the ultimate distribution of mutations determined by a balance between high-fidelity and error-prone DNA repair.

The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans.

Choanoflagellates are the closest known relatives of metazoans. To discover potential molecular mechanisms underlying the evolution of metazoan multicellularity, we sequenced and analysed the genome of the unicellular choanoflagellate Monosiga brevicollis. The genome contains approximately 9,200 intron-rich genes, including a number that encode cell adhesion and signalling protein domains that are otherwise restricted to metazoans. Here we show that the physical linkages among protein domains often differ between M. brevicollis and metazoans, suggesting that abundant domain shuffling followed the separation of the choanoflagellate and metazoan lineages. The completion of the M. brevicollis genome allows us to reconstruct with increasing resolution the genomic changes that accompanied the origin of metazoans.

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes.

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

Dollo Parsimony Overestimates Ancestral Gene Content Reconstructions.

Ancestral reconstruction is a widely used technique that has been applied to understand the evolutionary history of gain and loss of gene families. Ancestral gene content can be reconstructed via different phylogenetic methods, but many current and previous studies employ Dollo parsimony. We hypothesize that Dollo parsimony is not appropriate for ancestral gene content reconstruction inferences based on sequence homology, as Dollo parsimony is derived from the assumption that a complex character cannot be regained. This premise does not accurately model molecular sequence evolution, in which false orthology can result from sequence convergence or lateral gene transfer. The aim of this study is to test Dollo parsimony's suitability for ancestral gene content reconstruction and to compare its inferences with a maximum likelihood-based approach that allows a gene family to be gained more than once within a tree. We first compared the performance of the two approaches on a series of artificial data sets each of 5,000 genes that were simulated according to a spectrum of evolutionary rates without gene gain or loss, so that inferred deviations from the true gene count would arise only from errors in orthology inference and ancestral reconstruction. Next, we reconstructed protein domain evolution on a phylogeny representing known eukaryotic diversity. We observed that Dollo parsimony produced numerous ancestral gene content overestimations, especially at nodes closer to the root of the tree. These observations led us to the conclusion that, confirming our hypothesis, Dollo parsimony is not an appropriate method for ancestral reconstruction studies based on sequence homology.

The protist cultural renaissance.

Protists are key players in the biosphere. Here, we provide a perspective on integrating protist culturing with omics approaches, imaging, and high-throughput single-cell manipulation strategies, concluding with actions required for a successful return of the golden age of protist culturing.

Human genome sequence variation and the influence of gene history, mutation and recombination.

Variation in the human genome sequence is key to understanding susceptibility to disease in modern populations and the history of ancestral populations. Unlocking this information requires knowledge of the patterns and underlying causes of human sequence diversity. By applying a new population-genetic framework to two genome-wide polymorphism surveys, we find that the human genome contains sizeable regions (stretching over tens of thousands of base pairs) that have intrinsically high and low rates of sequence variation. We show that the primary determinant of these patterns is shared genealogical history. Only a fraction of the variation (at most 25%) is due to the local mutation rate. By measuring the average distance over which genealogical histories are typically preserved, these data provide the first genome-wide estimate of the average extent of correlation among variants (linkage disequilibrium). The results are best explained by extreme variability in the recombination rate at a fine scale, and provide the first empirical evidence that such recombination 'hot spots' are a general feature of the human genome and have a principal role in shaping genetic variation in the human population.

Taxonomic composition, community structure and molecular novelty of microeukaryotes in a temperate oligomesotrophic lake as revealed by metabarcoding.

Microbial eukaryotes are diverse and ecologically important organisms, yet sampling constraints have hindered the understanding of their distribution and diversity in freshwater ecosystems. Metabarcoding has provided a powerful complement to traditional limnological studies, revealing an unprecedented diversity of protists in freshwater environments. Here, we aim to expand our knowledge of the ecology and diversity of protists in lacustrine ecosystems by targeting the V4 hypervariable region of the 18S rRNA gene in water column, sediment and biofilm samples collected from Sanabria Lake (Spain) and surrounding freshwater ecosystems. Sanabria is a temperate lake, which are relatively understudied by metabarcoding in comparison to alpine and polar lakes. The phylogenetic diversity of microbial eukaryotes detected in Sanabria spans all currently recognized eukaryotic supergroups, with Stramenopiles being the most abundant and diverse supergroup in all sampling sites. Parasitic microeukaryotes account for 21% of the total protist ASVs identified in our study and were dominated by Chytridiomycota, both in terms of richness and abundance, in all sampling sites. Sediments, biofilms and water column samples harbour distinct microbial communities. Phylogenetic placement of poorly assigned and abundant ASVs indicates molecular novelty inside Rhodophyta, Bigyra, early-branching Nucletmycea and Apusomonadida. In addition, we report the first freshwater incidence of the previously exclusively marine genera Abeoforma and Sphaeroforma. Our results contribute to a deeper understanding of microeukaryotic communities in freshwater ecosystems, and provide the first molecular reference for future biomonitoring surveys in Sanabria Lake.

Genetic evidence for complex speciation of humans and chimpanzees.

The genetic divergence time between two species varies substantially across the genome, conveying important information about the timing and process of speciation. Here we develop a framework for studying this variation and apply it to about 20 million base pairs of aligned sequence from humans, chimpanzees, gorillas and more distantly related primates. Human-chimpanzee genetic divergence varies from less than 84% to more than 147% of the average, a range of more than 4 million years. Our analysis also shows that human-chimpanzee speciation occurred less than 6.3 million years ago and probably more recently, conflicting with some interpretations of ancient fossils. Most strikingly, chromosome X shows an extremely young genetic divergence time, close to the genome minimum along nearly its entire length. These unexpected features would be explained if the human and chimpanzee lineages initially diverged, then later exchanged genes before separating permanently.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems.

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the 'seascape') by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Oceans are brimming with life invisible to our eyes, a myriad of species of bacteria, viruses and other microscopic organisms essential for the health of the planet. These 'marine plankton' are unable to swim against currents and should therefore be constantly on the move, yet previous studies have suggested that distinct species of plankton may in fact inhabit different oceanic regions. However, proving this theory has been challenging; collecting plankton is logistically difficult, and it is often impossible to distinguish between species simply by examining them under a microscope. However, within the last decade, a research schooner called Tara has travelled the globe to gather thousands of plankton samples. At the same time, advances in genomics have made it possible to identify species based only on fragments of their DNA sequence. To understand the hidden geography of plankton communities in Earth's oceans, Richter et al. pored over DNA from the Tara Oceans expedition. This revealed that, despite being unable to resist the flow of water, various planktonic species which live close to the surface manage to occupy distinct, stable provinces shaped by currents. Different sizes of plankton are distributed in different sized provinces, with the smallest organisms tending to inhabit the smallest areas. Comparing DNA similarities and speeds of currents at the ocean surface revealed how these might stretch and mix plankton communities. Plankton play a critical role in the health of the ocean and the chemical cycles of planet Earth. These results could allow deeper investigation by marine modellers, ecologists, and evolutionary biologists. Meanwhile, work is already underway to investigate how climate change might impact this hidden geography.

Analysis of chimpanzee history based on genome sequence alignments.

Population geneticists often study small numbers of carefully chosen loci, but it has become possible to obtain orders of magnitude for more data from overlaps of genome sequences. Here, we generate tens of millions of base pairs of multiple sequence alignments from combinations of three western chimpanzees, three central chimpanzees, an eastern chimpanzee, a bonobo, a human, an orangutan, and a macaque. Analysis provides a more precise understanding of demographic history than was previously available. We show that bonobos and common chimpanzees were separated approximately 1,290,000 years ago, western and other common chimpanzees approximately 510,000 years ago, and eastern and central chimpanzees at least 50,000 years ago. We infer that the central chimpanzee population size increased by at least a factor of 4 since its separation from western chimpanzees, while the western chimpanzee effective population size decreased. Surprisingly, in about one percent of the genome, the genetic relationships between humans, chimpanzees, and bonobos appear to be different from the species relationships. We used PCR-based resequencing to confirm 11 regions where chimpanzees and bonobos are not most closely related. Study of such loci should provide information about the period of time 5-7 million years ago when the ancestors of humans separated from those of the chimpanzees.

Drug-sensitive FGFR2 mutations in endometrial carcinoma.

Oncogenic activation of tyrosine kinases is a common mechanism of carcinogenesis and, given the druggable nature of these enzymes, an attractive target for anticancer therapy. Here, we show that somatic mutations of the fibroblast growth factor receptor 2 (FGFR2) tyrosine kinase gene, FGFR2, are present in 12% of endometrial carcinomas, with additional instances found in lung squamous cell carcinoma and cervical carcinoma. These FGFR2 mutations, many of which are identical to mutations associated with congenital craniofacial developmental disorders, are constitutively activated and oncogenic when ectopically expressed in NIH 3T3 cells. Inhibition of FGFR2 kinase activity in endometrial carcinoma cell lines bearing such FGFR2 mutations inhibits transformation and survival, implicating FGFR2 as a novel therapeutic target in endometrial carcinoma.

The origin and evolution of cell-intrinsic antibacterial defenses in eukaryotes.

To survive in a world dominated by bacteria, eukaryotes have evolved numerous self-defense strategies. While some defenses are recent evolutionary innovations, others are ancient, with roots early in eukaryotic history. With a focus on antibacterial immunity, we highlight the evolution of pattern recognition receptors that detect bacteria, where diverse functional classes have been formed from the repeated use and reuse of a small set of protein domains. Next, we discuss core microbicidal strategies shared across eukaryotes, and how these systems may have been co-opted from ancient cellular mechanisms. We propose that studying antibacterial responses across diverse eukaryotes can reveal novel modes of defense, while highlighting the critical innovations that occurred early in the evolution of our own immune systems.

Gene family innovation, conservation and loss on the animal stem lineage.

Choanoflagellates, the closest living relatives of animals, can provide unique insights into the changes in gene content that preceded the origin of animals. However, only two choanoflagellate genomes are currently available, providing poor coverage of their diversity. We sequenced transcriptomes of 19 additional choanoflagellate species to produce a comprehensive reconstruction of the gains and losses that shaped the ancestral animal gene repertoire. We identified ~1944 gene families that originated on the animal stem lineage, of which only 39 are conserved across all animals in our study. In addition, ~372 gene families previously thought to be animal-specific, including Notch, Delta, and homologs of the animal Toll-like receptor genes, instead evolved prior to the animal-choanoflagellate divergence. Our findings contribute to an increasingly detailed portrait of the gene families that defined the biology of the Urmetazoan and that may underpin core features of extant animals.

A global ocean atlas of eukaryotic genes.

While our knowledge about the roles of microbes and viruses in the ocean has increased tremendously due to recent advances in genomics and metagenomics, research on marine microbial eukaryotes and zooplankton has benefited much less from these new technologies because of their larger genomes, their enormous diversity, and largely unexplored physiologies. Here, we use a metatranscriptomics approach to capture expressed genes in open ocean Tara Oceans stations across four organismal size fractions. The individual sequence reads cluster into 116 million unigenes representing the largest reference collection of eukaryotic transcripts from any single biome. The catalog is used to unveil functions expressed by eukaryotic marine plankton, and to assess their functional biogeography. Almost half of the sequences have no similarity with known proteins, and a great number belong to new gene families with a restricted distribution in the ocean. Overall, the resource provides the foundations for exploring the roles of marine eukaryotes in ocean ecology and biogeochemistry.

Bridging the gap between morphological species and molecular barcodes - Exemplified by loricate choanoflagellates.

Translating the vast amounts of molecular barcodes from global surveys of microbial eukaryotes into ecological insight depends critically on a well-curated reference database with adequate taxonomic coverage. In this respect, the choanoflagellates resemble other eukaryotic lineages: reasonable coverage at higher taxonomic levels, but missing diversity at the species level. The acanthoecid (loricate) choanoflagellates are well-characterized morphologically, with over 115 species described, but less than 10% with any sequence data. Because lorica shape is species-specific, the acanthoecids represent an opportunity to link morphological with molecular data within a lineage of eukaryotes. To match morphospecies to sequences, we sampled the Kattegat and the Isefjord in Denmark in September 2014 and February 2015. We identified 45 morphospecies and sequenced ribosomal DNA of nine previously unsequenced species, roughly doubling the number of acanthoecid species with sequence data, including the first data representing five genera: Bicosta, Calliacantha, Cosmoeca, Crinolina and Pleurasiga. Our phylogenetic analysis is mainly congruent with morphology-based systematics. Five of the newly sequenced species match a previously unidentified barcode from Tara Oceans, providing access to the global distribution of species isolated from Danish waters. One species, Calliacantha natans, is the second most globally abundant choanoflagellate present in Tara Oceans. Our project translating new ribosomal DNA sequences to distributions of described species on a global scale supports the approach linking morphology to molecular barcodes for microbial eukaryote ecology.

Transcriptome sequencing and delimitation of sympatric Oscarella species (O. carmela and O. pearsei sp. nov) from California, USA.

The homoscleromorph sponge Oscarella carmela, first described from central California, USA is shown to represent two superficially similar but both morphologically and phylogenetically distinct species that are co-distributed. We here describe a new species as Oscarella pearsei, sp. nov. and re-describe Oscarella carmela; the original description was based upon material from both species. Further, we correct the identification of published genomic/transcriptomic resources that were originally attributed to O. carmela, and present new Illumina-sequenced transcriptome assemblies for each of these species, and the mitochondrial genome sequence for O. pearsei sp. nov. Using SSU and LSU ribosomal DNA and the mitochondrial genome, we report the phylogenetic relationships of these species relative to other Oscarella species, and find strong support for the placement of O. pearsei sp. nov. in a distinct clade within genus Oscarella defined by the presence of spherulous cells that contain paracrystalline inclusions; O. carmela lacks this cell type. Oscarella pearsei sp. nov and O. carmela can be tentatively distinguished based upon gross morphological differences such as color, surface texture and extent of mucus production, but can be more reliably identified using mitochondrial and nuclear barcode sequencing, ultrastructural characteristics of cells in the mesohyl, and the morphology of the follicle epithelium which surrounds the developing embryo in reproductively active individuals.