Identification of a new gregarine parasite associated with mass mortality events of freshwater pearl mussels (Margaritifera margaritifera) in Sweden.

Freshwater bivalves play key ecological roles in lakes and rivers, largely contributing to healthy ecosystems. The freshwater pearl mussel, Margaritifera margaritifera, is found in Europe and on the East coast of North America. Once common in oxygenated streams, M. margaritifera is rapidly declining and consequently assessed as a threatened species worldwide. Deterioration of water quality has been considered the main factor for the mass mortality events affecting this species. Yet, the role of parasitic infections has not been investigated. Here, we report the discovery of three novel protist lineages found in Swedish populations of M. margaritifera belonging to one of the terrestrial groups of gregarines (Eugregarinorida, Apicomplexa). These lineages are closely related-but clearly separated-from the tadpole parasite Nematopsis temporariae. In one lineage, which is specifically associated with mortality events of M. margaritifera, we found cysts containing single vermiform zoites in the gills and other organs of diseased individuals using microscopy and in situ hybridization. This represents the first report of a parasitic infection in M. margaritifera that may be linked to the decline of this mussel species. We propose a tentative life cycle with the distribution of different developmental stages and potential exit from the host into the environment.

Unexpected cryptic species among streptophyte algae most distant to land plants.

Streptophytes are one of the major groups of the green lineage (Chloroplastida or Viridiplantae). During one billion years of evolution, streptophytes have radiated into an astounding diversity of uni- and multicellular green algae as well as land plants. Most divergent from land plants is a clade formed by Mesostigmatophyceae, Spirotaenia spp. and Chlorokybophyceae. All three lineages are species-poor and the Chlorokybophyceae consist of a single described species, Chlorokybus atmophyticus. In this study, we used phylogenomic analyses to shed light into the diversity within Chlorokybus using a sampling of isolates across its known distribution. We uncovered a consistent deep genetic structure within the Chlorokybus isolates, which prompted us to formally extend the Chlorokybophyceae by describing four new species. Gene expression differences among Chlorokybus species suggest certain constitutive variability that might influence their response to environmental factors. Failure to account for this diversity can hamper comparative genomic studies aiming to understand the evolution of stress response across streptophytes. Our data highlight that future studies on the evolution of plant form and function can tap into an unknown diversity at key deep branches of the streptophytes.

New Phylogenomic Analysis of the Enigmatic Phylum Telonemia Further Resolves the Eukaryote Tree of Life.

The resolution of the broad-scale tree of eukaryotes is constantly improving, but the evolutionary origin of several major groups remains unknown. Resolving the phylogenetic position of these "orphan" groups is important, especially those that originated early in evolution, because they represent missing evolutionary links between established groups. Telonemia is one such orphan taxon for which little is known. The group is composed of molecularly diverse biflagellated protists, often prevalent although not abundant in aquatic environments. Telonemia has been hypothesized to represent a deeply diverging eukaryotic phylum but no consensus exists as to where it is placed in the tree. Here, we established cultures and report the phylogenomic analyses of three new transcriptome data sets for divergent telonemid lineages. All our phylogenetic reconstructions, based on 248 genes and using site-heterogeneous mixture models, robustly resolve the evolutionary origin of Telonemia as sister to the Sar supergroup. This grouping remains well supported when as few as 60% of the genes are randomly subsampled, thus is not sensitive to the sets of genes used but requires a minimal alignment length to recover enough phylogenetic signal. Telonemia occupies a crucial position in the tree to examine the origin of Sar, one of the most lineage-rich eukaryote supergroups. We propose the moniker "TSAR" to accommodate this new mega-assemblage in the phylogeny of eukaryotes.

Benchmarking long-read sequencing strategies for obtaining ASV-resolved rRNA operons from environmental microeukaryotes.

The use of short-read metabarcoding for classifying microeukaryotes is challenged by the lack of comprehensive 18S rRNA reference databases. While recent advances in high-throughput long-read sequencing provide the potential to greatly increase the phylogenetic coverage of these databases, the performance of different sequencing technologies and subsequent bioinformatics processing remain to be evaluated, primarily because of the absence of well-defined eukaryotic mock communities. To address this challenge, we created a eukaryotic rRNA operon clone-library and turned it into a precisely defined synthetic eukaryotic mock community. This mock community was then used to evaluate the performance of three long-read sequencing strategies (PacBio circular consensus sequencing and two Nanopore approaches using unique molecular identifiers) and three tools for resolving amplicons sequence variants (ASVs) (USEARCH, VSEARCH, and DADA2). We investigated the sensitivity of the sequencing techniques based on the number of detected mock taxa, and the accuracy of the different ASV-calling tools with a specific focus on the presence of chimera among the final rRNA operon ASVs. Based on our findings, we provide recommendations and best practice protocols for how to cost-effectively obtain essentially error-free rRNA operons in high-throughput. An agricultural soil sample was used to demonstrate that the sequencing and bioinformatic results from the mock community also translates to highly diverse natural samples, which enables us to identify previously undescribed microeukaryotic lineages.

Global patterns and rates of habitat transitions across the eukaryotic tree of life.

The successful colonization of new habitats has played a fundamental role during the evolution of life. Salinity is one of the strongest barriers for organisms to cross, which has resulted in the evolution of distinct marine and non-marine (including both freshwater and soil) communities. Although microbes represent by far the vast majority of eukaryote diversity, the role of the salt barrier in shaping the diversity across the eukaryotic tree is poorly known. Traditional views suggest rare and ancient marine/non-marine transitions but this view is being challenged by the discovery of several recently transitioned lineages. Here, we investigate habitat evolution across the tree of eukaryotes using a unique set of taxon-rich phylogenies inferred from a combination of long-read and short-read environmental metabarcoding data spanning the ribosomal DNA operon. Our results show that, overall, marine and non-marine microbial communities are phylogenetically distinct but transitions have occurred in both directions in almost all major eukaryotic lineages, with hundreds of transition events detected. Some groups have experienced relatively high rates of transitions, most notably fungi for which crossing the salt barrier has probably been an important aspect of their successful diversification. At the deepest phylogenetic levels, ancestral habitat reconstruction analyses suggest that eukaryotes may have first evolved in non-marine habitats and that the two largest known eukaryotic assemblages (TSAR and Amorphea) arose in different habitats. Overall, our findings indicate that the salt barrier has played an important role during eukaryote evolution and provide a global perspective on habitat transitions in this domain of life.

Unravelling the Phylogeny of a Common Intestinal Protist: Intrageneric Diversity of Endolimax.

Endolimax nana is a common endobiont of the human intestine, but members of the genus have also been reported in non-human hosts and in non-intestinal organs. Limited information is available regarding the genetic diversity of Endolimax, which is necessary to delineate species, host specificity and potential differences in clinical impact on the host. Here, we used cloning of PCR products followed by Sanger sequencing and next-generation PacBio Sequencing to obtain Endolimax-related nuclear ribosomal gene sequences and undertook a phylogenetic analysis to gain additional insight into the taxonomy of Endolimax and related organisms. The new sequences confirmed that E. nana forms a discrete clade within the Archamoebae and is related to Endolimax piscium and Iodamoeba. However, we identified substantial sequence divergence within E. nana and evidence for two distinct clades, which we propose to name E. nana ribosomal lineage 1 and E. nana ribosomal lineage 2. Both of the sequencing approaches applied in the study helped us to improve our understanding of genetic diversity across Endolimax, and it is likely that wider application of next-generation sequencing technologies will facilitate the generation of Endolimax-related DNA sequence data and help complete our understanding of its phylogenetic position and intrageneric diversity.

metaPR2 : A database of eukaryotic 18S rRNA metabarcodes with an emphasis on protists.

In recent years, metabarcoding has become the method of choice for investigating the composition and assembly of microbial eukaryotic communities. The number of environmental data sets published has increased very rapidly. Although unprocessed sequence files are often publicly available, processed data, in particular clustered sequences, are rarely available in a usable format. Clustered sequences are reported as operational taxonomic units (OTUs) with different similarity levels or more recently as amplicon sequence variants (ASVs). This hampers comparative studies between different environments and data sets, for example examining the biogeographical patterns of specific groups/species, as well analysing the genetic microdiversity within these groups. Here, we present a newly-assembled database of processed 18S rRNA metabarcodes that are annotated with the PR2 reference sequence database. This database, called metaPR2 , contains 41 data sets corresponding to more than 4000 samples and 90,000 ASVs. The database, which is accessible through both a web-based interface (https://shiny.metapr2.org) and an R package, should prove very useful to all researchers working on protist diversity in a variety of systems.

On the origin of TSAR: morphology, diversity and phylogeny of Telonemia.

Telonemia is a poorly known major phylum of flagellated eukaryotes with a unique combination of morphological traits. Phylogenomics recently revealed the phylogenetic position of telonemids as sister to SAR, one of the largest groups of eukaryotes, comprising Stramenopiles, Alveolata and Rhizaria. Due to this key evolutionary position, investigations of telonemids are of critical importance for elucidating the origin and diversification of an astounding diversity of eukaryotic forms and life strategies. To date, however, only two species have been morphologically characterized from Telonemia, which do not represent this genetically very diverse group. In this study, we established cultures for six new telonemid strains, including the description of five new species and a new genus. We used these cultures to update the phylogeny of Telonemia and provide a detailed morphological and ultrastructural investigation. Our data elucidate the origin of TSAR from flagellates with complex morphology and reconstruction of the ancestral structure of stramenopiles, alveolates and rhizarians, and their main synapomorphic characters. Since telonemids are a common component of aquatic environments, the features of their feeding, behaviour and ecological preferences observed in clonal cultures and the results of global metabarcoding analysis contribute to a deeper understanding of organization of microbial food webs.

Diversity and ecology of protists revealed by metabarcoding.

Protists are the dominant eukaryotes in the biosphere where they play key functional roles. While protists have been studied for over a century, it is the high-throughput sequencing of molecular markers from environmental samples - the approach of metabarcoding - that has revealed just how diverse, and abundant, these small organisms are. Metabarcoding is now routine to survey environmental diversity, so data have rapidly accumulated from a multitude of environments and at different sampling scales. This mass of data has provided unprecedented opportunities to study the taxonomic and functional diversity of protists, and how this diversity is organised in space and time. Here, we use metabarcoding as a common thread to discuss the state of knowledge in protist diversity research, from technical considerations of the approach to important insights gained on diversity patterns and the processes that might have structured this diversity. In addition to these insights, we conclude that metabarcoding is on the verge of an exciting added dimension thanks to the maturation of high-throughput long-read sequencing, so that a robust eco-evolutionary framework of protist diversity is within reach.

Long-read metabarcoding of the eukaryotic rDNA operon to phylogenetically and taxonomically resolve environmental diversity.

High-throughput DNA metabarcoding of amplicon sizes below 500 bp has revolutionized the analysis of environmental microbial diversity. However, these short regions contain limited phylogenetic signal, which makes it impractical to use environmental DNA in full phylogenetic inferences. This lesser phylogenetic resolution of short amplicons may be overcome by new long-read sequencing technologies. To test this idea, we amplified soil DNA and used PacBio Circular Consensus Sequencing (CCS) to obtain an ~4500-bp region spanning most of the eukaryotic small subunit (18S) and large subunit (28S) ribosomal DNA genes. We first treated the CCS reads with a novel curation workflow, generating 650 high-quality operational taxonomic units (OTUs) containing the physically linked 18S and 28S regions. To assign taxonomy to these OTUs, we developed a phylogeny-aware approach based on the 18S region that showed greater accuracy and sensitivity than similarity-based methods. The taxonomically annotated OTUs were then combined with available 18S and 28S reference sequences to infer a well-resolved phylogeny spanning all major groups of eukaryotes, allowing us to accurately derive the evolutionary origin of environmental diversity. A total of 1,019 sequences were included, of which a majority (58%) corresponded to the new long environmental OTUs. The long reads also allowed us to directly investigate the relationships among environmental sequences themselves, which represents a key advantage over the placement of short reads on a reference phylogeny. Together, our results show that long amplicons can be treated in a full phylogenetic framework to provide greater taxonomic resolution and a robust evolutionary perspective to environmental DNA.

RNA Sequencing of Stentor Cell Fragments Reveals Transcriptional Changes during Cellular Regeneration.

While ciliates of the genus Stentor are known for their ability to regenerate when their cells are damaged or even fragmented, the physical and molecular mechanisms underlying this process are poorly understood. To identify genes involved in the regenerative capability of Stentor cells, RNA sequencing of individual Stentor polymorphus cell fragments was performed. After splitting a cell over the anterior-posterior axis, the posterior fragment has to regenerate the oral apparatus, while the anterior part needs to regenerate the hold fast. Altogether, differential expression analysis of both posterior and anterior S. polymorphus cell fragments for four different post-split time points revealed over 10,000 upregulated genes throughout the regeneration process. Among these, genes involved in cell signaling, microtubule-based movement, and cell cycle regulation seemed to be particularly important during cellular regeneration. We identified roughly nine times as many upregulated genes in regenerating S. polymorphus posterior fragments as compared to anterior fragments, indicating that regeneration of the anterior oral apparatus is a complex process that involves many genes. Our analyses identified several expanded groups of genes, such as dual-specific tyrosine-(Y)-phosphorylation-regulated kinases and MORN domain-containing proteins that seemingly act as key regulators of cellular regeneration. In agreement with earlier morphological and cell biological studies [1, 2], our differential expression analyses indicate that cellular regeneration and vegetative division share many similarities.

Genetic Diversity in the UV Sex Chromosomes of the Brown Alga Ectocarpus.

Three types of sex chromosome system exist in nature: diploid XY and ZW systems and haploid UV systems. For many years, research has focused exclusively on XY and ZW systems, leaving UV chromosomes and haploid sex determination largely neglected. Here, we perform a detailed analysis of DNA sequence neutral diversity levels across the U and V sex chromosomes of the model brown alga Ectocarpus using a large population dataset. We show that the U and V non-recombining regions of the sex chromosomes (SDR) exhibit about half as much neutral diversity as the autosomes. This difference is consistent with the reduced effective population size of these regions compared with the rest of the genome, suggesting that the influence of additional factors such as background selection or selective sweeps is minimal. The pseudoautosomal region (PAR) of this UV system, in contrast, exhibited surprisingly high neutral diversity and there were several indications that genes in this region may be under balancing selection. The PAR of Ectocarpus is known to exhibit unusual genomic features and our results lay the foundation for further work aimed at understanding whether, and to what extent, these structural features underlie the high level of genetic diversity. Overall, this study fills a gap between available information on genetic diversity in XY/ZW systems and UV systems and significantly contributes to advancing our knowledge of the evolution of UV sex chromosomes.