CBOL protist working group: barcoding eukaryotic richness beyond the animal, plant, and fungal kingdoms

DeepMicroClass sorts metagenomic contigs into prokaryotes, eukaryotes and viruses

Sequence classification facilitates a fundamental understanding of the structure of microbial communities. Binary metagenomic sequence classifiers are insufficient because environmental metagenomes are typically derived from multiple sequence sources. Here we introduce a deep-learning based sequence classifier, DeepMicroClass, that classifies metagenomic contigs into five sequence classes, i.e. viruses infecting prokaryotic or eukaryotic hosts, eukaryotic or prokaryotic chromosomes, and prokaryotic plasmids. DeepMicroClass achieved high performance for all sequence classes at various tested sequence lengths ranging from 500 bp to 100 kbps. By benchmarking on a synthetic dataset with variable sequence class composition, we showed that DeepMicroClass obtained better performance for eukaryotic, plasmid and viral contig classification than other state-of-the-art predictors. DeepMicroClass achieved comparable performance on viral sequence classification with geNomad and VirSorter2 when benchmarked on the CAMI II marine dataset. Using a coastal daily time-series metagenomic dataset as a case study, we showed that microbial eukaryotes and prokaryotic viruses are integral to microbial communities. By analyzing monthly metagenomes collected at HOT and BATS, we found relatively higher viral read proportions in the subsurface layer in late summer, consistent with the seasonal viral infection patterns prevalent in these areas. We expect DeepMicroClass will promote metagenomic studies of under-appreciated sequence types.

Habitat conditions and not moss composition mediate microbial community structure in Swiss peatlands

Peatlands, one of the oldest ecosystems, globally store significant amounts of carbon and freshwater. However, they are under severe threat from human activities, leading to changes in water, nutrient and temperature regimes in these delicate systems. Such shifts can trigger a substantial carbon flux into the atmosphere and diminish the water-holding capacity of peatlands. Microbes associated with moss in peatlands play a crucial role in providing these ecosystem services, which are at risk due to global change. Therefore, understanding the factors influencing microbial composition and function is vital. Our study focused on five peatlands along an altitudinal gradient in Switzerland, where we sampled moss on hummocks containing Sarracenia purpurea. Structural equation modelling revealed that habitat condition was the primary predictor of community structure and directly influenced other environmental variables. Interestingly, the microbial composition was not linked to the local moss species identity. Instead, microbial communities varied significantly between sites due to differences in acidity levels and nitrogen availability. This finding was also mirrored in a co-occurrence network analysis, which displayed a distinct distribution of indicator species for acidity and nitrogen availability. Therefore, peatland conservation should take into account the critical habitat characteristics of moss-associated microbial communities.

DNA metabarcoding reveals the impact of Cu2+ on soil cercozoan diversity

Although copper (Cu2+) is a micronutrient, the metal may be toxic if present in high concentrations in soil ecosystems and subsequently affect various organisms, ranging from microorganisms to earthworms. We performed a microcosm study with an array of Cu2+ concentrations, with a specific focus on Cercozoa, an important protozoan group in most soil food webs. Research on Cercozoa is still scarce in terms of both diversity and ecology; hence, to explore this group in more depth, we used high-throughput sequencing to detect Cu2+ induced community changes. Increased levels of Cu2+ caused a shift in the cercozoan community, and we observed decreased cercozoan relative abundance across the majority of orders, families and genera. Due to their key role in soil food webs, especially as bacterial predators and providers of nutrients to plants, the reduction of cercozoan abundance and diversity may seriously affect soil functionality. Our results indicate that the increase of Cu2+ concentrations in the soil could potentially have this effect and the consequences need exploration.

eDNA metabarcoding reveals biodiversity and depth stratification patterns of dinoflagellate assemblages within the epipelagic zone of the western Coral Sea

BACKGROUND

Dinoflagellates play critical roles in the functioning of marine ecosystems but also may pose a hazard to human and ecosystem health by causing harmful algal blooms (HABs). The Coral Sea is a biodiversity hotspot, but its dinoflagellate assemblages in pelagic waters have not been studied by modern sequencing methods. We used metabarcoding of the 18 S rRNA V4 amplicon to assess the diversity and structure of dinoflagellate assemblages throughout the water column to a depth of 150 m at three stations in the Western Coral Sea. Additionally, at one station we compared metabarcoding with morphological methods to optimise identification and detection of dinoflagellates.

RESULTS

Stratification of dinoflagellate assemblages was evident in depth-specific relative abundances of taxonomic groups; the greatest difference was between the 5-30 m assemblages and the 130-150 m assemblages. The relative abundance of Dinophyceae (photosynthetic and heterotrophic) decreased with increasing depth, whereas that of Syndiniales (parasitic) increased with increasing depth. The composition of major taxonomic groups was similar among stations. Taxonomic richness and diversity of amplicon sequence variants (ASVs) were similar among depths and stations; however, the abundance of dominant taxa was highest within 0-30 m, and the abundance of rare taxa was highest within 130-150 m, indicating adaptations to specific depth strata. The number of unclassified ASVs at the family and species levels was very high, particularly for Syndinian representatives.

CONCLUSIONS

Dinoflagellate assemblages in open water of the Coral Sea are highly diverse and taxonomically stratified by depth; patterns of relative abundance along the depth gradient reflect environmental factors and ecological processes. Metabarcoding detects more species richness than does traditional microscopical methods of sample analysis, yet the methods are complementary, with morphological analysis revealing additional richness. The large number of unclassified dinoflagellate-ASVs indicates a need for improved taxonomic reference databases and suggests presence of dinoflagellate-crypto and-morphospecies.

The second most abundant dinophyte in the ponds of a botanical garden is a species new to science

In the microscopy realm, a large body of dark biodiversity still awaits to be uncovered. Unarmoured dinophytes are particularly neglected here, as they only present inconspicuous traits. In a remote German locality, we collected cells, from which a monoclonal strain was established, to study morphology using light and electron microscopy and to gain DNA sequences from the rRNA operon. In parallel, we detected unicellular eukaryotes in ponds of the Botanical Garden Munich-Nymphenburg by DNA-metabarcoding (V4 region of the 18S rRNA gene), weekly sampled over the course of a year. Strain GeoK*077 turned out to be a new species of Borghiella with a distinct position in molecular phylogenetics and characteristic coccoid cells of ovoid shape as the most important diagnostic trait. Borghiella ovum, sp. nov., was also present in artificial ponds of the Botanical Garden and was the second most abundant dinophyte detected in the samples. More specifically, Borghiella ovum, sp. nov., shows a clear seasonality, with high frequency during winter months and complete absence during summer months. The study underlines the necessity to assess the biodiversity, particularly of the microscopy realm more ambitiously, if even common species such as formerly Borghiella ovum are yet unknown to science.

How good are global DNA-based environmental surveys for detecting all protist diversity? Arcellinida as an example of biased representation

Metabarcoding approaches targeting microeukaryotes have deeply changed our vision of protist environmental diversity. The public repository EukBank consists of 18S v4 metabarcodes from 12,672 samples worldwide. To estimate how far this database provides a reasonable overview of all eukaryotic diversity, we used Arcellinida (lobose testate amoebae) as a case study. We hypothesised that (1) this approach would allow the discovery of unexpected diversity, but also that (2) some groups would be underrepresented because of primer/sequencing biases. Most of the Arcellinida sequences appeared in freshwater and soil, but their abundance and diversity appeared underrepresented. Moreover, 84% of ASVs belonged to the suborder Phryganellina, a supposedly species-poor clade, whereas the best-documented suborder (Glutinoconcha, 600 described species) was only marginally represented. We explored some possible causes of these biases. Mismatches in the primer-binding site seem to play a minor role. Excessive length of the target region could explain some of these biases, but not all. There must be some other unknown factors involved. Altogether, while metabarcoding based on ribosomal genes remains a good first approach to document microbial eukaryotic clades, alternative approaches based on other genes or sequencing techniques must be considered for an unbiased picture of the diversity of some groups.

MuDoGeR: Multi-Domain Genome recovery from metagenomes made easy

Several computational frameworks and workflows that recover genomes from prokaryotes, eukaryotes and viruses from metagenomes exist. Yet, it is difficult for scientists with little bioinformatics experience to evaluate quality, annotate genes, dereplicate, assign taxonomy and calculate relative abundance and coverage of genomes belonging to different domains. MuDoGeR is a user-friendly tool tailored for those familiar with Unix command-line environment that makes it easy to recover genomes of prokaryotes, eukaryotes and viruses from metagenomes, either alone or in combination. We tested MuDoGeR using 24 individual-isolated genomes and 574 metagenomes, demonstrating the applicability for a few samples and high throughput. While MuDoGeR can recover eukaryotic viral sequences, its characterization is predominantly skewed towards bacterial and archaeal viruses, reflecting the field's current state. However, acting as a dynamic wrapper, the MuDoGeR is designed to constantly incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field. MuDoGeR is open-source software available at https://github.com/mdsufz/MuDoGeR. Additionally, MuDoGeR is also available as a Singularity container.

Subtropical coastal microbiome variations due to massive river runoff after a cyclonic event

BACKGROUND

Coastal ecosystem variability at tropical latitudes is dependent on climatic conditions. During the wet, rainy season, extreme climatic events such as cyclones, precipitation, and winds can be intense over a short period and may have a significant impact on the entire land‒sea continuum. This study focused on the effect of river runoff across the southwest coral lagoon ecosystem of Grand Terre Island of New Caledonia (South Pacific) after a cyclonic event, which is considered a pulse disturbance at our study site. The variability of coastal microbiomes, studied by the metabarcoding of V4 18S (protists) and V4-V5 16S (bacteria) rDNA genes, after the cyclone passage was associated with key environmental parameters describing the runoff impact (salinity, organic matter proxies, terrestrial rock origin metals) and compared to community structures observed during the dry season.

RESULTS

Microbiome biodiversity patterns of the dry season were destructured because of the runoff impact, and land-origin taxa were observed in the coastal areas. After the rainy event, different daily community dynamics were observed locally, with specific microbial taxa explaining these variabilities. Plume dispersal modeling revealed the extent of low salinity areas up to the coral reef area (16 km offshore), but a rapid (< 6 days) recovery to typical steady conditions of the lagoon's hydrology was observed. Conversely, during the same time, some biological components (microbial communities, Chl a) and biogeochemical components (particulate nickel, terrigenous organic matter) of the ecosystem did not recover to values observed during the dry season conditions.

CONCLUSION

The ecosystem resilience of subtropical ecosystems must be evaluated from a multidisciplinary, holistic perspective and over the long term. This allows evaluating the risk associated with a potential continued and long-term disequilibrium of the ecosystem, triggered by the change in the frequency and intensity of extreme climatic events in the era of planetary climatic changes.

Mitochondrial genomes revisited: why do different lineages retain different genes?

The mitochondria contain their own genome derived from an alphaproteobacterial endosymbiont. From thousands of protein-coding genes originally encoded by their ancestor, only between 1 and about 70 are encoded on extant mitochondrial genomes (mitogenomes). Thanks to a dramatically increasing number of sequenced and annotated mitogenomes a coherent picture of why some genes were lost, or relocated to the nucleus, is emerging. In this review, we describe the characteristics of mitochondria-to-nucleus gene transfer and the resulting varied content of mitogenomes across eukaryotes. We introduce a 'burst-upon-drift' model to best explain nuclear-mitochondrial population genetics with flares of transfer due to genetic drift.

Meteora sporadica, a protist with incredible cell architecture, is related to Hemimastigophora

"Kingdom-level" branches are being added to the tree of eukaryotes at a rate approaching one per year, with no signs of slowing down.1,2,3,4 Some are completely new discoveries, whereas others are morphologically unusual protists that were previously described but lacked molecular data. For example, Hemimastigophora are predatory protists with two rows of flagella that were known since the 19th century but proved to represent a new deep-branching eukaryote lineage when phylogenomic analyses were conducted.2Meteora sporadica5 is a protist with a unique morphology; cells glide over substrates along a long axis of anterior and posterior projections while a pair of lateral "arms" swing back and forth, a motility system without any obvious parallels. Originally, Meteora was described by light microscopy only, from a short-term enrichment of deep-sea sediment. A small subunit ribosomal RNA (SSU rRNA) sequence was reported recently, but the phylogenetic placement of Meteora remained unresolved.6 Here, we investigated two cultivated Meteora sporadica isolates in detail. Transmission electron microscopy showed that both the anterior-posterior projections and the arms are supported by microtubules originating from a cluster of subnuclear microtubule organizing centers (MTOCs). Neither have a flagellar axoneme-like structure. Sequencing the mitochondrial genome showed this to be among the most gene-rich known, outside jakobids. Remarkably, phylogenomic analyses of 254 nuclear protein-coding genes robustly support a close relationship with Hemimastigophora. Our study suggests that Meteora and Hemimastigophora together represent a morphologically diverse "supergroup" and thus are important for resolving the tree of eukaryote life and early eukaryote evolution.

Evolution: No end in sight for novel incredible (heterotrophic) protists

The microbial eukaryotes known as protists are of immense importance for our understanding of eukaryotic biology. Although it is often difficult to convince funding bodies to sponsor research projects aimed at finding new protist lineages, such discoveries usually provide new and fundamental insights into cell and evolutionary biology, and ecology.

Genetic variation of the freshwater snail Indoplanorbis exustus (Gastropoda: Planorbidae) in Thailand, inferred from 18S and 28S rDNA sequences

Indoplanorbis exustus, a freshwater pulmonate snail, is widely distributed in tropical and subtropical zones and plays a significant role as an intermediate host for trematode parasites. Various genetic markers have been used for species identification and phylogenetic studies of this snail. However, there are limited studies about their molecular genetics based on nuclear ribosomal DNA (rDNA) genes. A genetic analysis of I. exustus in Thailand was conducted based on the nuclear 18S rDNA (339 bp) and 28S rDNA (1036 bp) genes. Indoplanorbis snails were collected from 29 localities in 21 provinces covering six regions of Thailand. Nucleotide sequences from 44 snails together with sequences from the GenBank database were examined for phylogenetic relationships and genetic diversity. All sequences of the selected nucleotide regions exhibited a high level of similarity (99%) to the sequences of I. exustus in the GenBank database. The maximum likelihood tree based on the 18S and 28S rDNA fragment sequences of I. exustus in Thailand revealed only one group with clear separation from another genus in the family Planorbidae. The I. exustus 28S rDNA sequences showed intraspecific genetic divergence ranging from 0 to 0.78% and were classified into 8 different haplotypes. Conversely, the 18S rDNA data showed lower variation than the 28S rDNA data and revealed a single haplotype and intraspecific distances of zero among all sampled individuals. The haplotype network of 28S rDNA sequences of I. exustus in Thailand revealed six unique haplotypes and two haplotypes shared by at least two regions. Overall, both markers were successful in the identification of I. exustus. However, these markers, particularly the 18S rDNA, may not be suitable for genetic analysis within the species, particularly for population genetic studies, due to their limited variation as seen in this study. In summary, this study not only enhances understanding of genetic variation in I. exustus but is also useful for the selection of molecular markers in future genetic research.

DNA Sequencing Technologies and DNA Barcoding

DNA barcodes are short, standardized DNA segments that geneticists can use to identify all living taxa. On the other hand, DNA barcoding identifies species by analyzing these specific regions against a DNA barcode reference library. In its initial years, DNA barcodes sequenced by Sanger's method were extensively used by taxonomists for the characterization and identification of species. But in recent years, DNA barcoding by next-generation sequencing (NGS) has found broader applications, such as quality control, biomonitoring of protected species, and biodiversity assessment. Technological advancements have also paved the way to metabarcoding, which has enabled massive parallel sequ.encing of complex bulk samples using high-throughput sequencing techniques. In future, DNA barcoding along with high-throughput techniques will show stupendous progress in taxonomic classification with reference to available sequence data.

Metagenomics and metatranscriptomics as potential driving forces for the exploration of diversity and functions of micro-eukaryotes in soil

Micro-eukaryotes are ubiquitous and play vital roles in diverse ecological systems, yet their diversity and functions are scarcely known. This may be due to the limitations of formerly used conventional culture-based methods. Metagenomics and metatranscriptomics are enabling to unravel the genomic, metabolic, and phylogenetic diversity of micro-eukaryotes inhabiting in different ecosystems in a more comprehensive manner. The in-depth study of structural and functional characteristics of micro-eukaryote community residing in soil is crucial for the complete understanding of this major ecosystem. This review provides a deep insight into the methodologies employed under these approaches to study soil micro-eukaryotic organisms. Furthermore, the review describes available computational tools, pipelines, and database sources and their manipulation for the analysis of sequence data of micro-eukaryotic origin. The challenges and limitations of these approaches are also discussed in detail. In addition, this review summarizes the key findings of metagenomic and metatranscriptomic studies on soil micro-eukaryotes. It also highlights the exploitation of these methods to study the structural as well as functional profiles of soil micro-eukaryotic community and to screen functional eukaryotic protein coding genes for biotechnological applications along with the future perspectives in the field.

Comparing Microsporidia-targeting primers for environmental DNA sequencing

Metabarcoding is a powerful tool to detect classical, and well-known "long-branch" Microsporidia in environmental samples. Several primer pairs were developed to target these unique microbial parasites, the majority of which remain undetected when using general metabarcoding primers. Most of these Microsporidia-targeting primer pairs amplify fragments of different length of the small subunit ribosomal RNA (SSU-rRNA) gene. However, we lack a broad comparison of the efficacy of those primers. Here, we conducted in silico PCRs with three short-read (which amplify a few-hundred base pairs) and two long-read (which amplify over a thousand base pairs) metabarcoding primer pairs on a variety of publicly available Microsporidia sensu lato SSU-rRNA gene sequences to test which primers capture most of the Microsporidia diversity. Our results indicate that the primer pairs do result in slight differences in inferred richness. Furthermore, some of the reverse primers are also able to bind to microsporidian subtaxa beyond the classical Microsporidia, which include the metchnikovellidan Amphiamblys spp., the chytridiopsid Chytridiopsis typographi and the "short-branch" microsporidian Mitosporidium daphniae.

Soil community composition in dynamic stages of semi-natural calcareous grassland

European dry thin-soil calcareous grasslands (alvars) are species-rich semi-natural habitats. Cessation of traditional management, such as mowing and grazing, leads to shrub and tree encroachment and the local extinction of characteristic alvar species. While soil microbes are known to play a critical role in driving vegetation and ecosystem dynamics, more information is needed about their composition and function in grasslands of different dynamic stages. Here we assess the composition of soil fungal, prokaryotic, and plant communities using soil environmental DNA from restored alvar grasslands in Estonia. The study areas included grasslands that had experienced different degrees of woody encroachment prior to restoration (woody plant removal and grazing), as well as unmanaged open grasslands. We found that, in general, different taxonomic groups exhibited correlated patterns of between-community variation. Previous forest sites, which had prior to restoration experienced a high degree of woody encroachment by ectomycorrhizal Scots pine, were compositionally most distinct from managed open grasslands, which had little woody vegetation even prior to restoration. The functional structure of plant and fungal communities varied in ways that were consistent with the representation of mycorrhizal types in the ecosystems prior to restoration. Compositional differences between managed and unmanaged open grasslands reflecting the implementation of grazing without further management interventions were clearer among fungal, and to an extent prokaryotic, communities than among plant communities. While previous studies have shown that during woody encroachment of alvar grassland, plant communities change first and fungal communities follow, our DNA-based results suggest that microbial communities reacted faster than plant communities during the restoration of grazing management in alvar grassland. We conclude that while the plant community responds faster to cessation of management, the fungal community responds faster to restoration of management. This may indicate hysteresis, where the eventual pathway back to the original state (grazed ecosystem) differs from the pathway taken towards the alternative state (abandoned semi-natural grassland ecosystem).

Protist taxonomic and functional diversity in aquatic ecosystems of the Brazilian Atlantic Forest

The Brazilian Atlantic Forest and its associated ecosystems are highly biodiverse but still understudied, especially with respect to eukaryotic microbes. Protists represent the largest proportion of eukaryotic diversity and play important roles in nutrient cycling and maintenance of the ecosystems in which they occur. However, much of protist diversity remains unknown, particularly in the Neotropics. Understanding the taxonomic and functional diversity of these organisms is urgently needed, not only to fill this gap in our knowledge, but also to enable the development of public policies for biological conservation. This is the first study to investigate the taxonomic and trophic diversity of the major protist groups in freshwater systems and brackish coastal lagoons located in fragments of the Brazilian Atlantic Forest by DNA metabarcoding, using high-throughput sequencing of the gene coding for the V4 region of the 18S rRNA gene. We compared α and β diversity for all protist communities and assessed the relative abundance of phototrophic, consumer, and parasitic taxa. We found that the protist communities of coastal lagoons are as diverse as the freshwater systems studied in terms of α diversity, although differed significantly in terms of taxonomic composition. Our results still showed a notable functional homogeneity between the trophic groups in freshwater environments. Beta diversity was higher among freshwater samples, suggesting a greater level of heterogeneity within this group of samples concerning the composition and abundance of OTUs.Ciliophora was the most represented group in freshwater, while Diatomea dominated diversity in coastal lagoons.

Genetic Markers for Metabarcoding of Freshwater Microalgae: Review

The metabarcoding methods for studying the diversity of freshwater microalgae and routine biomonitoring are actively used in modern research. A lot of experience has been accumulated already, and many methodological questions have been solved (such as the influence of the methods and time of sample conservation, DNA extraction and bioinformatical processing). The reproducibility of the method has been tested and confirmed. However, one of the main problems-choosing a genetic marker for the study-still lacks a clear answer. We analyzed 70 publications and found out that studies on eukaryotic freshwater microalgae use 12 markers (different nuclear regions 18S and ITS and plastids rbcL, 23S and 16S). Each marker has its peculiarities; they amplify differently and have various levels of efficiency (variability) in different groups of algae. The V4 and V9 18S and rbcL regions are used most often. We concentrated especially on the studies that compare the results of using different markers and microscopy. We summarize the data on the primers for each region and on how the choice of a marker affects the taxonomic composition of a community.

Community dynamics and co-occurrence relationships of pelagic ciliates and their potential prey at a coastal and an offshore station in the ultra-oligotrophic Eastern Mediterranean Sea

Ciliates have been recognized as one of the major components of the microbial food web, especially in ultra-oligotrophic waters, such as the Eastern Mediterranean Sea, where nutrients are scarce and the microbial community is dominated by pico- and nano-sized organisms. For this reason, ciliates play an important role in these ecosystems since they are the main planktonic grazers. Regardless the importance of these organisms, little is known about the community structure of heterotrophic and mixotrophic ciliates and how they are associated to their potential prey. In this study, we used 18S V4 rRNA gene metabarcoding to analyze ciliate community dynamics and how the relationship with potential prey changes according to different seasons and depths. Samples were collected seasonally at two stations of the Eastern Mediterranean Sea (HCB: coastal, M3A: offshore) from the surface and deep chlorophyll maximum (DCM) layers. The ciliate community structure varied across depths in HCB and across seasons in M3A, and the network analysis showed that in both stations, mixotrophic oligotrichs were positively associated with diatoms and showed few negative associations with ASVs annotated as marine Stramenopiles (MAST). On the other hand, heterotrophic tintinnids showed negative relationships in both HCB and M3A stations, mostly with Ochrophyta and Chlorophyta. These results showed, in first place that, although the two stations are close to each other, the ciliate dynamics differed between them. Moreover, mixotrophic and heterotrophic ciliates may have different ecological niches since mixotrophic ciliates may be more selective compared to heterotrophic species regarding their prey. These findings are the first glimpse into an understanding of the dynamics between heterotrophic and mixotrophic ciliates and their role in microbial assemblages and dynamics of ultra-oligotrophic environments.

Looking outside the box: a comparative cross-kingdom view on the cell biology of the three major lineages of eukaryotic multicellular life

Many cell biological facts that can be found in dedicated scientific textbooks are based on findings originally made in humans and/or other mammals, including respective tissue culture systems. They are often presented as if they were universally valid, neglecting that many aspects differ-in part considerably-between the three major kingdoms of multicellular eukaryotic life, comprising animals, plants and fungi. Here, we provide a comparative cross-kingdom view on the basic cell biology across these lineages, highlighting in particular essential differences in cellular structures and processes between phyla. We focus on key dissimilarities in cellular organization, e.g. regarding cell size and shape, the composition of the extracellular matrix, the types of cell-cell junctions, the presence of specific membrane-bound organelles and the organization of the cytoskeleton. We further highlight essential disparities in important cellular processes such as signal transduction, intracellular transport, cell cycle regulation, apoptosis and cytokinesis. Our comprehensive cross-kingdom comparison emphasizes overlaps but also marked differences between the major lineages of the three kingdoms and, thus, adds to a more holistic view of multicellular eukaryotic cell biology.

Extensive Cryptic Diversity and Ecological Associations Uncovered among Mexican and Global Collections of Naegleria and Vermamoeba Species by 18S Ribosomal DNA, Internal Transcribed Spacer, and Cytochrome Oxidase Subunit I Sequence Analysis

Free-living amoebae (FLA) are phagocytic protists that play crucial roles in microbial communities as significant microbial grazers. However, our current knowledge of their diversity, ecology, and population genetic structures is marginal due to the shallow and biased sampling of ecosystems and the use of few, poorly resolving molecular markers. Thirty-two FLA were isolated from soil and water samples collected across representative ecosystems of the State of Morelos in Central Mexico, including the drinking water distribution system (DWDS) from the state capital. We classified our isolates as members of Acanthamoeba, Vermamoeba, Naegleria, and Tetramitus by 18S ribosomal DNA (rDNA) sequencing. Vermamoeba isolates were recovered exclusively from the DWDS samples. In contrast, Naegleria strains displayed a broad distribution in soil and water samples across the natural ecosystems. We used a combination of phylogenetic and population genetic analyses of internal transcribed spacer (ITS) and cytochrome oxidase subunit I (COI) sequences from our isolates and a comprehensive set of reference sequences to analyze the currently known diversity of Naegleria spp. Significant associations were uncovered between the most prevalent lineages of Naegleria and Vermamoeba and broad ecological and geographical variables at regional and global levels. The population structure and cryptic diversity within the Naegleria galeacystis-Naegleria americana and Vermamoeba vermiformis species complexes were thoroughly analyzed. Our results prove that the genus Vermamoeba, which was previously thought to consist of only one species, actually encompasses at least seven widely distributed species, as indicated by consistent evidence from Bayesian phylogenetics, two species-delimitation programs, and population genetics analyses. IMPORTANCE Our study sheds new light on the population genetic structure of V. vermiformis and diverse Naegleria species. Using improved molecular markers and advanced analytical approaches, we discovered that N. americana, previously considered a single species, actually contains multiple distinct lineages, as revealed by COI sequencing. These lineages are highly differentiated, with little gene flow between them. Our findings demonstrate that the genus Vermamoeba holds multiple cryptic species, requiring a significant taxonomic revision in light of multilocus sequence analyses. These results advance our understanding of the ecology, molecular systematics, and biogeography of these genera and species complexes at both regional and global scales. This study has significant implications for diagnosing amoebal infections and evaluating health risks associated with FLA in domestic and recreational waters.

Beyond the limits of the unassigned protist microbiome: inferring large-scale spatio-temporal patterns of Syndiniales marine parasites

Marine protists are major components of the oceanic microbiome that remain largely unrepresented in culture collections and genomic reference databases. The exploration of this uncharted protist diversity in oceanic communities relies essentially on studying genetic markers from the environment as taxonomic barcodes. Here we report that across 6 large scale spatio-temporal planktonic surveys, half of the genetic barcodes remain taxonomically unassigned at the genus level, preventing a fine ecological understanding for numerous protist lineages. Among them, parasitic Syndiniales (Dinoflagellata) appear as the least described protist group. We have developed a computational workflow, integrating diverse 18S rDNA gene metabarcoding datasets, in order to infer large-scale ecological patterns at 100% similarity of the genetic marker, overcoming the limitation of taxonomic assignment. From a spatial perspective, we identified 2171 unassigned clusters, i.e., Syndiniales sequences with 100% similarity, exclusively shared between the Tropical/Subtropical Ocean and the Mediterranean Sea among all Syndiniales orders and 25 ubiquitous clusters shared within all the studied marine regions. From a temporal perspective, over 3 time-series, we highlighted 39 unassigned clusters that follow rhythmic patterns of recurrence and are the best indicators of parasite community's variation. These clusters withhold potential as ecosystem change indicators, mirroring their associated host community responses. Our results underline the importance of Syndiniales in structuring planktonic communities through space and time, raising questions regarding host-parasite association specificity and the trophic mode of persistent Syndiniales, while providing an innovative framework for prioritizing unassigned protist taxa for further description.

Recovery of 197 eukaryotic bins reveals major challenges for eukaryote genome reconstruction from terrestrial metagenomes

As most eukaryotic genomes are yet to be sequenced, the mechanisms underlying their contribution to different ecosystem processes remain untapped. Although approaches to recovering Prokaryotic genomes have become common in genome biology, few studies have tackled the recovery of eukaryotic genomes from metagenomes. This study assessed the reconstruction of microbial eukaryotic genomes using 6000 metagenomes from terrestrial and some transition environments using the EukRep pipeline. Only 215 metagenomic libraries yielded eukaryotic bins. From a total of 447 eukaryotic bins recovered 197 were classified at the phylum level. Streptophytes and fungi were the most represented clades with 83 and 73 bins, respectively. More than 78% of the obtained eukaryotic bins were recovered from samples whose biomes were classified as host-associated, aquatic, and anthropogenic terrestrial. However, only 93 bins were taxonomically assigned at the genus level and 17 bins at the species level. Completeness and contamination estimates were obtained for a total of 193 bins and consisted of 44.64% (σ = 27.41%) and 3.97% (σ = 6.53%), respectively. Micromonas commoda was the most frequent taxon found while Saccharomyces cerevisiae presented the highest completeness, probably because more reference genomes are available. Current measures of completeness are based on the presence of single-copy genes. However, mapping of the contigs from the recovered eukaryotic bins to the chromosomes of the reference genomes showed many gaps, suggesting that completeness measures should also include chromosome coverage. Recovering eukaryotic genomes will benefit significantly from long-read sequencing, development of tools for dealing with repeat-rich genomes, and improved reference genomes databases.

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.

Biogeographical distribution of Hartaetosiga strains (Choanoflagellatea, Craspedida, Salpingoecidae) including morphological and transcriptomic data from a transect across the Atlantic Ocean

The genus Hartaetosiga Carr, Richter and Nitsche, 2017 comprised up to now only three species, H. gracilis (Kent) Carr, Richter, Nitsche, 2017, H. balthica (Wylezich and Karpov) Carr, Richter and Nitsche, 2017 and H. minima (Wylezich and Karpov) Carr, Richter and Nitsche, 2017. Based on distinct molecular data these species were relocated from the strictly freshwater genus Codosiga (Ehrenberg) Bütschli, 1878 to a new genus comprising brackish and marine species. During the cruise MSM82/2 across the Atlantic Ocean in 2019, surface water samples were taken from 15 stations along a transect ranging from 35°S to 23°N. We were able to isolate and cultivate 14 strains of the genus Hartaetosiga. Morphometric data showed no distinct morphological traits allowing for a species delineation, indicating a cryptic species complex within the genus. Based on cultivation, morphological data, and molecular analyses, we recorded H. gracilis for the first time from off-shelf waters of the Atlantic Ocean and could describe a new species, H. australis n. sp. This new species was recorded from sampling stations in the Southern Hemisphere only, which may indicate a potential biogeographic restriction likely caused by the Equatorial Counter Current (ECC), dividing the northern and southern surface waters.

DNA barcoding, an effective tool for species identification: a review

DNA barcoding is a powerful taxonomic tool to identify and discover species. DNA barcoding utilizes one or more standardized short DNA regions for taxon identification. With the emergence of new sequencing techniques, such as Next-generation sequencing (NGS), ONT MinION nanopore sequencing, and Pac Bio sequencing, DNA barcoding has become more accurate, fast, and reliable. Rapid species identification by DNA barcodes has been used in a variety of fields, including forensic science, control of the food supply chain, and disease understanding. The Consortium for Barcode of Life (CBOL) presents various working groups to identify the universal barcode gene, such as COI in metazoans; rbcL, matK, and ITS in plants; ITS in fungi; 16S rRNA gene in bacteria and archaea, and creating a reference DNA barcode library. In this article, an attempt has been made to analyze the various proposed DNA barcode for different organisms, strengths & limitations, recent advancements in DNA barcoding, and methods to speed up the DNA barcode reference library construction. This study concludes that constructing a reference library with high species coverage would be a major step toward identifying species by DNA barcodes. This can be achieved in a short period of time by using advanced sequencing and data analysis methods.

Coinfection patterns of two marine apicomplexans are not associated with genetic diversity of their polychaete host

Coinfections of two or more parasites within one host are more of a rule than an exception in nature. Interactions between coinfecting parasites can greatly affect their abundance and prevalence. Characteristics of the host, such as genetic diversity, can also affect the infection dynamics of coinfecting parasites. Here, we investigate for the first time the association of coinfection patterns of two marine apicomplexans, Rhytidocystis sp. and Selenidium pygospionis, with the genetic diversity of their host, the polychaete Pygospio elegans, from natural populations. Host genetic diversity was determined with seven microsatellite loci and summarized as allelic richness, inbreeding coefficient, and individual heterozygosity. We detected nonsignificant correlations between infection loads and both individual host heterozygosity and population genetic diversity. Prevalence and infection load of Rhytidocystis sp. were higher than those of S. pygospionis, and both varied spatially. Coinfections were common, and almost all hosts infected by S. pygospionis were also infected by Rhytidocystis sp. Rhytidocystis sp. infection load was significantly higher in dual infections. Our results suggest that factors other than host genetic diversity might be more important in marine apicomplexan infection patterns and experimental approaches would be needed to further determine how interactions between the apicomplexans and their host influence infection.

Ciliate Morpho-Taxonomy and Practical Considerations before Deploying Metabarcoding to Ciliate Community Diversity Surveys in Urban Receiving Waters

Disentangling biodiversity and community assembly effects on ecosystem function has always been an important topic in ecological research. The development and application of a DNA metabarcoding method has fundamentally changed the way we describe prokaryotic communities and estimate biodiversity. Compared to prokaryotes (bacteria and archaea), the eukaryotic microbes (unicellular eukaryotes) also fulfill extremely important ecological functions in different ecosystems regarding their intermediate trophic positions. For instance, ciliated microbes (accounting for a substantial portion of the diversity of unicellular eukaryotes) perform pivotal roles in microbial loops and are essential components in different ecosystems, especially in water purification processes. Therefore, the community composition of ciliated species has been widely utilized as a proxy for water quality and biological assessment in urban river ecosystems and WWTPs (wastewater treatment plants). Unfortunately, investigating the dynamic changes and compositions in ciliate communities relies heavily on existing morpho-taxonomical descriptions, which is limited by traditional microscopic approaches. To deal with this dilemma, we discuss the DNA-based taxonomy of ciliates, the relative merits and challenges of deploying its application using DNA metabarcoding for surveys of ciliate community diversity in urban waterbodies, and provide suggestions for minimizing relevant sources of biases in its implementation. We expect that DNA metabarcoding could untangle relationships between community assembly and environmental changes affecting ciliate communities. These analyses and discussions could offer a replicable method in support of the application of evaluating communities of ciliated protozoa as indicators of urban freshwater ecosystems.

Long-Read Sequencing Analysis Revealed the Impact of Forest Conversion on Soil Fungal Diversity in Limu Mountain, Hainan

Soil fungi are essential to soil microorganisms that play an important role in the ecosystem's soil carbon cycle and mineral nutrient transformation. Understanding the structural characteristics and diversity of soil fungal communities helps understand the health of forest ecosystems. The transition from tropical rainforest to artificial forest greatly impacts the composition and diversity of fungal communities. Hainan Limushan tropical rainforest National Park has a large area of artificial forests. Ecologists have conducted in-depth studies on the succession of animals and plants to regenerate tropical rainforests. There are few reports on the diversity of soil fungi and its influencing factors in the succession of tropical rainforests in Limu Mountain. In this study, 44 soil samples from five different stands were collected in the tropical rainforest of Limushan, Hainan. High-throughput sequencing of rDNA in its region was used to analyze fungal communities and study their α and β diversity. Analysis of variance and multiple regression models was used to analyze soil variables and fungal functional groups to determine the effects of interaction between fungi and environmental factors. A total of 273,996 reads and 1290 operational taxonomic units (OTUs) were obtained, belonging to 418 species, 325 genera, 159 families, eight phyla, 30 classes, and 73 orders. The results showed that the composition of soil fungal communities in the five stands was similar, with ascomycetes accounting for 70.5% and basidiomycetes accounting for 14.7%. α and β diversity analysis showed that soil fungi in Limushan tropical rainforest had high abundance and diversity. Multiple regression analysis between soil variables and functional groups showed that organic matter, TN, TP, TK, and AK were excellent predictors for soil fungi. TP was the strongest predictor in all functional groups except soil saprotroph. Organic matter and total nitrogen were the strongest predictors of soil rot. The transformation from tropical rainforest to artificial forest in Limushan did not change the soil fungal community structure, but the richness and diversity of soil fungi changed. The forest transformation did not lead to decreased soil fungal abundance and diversity. Different vegetation types and soil properties affect the diversity of soil fungal communities. We found that Caribbean pine plantations can improve soil fungal diversity, while long-term Eucalyptus spp. plantations may reduce soil fungal 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 PR² reference sequence database. This database, called metaPR² , 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.

Discrepancies between prokaryotes and eukaryotes need to be considered in soil DNA-based studies

Metabarcoding approaches are exponentially increasing our understanding of soil biodiversity, with a major focus on the bacterial part of the microbiome. Part of the soil diversity are also eukaryotes that include fungi, algae, protists and Metazoa. Nowadays, soil eukaryotes are targeted with the same approaches developed for bacteria and archaea (prokaryotes). However, fundamental differences exist between domains. After providing a short historical overview of the developments of metabarcoding applied to environmental microbiology, we compile the most important differences between domains that prevent direct method transfers between prokaryotic and eukaryotic soil metabarcoding approaches, currently dominated by short-read sequencing. These include the existence of divergent diversity concepts and the variations in eukaryotic morphology that affect sampling and DNA extraction. Furthermore, eukaryotes experienced much more variable evolutionary rates than prokaryotes, which prevent capturing the entire eukaryotic diversity in a soil with a single amplification protocol fit for short-read sequencing. In the final part we focus on future potentials for optimization of eukaryotic metabarcoding that include superior possibility of functionally characterizing eukaryotes and to extend the current information obtained, such as by adding a real quantitative component. This review should optimize future metabarcoding approaches targeting soil eukaryotes and kickstart this promising research direction.

Metagenomics versus total RNA sequencing: most accurate data-processing tools, microbial identification accuracy and perspectives for ecological assessments

Metagenomics and total RNA sequencing (total RNA-Seq) have the potential to improve the taxonomic identification of diverse microbial communities, which could allow for the incorporation of microbes into routine ecological assessments. However, these target-PCR-free techniques require more testing and optimization. In this study, we processed metagenomics and total RNA-Seq data from a commercially available microbial mock community using 672 data-processing workflows, identified the most accurate data-processing tools, and compared their microbial identification accuracy at equal and increasing sequencing depths. The accuracy of data-processing tools substantially varied among replicates. Total RNA-Seq was more accurate than metagenomics at equal sequencing depths and even at sequencing depths almost one order of magnitude lower than those of metagenomics. We show that while data-processing tools require further exploration, total RNA-Seq might be a favorable alternative to metagenomics for target-PCR-free taxonomic identifications of microbial communities and might enable a substantial reduction in sequencing costs while maintaining accuracy. This could be particularly an advantage for routine ecological assessments, which require cost-effective yet accurate methods, and might allow for the incorporation of microbes into ecological assessments.

The Toxoplasma glucan phosphatase TgLaforin utilizes a distinct functional mechanism that can be exploited by therapeutic inhibitors

Toxoplasma gondii is an intracellular parasite that generates amylopectin granules (AGs), a polysaccharide associated with bradyzoites that define chronic T. gondii infection. AGs are postulated to act as an essential energy storage molecule that enable bradyzoite persistence, transmission, and reactivation. Importantly, reactivation can result in the life-threatening symptoms of toxoplasmosis. T. gondii encodes glucan dikinase and glucan phosphatase enzymes that are homologous to the plant and animal enzymes involved in reversible glucan phosphorylation and which are required for efficient polysaccharide degradation and utilization. However, the structural determinants that regulate reversible glucan phosphorylation in T. gondii are unclear. Herein, we define key functional aspects of the T. gondii glucan phosphatase TgLaforin (TGME49_205290). We demonstrate that TgLaforin possesses an atypical split carbohydrate-binding-module domain. AlphaFold2 modeling combined with hydrogen-deuterium exchange mass spectrometry and differential scanning fluorimetry also demonstrate the unique structural dynamics of TgLaforin with regard to glucan binding. Moreover, we show that TgLaforin forms a dual specificity phosphatase domain-mediated dimer. Finally, the distinct properties of the glucan phosphatase catalytic domain were exploited to identify a small molecule inhibitor of TgLaforin catalytic activity. Together, these studies define a distinct mechanism of TgLaforin activity, opening up a new avenue of T. gondii bradyzoite biology as a therapeutic target.

Study on Analysis of Several Molecular Identification Methods for Ciliates of Colpodea (Protista, Ciliophora)

The application of molecular techniques to accurately identify protozoan species can correct previous misidentifications based on traditional morphological identification. Colpodea ciliates have many toxicological and cytological applications, but their subtle morphological differences and small body size hinder species delineation. Herein, we used Cox I and β-tubulin genes, alongside fluorescence in situ hybridization (FISH), to evaluate each method in delineating Colpodea species. For this analysis, Colpoda harbinensis n. sp., C. reniformis, two populations of C. inflata, Colpoda compare grandis, and five populations of Paracolpoda steinii, from the soil in northeastern China, were used. We determined that (1) the Cox I gene was more suitable than the β-tubulin gene as a molecular marker for defining intra- and interspecific level relationships of Colpoda. (2) FISH probes designed for Colpoda sp., C. inflata, Colpoda compare grandis, and Paracolpoda steinii, provided rapid interspecific differentiation of Colpodea species. (3) Colpoda harbinensis n. sp. was established and mainly characterized by its size in vivo (approximately $80\times 60\mu \text{m}$ ), a reniform body in outline, one macronucleus, its spherical shape, a sometimes nonexistent micronucleus, 11–15 somatic kineties, and five or six postoral kineties. In conclusion, combining oligonucleotide probes, DNA barcoding, and morphology for the first time, we have greatly improved the delineation of Colpodea and confirmed that Cox I gene was a promising DNA barcoding marker for species of Colpodea, and FISH could provide useful morphological information as complementing traditional techniques such as silver carbonate.

Microbial Communities in Biocrusts Are Recruited From the Neighboring Sand at Coastal Dunes Along the Baltic Sea

Biological soil crusts occur worldwide as pioneer communities stabilizing the soil surface. In coastal primary sand dunes, vascular plants cannot sustain due to scarce nutrients and the low-water-holding capacity of the sand sediment. Thus, besides planted dune grass, biocrusts are the only vegetation there. Although biocrusts can reach high coverage rates in coastal sand dunes, studies about their biodiversity are rare. Here, we present a comprehensive overview of the biodiversity of microorganisms in such biocrusts and the neighboring sand from sampling sites along the Baltic Sea coast. The biodiversity of Bacteria, Cyanobacteria, Fungi, and other microbial Eukaryota were assessed using high-throughput sequencing (HTS) with a mixture of universal and group-specific primers. The results showed that the biocrusts recruit their microorganisms mainly from the neighboring sand rather than supporting a universal biocrust microbiome. Although in biocrusts the taxa richness was lower than in sand, five times more co-occurrences were identified using network analysis. This study showed that by comparing neighboring bare surface substrates with biocrusts holds the potential to better understand biocrust development. In addition, the target sequencing approach helps outline potential biotic interactions between different microorganisms groups and identify key players during biocrust development.

DNA barcoding of fungal specimens using PacBio long-read high-throughput sequencing

Molecular methods are increasingly used to identify species that lack conspicuous macro- or micromorphological characters. Taxonomic and ecological research teams barcode large numbers of collected voucher specimens annually. In this study we assessed the efficiency of long-read high throughput sequencing (HTS) as opposed to the traditionally used Sanger method for taxonomic identification of multiple vouchered fungal specimens. We also evaluated whether this method can provide reference information about intra-individual gene polymorphism. We developed a workflow based on a test set of 423 basidiomycete specimens (representing 195 species), the PacBio HTS method, and ribosomal rRNA operon internal transcribed spacer (ITS) and 28S rRNA gene (LSU) markers. The PacBio HTS had a higher success rate than Sanger sequencing at a comparable cost. Species identification based on PacBio reads was usually straightforward, because the dominant operational taxonomic unit (OTU) typically represented the targeted organism. The PacBio HTS also enabled us to detect widespread polymorphism within the ITS marker. We conclude that multiplex DNA barcoding of the fungal ITS and LSU markers using PacBio HTS is a useful tool for taxonomic identification of large amounts of collected voucher specimens at a competitive price. Furthermore, PacBio HTS accurately recovers various alleles and paralogs, which can provide crucial information for species delimitation and population-level studies.

Intra-genomic rRNA gene variability of Nassellaria and Spumellaria (Rhizaria, Radiolaria) assessed by Sanger, MinION and Illumina sequencing

Ribosomal RNA (rRNA) genes are known to be valuable markers for the barcoding of eukaryotic life and its phylogenetic classification at various taxonomic levels. The large-scale exploration of environmental microbial diversity through metabarcoding approaches have been focused mainly on the V4 and V9 regions of the 18S rRNA gene. The accurate interpretation of such environmental surveys is hampered by technical (e.g., PCR and sequencing errors) and biological biases (e.g., intra-genomic variability). Here we explored the intra-genomic diversity of Nassellaria and Spumellaria specimens (Radiolaria) by comparing Sanger sequencing with Illumina and Oxford Nanopore Technologies (MinION). Our analysis determined that intra-genomic variability of Nassellaria and Spumellaria is generally low, yet some Spumellaria specimens showed two different copies of the V4 with <97% similarity. From the different sequencing methods, Illumina showed the highest number of contaminations (i.e., environmental DNA, cross-contamination, tag-jumping), revealed by its high sequencing depth; and MinION showed the highest sequencing rate error (~14%). Yet the long reads produced by MinION (~2900 bp) allowed accurate phylogenetic reconstruction studies. These results highlight the requirement for a careful interpretation of Illumina based metabarcoding studies, in particular regarding low abundant amplicons, and open future perspectives towards full-length rDNA environmental metabarcoding surveys.

Eukaryotic cellular intricacies shape mitochondrial proteomic complexity

Mitochondria have been fundamental to the eco-physiological success of eukaryotes since the last eukaryotic common ancestor (LECA). They contribute essential functions to eukaryotic cells, above and beyond classical respiration. Mitochondria interact with, and complement, metabolic pathways occurring in other organelles, notably diversifying the chloroplast metabolism of photosynthetic organisms. Here, we integrate existing literature to investigate how mitochondrial metabolism varies across the landscape of eukaryotic evolution. We illustrate the mitochondrial remodelling and proteomic changes undergone in conjunction with major evolutionary transitions. We explore how the mitochondrial complexity of the LECA has been remodelled in specific groups to support subsequent evolutionary transitions, such as the acquisition of chloroplasts in photosynthetic species and the emergence of multicellularity. We highlight the versatile and crucial roles played by mitochondria during eukaryotic evolution, extending from its huge contribution to the development of the LECA itself to the dynamic evolution of individual eukaryote groups, reflecting both their current ecologies and evolutionary histories.

Parasitic infections by Group II Syndiniales target selected dinoflagellate host populations within diverse protist assemblages in a model coastal pond

Protists are integral to marine food webs and biogeochemical cycles; however, there is a paucity of data describing specific ecological niches for some of the most abundant taxa in marker gene libraries. Syndiniales are one such group, often representing the majority of sequence reads recovered from picoplankton samples across the global ocean. However, the prevalence and impacts of syndinian parasitism in marine environments remains unclear. We began to address these critical knowledge gaps by generating a high-resolution time series (March-October 2018) in a productive coastal pond. Seasonal shifts in protist populations, including parasitic Syndiniales, were documented during periods of higher primary productivity and increased summer temperature-driven stratification. Elevated concentrations of infected hosts and free-living parasite spores occurred at nearly monthly intervals in July, August, and September. We suggest intensifying stratification during this period correlated with the increased prevalence of dinoflagellates that were parasitized by Group II Syndiniales. Infections in some protist populations were comparable to previously reported large single-taxon dinoflagellate blooms. Infection dynamics in Salt Pond demonstrated the propagation of syndinian parasites through mixed protist assemblages and highlighted patterns of host/parasite interactions that better reflect many other marine environments where single taxon blooms are uncommon.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

18S rRNA gene amplicon sequencing combined with culture-based surveys of maize rhizosphere protists reveal dominant, plant-enriched and culturable community members

Protists play important roles in shaping the microbial community of the rhizosphere and defining these roles will require the study of protist isolates. However, there is still a limited understanding of how well protist isolation efforts can capture the diversity and composition of rhizosphere protistan communities. Here, we report a simultaneous isolation and 18S rRNA gene amplicon sequencing survey describing the protist diversity of maize rhizospheres in two climatically and pedologically distinct sites. We demonstrated that the maize rhizosphere exerted significant and site-dependent effects on the protistan community structure and defined a set of core and rhizosphere-enriched protists. From the same root samples, we generated a library of 103 protist isolates representing 46 18S rRNA gene sequence variants from six eukaryotic supergroups. While cultured isolates represented a small proportion of total protist diversity recovered by sequencing, they included taxa enriched in rhizosphere soils across all samples, encompassing 9% of all core sequence variants. The isolation approach also captured 17 protists not detected through 18S rRNA gene amplicon sequencing. This study demonstrated that maize roots select for distinct protistan communities, and established a diverse protist culture collection that can be used for future research linking protists to rhizosphere status and plant health.

Unearthing soil arthropod diversity through DNA metabarcoding

DNA metabarcoding has the potential to greatly advance understanding of soil biodiversity, but this approach has seen limited application for the most abundant and species-rich group of soil fauna-the arthropods. This study begins to address this gap by comparing information on species composition recovered from metabarcoding two types of bulk samples (specimens, soil) from a temperate zone site and from bulk soil samples collected at eight sites in the Arctic. Analysis of 22 samples (3 specimen, 19 soil) revealed 410 arthropod OTUs belonging to 112 families, 25 orders, and nine classes. Studies at the temperate zone site revealed little overlap in species composition between soil and specimen samples, but more overlap at higher taxonomic levels (families, orders) and congruent patterns of α- and β-diversity. Expansion of soil analyses to the Arctic revealed locally rich, highly dissimilar, and spatially structured assemblages compatible with dispersal limited and environmentally driven assembly. The current study demonstrates that DNA metabarcoding of bulk soil enables rapid, large-scale assessments of soil arthropod diversity. However, deep sequence coverage is required to adequately capture the species present in these samples, and expansion of the DNA barcode reference library is necessary to improve taxonomic resolution of the sequences recovered through this approach.

Horizontal and vertical distribution of Gambierdiscus spp. (Dinophyceae) including novel phylotypes in Japan identified by 18S rDNA metabarcoding

The genus Gambierdiscus is a marine benthic/epiphytic dinoflagellate considered the causative agent of ciguatera poisoning (CP). Clarifying the geographical distribution of this genus to understand the potential risk of CP is important. Many studies have focused only on the species/phylotype composition of Gambierdiscus in shallow waters, but no study has investigated the species/phylotype composition of the genus in deep waters. In the present study, the distributions of Gambierdiscus species/phylotypes at two depths (2-8 and 30 m) and two sampling sites (temperate and subtropical) in Japan was investigated using high throughput sequencing (HTS) with a newly developed primer set that preferentially amplifies the 18S rDNA V8-V9 region of Alveolata. A phylogenetic analysis using 89 samples collected over three years revealed of ten Gambierdiscus species/phylotypes including not only two species that have not been reported in Japan (G. caribaeus and G. silvae) but also four novel phylotypes (Gambierdiscus spp. Clade II_1, Clade II_2, Clade II_3, and Clade VI_1). Uncorrected genetic distances also supported that these new phylotypes clearly diverged from other Gambierdiscus species. All four new phylotypes, G. caribaeus, and G. silvae were distributed in the subtropical region. Among them, Clade II_2, Clade VI_1, and G. silvae were also distributed in the temperate region. Four species/phylotypes previously reported from Japan showed a similar distribution as reported previously. Among the ten species/phylotypes, Gambierdiscus sp. type 3 and Clade VI_1 were found only in deep waters, whereas five species/phylotypes were observed only in shallow waters. The other three species/phylotypes were found in both deep and shallow waters. The results of the horizontal and vertical distribution suggest that the growth characteristics of each species/phylotypes found in Japan might adapt to the ambient environmental conditions. This study revealed an inclusive assemblage of Gambierdiscus species/phylotypes in Japan through metabarcoding using the Alveolata primer set. In the future, the abundance and toxicities/toxin productions of the newly reported species/phylotypes need to be clarified to understand the mechanism of CP outbreaks in Japan.

A high-throughput amplicon sequencing approach for population-wide species diversity and composition survey

Management of agricultural pests requires an understanding of pest species diversity, their interactions with beneficial insects and spatial-temporal patterns of pest abundance. Invasive and agriculturally important insect pests can build up very high populations, especially in cropping landscapes. Traditionally, sampling effort for species identification involves small sample sizes and is labour intensive. Here, we describe a multi-primer high throughput sequencing (HTS) metabarcoding method and associated analytical workflow for a rapid, intensive, high-volume survey of pest species compositions. We demonstrate our method using the taxonomically challenging Bemisia pest cryptic species complex as examples. The whiteflies Bemisia including the 'tabaci' species are agriculturally important capable of vectoring diverse plant viruses that cause diseases and crop losses. Our multi-primer metabarcoding HTS amplicon approach simultaneously process high volumes of whitefly individuals, with efficiency to detect rare (i.e., 1%) test-species, while our improved whitefly primers for metabarcoding also detected beneficial hymenopteran parasitoid species from whitefly nymphs. Field-testing our redesigned Bemisia metabarcoding primer sets across the Tanzania, Uganda and Malawi cassava cultivation landscapes, we identified the sub-Saharan Africa 1 Bemisia putative species as the dominant pest species, with other cryptic Bemisia species being detected at various abundances. We also provide evidence that Bemisia species compositions can be affected by host crops and sampling techniques that target either nymphs or adults. Our multi-primer HTS metabarcoding method incorporated two over-lapping amplicons of 472bp and 518bp that spanned the entire 657bp 3' barcoding region for Bemisia, and is particularly suitable to molecular diagnostic surveys of this highly cryptic insect pest species complex that also typically exhibited high population densities in heavy crop infestation episodes. Our approach can be adopted to understand species biodiversity across landscapes, with broad implications for improving trans-boundary biosecurity preparedness, thus contributing to molecular ecological knowledge and the development of control strategies for high-density, cryptic, pest-species complexes.

Tubulin post-translational modifications in protists - Tiny models for solving big questions

Protists are an exceptionally diverse group of mostly single-celled eukaryotes. The organization of the microtubular cytoskeleton in protists from various evolutionary lineages has different levels of sophistication, from a network of microtubules (MTs) supporting intracellular trafficking as in Dictyostelium, to complex structures such as basal bodies and cilia/flagella enabling cell motility, and lineage-specific adaptations such as the ventral disc in Giardia. MTs building these diverse structures have specific properties partly due to the presence of tubulin post-translational modifications (PTMs). Among them there are highly evolutionarily conserved PTMs: acetylation, detyrosination, (poly)glutamylation and (poly)glycylation. In some protists also less common tubulin PTMs were identified, including phosphorylation, methylation, Δ2-, Δ5- of α-tubulin, polyubiquitination, sumoylation, or S-palmitoylation. Not surprisingly, several single-celled organisms become models to study tubulin PTMs, including their effect on MT properties and discovery of the modifying enzymes. Here, we briefly summarize the current knowledge on tubulin PTMs in unicellular eukaryotes and highlight key findings in protists as model organisms.

Fine-scale genetic diversity and putative ecotypes of oxymonad protists coinhabiting the hindgut of Reticulitermes speratus

The hindgut of lower termites is generally coinhabited by multiple morphologically identifiable protist species. However, it is unclear how many protist species truly coexist in this miniaturized environment, and moreover, it is difficult to define the fundamental unit of protist diversity. Species delineation of termite gut protists has therefore been guided without a theory-based concept of species. Here, we focused on the hindgut of the termite Reticulitermes speratus, where 10 or 11 morphologically distinct oxymonad cell types, that is, morphospecies, coexist. We elucidated the phylogenetic structure of all co-occurring oxymonads and addressed whether their diversity can be explained by the "ecotype" hypothesis. Oxymonad-specific 18S rRNA gene amplicon sequencing analyses of whole-gut samples, combined with single-cell 18S rRNA sequencing of the oxymonad morphospecies, identified 210 one-nucleotide-level variants. The phylogenetic analysis of these variants revealed the presence of microdiverse clusters typically within 1% sequence divergence. Each known oxymonad morphospecies comprised one to several monophyletic or paraphyletic microdiverse clusters. Using these sequence data sets, we conducted computational simulation to predict the rates of ecotype formation and periodic selection, and to demarcate putative ecotypes. Our simulations suggested that the oxymonad genetic divergence is constrained primarily by strong selection, in spite of limited population size and possible bottlenecks during intergenerational transmission. A total of 33 oxymonad ecotypes were predicted, and most of the putative ecotypes were consistently detected among different colonies and host individuals. These findings provide a possible theoretical basis for species diversity and underlying mechanisms of coexistence of termite gut protists.

Droplet digital PCR as a tool for investigating dynamics of cryptic symbionts

Interactions among symbiotic organisms and their hosts are major drivers of ecological and evolutionary processes. Monitoring the infection patterns among natural populations and identifying factors affecting these interactions are critical for understanding symbiont-host relationships. However, many of these interactions remain understudied since the knowledge about the symbiont species is lacking, which hinders the development of appropriate tools. In this study, we developed a digital droplet PCR (ddPCR) assay based on apicomplexan COX1 gene to detect an undescribed agamococcidian symbiont. We show that the method gives precise and reproducible results and enables detecting cryptic symbionts in low target concentration. We further exemplify the assay's use to survey seasonally sampled natural host (Pygospio elegans) populations for symbiont infection dynamics. We found that symbiont prevalence differs spatially but does not show seasonal changes. Infection load differed between populations and was low in spring and significantly increased towards fall in all populations. We also found that the symbiont prevalence is affected by host length and population density. Larger hosts were more likely to be infected, and high host densities were found to have a lower probability of infection. The observed variations could be due to characteristics of both symbiont and host biology, especially the seasonal variation in encounter rates. Our findings show that the developed ddPCR assay is a robust tool for detecting undescribed symbionts that are otherwise difficult to quantify, enabling further insight into the impact cryptic symbionts have on their hosts.