Global distribution, diversity, and ecological niche of Picozoa, a widespread and enigmatic marine protist lineage.

BACKGROUND: The backbone of the eukaryotic tree of life contains taxa only found in molecular surveys, of which we still have a limited understanding. Such is the case of Picozoa, an enigmatic lineage of heterotrophic picoeukaryotes within the supergroup Archaeplastida, which has emerged as a significant component of marine microbial planktonic communities. To enhance our understanding of the diversity, distribution, and ecology of Picozoa, we conduct a comprehensive assessment at different levels, from assemblages to taxa, employing phylogenetic analysis, species distribution modeling, and ecological niche characterization. RESULTS: Picozoa was among the ten most abundant eukaryotic groups, found almost exclusively in marine environments. The phylum was represented by 179 Picozoa's OTU (pOTUs) placed in five phylogenetic clades. Picozoa community structure had a clear latitudinal pattern, with polar assemblages tending to cluster separately from non-polar ones. Based on the abundance and occupancy pattern, the pOTUs were classified into four categories: Low-abundant, Widespread, Polar, and Non-polar. We calculated the ecological niche of each of these categories. Notably, pOTUs sharing similar ecological niches were not closely related species, indicating a phylogenetic overdispersion in Picozoa communities. This could be attributed to competitive exclusion and the strong influence of the seasonal amplitude of variations in environmental factors, such as temperature, shaping physiological and ecological traits. CONCLUSIONS: Overall, this work advances our understanding of uncharted protists' evolutionary dynamics and ecological strategies. Our results highlight the importance of understanding the species-level ecology of marine heteroflagellates like Picozoa. The observed phylogenetic overdispersion challenges the concept of phylogenetic niche conservatism in protist communities, suggesting that closely related species do not necessarily share similar ecological niches. Video Abstract.

Protist diversity and co-occurrence patterns obtained by metabarcoding of terricolous lichens, coastal cliffs and a microbial mat in the Atacama Desert, northern Chile.

Protists can endure challenging environments sustaining key ecosystem processes of the microbial food webs even under aridic or hypersaline conditions. We studied the diversity of protists at different latitudes of the Atacama Desert by massive sequencing of the hypervariable region V9 of the 18S rRNA gene from soils and microbial mats collected in the Andes. The main protist groups in soils detected in active stage through cDNA were cercozoans, ciliates, and kinetoplastids, while the diversity of protists was higher including diatoms and amoebae in the microbial mat detected solely through DNA. Co-occurrence networks from soils indicated similar assemblages dominated by amplicon sequence variants (ASVs) identified as Rhogostoma, Euplotes, and Neobodo. Microbial mat networks, on the other hand, were structured by ASVs classified as raphid-pennate diatoms and amoebae from the genera Hartmannella and Vannella, mostly negatively correlated to flagellates and microalgae. Additionally, our phylogenetic inferences of ASVs classified as Euplotes, Neobodo, and Rhogostoma were supported by sequence data of strains isolated during this study. Our results represent the first snapshot of the diversity patterns of culturable and unculturable protists and putative keystone taxa detected at remote habitats from the Atacama Desert.

Unveiling potential PET degrading eukaryotes through in silico bioprospecting of PETases.

This study addresses the environmental problem of PET plastic through in silico bioprospecting for the identification and experimental validation of novel PET degrading eukaryotes through the in silico bioprospectingI of PETases, employing a methodology that combines Hidden Markov Models (HMMs), clustering techniques, molecular docking, and dynamic simulations. A total of 424 putative PETase sequences were identified from 219 eukaryotic organisms, highlighting six sequences with low affinity energies. The Aspergillus luchuensis sequence showed the lowest Gibbs free energy and exhibited stability at different temperatures in molecular dynamics assays. Experimental validation, through a plate clearance assay and HPLC, confirmed PETase activity in three wild-type fungal strains, with A. luchuensis showing the highest efficiency. The results obtained demonstrate the effectiveness of combining computational and experimental approaches as proof of concept to discover and validate eukaryotes with PET-degrading capabilities opening new perspectives for the sustainable management of this type of waste and contributing to its environmental mitigation.

Amoebozoan testate amoebae illuminate the diversity of heterotrophs and the complexity of ecosystems throughout geological time.

Heterotrophic protists are vital in Earth's ecosystems, influencing carbon and nutrient cycles and occupying key positions in food webs as microbial predators. Fossils and molecular data suggest the emergence of predatory microeukaryotes and the transition to a eukaryote-rich marine environment by 800 million years ago (Ma). Neoproterozoic vase-shaped microfossils (VSMs) linked to Arcellinida testate amoebae represent the oldest evidence of heterotrophic microeukaryotes. This study explores the phylogenetic relationship and divergence times of modern Arcellinida and related taxa using a relaxed molecular clock approach. We estimate the origin of nodes leading to extant members of the Arcellinida Order to have happened during the latest Mesoproterozoic and Neoproterozoic (1054 to 661 Ma), while the divergence of extant infraorders postdates the Silurian. Our results demonstrate that at least one major heterotrophic eukaryote lineage originated during the Neoproterozoic. A putative radiation of eukaryotic groups (e.g., Arcellinida) during the early-Neoproterozoic sustained by favorable ecological and environmental conditions may have contributed to eukaryotic life endurance during the Cryogenian severe ice ages. Moreover, we infer that Arcellinida most likely already inhabited terrestrial habitats during the Neoproterozoic, coexisting with terrestrial Fungi and green algae, before land plant radiation. The most recent extant Arcellinida groups diverged during the Silurian Period, alongside other taxa within Fungi and flowering plants. These findings shed light on heterotrophic microeukaryotes' evolutionary history and ecological significance in Earth's ecosystems, using testate amoebae as a proxy.

A sedimentary DNA record of the Atacama Trench reveals biodiversity changes in the most productive marine ecosystem.

The hadopelagic environment remains highly understudied due to the inherent difficulties in sampling at these depths. The use of sediment environmental DNA (eDNA) can overcome some of these restrictions as settled and preserved DNA represent an archive of the biological communities. We use sediment eDNA to assess changes in the community within one of the world's most productive open-ocean ecosystems: the Atacama Trench. The ecosystems around the Atacama Trench have been intensively fished and are affected by climate oscillations, but the understanding of potential impacts on the marine community is limited. We sampled five sites using sediment cores at water depths from 2400 to ~8000 m. The chronologies of the sedimentary record were determined using 210Pbex. Environmental DNA was extracted from core slices and metabarcoding was used to identify the eukaryote community using two separate primer pairs for different sections of the 18S rRNA gene (V9 and V7) effectively targeting pelagic taxa. The reconstructed communities were similar among markers and mainly composed of chordates and members of the Chromista kingdom. Alpha diversity was estimated for all sites in intervals of 15 years (from 1842 to 2018), showing a severe drop in biodiversity from 1970 to 1985 that aligns with one of the strongest known El Niño events and extensive fishing efforts during the time. We find a direct impact of sea surface temperature on the community composition over time. Fish and cnidarian read abundance was examined separately to determine whether fishing had a direct impact, but no direct relation was found. These results demonstrate that sediment eDNA can be a valuable emerging tool providing insight in historical perspectives on ecosystem developments. This study constitutes an important step toward an improved understanding of the importance of environmental and anthropogenic drivers in affecting open and deep ocean communities.

DeepReg: a deep learning hybrid model for predicting transcription factors in eukaryotic and prokaryotic genomes.

Deep learning models (DLMs) have gained importance in predicting, detecting, translating, and classifying a diversity of inputs. In bioinformatics, DLMs have been used to predict protein structures, transcription factor-binding sites, and promoters. In this work, we propose a hybrid model to identify transcription factors (TFs) among prokaryotic and eukaryotic protein sequences, named Deep Regulation (DeepReg) model. Two architectures were used in the DL model: a convolutional neural network (CNN), and a bidirectional long-short-term memory (BiLSTM). DeepReg reached a precision of 0.99, a recall of 0.97, and an F1-score of 0.98. The quality of our predictions, the bias-variance trade-off approach, and the characterization of new TF predictions were evaluated and compared against those produced by DeepTFactor, as well as against experimental data from three model organisms. Predictions based on our DLM tended to exhibit less variance and bias than those from DeepTFactor, thus increasing reliability and decreasing overfitting.

Contribution of microscopy to a better understanding of the anatomy of pathogenic protists.

In this Inaugural Article the author briefly revises its scientific career and how he starts to work with parasitic protozoa. Emphasis is given to his contribution to topics such as a) the structural organization of the surface of protozoa using freeze-fracture and deep-etching; b) the cytoskeleton of protozoa, especially structures such as the subpellicular microtubules of trypanosomatids, the conoid of Toxoplasma gondii, microtubules and inner membrane complex of this protozoan, and the costa of Tritrichomonas foetus; c) the flagellulm of trypanosomatids, that in addition to the axoneme contains a complex network of filaments that constitute the paraflagellar rod; d) special organelles such as the acidocalcisome, hydrogenosome, and glycosome; and e) the highly polarized endocytic pathway found in epimastigote forms of Trypanosoma cruzi.

The consequences of viral infection on protists.

Protists encompass a vast widely distributed group of organisms, surpassing the diversity observed in metazoans. Their diverse ecological niches and life forms are intriguing characteristics that render them valuable subjects for in-depth cell biology studies. Throughout history, viruses have played a pivotal role in elucidating complex cellular processes, particularly in the context of cellular responses to viral infections. In this comprehensive review, we provide an overview of the cellular alterations that are triggered in specific hosts following different viral infections and explore intricate biological interactions observed in experimental conditions using different host-pathogen groups.

Functional analysis of the AUG initiator codon context reveals novel conserved sequences that disfavor mRNA translation in eukaryotes.

mRNA translation is a fundamental process for life. Selection of the translation initiation site (TIS) is crucial, as it establishes the correct open reading frame for mRNA decoding. Studies in vertebrate mRNAs discovered that a purine at -3 and a G at +4 (where A of the AUG initiator codon is numbered + 1), promote TIS recognition. However, the TIS context in other eukaryotes has been poorly experimentally analyzed. We analyzed in vitro the influence of the -3, -2, -1 and + 4 positions of the TIS context in rabbit, Drosophila, wheat, and yeast. We observed that -3A conferred the best translational efficiency across these species. However, we found variability at the + 4 position for optimal translation. In addition, the Kozak motif that was defined from mammalian cells was only weakly predictive for wheat and essentially non-predictive for yeast. We discovered eight conserved sequences that significantly disfavored translation. Due to the big differences in translational efficiency observed among weak TIS context sequences, we define a novel category that we termed 'barren AUG context sequences (BACS)', which represent sequences disfavoring translation. Analysis of mRNA-ribosomal complexes structures provided insights into the function of BACS. The gene ontology of the BACS-containing mRNAs is presented.

Global biogeography of the smallest plankton across ocean depths.

Tiny ocean plankton (picoplankton) are fundamental for the functioning of the biosphere, but the ecological mechanisms shaping their biogeography were partially understood. Comprehending whether these microorganisms are structured by niche versus neutral processes is relevant in the context of global change. We investigate the ecological processes (selection, dispersal, and drift) structuring global-ocean picoplanktonic communities inhabiting the epipelagic (0 to 200 meters), mesopelagic (200 to 1000 meters), and bathypelagic (1000 to 4000 meters) zones. We found that selection decreased, while dispersal limitation increased with depth, possibly due to differences in habitat heterogeneity and dispersal barriers such as water masses and bottom topography. Picoplankton ß-diversity positively correlated with environmental heterogeneity and water mass variability, but this relationship tended to be weaker for eukaryotes than for prokaryotes. Community patterns were more pronounced in the Mediterranean Sea, probably because of its cross-basin environmental heterogeneity and deep-water isolation. We conclude that different combinations of ecological mechanisms shape the biogeography of the ocean microbiome across depths.

Conserved and divergent signals in 5' splice site sequences across fungi, metazoa and plants.

In eukaryotic organisms the ensemble of 5' splice site sequences reflects the balance between natural nucleotide variability and minimal molecular constraints necessary to ensure splicing fidelity. This compromise shapes the underlying statistical patterns in the composition of donor splice site sequences. The scope of this study was to mine conserved and divergent signals in the composition of 5' splice site sequences. Because 5' donor sequences are a major cue for proper recognition of splice sites, we reasoned that statistical regularities in their composition could reflect the biological functionality and evolutionary history associated with splicing mechanisms. Results: We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequence datasets corresponding to 30 different eukaryotes. For each analyzed species, we identified minimal sets of two-site coupling patterns that were able to replicate, at a given regularization level, the observed one-site and two-site frequencies in donor sequences. By performing a systematic and comparative analysis of 5'splice sites we showed that lineage information could be traced from joint di-nucleotide probabilities. We were able to identify characteristic two-site coupling patterns for plants and animals, and propose that they may echo differences in splicing regulation previously reported between these groups.

Current knowledge of the Southern Hemisphere marine microbiome in eukaryotic hosts and the Strait of Magellan surface microbiome project.

Host-microbe interactions are ubiquitous and play important roles in host biology, ecology, and evolution. Yet, host-microbe research has focused on inland species, whereas marine hosts and their associated microbes remain largely unexplored, especially in developing countries in the Southern Hemisphere. Here, we review the current knowledge of marine host microbiomes in the Southern Hemisphere. Our results revealed important biases in marine host species sampling for studies conducted in the Southern Hemisphere, where sponges and marine mammals have received the greatest attention. Sponge-associated microbes vary greatly across geographic regions and species. Nevertheless, besides taxonomic heterogeneity, sponge microbiomes have functional consistency, whereas geography and aging are important drivers of marine mammal microbiomes. Seabird and macroalgal microbiomes in the Southern Hemisphere were also common. Most seabird microbiome has focused on feces, whereas macroalgal microbiome has focused on the epibiotic community. Important drivers of seabird fecal microbiome are aging, sex, and species-specific factors. In contrast, host-derived deterministic factors drive the macroalgal epibiotic microbiome, in a process known as "microbial gardening". In turn, marine invertebrates (especially crustaceans) and fish microbiomes have received less attention in the Southern Hemisphere. In general, the predominant approach to study host marine microbiomes has been the sequencing of the 16S rRNA gene. Interestingly, there are some marine holobiont studies (i.e., studies that simultaneously analyze host (e.g., genomics, transcriptomics) and microbiome (e.g., 16S rRNA gene, metagenome) traits), but only in some marine invertebrates and macroalgae from Africa and Australia. Finally, we introduce an ongoing project on the surface microbiome of key species in the Strait of Magellan. This is an international project that will provide novel microbiome information of several species in the Strait of Magellan. In the short-term, the project will improve our knowledge about microbial diversity in the region, while long-term potential benefits include the use of these data to assess host-microbial responses to the Anthropocene derived climate change.

Antarctic snow algae: unraveling the processes underlying microbial community assembly during blooms formation.

BACKGROUND AND AIMS: At the West Antarctic Peninsula, snow algae blooms are composed of complex microbial communities dominated by green microalgae and bacteria. During their progression, the assembly of these microbial communities occurs under harsh environmental conditions and variable nutrient content due to fast snow melting. To date, it is still unclear what are the ecological mechanisms governing the composition and abundance of microorganisms during the formation of snow algae blooms. In this study, we aim to examine the main ecological mechanisms governing the assembly of snow algae blooms from early stages to colorful stages blooms. METHODS: The composition of the microbial communities within snow algae blooms was recorded in the West Antarctic Peninsula (Isabel Riquelme Islet) during a 35-day period using 16S rRNA and 18S rRNA metabarcoding. In addition, the contribution of different ecological processes to the assembly of the microbial community was quantified using phylogenetic bin-based null model analysis. RESULTS: Our results showed that alpha diversity indices of the eukaryotic communities displayed a higher variation during the formation of the algae bloom compared with the bacterial community. Additionally, in a macronutrients rich environment, the content of nitrate, ammonium, phosphate, and organic carbon did not play a major role in structuring the community. The quantification of ecological processes showed that the bacterial community assembly was governed by selective processes such as homogenous selection. In contrast, stochastic processes such as dispersal limitation and drift, and to a lesser extent, homogenous selection, regulate the eukaryotic community. CONCLUSIONS: Overall, our study highlights the differences in the microbial assembly between bacteria and eukaryotes in snow algae blooms and proposes a model to integrate both assembly processes. Video Abstract.

Genomic object detection: An improved approach for transposable elements detection and classification using convolutional neural networks.

Analysis of eukaryotic genomes requires the detection and classification of transposable elements (TEs), a crucial but complex and time-consuming task. To improve the performance of tools that accomplish these tasks, Machine Learning approaches (ML) that leverage computer resources, such as GPUs (Graphical Processing Unit) and multiple CPU (Central Processing Unit) cores, have been adopted. However, until now, the use of ML techniques has mostly been limited to classification of TEs. Herein, a detection-classification strategy (named YORO) based on convolutional neural networks is adapted from computer vision (YOLO) to genomics. This approach enables the detection of genomic objects through the prediction of the position, length, and classification in large DNA sequences such as fully sequenced genomes. As a proof of concept, the internal protein-coding domains of LTR-retrotransposons are used to train the proposed neural network. Precision, recall, accuracy, F1-score, execution times and time ratios, as well as several graphical representations were used as metrics to measure performance. These promising results open the door for a new generation of Deep Learning tools for genomics. YORO architecture is available at https://github.com/simonorozcoarias/YORO.

Response of protists to nitrogen addition, arbuscular mycorrhizal fungi manipulation, and mesofauna reduction in a tropical montane rainforest in southern Ecuador.

The tropical Andes are a species-rich and nitrogen-limited system, susceptible to increased nitrogen (N) inputs from the atmosphere. However, our understanding of the impacts of increased N input on belowground systems, in particular on protists and their role in nutrient cycling, remains limited. We explored how increased N affects protists in tropical montane rainforests in Ecuador using high-throughput sequencing (HTS) of environmental DNA from two litter layers. In addition, we manipulated the amount of arbuscular mycorrhizal fungi (AMF) and mesofauna, both playing a significant role in N cycling and interacting in complex ways with protist communities. We found that N strongly affected protist community composition in both layers, while mesofauna reduction had a stronger effect on the lower layer. Changes in concentration of the AMF marker lipid had little effect on protists. In both layers, the addition of N increased phagotrophs and animal parasites and decreased plant parasites, while mixotrophs decreased in the upper layer but increased in the lower layer. In the upper layer with higher AMF concentration, mixotrophs decreased, while in the lower layer, photoautotrophs increased and plant parasites decreased. With reduced mesofauna, phagotrophs increased and animal parasites decreased in both layers, while plant parasites increased only in the upper layer. The findings indicate that to understand the intricate response of protist communities to environmental changes, it is critical to thoroughly analyze these communities across litter and soil layers, and to include HTS.

Diversity of eukaryotic and prokaryotic microbiota revealed by metabarcoding in Neotropical floodplain lakes.

The diversity of eukaryotic and prokaryotic communities has been assessed by morphological and genetic approaches, which are used to characterize the microbiota in different environments. Here, planktonic prokaryotic and eukaryotic communities of the Araguaia River, located in the Central region of Brazil, were analyzed based on metabarcoding analysis of rRNA genes to evaluate the diversity of these groups in tropical floodplain lakes. Also, we tested their spatial concordance throughout the Araguaia river. Water samples were collected from 8 floodplain lakes in Araguaia River. The 16S and 18S rRNA genes were amplified and sequenced using Illumina MiSeq. For eukaryotes, 34,242 merged reads were obtained and 225 distinct OTUs were delineated, of which 106 OTUs were taxonomically classified. For prokaryotes, 26,426 sequences were obtained and 351 OTUs were detected. Of them, 231 were classified in at least one taxonomic category. The most representative eukaryotes belonged to Ciliophora, Chlorophyta and Charophyta. The prokaryotic phylum with the most OTUs classified were Proteobacteria, Actinobacteria and Bacteroidetes. The lakes did not show spatial concordance when comparing the similarity between their microbiota. The knowledge of freshwater biodiversity using DNA sequencing for important rivers, such as Araguaia River, can improve microbiota inventories of tropical biodiversity hotspots.

Pangenomic Analysis of Nucleo-Cytoplasmic Large DNA Viruses. I: The Phylogenetic Distribution of Conserved Oxygen-Dependent Enzymes Reveals a Capture-Gene Process.

The Nucleo-Cytoplasmic Large DNA Viruses (NCLDVs) infect a wide range of eukaryotic species, including amoeba, algae, fish, amphibia, arthropods, birds, and mammals. This group of viruses has linear or circular double-stranded DNA genomes whose size spans approximately one order of magnitude, from 100 to 2500 kbp. The ultimate origin of this peculiar group of viruses remains an open issue. Some have argued that NCLDVs' origin may lie in a bacteriophage ancestor that increased its genome size by subsequent recruitment of eukaryotic and bacterial genes. Others have suggested that NCLDVs families originated from cells that underwent an irreversible process of genome reduction. However, the hypothesis that a number of NCLDVs sequences have been recruited from the host genomes has been largely ignored. In the present work, we have performed pangenomic analyses of each of the seven known NCLDVs families. We show that these families' core- and shell genes have cellular homologs, supporting possible escaping-gene events as part of its evolution. Furthermore, the detection of sequences that belong to two protein families (small chain ribonucleotide reductase and Erv1/Air) and to one superfamily [2OG-Fe(II) oxygenases] that are for distribution in all NCLDVs core and shell clusters encoding for oxygen-dependent enzymes suggests that the highly conserved core these viruses originated after the Proterozoic Great Oxidation Event that transformed the terrestrial atmosphere 2.4-2.3 Ga ago.

Evaluation of amino acid profile by targeted metabolomics in the eukaryotic model under exposure of benzo[a]pyrene as the exclusive stressor.

Amino acids (AAs) are a class of important metabolites in metabolomics methodology that investigates metabolite changes in a cell, tissue, or organism for early diagnosis of diseases. Benzo[a]pyrene (BaP) is considered a priority contaminant by different environmental control agencies because it is a proven carcinogenic compound for humans. Therefore, it is important to evaluate the BaP interference in the metabolism of amino acids. In this work, a new amino acid extraction procedure (derivatized with propyl chloroformate/propanol) using functionalized magnetic carbon nanotubes was developed and optimized. A hybrid nanotube was used followed by desorption without heating, and excellent extraction of analytes was obtained. After exposure of Saccharomyces cerevisiae, the BaP concentration of 25.0 µmol L-1 caused changes in cell viability, indicating metabolic changes. A fast and efficient GC/MS method using a Phenomenex ZB-AAA column was optimized, enabling the determination of 16 AAs in yeasts exposed or not to BaP. A comparison of AA concentrations obtained in the two experimental groups showed that glycine (Gly), serine (Ser), phenylalanine (Phe), proline (Pro), asparagine (Asn), aspartic acid (Asp), glutamic acid (Glu), tyrosine (Tyr), and leucine (Leu) statistically differentiated, after subsequent application of ANOVA with Bonferroni post-hoc test, with a confidence level of 95%. This amino acid pathway analysis confirmed previous studies that revealed the potential of these AAs as toxicity biomarker candidates.