Four high-quality draft genome assemblies of the marine heterotrophic nanoflagellate Cafeteria roenbergensis

Endogenous virophages are active and mitigate giant virus infection in the marine protist Cafeteria burkhardae

Endogenous viral elements (EVEs) are common genetic passengers in various protists. Some EVEs represent viral fossils, whereas others are still active. The marine heterotrophic flagellate Cafeteria burkhardae contains several EVE types related to the virophage mavirus, a small DNA virus that parasitizes the lytic giant virus CroV. We hypothesized that endogenous virophages may act as an antiviral defense system in protists, but no protective effect of virophages in wild host populations has been shown so far. Here, we tested the activity of virophage EVEs and studied their impact on giant virus replication. We found that endogenous mavirus-like elements (EMALEs) from globally distributed Cafeteria populations produced infectious virus particles specifically in response to CroV infection. However, reactivation was stochastic, often inefficient, and poorly reproducible. Interestingly, only one of eight EMALE types responded to CroV infection, implying that other EMALEs may be linked to different giant viruses. We isolated and cloned several reactivated virophages and characterized their particles, genomes, and infection dynamics. All tested virophages inhibited the production of CroV during coinfection, thereby preventing lysis of the host cultures in a dose-dependent manner. Comparative genomics of different C. burkhardae strains revealed that inducible EMALEs are common and are not linked to specific geographic locations. We demonstrate that naturally occurring virophage EVEs reactivate upon giant virus infection, thus providing a striking example that eukaryotic EVEs can become active under specific conditions. Moreover, our results support the hypothesis that virophages can act as an adaptive antiviral defense system in protists.

Membrane-anchored calpains - hidden regulators of growth and development beyond plants?

Calpains are modulatory proteases that modify diverse cellular substrates and play essential roles in eukaryots. The best studied are animal cytosolic calpains. Here, we focus on enigmatic membrane-anchored calpains, their structural and functional features as well as phylogenetic distribution. Based on domain composition, we identified four types of membrane-anchored calpains. Type 1 and 2 show broad phylogenetic distribution among unicellular protists and streptophytes suggesting their ancient evolutionary origin. Type 3 and 4 diversified early and are present in brown algae and oomycetes. The plant DEK1 protein is the only representative of membrane-anchored calpains that has been functionally studied. Here, we present up to date knowledge about its structural features, putative regulation, posttranslational modifications, and biological role. Finally, we discuss potential model organisms and available tools for functional studies of membrane-anchored calpains with yet unknown biological role. Mechanistic understanding of membrane-anchored calpains may provide important insights into fundamental principles of cell polarization, cell fate control, and morphogenesis beyond plants.

Hybrid assemblies of microbiome Blastocystis protists reveal evolutionary diversification reflecting host ecology

The most prevalent microbial eukaryote in the human gut is Blastocystis , an obligate commensal protist also common in many other vertebrates. Blastocystis is descended from free-living stramenopile ancestors; how it has adapted to thrive within humans and a wide range of hosts is unclear. Here, we cultivated six Blastocystis strains spanning the diversity of the genus and generated highly contiguous, annotated genomes with long-read DNA-seq, Hi-C, and RNA-seq. Comparative genomics between these strains and two closely related stramenopiles with different lifestyles, the lizard gut symbiont Proteromonas lacertae and the free-living marine flagellate Cafeteria burkhardae , reveal the evolutionary history of the Blastocystis genus. We find substantial gene content variability between Blastocystis strains. Blastocystis isolated from an herbivorous tortoise has many plant carbohydrate metabolizing enzymes, some horizontally acquired from bacteria, likely reflecting fermentation within the host gut. In contrast, human- isolated Blastocystis have gained many heat shock proteins, and we find numerous subtype- specific expansions of host-interfacing genes, including cell adhesion and cell surface glycan genes. In addition, we observe that human-isolated Blastocystis have substantial changes in gene structure, including shortened introns and intergenic regions, as well as genes lacking canonical termination codons. Finally, our data indicate that the common ancestor of Blastocystis lost nearly all ancestral genes for heterokont flagella morphology, including cilia proteins, microtubule motor proteins, and ion channel proteins. Together, these findings underscore the huge functional variability within the Blastocystis genus and provide candidate genes for the adaptations these lineages have undergone to thrive in the gut microbiomes of diverse vertebrates.

Updated Virophage Taxonomy and Distinction from Polinton-like Viruses

Virophages are small dsDNA viruses that hijack the machinery of giant viruses during the co-infection of a protist (i.e., microeukaryotic) host and represent an exceptional case of "hyperparasitism" in the viral world. While only a handful of virophages have been isolated, a vast diversity of virophage-like sequences have been uncovered from diverse metagenomes. Their wide ecological distribution, idiosyncratic infection and replication strategy, ability to integrate into protist and giant virus genomes and potential role in antiviral defense have made virophages a topic of broad interest. However, one limitation for further studies is the lack of clarity regarding the nomenclature and taxonomy of this group of viruses. Specifically, virophages have been linked in the literature to other "virophage-like" mobile genetic elements and viruses, including polinton-like viruses (PLVs), but there are no formal demarcation criteria and proper nomenclature for either group, i.e., virophage or PLVs. Here, as part of the ICTV Virophage Study Group, we leverage a large set of genomes gathered from published datasets as well as newly generated protist genomes to propose delineation criteria and classification methods at multiple taxonomic ranks for virophages 'sensu stricto', i.e., genomes related to the prototype isolates Sputnik and mavirus. Based on a combination of comparative genomics and phylogenetic analyses, we show that this group of virophages forms a cohesive taxon that we propose to establish at the class level and suggest a subdivision into four orders and seven families with distinctive ecogenomic features. Finally, to illustrate how the proposed delineation criteria and classification method would be used, we apply these to two recently published datasets, which we show include both virophages and other virophage-related elements. Overall, we see this proposed classification as a necessary first step to provide a robust taxonomic framework in this area of the virosphere, which will need to be expanded in the future to cover other virophage-related viruses such as PLVs.

Crystal structures of FNIP/FGxxFN motif-containing leucine-rich repeat proteins

The Cafeteria roenbergensis virus (Crov), Dictyostelium, and other species encode a large family of leucine-rich repeat (LRR) proteins with FGxxFN motifs. We determined the structures of two of them and observed several unique structural features that set them aside from previously characterized LRR family members. Crov588 comprises 25 regular repeats with a LxxLxFGxxFNQxIxENVLPxx consensus, forming a unique closed circular repeat structure. Novel features include a repositioning of a conserved asparagine at the middle of the repeat, a double phenylalanine spine that generates an alternate core packing arrangement, and a histidine/tyrosine ladder on the concave surface. Crov539 is smaller, comprising 12 repeats of a similar LxxLxFGxxFNQPIExVxW/LPxx consensus and forming an unusual cap-swapped dimer structure. The phenylalanine spine of Crov539 is supplemented with a tryptophan spine, while a hydrophobic isoleucine-rich patch is found on the central concave surface. We present a detailed analysis of the structures of Crov588 and Crov539 and compare them to related repeat proteins and other LRR classes.

Draft Genome Sequence of the Freshwater Diatom Fragilaria crotonensis SAG 28.96

Here, we report the assembled and annotated genome of the freshwater diatom Fragilaria crotonensis SAG 28.96. The 61.85-Mb nuclear genome was assembled into 879 contigs, has a GC content of 47.40%, contains 26,015 predicted genes, and shows completeness of 81%.

Aureococcus anophagefferens (Pelagophyceae) genomes improve evaluation of nutrient acquisition strategies involved in brown tide dynamics

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment.

Emerging Diversity of Channelrhodopsins and Their Structure-Function Relationships

Cation and anion channelrhodopsins (CCRs and ACRs, respectively) from phototactic algae have become widely used as genetically encoded molecular tools to control cell membrane potential with light. Recent advances in polynucleotide sequencing, especially in environmental samples, have led to identification of hundreds of channelrhodopsin homologs in many phylogenetic lineages, including non-photosynthetic protists. Only a few CCRs and ACRs have been characterized in detail, but there are indications that ion channel function has evolved within the rhodopsin superfamily by convergent routes. The diversity of channelrhodopsins provides an exceptional platform for the study of structure-function evolution in membrane proteins. Here we review the current state of channelrhodopsin research and outline perspectives for its further development.

MOSGA 2: Comparative genomics and validation tools

Due to the highly growing number of available genomic information, the need for accessible and easy-to-use analysis tools is increasing. To facilitate eukaryotic genome annotations, we created MOSGA. In this work, we show how MOSGA 2 is developed by including several advanced analyses for genomic data. Since the genomic data quality greatly impacts the annotation quality, we included multiple tools to validate and ensure high-quality user-submitted genome assemblies. Moreover, thanks to the integration of comparative genomics methods, users can benefit from a broader genomic view by analyzing multiple genomic data sets simultaneously. Further, we demonstrate the new functionalities of MOSGA 2 by different use-cases and practical examples. MOSGA 2 extends the already established application to the quality control of the genomic data and integrates and analyzes multiple genomes in a larger context, e.g., by phylogenetics.

Virophages and retrotransposons colonize the genomes of a heterotrophic flagellate

Virophages can parasitize giant DNA viruses and may provide adaptive anti-giant virus defense in unicellular eukaryotes. Under laboratory conditions, the virophage mavirus integrates into the nuclear genome of the marine flagellate Cafeteria burkhardae and reactivates upon superinfection with the giant virus CroV. In natural systems, however, the prevalence and diversity of host-virophage associations has not been systematically explored. Here, we report dozens of integrated virophages in four globally sampled C. burkhardae strains that constitute up to 2% of their host genomes. These endogenous mavirus-like elements (EMALEs) separated into eight types based on GC-content, nucleotide similarity, and coding potential and carried diverse promoter motifs implicating interactions with different giant viruses. Between host strains, some EMALE insertion loci were conserved indicating ancient integration events, whereas the majority of insertion sites were unique to a given host strain suggesting that EMALEs are active and mobile. Furthermore, we uncovered a unique association between EMALEs and a group of tyrosine recombinase retrotransposons, revealing yet another layer of parasitism in this nested microbial system. Our findings show that virophages are widespread and dynamic in wild Cafeteria populations, supporting their potential role in antiviral defense in protists.

Viruses exist in all ecosystems in vast numbers and infect many organisms. Some of them are harmful but others can protect the organisms they infect. For example, a group of viruses called virophages protect microscopic sea creatures called plankton from deadly infections by so-called giant viruses. In fact, virophages need plankton infected with giant viruses to survive because they use enzymes from the giant viruses to turn on their own genes.

A virophage called mavirus integrates its genes into the DNA of a type of plankton called Cafeteria. It lays dormant in the DNA until a giant virus called CroV infects the plankton. This suggests that the mavirus may be a built-in defense against CroV infections and laboratory studies seem to confirm this. But whether wild Cafeteria also use these defenses is unknown.

Hackl et al. show that virophages are common in the DNA of wild Cafeteria and that the two appear to have a mutually beneficial relationship. In the experiments, the researchers sequenced the genomes of four Cafeteria populations from the Atlantic and Pacific Oceans and looked for virophages in their DNA. Each of the four Cafeteria genomes contained dozens of virophages, which suggests that virophages are important to these plankton. This included several relatives of the mavirus and seven new virophages. Virophage genes were often interrupted by so called jumping genes, which may take advantage of the virophages the way the virophages use giant viruses to meet their own needs.

The experiments show that virophages often co-exist with marine plankton from around the world and these relationships are likely beneficial. In fact, the experiments suggest that the virophages may have played an important role in the evolution of these plankton. Further studies may help scientists learn more about virus ecology and how viruses have shaped the evolution of other creatures.

Niche adaptation promoted the evolutionary diversification of tiny ocean predators

Unicellular eukaryotic predators play a crucial role in the functioning of the ocean ecosystem by recycling nutrients and energy that are channeled to upper trophic levels. Traditionally, these evolutionarily diverse organisms have been combined into a single functional group (heterotrophic flagellates), overlooking their organismal differences. Here, we investigated four evolutionarily related species belonging to one cosmopolitan group of uncultured marine picoeukaryotic predators: marine stramenopiles (MAST)-4 (species A, B, C, and E). Co-occurrence and distribution analyses in the global surface ocean indicated contrasting patterns in MAST-4A and C, suggesting adaptation to different temperatures. We then investigated whether these spatial distribution patterns were mirrored by MAST-4 genomic content using single-cell genomics. Analyses of 69 single cells recovered 66 to 83% of the MAST-4A/B/C/E genomes, which displayed substantial interspecies divergence. MAST-4 genomes were similar in terms of broad gene functional categories, but they differed in enzymes of ecological relevance, such as glycoside hydrolases (GHs), which are part of the food degradation machinery in MAST-4. Interestingly, MAST-4 species featuring a similar GH composition (A and C) coexcluded each other in the surface global ocean, while species with a different set of GHs (B and C) appeared to be able to coexist, suggesting further niche diversification associated with prey digestion. We propose that differential niche adaptation to temperature and prey type has promoted adaptive evolutionary diversification in MAST-4. We show that minute ocean predators from the same phylogenetic group may have different biogeography and genomic content, which needs to be accounted for to better comprehend marine food webs.

MOSGA: Modular Open-Source Genome Annotator

MOTIVATION

The generation of high-quality assemblies, even for large eukaryotic genomes, has become a routine task for many biologists thanks to recent advances in sequencing technologies. However, the annotation of these assemblies-a crucial step toward unlocking the biology of the organism of interest-has remained a complex challenge that often requires advanced bioinformatics expertise.

RESULTS

Here, we present MOSGA (Modular Open-Source Genome Annotator), a genome annotation framework for eukaryotic genomes with a user-friendly web-interface that generates and integrates annotations from various tools. The aggregated results can be analyzed with a fully integrated genome browser and are provided in a format ready for submission to NCBI. MOSGA is built on a portable, customizable and easily extendible Snakemake backend, and thus, can be tailored to a wide range of users and projects.

AVAILABILITY AND IMPLEMENTATION

We provide MOSGA as a web service at https://mosga.mathematik.uni-marburg.de and as a docker container at registry.gitlab.com/mosga/mosga: latest. Source code can be found at https://gitlab.com/mosga/mosga.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Gene expression during bacterivorous growth of a widespread marine heterotrophic flagellate

Phagocytosis is a fundamental process in marine ecosystems by which prey organisms are consumed and their biomass incorporated in food webs or remineralized. However, studies searching for the genes underlying this key ecological process in free-living phagocytizing protists are still scarce, in part due to the lack of appropriate ecological models. Our reanalysis of recent molecular datasets revealed that the cultured heterotrophic flagellate Cafeteria burkhardae is widespread in the global oceans, which prompted us to design a transcriptomics study with this species, grown with the cultured flavobacterium Dokdonia sp. We compared the gene expression between exponential and stationary phases, which were complemented with three starvation by dilution phases that appeared as intermediate states. We found distinct expression profiles in each condition and identified 2056 differentially expressed genes between exponential and stationary samples. Upregulated genes at the exponential phase were related to DNA duplication, transcription and translational machinery, protein remodeling, respiration and phagocytosis, whereas upregulated genes in the stationary phase were involved in signal transduction, cell adhesion, and lipid metabolism. We identified a few highly expressed phagocytosis genes, like peptidases and proton pumps, which could be used to target this ecologically relevant process in marine ecosystems.

Estimating the quality of eukaryotic genomes recovered from metagenomic analysis with EukCC

Microbial eukaryotes constitute a significant fraction of biodiversity and have recently gained more attention, but the recovery of high-quality metagenomic assembled eukaryotic genomes is limited by the current availability of tools. To help address this, we have developed EukCC, a tool for estimating the quality of eukaryotic genomes based on the automated dynamic selection of single copy marker gene sets. We demonstrate that our method outperforms current genome quality estimators, particularly for estimating contamination, and have applied EukCC to datasets derived from two different environments to enable the identification of novel eukaryote genomes, including one from the human skin.

Polycystic Kidney Disease Ryanodine Receptor Domain (PKDRR) Proteins in Oomycetes

In eukaryotes, two sources of Ca²⁺ are accessed to allow rapid changes in the cytosolic levels of this second messenger: the extracellular medium and intracellular Ca²⁺ stores, such as the endoplasmic reticulum. One class of channel that permits Ca²⁺ entry is the transient receptor potential (TRP) superfamily, including the polycystic kidney disease (PKD) proteins, or polycystins. Channels that release Ca²⁺ from intracellular stores include the inositol 1,4,5-trisphosphate/ryanodine receptor (ITPR/RyR) superfamily. Here, we characterise a family of proteins that are only encoded by oomycete genomes, that we have named PKDRR, since they share domains with both PKD and RyR channels. We provide evidence that these proteins belong to the TRP superfamily and are distinct from the ITPR/RyR superfamily in terms of their evolutionary relationships, protein domain architectures and predicted ion channel structures. We also demonstrate that a hypothetical PKDRR protein from Phytophthora infestans is produced by this organism, is located in the cell-surface membrane and forms multimeric protein complexes. Efforts to functionally characterise this protein in a heterologous expression system were unsuccessful but support a cell-surface localisation. These PKDRR proteins represent potential targets for the development of new "fungicides", since they are of a distinctive structure that is only found in oomycetes and not in any other cellular organisms.