Evolution of Rhizaria: new insights from phylogenomic analysis of uncultivated protists

Antiviral Mx proteins have an ancient origin and widespread distribution among eukaryotes

First identified in mammals, Mx proteins are potent antivirals against a broad swathe of viruses. Mx proteins arose within the Dynamin superfamily of proteins (DSP), mediating critical cellular processes, such as endocytosis and mitochondrial, plastid, and peroxisomal dynamics. And yet, the evolutionary origins of Mx proteins are poorly understood. Using a series of phylogenomic analyses with stepwise increments in taxonomic coverage, we show that Mx proteins predate the interferon signaling system in vertebrates. Our analyses find an ancient monophyletic DSP lineage in eukaryotes that groups vertebrate and invertebrate Mx proteins with previously undescribed fungal MxF proteins, the relatively uncharacterized plant and algal Dynamin 4A/4C proteins, and representatives from several early-branching eukaryotic lineages. Thus, Mx-like proteins date back close to the origin of Eukarya. Our phylogenetic analyses also reveal that host-encoded and NCLDV (nucleocytoplasmic large DNA viruses)-encoded DSPs are interspersed in four distinct DSP lineages, indicating recurrent viral theft of host DSPs. Our analyses thus reveal an ancient history of viral and antiviral functions encoded by the Dynamin superfamily in eukaryotes.

Comparative genomics of Ascetosporea gives new insight into the evolutionary basis for animal parasitism in Rhizaria

BACKGROUND

Ascetosporea (Endomyxa, Rhizaria) is a group of unicellular parasites infecting aquatic invertebrates. They are increasingly being recognized as widespread and important in marine environments, causing large annual losses in invertebrate aquaculture. Despite their importance, little molecular data of Ascetosporea exist, with only two genome assemblies published to date. Accordingly, the evolutionary origin of these parasites is unclear, including their phylogenetic position and the genomic adaptations that accompanied the transition from a free-living lifestyle to parasitism. Here, we sequenced and assembled three new ascetosporean genomes, as well as the genome of a closely related amphizoic species, to investigate the phylogeny, origin, and genomic adaptations to parasitism in Ascetosporea.

RESULTS

Using a phylogenomic approach, we confirm the monophyly of Ascetosporea and show that Paramyxida group with Mikrocytida, with Haplosporida being sister to both groups. We report that the genomes of these parasites are relatively small (12-36 Mb) and gene-sparse (~ 2300-5200 genes), while containing surprisingly high amounts of non-coding sequence (~ 70-90% of the genomes). Performing gene-tree aware ancestral reconstruction of gene families, we demonstrate extensive gene losses at the origin of parasitism in Ascetosporea, primarily of metabolic functions, and little gene gain except on terminal branches. Finally, we highlight some functional gene classes that have undergone expansions during evolution of the group.

CONCLUSIONS

We present important new genomic information from a lineage of enigmatic but important parasites of invertebrates and illuminate some of the genomic innovations accompanying the evolutionary transition to parasitism in this lineage. Our results and data provide a genetic basis for the development of control measures against these parasites.

Starch Plays a Key Role in Sporosorus Formation by the Powdery Scab Pathogen Spongospora subterranea

Powdery scab disease, caused by the soilborne protist Spongospora subterranea f. sp. subterranea, poses a major constraint to potato production worldwide. Disease symptoms include damage to the tuber skin and the formation of root galls. This study aimed to investigate the potential mechanism behind the formation of sporosori, which are aggregates of resting spores, within root galls. Scanning electron microscopy analysis revealed that the early stage of gall formation, characterized by a white color, involved the accumulation of starch grains, which later disappeared as the gall matured and turned brown. The mature brown galls were found to contain fully formed sporosori. Light microscopy examination of ultramicrotome sections of the root galls showed that the high-amylopectin starches were surrounded by a plasmodium, a precursor to sporosorus. These findings suggest that starch grains contribute to the formation of a sponge-like structure within the sporosori. A significant reduction in total starch levels in both the root galls and their associated roots was observed compared with healthy roots. These findings indicate starch consumption by sporosori during the maturation of root galls. Interestingly, analysis of the transcript levels of starch-related genes showed downregulation of genes encoding starch degrading enzymes and an amylopectin-debranching enzyme, whereas genes encoding a starch synthase and a protein facilitating starch synthesis were upregulated in the infected roots. Overall, our results demonstrate that starch is consumed during sporosorus formation, and the pathogen likely manipulates starch homeostasis to its advantage for sporosorus development within the root galls.

Plasmodiophora brassicae affects host gene expression by secreting the transcription factor-type effector PbZFE1

The protist pathogen Plasmodiophora brassicae hijacks the metabolism and development of host cruciferous plants and induces clubroot formation, but little is known about its regulatory mechanisms. Previously, the Pnit2int2 sequence, a sequence around the second intron of the nitrilase gene (BrNIT2) involved in auxin biosynthesis in Brassica rapa ssp. pekinensis, was identified as a specific promoter activated during clubroot formation. In this study, we hypothesized that analysis of the transcriptional regulation of Pnit2int2 could reveal how P. brassicae affects the host gene regulatory system during clubroot development. By yeast one-hybrid screening, the pathogen zinc finger protein PbZFE1 was identified to specifically bind to Pnit2int2. Specific binding of PbZFE1 to Pnit2int2 was also confirmed by electrophoretic mobility shift assay. The binding site of PbZFE1 is essential for promoter activity of Pnit2int2 in clubbed roots of transgenic Arabidopsis thaliana (Pnit2int2-2::GUS), indicating that PbZFE1 is secreted from P. brassicae and functions within plant cells. Ectopic expression of PbZEF1 in A. thaliana delayed growth and flowering time, suggesting that PbZFE1 has significant impacts on host development and metabolic systems. Thus, P. brassicae appears to secrete PbZFE1 into host cells as a transcription factor-type effector during pathogenesis.

Lipid bilayer properties potentially contributed to the evolutionary disappearance of betaine lipids in seed plants

BACKGROUND

Many organisms rely on mineral nutrients taken directly from the soil or aquatic environment, and therefore, developed mechanisms to cope with the limitation of a given essential nutrient. For example, photosynthetic cells have well-defined responses to phosphate limitation, including the replacement of cellular membrane phospholipids with non-phosphorous lipids. Under phosphate starvation, phospholipids in extraplastidial membranes are replaced by betaine lipids in microalgae. In higher plants, the synthesis of betaine lipid is lost, driving plants to other strategies to cope with phosphate starvation where they replace their phospholipids by glycolipids.

RESULTS

The aim of this work was to evaluate to what extent betaine lipids and PC lipids share physicochemical properties and could substitute for each other. By neutron diffraction experiments and dynamic molecular simulation of two synthetic lipids, the dipalmitoylphosphatidylcholine (DPPC) and the dipalmitoyl-diacylglyceryl-N,N,N-trimethylhomoserine (DP-DGTS), we found that DP-DGTS bilayers are thicker than DPPC bilayers and therefore are more rigid. Furthermore, DP-DGTS bilayers are more repulsive, especially at long range, maybe due to unexpected unscreened electrostatic contribution. Finally, DP-DGTS bilayers could coexist in the gel and fluid phases.

CONCLUSION

The different properties and hydration responses of PC and DGTS provide an explanation for the diversity of betaine lipids observed in marine organisms and for their disappearance in seed plants.

Comparative spatial proteomics of Plasmodium-infected erythrocytes

Plasmodium parasites contribute to one of the highest global infectious disease burdens. To achieve this success, the parasite has evolved a range of specialized subcellular compartments to extensively remodel the host cell for its survival. The information to fully understand these compartments is likely hidden in the so far poorly characterized Plasmodium species spatial proteome. To address this question, we determined the steady-state subcellular location of more than 12,000 parasite proteins across five different species by extensive subcellular fractionation of erythrocytes infected by Plasmodium falciparum, Plasmodium knowlesi, Plasmodium yoelii, Plasmodium berghei, and Plasmodium chabaudi. This comparison of the pan-species spatial proteomes and their expression patterns indicates increasing species-specific proteins associated with the more external compartments, supporting host adaptations and post-transcriptional regulation. The spatial proteome offers comprehensive insight into the different human, simian, and rodent Plasmodium species, establishing a powerful resource for understanding species-specific host adaptation processes in the parasite.

Enzymatic Investigation of Spongospora subterranea Zoospore Attachment to Roots of Potato Cultivars Resistant or Susceptible to Powdery Scab Disease

For potato crops, host resistance is currently the most effective and sustainable tool to manage diseases caused by the plasmodiophorid Spongospora subterranea. Arguably, zoospore root attachment is the most critical phase of infection; however, the underlying mechanisms remain unknown. This study investigated the potential role of root-surface cell-wall polysaccharides and proteins in cultivars resistant/susceptible to zoospore attachment. We first compared the effects of enzymatic removal of root cell-wall proteins, N-linked glycans and polysaccharides on S. subterranea attachment. Subsequent analysis of peptides released by trypsin shaving (TS) of root segments identified 262 proteins that were differentially abundant between cultivars. These were enriched in root-surface-derived peptides but also included intracellular proteins, e.g., proteins associated with glutathione metabolism and lignin biosynthesis, which were more abundant in the resistant cultivar. Comparison with whole-root proteomic analysis of the same cultivars identified 226 proteins specific to the TS dataset, of which 188 were significantly different. Among these, the pathogen-defence-related cell-wall protein stem 28 kDa glycoprotein and two major latex proteins were significantly less abundant in the resistant cultivar. A further major latex protein was reduced in the resistant cultivar in both the TS and whole-root datasets. In contrast, three glutathione S-transferase proteins were more abundant in the resistant cultivar (TS-specific), while the protein glucan endo-1,3-beta-glucosidase was increased in both datasets. These results imply a particular role for major latex proteins and glucan endo-1,3-beta-glucosidase in regulating zoospore binding to potato roots and susceptibility to S. subterranea.

Microbial food webs share similar biogeographic patterns and driving mechanisms with depths in oligotrophic tropical western Pacific Ocean

Microbial food web (MFW) dominates the energy flow in oligotrophic tropical open ocean pelagic ecosystems. Understanding biogeographic patterns and driving mechanisms of key components of the MFW is one of the central topics in current marine ecology. Investigations were conducted along an 1,100-km horizontal gradient and in the full-water column vertical gradient of the oligotrophic tropical western Pacific Ocean. High-throughput sequencing and association networking methods were used to analyze the community structure and interspecies interactions of MFW. The structure of MFW significantly differed with depths, but not across horizontal gradients. Bacteria and microeukaryotes were interconnected and had more predominantly positive and negative linkages in the aphotic layers. Key components of MFW exhibited similar biogeographic patterns and driving mechanisms. Geographic distance exerted minimal effects on the distribution patterns of the microbial food web, while environmental factors played more important roles, especially for temperature and inorganic nutrients. Stochastic processes were more important in the microbial food webs of the 5-200  m layer than the >500  m layer, and drift explained the majority of stochastic processes. Moreover, only a weak but not significant driving force for North Equatorial Current on the east-west connectivity of the microbial food web was found in the upper layers. This knowledge is a critical fundamental data for future planning of marine protected areas targeting the protection of tuna fishing in the western Pacific Ocean.

Metabarcoding reveals high diversity of benthic foraminifera linked to water masses circulation at coastal Svalbard

Arctic marine biodiversity is undergoing rapid changes due to global warming and modifications of oceanic water masses circulation. These changes have been demonstrated in the case of mega- and macrofauna, but much less is known about their impact on the biodiversity of smaller size organisms, such as foraminifera that represent a main component of meiofauna in the Arctic. Several studies analyzed the distribution and diversity of Arctic foraminifera. However, all these studies are based exclusively on the morphological identification of specimens sorted from sediment samples. Here, we present the first assessment of Arctic foraminifera diversity based on metabarcoding of sediment DNA samples collected in fjords and open sea areas in the Svalbard Archipelago. We obtained a total of 5,968,786 reads that represented 1384 amplicon sequence variants (ASVs). More than half of the ASVs (51.7%) could not be assigned to any group in the reference database suggesting a high genetic novelty of Svalbard foraminifera. The sieved and unsieved samples resolved comparable communities, sharing 1023 ASVs, comprising over 97% of reads. Our analyses show that the foraminiferal assemblage differs between the localities, with communities distinctly separated between fjord and open sea stations. Each locality was characterized by a specific assemblage, with only a small overlap in the case of open sea areas. Our study demonstrates a clear pattern of the influence of water masses on the structure of foraminiferal communities. The stations situated on the western coast of Svalbard that are strongly influenced by warm and salty Atlantic water (AW) are characterized by much higher diversity than stations in the northern and eastern part, where the impact of AW is less pronounced. This high diversity and specificity of Svalbard foraminifera associated with water mass distribution indicate that the foraminiferal metabarcoding data can be very useful for inferring present and past environmental conditions in the Arctic.

Evolutionary diversification of the autophagy-related ubiquitin-like conjugation systems

Two autophagy-related (ATG) ubiquitin-like conjugation systems, the ATG12 and ATG8 systems, play important roles in macroautophagy. While multiple duplications and losses of the ATG conjugation system proteins are found in different lineages, the extent to which the underlying systems diversified across eukaryotes is not fully understood. Here, in order to understand the evolution of the ATG conjugation systems, we constructed a transcriptome database consisting of 94 eukaryotic species covering major eukaryotic clades and systematically identified ATG conjugation system components. Both ATG10 and the C-terminal glycine of ATG12 are essential for the canonical ubiquitin-like conjugation of ATG12 and ATG5. However, loss of ATG10 or the C-terminal glycine of ATG12 occurred at least 16 times in a wide range of lineages, suggesting that possible covalent-to-non-covalent transition is not limited to the species that we previously reported such as Alveolata and some yeast species. Some species have only the ATG8 system (with conjugation enzymes) or only ATG8 (without conjugation enzymes). More than 10 species have ATG8 homologs without the conserved C-terminal glycine, and Tetrahymena has an ATG8 homolog with a predicted transmembrane domain, which may be able to anchor to the membrane independent of the ATG conjugation systems. We discuss the possibility that the ancestor of the ATG12 and ATG8 systems is more similar to ATG8. Overall, our study offers a whole picture of the evolution and diversity of the ATG conjugation systems among eukaryotes, and provides evidence that functional diversifications of the systems are more common than previously thought.Abbreviations: APEAR: ATG8-PE association region; ATG: autophagy-related; LIR: LC3-interacting region; NEDD8: neural precursor cell expressed, developmentally down-regulated gene 8; PE: phosphatidylethanolamine; SAMP: small archaeal modifier protein; SAR: Stramenopiles, Alveolata, and Rhizaria; SMC: structural maintenance of chromosomes; SUMO: small ubiquitin like modifier; TACK: Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota; UBA: ubiquitin like modifier activating enzyme; UFM: ubiquitin fold modifier; URM: ubiquitin related modifier.

Mitochondrial cytochrome c oxidase subunit I (COI) metabarcoding of Foraminifera communities using taxon-specific primers

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in most marine environments. Molecular methods such as metabarcoding have revealed a high, yet undescribed diversity of Foraminifera. However, so far only one molecular marker, the 18S ribosomal RNA, was available for metabarcoding studies on Foraminifera. Primers that allow amplification of foraminiferal mitochondrial cytochrome oxidase I (COI) and identification of Foraminifera species were recently published. Here we test the performance of these primers for the amplification of whole foraminiferal communities, and compare their performance to that of the highly degenerate LerayXT primers, which amplify the same COI region in a wide range of eukaryotes. We applied metabarcoding to 48 samples taken along three transects spanning a North Sea beach in the Netherlands from dunes to the low tide level, and analysed both sediment samples and meiofauna samples, which contained taxa between 42 µm and 1 mm in body size obtained by decantation from sand samples. We used single-cell metabarcoding (Girard et al., 2022) to generate a COI reference library containing 32 species of Foraminifera, and used this to taxonomically annotate our community metabarcoding data. Our analyses show that the highly degenerate LerayXT primers do not amplify Foraminifera, while the Foraminifera primers are highly Foraminifera- specific, with about 90% of reads assigned to Foraminifera and amplifying taxa from all major groups, i.e., monothalamids, Globothalamea, and Tubothalamea. We identified 176 Foraminifera ASVs and found a change in Foraminifera community composition along the beach transects from high tide to low tide level, and a dominance of single-chambered monothalamid Foraminifera. Our results highlight that COI metabarcoding can be a powerful tool for assessing Foraminiferal communities.

Taxon-rich transcriptomics supports higher-level phylogeny and major evolutionary trends in Foraminifera

Foraminifera, classified in the supergroup Rhizaria, are a common and highly diverse group of mainly marine protists. Despite their evolutionary and ecological importance, only limited genomic data (one partial genome and nine transcriptomic datasets) have been published for this group. Foraminiferal molecular phylogeny is largely based on 18S rRNA gene sequence analysis. However, due to highly variable evolutionary rates of substitution in ribosomal genes plus the existence of intragenomic variation at this locus, the relationships between and within foraminiferal classes remain uncertain. We analyze transcriptomic data from 28 species, adding 19 new species to the previously published dataset, including members of the strongly under-represented class Monothalamea. A phylogenomic reconstruction of Rhizaria, rooted with alveolates and stramenopiles, based on 199 genes and 68 species supports the monophyly of Foraminifera and their sister relationship to Polycystinea. The phylogenomic tree of Foraminifera is very similar to the 18S rRNA tree, with the paraphyletic single-chambered monothalamids giving rise to the multi-chambered Tubothalamea and Globothalamea. Within the Monothalamea, our analyses confirm the monophyly of the giant, deep-sea xenophyophores that branch within clade C and indicate the basal position of monothalamous clades D and E. The multi-chambered Globothalamea are monophyletic and comprise the paraphyletic Textulariida and monophyletic Rotaliida. Our phylogenomic analyses support major evolutionary trends of Foraminifera revealed by ribosomal phylogenies and reinforce their current higher-level classification.

What Can We Learn from -Omics Approaches to Understand Clubroot Disease?

Clubroot is one of the most economically significant diseases worldwide. As a result, many investigations focus on both curing the disease and in-depth molecular studies. Although the first transcriptome dataset for the clubroot disease describing the clubroot disease was published in 2006, many different pathogen-host plant combinations have only recently been investigated and published. Articles presenting -omics data and the clubroot pathogen Plasmodiophora brassicae as well as different host plants were analyzed to summarize the findings in the richness of these datasets. Although genome data for the protist have only recently become available, many effector candidates have been identified, but their functional characterization is incomplete. A better understanding of the life cycle is clearly required to comprehend its function. While only a few proteome studies and metabolome analyses were performed, the majority of studies used microarrays and RNAseq approaches to study transcriptomes. Metabolites, comprising chemical groups like hormones were generally studied in a more targeted manner. Furthermore, functional approaches based on such datasets have been carried out employing mutants, transgenic lines, or ecotypes/cultivars of either Arabidopsis thaliana or other economically important host plants of the Brassica family. This has led to new discoveries of potential genes involved in disease development or in (partial) resistance or tolerance to P. brassicae. The overall contribution of individual experimental setups to a larger picture will be discussed in this review.

Sugar Transporters in Plasmodiophora brassicae: Genome-Wide Identification and Functional Verification

Plasmodiophora brassicae, an obligate intracellular pathogen, can hijack the host’s carbohydrates for survival. When the host plant is infected by P. brassicae, a large amount of soluble sugar accumulates in the roots, especially glucose, which probably facilitates the development of this pathogen. Although a complete glycolytic and tricarboxylic acid cycle (TCA) cycle existed in P. brassicae, very little information about the hexose transport system has been reported. In this study, we screened 17 putative sugar transporters based on information about their typical domains. The structure of these transporters showed a lot of variation compared with that of other organisms, especially the number of transmembrane helices (TMHs). Phylogenetic analysis indicated that these sugar transporters were far from the evolutionary relationship of other organisms and were unique in P. brassicae. The hexose transport activity assay indicated that eight transporters transported glucose or fructose and could restore the growth of yeast strain EBY.VW4000, which was deficient in hexose transport. The expression level of these glucose transporters was significantly upregulated at the late inoculation time when resting spores and galls were developing and a large amount of energy was needed. Our study provides new insights into the mechanism of P. brassicae survival in host cells by hijacking and utilizing the carbohydrates of the host.

Control Strategies of Clubroot Disease Caused by Plasmodiophora brassicae

The clubroot disease caused by the soil-borne pathogen Plasmodiophora brassicae is one of the most important diseases of cruciferous crops worldwide. As with many plant pathogens, the spread is closely related to the cultivation of suitable host plants. In addition, temperature and water availability are crucial determinants for the occurrence and reproduction of clubroot disease. Current global changes are contributing to the widespread incidence of clubroot disease. On the one hand, global trade and high prices are leading to an increase in the cultivation of the host plant rapeseed worldwide. On the other hand, climate change is improving the living conditions of the pathogen P. brassicae in temperate climates and leading to its increased occurrence. Well-known ways to control efficiently this disease include arable farming strategies: growing host plants in wide crop rotations, liming the contaminated soils, and using resistant host plants. Since chemical control of the clubroot disease is not possible or not ecologically compatible, more and more alternative control options are being investigated. In this review, we address the challenges for its control, with a focus on biological control options.

Metagenomics approach for Polymyxa betae genome assembly enables comparative analysis towards deciphering the intracellular parasitic lifestyle of the plasmodiophorids

Genomic knowledge of the tree of life is biased to specific groups of organisms. For example, only six full genomes are currently available in the rhizaria clade. Here, we have applied metagenomic techniques enabling the assembly of the genome of Polymyxa betae (Rhizaria, Plasmodiophorida) RES F41 isolate from unpurified zoospore holobiont and comparison with the A26-41 isolate. Furthermore, the first P. betae mitochondrial genome was assembled. The two P. betae nuclear genomes were highly similar, each with just ~10.2 k predicted protein coding genes, ~3% of which were unique to each isolate. Extending genomic comparisons revealed a greater overlap with Spongospora subterranea than with Plasmodiophora brassicae, including orthologs of the mammalian cation channel sperm-associated proteins, raising some intriguing questions about zoospore physiology. This work validates our metagenomics pipeline for eukaryote genome assembly from unpurified samples and enriches plasmodiophorid genomics; providing the first full annotation of the P. betae genome.

First report of mitochondrial COI in foraminifera and implications for DNA barcoding

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in many marine environments and play a major role in global carbon cycling. Species recognition in Foraminifera is mainly based on morphological characters and nuclear 18S ribosomal RNA barcoding. The 18S rRNA contains variable sequence regions that allow for the identification of most foraminiferal species. Still, some species show limited variability, while others contain high levels of intragenomic polymorphisms, thereby complicating species identification. The use of additional, easily obtainable molecular markers other than 18S rRNA will enable more detailed investigation of evolutionary history, population genetics and speciation in Foraminifera. Here we present the first mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences ("barcodes") of Foraminifera. We applied shotgun sequencing to single foraminiferal specimens, assembled COI, and developed primers that allow amplification of COI in a wide range of foraminiferal species. We obtained COI sequences of 49 specimens from 17 species from the orders Rotaliida and Miliolida. Phylogenetic analysis showed that the COI tree is largely congruent with previously published 18S rRNA phylogenies. Furthermore, species delimitation with ASAP and ABGD algorithms showed that foraminiferal species can be identified based on COI barcodes.

A clubroot pathogen effector targets cruciferous cysteine proteases to suppress plant immunity

Plant pathogen effector proteins are key to pathogen virulence. In susceptible host Brassicas, the clubroot pathogen, Plasmodiophora brassicae, induces the production of nutrient-sink root galls, at the site of infection. Among a list of 32 P. brassiae effector candidates previously reported by our group, we identified SSPbP53 as a putative apoplastic cystatin-like protein highly expressed during the secondary infection. Here we found that SSPbP53 encoding gene is conserved among several P. brassicae pathotypes and that SSPbP53 is an apoplastic protein able to directly interact with and inhibit cruciferous papain-like cysteine proteases (PLCPs), specifically Arabidopsis XYLEM CYSTEINE PEPTIDASE 1 (AtXCP1). The severity of clubroot disease is greatly reduced in the Arabidopsis xcp1 null mutant (AtΔxcp1) after infection with P. brassicae resting spores, indicating that the interaction of P. brassicae SSPbP53 with XCP1 is important to clubroot susceptibility. SSPbP53 is the first cystatin-like effector identified and characterized for a plant pathogenic protist.

Gene Duplications Trace Mitochondria to the Onset of Eukaryote Complexity

The last eukaryote common ancestor (LECA) possessed mitochondria and all key traits that make eukaryotic cells more complex than their prokaryotic ancestors, yet the timing of mitochondrial acquisition and the role of mitochondria in the origin of eukaryote complexity remain debated. Here, we report evidence from gene duplications in LECA indicating an early origin of mitochondria. Among 163,545 duplications in 24,571 gene trees spanning 150 sequenced eukaryotic genomes, we identify 713 gene duplication events that occurred in LECA. LECA's bacterial-derived genes include numerous mitochondrial functions and were duplicated significantly more often than archaeal-derived and eukaryote-specific genes. The surplus of bacterial-derived duplications in LECA most likely reflects the serial copying of genes from the mitochondrial endosymbiont to the archaeal host's chromosomes. Clustering, phylogenies and likelihood ratio tests for 22.4 million genes from 5,655 prokaryotic and 150 eukaryotic genomes reveal no evidence for lineage-specific gene acquisitions in eukaryotes, except from the plastid in the plant lineage. That finding, and the functions of bacterial genes duplicated in LECA, suggests that the bacterial genes in eukaryotes are acquisitions from the mitochondrion, followed by vertical gene evolution and differential loss across eukaryotic lineages, flanked by concomitant lateral gene transfer among prokaryotes. Overall, the data indicate that recurrent gene transfer via the copying of genes from a resident mitochondrial endosymbiont to archaeal host chromosomes preceded the onset of eukaryotic cellular complexity, favoring mitochondria-early over mitochondria-late hypotheses for eukaryote origin.

Non-coding RNAs as emerging targets for crop improvement

Food security is affected by climate change, population growth, as well as abiotic and biotic stresses. Conventional and molecular marker assisted breeding and genetic engineering techniques have been employed extensively for improving resistance to biotic stress in crop plants. Advances in next-generation sequencing technologies have permitted the exploration and identification of parts of the genome that extend beyond the regions with protein coding potential. These non-coding regions of the genome are transcribed to generate many types of non-coding RNAs (ncRNAs). These ncRNAs are involved in the regulation of growth, development, and response to stresses at transcriptional and translational levels. ncRNAs, including long ncRNAs (lncRNAs), small RNAs and circular RNAs have been recognized as important regulators of gene expression in plants and have been suggested to play important roles in plant immunity and adaptation to abiotic and biotic stresses. In this article, we have reviewed the current state of knowledge with respect to lncRNAs and their mechanism(s) of action as well as their regulatory functions, specifically within the context of biotic stresses. Additionally, we have provided insights into how our increased knowledge about lncRNAs may be used to improve crop tolerance to these devastating biotic stresses.

Transcriptome Analysis Identifies Plasmodiophora brassicae Secondary Infection Effector Candidates

Plasmodiophora brassicae (Wor.) is an obligate intracellular plant pathogen affecting Brassicas worldwide. Identification of effector proteins is key to understanding the interaction between P. brassicae and its susceptible host plants. To date, there is very little information available on putative effector proteins secreted by P. brassicae during a secondary infection of susceptible host plants, resulting in root gall production. A bioinformatics pipeline approach to RNA-Seq data from Arabidopsis thaliana (L.) Heynh. root tissues at 17, 20, and 24 d postinoculation (dpi) identified 32 small secreted P. brassicae proteins (SSPbPs) that were highly expressed over this secondary infection time frame. Functional signal peptides were confirmed for 31 of the SSPbPs, supporting the accuracy of the pipeline designed to identify secreted proteins. Expression profiles at 0, 2, 5, 7, 14, 21, and 28 dpi verified the involvement of some of the SSPbPs in secondary infection. For seven of the SSPbPs, a functional domain was identified using Blast2GO and 3D structure analysis and domain functionality was confirmed for SSPbP22, a kinase localized to the cytoplasm and nucleus.

Whole-genome DNA similarity and population structure of Plasmodiophora brassicae strains from Canada

Background

Clubroot is an important disease of brassica crops world-wide. The causal agent, Plasmodiophora brassicae, has been present in Canada for over a century but was first identified on canola (Brassica napus) in Alberta, Canada in 2003. Genetic resistance to clubroot in an adapted canola cultivar has been available since 2009, but resistance breakdown was detected in 2013 and new pathotypes are increasing rapidly. Information on genetic similarity among pathogen populations across Canada could be useful in estimating the genetic variation in pathogen populations, predicting the effect of subsequent selection pressure on changes in the pathogen population over time, and even in identifying the origin of the initial pathogen introduction to canola in Alberta.

Results

The genomic sequences of 43 strains (34 field collections, 9 single-spore isolates) of P. brassicae from Canada, the United States, and China clustered into five clades based on SNP similarity. The strains from Canada separated into four clades, with two containing mostly strains from the Prairies (provinces of Alberta, Saskatchewan, and Manitoba) and two that were mostly from the rest of Canada or the USA. Several strains from China formed a separate clade. More than one pathotype and host were present in all four Canadian clades. The initial pathotypes from canola on the Prairies clustered separately from the pathotypes on canola that could overcome resistance to the initial pathotypes. Similarly, at one site in central Canada where resistance had broken down, about half of the genes differed (based on SNPs) between strains before and after the breakdown.

Conclusion

Clustering based on genome-wide DNA sequencing demonstrated that the initial pathotypes on canola on the Prairies clustered separately from the new virulent pathotypes on the Prairies. Analysis indicated that these ‘new’ pathotypes were likely present in the pathogen population at very low frequency, maintained through balancing selection, and increased rapidly in response to selection from repeated exposure to host resistance.

The New Tree of Eukaryotes

For 15 years, the eukaryote Tree of Life (eToL) has been divided into five to eight major groupings, known as 'supergroups'. However, the tree has been profoundly rearranged during this time. The new eToL results from the widespread application of phylogenomics and numerous discoveries of major lineages of eukaryotes, mostly free-living heterotrophic protists. The evidence that supports the tree has transitioned from a synthesis of molecular phylogenetics and biological characters to purely molecular phylogenetics. Most current supergroups lack defining morphological or cell-biological characteristics, making the supergroup label even more arbitrary than before. Going forward, the combination of traditional culturing with maturing culture-free approaches and phylogenomics should accelerate the process of completing and resolving the eToL at its deepest levels.

Identification of lncRNAs Responsive to Infection by Plasmodiophora brassicae in Clubroot-Susceptible and -Resistant Brassica napus Lines Carrying Resistance Introgressed from Rutabaga

Clubroot disease, caused by Plasmodiophora brassicae Woronin, is a major threat to the production of Brassica' crops. Resistance to different P. brassicae pathotypes has been reported in the A genome, chromosome A08; however, the molecular mechanism of this resistance, especially the involvement of long noncoding RNAs (lncRNAs), is not understood. We have used a strand-specific lncRNA-Seq approach to catalog lncRNAs from the roots of clubroot-susceptible and -resistant Brassica napus lines. In total, 530 differentially expressed (DE) lncRNAs were identified, including 88% of long intergenic RNAs and 11% natural antisense transcripts. Sixteen lncRNAs were identified as target mimics of the microRNAs (miRNAs) and eight were identified as the precursors of miRNAs. KEGG pathway analysis of the DE lncRNAs showed that the cis-regulated target genes mostly belong to the phenylpropanoid biosynthetic pathway (15%) and plant-pathogen interactions (15%) while the transregulated target genes mostly belong to carbon (18%) and amino acid biosynthesis pathway (19%). In all, 24 DE lncRNAs were identified from chromosome A08, which is known to harbor a quantitative trait locus conferring resistance to different P. brassicae pathotypes; however, eight of these lncRNAs showed expression only in the resistant plants. These results could form the basis for future studies aimed at delineating the roles of lncRNAs in plant-microbe interactions.

Genome-wide identification of genes encoding putative secreted E3 ubiquitin ligases and functional characterization of PbRING1 in the biotrophic protist Plasmodiophora brassicae

The E3 ubiquitin ligases are key regulators of protein ubiquitination, which have been shown to be involved in a variety of cellular responses to both biotic and abiotic stresses in eukaryotes. However, the E3 ubiquitin ligase homologues in the soil-borne plant pathogen Plasmodiophora brassicae, the causal agent of clubroot disease of crucifer crops worldwide, remain largely unknown. In this study, we characterized secreted E3 ubiquitin ligases, a group of proteins known to be involved in virulence in many pathogens, in a plasmodiophorid P. brassicae. Genome-wide search in the P. brassicae genome retrieved 139 putative E3 ubiquitin ligases, comprising of 115 RING, 15 HECT, 1 HECT-like, and 8 U-box E3 ubiquitin ligases. Among these E3 ubiquitin ligases, 11 RING, 1 U-box, and 3 HECT were found to harbor signal peptide. Based on published RNA-seq data (Schwelm et al. in Sci Rep 5:11153, 2015), we found that these genes were differentially expressed in distinct life stages including germinating spores, maturing spores, and plasmodia. We characterized one potential secreted E3 ubiquitin ligase, PbRING1 (PBRA_000499). Yeast invertase assay showed that PbRING1 harbors a functional N-terminal signal peptide. PbRING1 also harbors a really interested new gene (RING) domain at its C terminus, which was found to display the E3 ligase activity in vitro. Collectively, this study provides a comprehensive insight into the reservoir of putative secreted E3 ligases in P. brassicae.

Coupling between taxonomic and functional diversity in protistan coastal communities

The study of protistan functional diversity is crucial to understand the dynamics of oceanic ecological processes. We combined the metabarcoding data of various coastal ecosystems and a newly developed trait-based approach to study the link between taxonomic and functional diversity across marine protistan communities of different size-classes. Environmental DNA was extracted and the V4 18S rDNA genomic region was amplified and sequenced. In parallel, we tried to annotate the operational taxonomic units (OTUs) from our metabarcoding dataset to 30 biological traits using published and accessible information on protists. We then developed a method to study trait correlations across protists (i.e. trade-offs) in order to build the best functional groups. Based on the annotated OTUs and our functional groups, we demonstrated that the functional diversity of marine protist communities varied in parallel with their taxonomic diversity. The coupling between functional and taxonomic diversity was conserved across different protist size classes. However, the smallest size-fraction was characterized by wider taxonomic and functional groups diversity, corroborating the idea that nanoplankton and picoplankton are part of a more stable ecological background on which larger protists and metazoans might develop.

Proto-oncogenes in a eukaryotic unicellular organism play essential roles in plasmodial growth in host cells

BACKGROUND

The eukaryotic unicellular protist Plasmodiophora brassicae is an endocellular parasite of cruciferous plants. In host cortical cells, this protist develops a unicellular structure that is termed the plasmodium. The plasmodium is actually a multinucleated cell, which subsequently splits and forms resting spores. The mechanism for the growth of this endocellular parasite in host cell is unclear.

RESULTS

Here, combining de novo genome sequence and transcriptome analysis of strain ZJ-1, we identified top five significant enriched KEGG pathways of differentially expressed genes (DEGs), namely translation, cell growth and death, cell communication, cell motility and cancers. We detected 171 proto-oncogenes from the genome of P. brassicae that were implicated in cancer-related pathways, of which 46 were differential expression genes. Three predicted proto-oncogenes (Pb-Raf1, Pb-Raf2, and Pb-MYB), which showed homology to the human proto-oncogenes Raf and MYB, were specifically activated during the plasmodial growth in host cortical cells, demonstrating their involvement in the multinucleate development stage of the unicellular protist organism. Gene networks involved in the tumorigenic-related signaling transduction pathways and the activation of 12 core genes were identified. Inhibition of phosphoinositol-3-kinase relieved the clubroot symptom and significantly suppressed the development process of plasmodia.

CONCLUSIONS

Proto-oncogene-related regulatory mechanisms play an important role in the plasmodial growth of P. brassicae.

Computational analysis of the Plasmodiophora brassicae genome: mitochondrial sequence description and metabolic pathway database design

Plasmodiophora brassicae is an obligate biotrophic pathogenic protist responsible for clubroot, a root gall disease of Brassicaceae species. In addition to the reference genome of the P. brassicae European e3 isolate and the draft genomes of Canadian or Chinese isolates, we present the genome of eH, a second European isolate. Refinement of the annotation of the eH genome led to the identification of the mitochondrial genome sequence, which was found to be bigger than that of Spongospora subterranea, another plant parasitic Plasmodiophorid phylogenetically related to P. brassicae. New pathways were also predicted, such as those for the synthesis of spermidine, a polyamine up-regulated in clubbed regions of roots. A P. brassicae pathway genome database was created to facilitate the functional study of metabolic pathways in transcriptomics approaches. These available tools can help in our understanding of the regulation of P. brassicae metabolism during infection and in response to diverse constraints.

Multigene phylogeny and cell evolution of chromist infrakingdom Rhizaria: contrasting cell organisation of sister phyla Cercozoa and Retaria

Infrakingdom Rhizaria is one of four major subgroups with distinct cell body plans that comprise eukaryotic kingdom Chromista. Unlike other chromists, Rhizaria are mostly heterotrophic flagellates, amoebae or amoeboflagellates, commonly with reticulose (net-like) or filose (thread-like) feeding pseudopodia; uniquely for eukaryotes, cilia have proximal ciliary transition-zone hub-lattices. They comprise predominantly flagellate phylum Cercozoa and reticulopodial phylum Retaria, whose exact phylogenetic relationship has been uncertain. Given even less clear relationships amongst cercozoan classes, we sequenced partial transcriptomes of seven Cercozoa representing five classes and endomyxan retarian Filoreta marina to establish 187-gene multiprotein phylogenies. Ectoreta (retarian infraphyla Foraminifera, Radiozoa) branch within classical Cercozoa as sister to reticulose Endomyxa. This supports recent transfer of subphylum Endomyxa from Cercozoa to Retaria alongside subphylum Ectoreta which embraces classical retarians where capsules or tests subdivide cells into organelle-containing endoplasm and anastomosing pseudopodial net-like ectoplasm. Cercozoa are more homogeneously filose, often with filose pseudopodia and/or posterior ciliary gliding motility: zooflagellate Helkesimastix and amoeboid Guttulinopsis form a strongly supported clade, order Helkesida. Cercomonads are polyphyletic (Cercomonadida sister to glissomonads; Paracercomonadida deeper). Thecofilosea are a clade, whereas Imbricatea may not be; Sarcomonadea may be paraphyletic. Helkesea and Metromonadea are successively deeper outgroups within cercozoan subphylum Monadofilosa; subphylum Reticulofilosa (paraphyletic on site-heterogeneous trees) branches earliest, Granofilosea before Chlorarachnea. Our multiprotein trees confirm that Rhizaria are sisters of infrakingdom Halvaria (Alveolata, Heterokonta) within chromist subkingdom Harosa (= SAR); they further support holophyly of chromist subkingdom Hacrobia, and are consistent with holophyly of Chromista as sister of kingdom Plantae. Site-heterogeneous rDNA trees group Kraken with environmental DNA clade 'eSarcomonad', not Paracercomonadida. Ectoretan fossil dates evidence ultrarapid episodic stem sequence evolution. We discuss early rhizarian cell evolution and multigene tree coevolutionary patterns, gene-paralogue evidence for chromist monophyly, and integrate this with fossil evidence for the age of Rhizaria and eukaryote cells, and revise rhizarian classification.

Marteilia refringens and Marteilia pararefringens sp. nov. are distinct parasites of bivalves and have different European distributions

Marteilia refringens causes marteiliosis in oysters, mussels and other bivalve molluscs. This parasite previously comprised two species, M. refringens and Marteilia maurini, which were synonymized in 2007 and subsequently referred to as M. refringens ‘O-type’ and ‘M-type’. O-type has caused mass mortalities of the flat oyster Ostrea edulis. We used high throughput sequencing and histology to intensively screen flat oysters and mussels (Mytilus edulis) from the UK, Sweden and Norway for infection by both types and to generate multi-gene datasets to clarify their genetic distinctiveness. Mussels from the UK, Norway and Sweden were more frequently polymerase chain reaction (PCR)-positive for M-type (75/849) than oysters (11/542). We did not detect O-type in any northern European samples, and no histology-confirmed Marteilia-infected oysters were found in the UK, Norway and Sweden, even where co-habiting mussels were infected by the M-type. The two genetic lineages within ‘M. refringens’ are robustly distinguishable at species level. We therefore formally define them as separate species: M. refringens (previously O-type) and Marteilia pararefringens sp. nov. (M-type). We designed and tested new Marteilia-specific PCR primers amplifying from the 3’ end of the 18S rRNA gene through to the 5.8S gene, which specifically amplified the target region from both tissue and environmental samples.

Genotyping of Plasmodiophora brassicae reveals the presence of distinct populations

Background

Plasmodiophora brassicae is a soilborne pathogen of the family Brassicaceae and the causal agent of clubroot disease. In Canada, P. brassicae is now one of the most important constraints to canola (Brassica napus) production, and is managed mainly by the deployment of resistant cultivars. In recent years, however, new strains of the pathogen have emerged that are capable of overcoming host resistance, posing new challenges for disease management. Despite its economic significance, molecular studies of P. brassicae are rare, mainly because this microorganism cannot be cultured outside of its host.

Results

Restriction site-associated DNA sequencing (RADseq) was used to examine the genetic diversity within P. brassicae single-spore and field isolates collected from across Canada. The isolates included individuals that were either capable or incapable of causing disease on clubroot resistant canola cultivars. Over 8750 variants were identified through RADseq. Population analysis indicated that most isolates belonged to one of two distinct populations, corresponding with the ability of isolates to cause disease on resistant cultivars. Within each population, there were low levels of genetic diversity. One thousand and fifty of the genetic variants that distinguished the two populations were nonsynonymous, altering the coding sequences of genes.

Conclusion

The application of RADseq revealed two distinct populations of P. brassicae in Canada, suggesting multiple introductions of the pathogen into the country. The genetic variation found here will be important for future research and monitoring of the pathogen.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4658-1) contains supplementary material, which is available to authorized users.

Not in your usual Top 10: protists that infect plants and algae

Fungi, nematodes and oomycetes belong to the most prominent eukaryotic plant pathogenic organisms. Unicellular organisms from other eukaryotic lineages, commonly addressed as protists, also infect plants. This review provides an introduction to plant pathogenic protists, including algae infecting oomycetes, and their current state of research.

Phytocercomonas venanatans, a New Species of Cercozoa Associated with Chlorotic Streak of Sugarcane

Chlorotic streak is a global disease of commercial sugarcane (Saccharum spp. hybrids). The disease is transmitted by wet soil, water, as well as in diseased planting material. Although first recognized almost 90 years ago and despite significant research effort, the identity of the causal agent has been elusive. Metagenomic high throughput sequencing (HTS) facilitated the discovery of novel protistan ribosomal and nuclear genes in chlorotic streak-infected sugarcane. These sequences suggest a possible causal agent belonging to the order Cercomonadida (Rhizaria, phylum Cercozoa). An organism with morphological features similar to cercomonads (=Cercomonadida) was isolated into pure axenic culture from internal stalk tissues of infected sugarcane. The isolated organism contained DNA sequences identical to those identified in infected plants by HTS. The DNA sequences and the morphology of the organism did not match any known species. Here we present a new genus and species, Phytocercomonas venanatans, which is associated with chlorotic streak of sugarcane. Amplicon sequencing also supports that P. venanatans is associated with this disease. This is the first reported member from Cercomonadida showing a probable pathogenic association with higher plants.

Transcriptome analysis of response to Plasmodiophora brassicae infection in the Arabidopsis shoot and root

BACKGROUND

Clubroot is an important disease caused by the obligate parasite Plasmodiophora brassicae that infects the Brassicaceae. As a soil-borne pathogen, P. brassicae induces the generation of abnormal tissue in the root, resulting in the formation of galls. Root infection negatively affects the uptake of water and nutrients in host plants, severely reducing their growth and productivity. Many studies have emphasized the molecular and physiological effects of the clubroot disease on root tissues. The aim of the present study is to better understand the effect of P. brassicae on the transcriptome of both shoot and root tissues of Arabidopsis thaliana.

RESULTS

Transcriptome profiling using RNA-seq was performed on both shoot and root tissues at 17, 20 and 24 days post inoculation (dpi) of A. thaliana, a model plant host for P. brassicae. The number of differentially expressed genes (DEGs) between infected and uninfected samples was larger in shoot than in root. In both shoot and root, more genes were differentially regulated at 24 dpi than the two earlier time points. Genes that were highly regulated in response to infection in both shoot and root primarily were involved in the metabolism of cell wall compounds, lipids, and shikimate pathway metabolites. Among hormone-related pathways, several jasmonic acid biosynthesis genes were upregulated in both shoot and root tissue. Genes encoding enzymes involved in cell wall modification, biosynthesis of sucrose and starch, and several classes of transcription factors were generally differently regulated in shoot and root.

CONCLUSIONS

These results highlight the similarities and differences in the transcriptomic response of above- and below-ground tissues of the model host Arabidopsis following P. brassicae infection. The main transcriptomic changes in root metabolism during clubroot disease progression were identified. An overview of DEGs in the shoot underlined the physiological changes in above-ground tissues following pathogen establishment and disease progression. This study provides insights into host tissue-specific molecular responses to clubroot development and may have applications in the development of clubroot markers for more effective breeding strategies.

Kingdom Chromista and its eight phyla: a new synthesis emphasising periplastid protein targeting, cytoskeletal and periplastid evolution, and ancient divergences

In 1981 I established kingdom Chromista, distinguished from Plantae because of its more complex chloroplast-associated membrane topology and rigid tubular multipartite ciliary hairs. Plantae originated by converting a cyanobacterium to chloroplasts with Toc/Tic translocons; most evolved cell walls early, thereby losing phagotrophy. Chromists originated by enslaving a phagocytosed red alga, surrounding plastids by two extra membranes, placing them within the endomembrane system, necessitating novel protein import machineries. Early chromists retained phagotrophy, remaining naked and repeatedly reverted to heterotrophy by losing chloroplasts. Therefore, Chromista include secondary phagoheterotrophs (notably ciliates, many dinoflagellates, Opalozoa, Rhizaria, heliozoans) or walled osmotrophs (Pseudofungi, Labyrinthulea), formerly considered protozoa or fungi respectively, plus endoparasites (e.g. Sporozoa) and all chromophyte algae (other dinoflagellates, chromeroids, ochrophytes, haptophytes, cryptophytes). I discuss their origin, evolutionary diversification, and reasons for making chromists one kingdom despite highly divergent cytoskeletons and trophic modes, including improved explanations for periplastid/chloroplast protein targeting, derlin evolution, and ciliary/cytoskeletal diversification. I conjecture that transit-peptide-receptor-mediated 'endocytosis' from periplastid membranes generates periplastid vesicles that fuse with the arguably derlin-translocon-containing periplastid reticulum (putative red algal trans-Golgi network homologue; present in all chromophytes except dinoflagellates). I explain chromist origin from ancestral corticates and neokaryotes, reappraising tertiary symbiogenesis; a chromist cytoskeletal synapomorphy, a bypassing microtubule band dextral to both centrioles, favoured multiple axopodial origins. I revise chromist higher classification by transferring rhizarian subphylum Endomyxa from Cercozoa to Retaria; establishing retarian subphylum Ectoreta for Foraminifera plus Radiozoa, apicomonad subclasses, new dinozoan classes Myzodinea (grouping Colpovora gen. n., Psammosa), Endodinea, Sulcodinea, and subclass Karlodinia; and ranking heterokont Gyrista as phylum not superphylum.

Identification of Plasmodiophora brassicae effectors - A challenging goal

Clubroot is an economically important disease affecting Brassica plants worldwide. Plasmodiophora brassicae is the protist pathogen associated with the disease, and its soil-borne obligate parasitic nature has impeded studies related to its biology and the mechanisms involved in its infection of the plant host. The identification of effector proteins is key to understanding how the pathogen manipulates the plant's immune response and the genes involved in resistance. After more than 140 years studying clubroot and P. brassicae, very little is known about the effectors playing key roles in the infection process and subsequent disease progression. Here we analyze the information available for identified effectors and suggest several features of effector genes that can be used in the search for others. Based on the information presented in this review, we propose a comprehensive bioinformatics pipeline for effector identification and provide a list of the bioinformatics tools available for such.

Single Cell Transcriptomics, Mega-Phylogeny, and the Genetic Basis of Morphological Innovations in Rhizaria

The innovation of the eukaryote cytoskeleton enabled phagocytosis, intracellular transport, and cytokinesis, and is largely responsible for the diversity of morphologies among eukaryotes. Still, the relationship between phenotypic innovations in the cytoskeleton and their underlying genotype is poorly understood. To explore the genetic mechanism of morphological evolution of the eukaryotic cytoskeleton, we provide the first single cell transcriptomes from uncultured, free-living unicellular eukaryotes: the polycystine radiolarian Lithomelissa setosa (Nassellaria) and Sticholonche zanclea (Taxopodida). A phylogenomic approach using 255 genes finds Radiolaria and Foraminifera as separate monophyletic groups (together as Retaria), while Cercozoa is shown to be paraphyletic where Endomyxa is sister to Retaria. Analysis of the genetic components of the cytoskeleton and mapping of the evolution of these on the revised phylogeny of Rhizaria reveal lineage-specific gene duplications and neofunctionalization of α and β tubulin in Retaria, actin in Retaria and Endomyxa, and Arp2/3 complex genes in Chlorarachniophyta. We show how genetic innovations have shaped cytoskeletal structures in Rhizaria, and how single cell transcriptomics can be applied for resolving deep phylogenies and studying gene evolution in uncultured protist species.

Sinks and Sources of Intracellular Nitrate in Gromiids

A substantial nitrate pool is stored within living cells in various benthic marine environments. The fate of this bioavailable nitrogen differs according to the organisms managing the intracellular nitrate (ICN). While some light has been shed on the nitrate carried by diatoms and foraminiferans, no study has so far followed the nitrate kept by gromiids. Gromiids are large protists and their ICN concentration can exceed 1000x the ambient nitrate concentration. In the present study we investigated gromiids from diverse habitats and showed that they contained nitrate at concentrations ranging from 1 to 370 mM. We used ¹⁵N tracer techniques to investigate the source of this ICN, and found that it originated both from active nitrate uptake from the environment and from intracellular production, most likely through bacterial nitrification. Microsensor measurements showed that part of the ICN was denitrified to N₂ when gromiids were exposed to anoxia. Denitrification seemed to be mediated by endobiotic bacteria because antibiotics inhibited denitrification activity. The active uptake of nitrate suggests that ICN plays a role in gromiid physiology and is not merely a consequence of the gromiid hosting a diverse bacterial community. Measurements of aerobic respiration rates and modeling of oxygen consumption by individual gromiid cells suggested that gromiids may occasionally turn anoxic by their own respiration activity and thus need strategies for coping with this self-inflicted anoxia.

The Large Subunit rDNA Sequence of Plasmodiophora brassicae Does not Contain Intra-species Polymorphism

Clubroot disease caused by Plasmodiophora brassicae is one of the most important diseases of cultivated brassicas. P. brassicae occurs in pathotypes which differ in the aggressiveness towards their Brassica host plants. To date no DNA based method to distinguish these pathotypes has been described. In 2011 polymorphism within the 28S rDNA of P. brassicae was reported which potentially could allow to distinguish pathotypes without the need of time-consuming bioassays. However, isolates of P. brassicae from around the world analysed in this study do not show polymorphism in their LSU rDNA sequences. The previously described polymorphism most likely derived from soil inhabiting Cercozoa more specifically Neoheteromita-like glissomonads. Here we correct the LSU rDNA sequence of P. brassicae. By using FISH we demonstrate that our newly generated sequence belongs to the causal agent of clubroot disease.

The compact genome of the plant pathogen Plasmodiophora brassicae is adapted to intracellular interactions with host Brassica spp

BACKGROUND

The protist Plasmodiophora brassicae is a soil-borne pathogen of cruciferous species and the causal agent of clubroot disease of Brassicas including agriculturally important crops such as canola/rapeseed (Brassica napus). P. brassicae has remained an enigmatic plant pathogen and is a rare example of an obligate biotroph that resides entirely inside the host plant cell. The pathogen is the cause of severe yield losses and can render infested fields unsuitable for Brassica crop growth due to the persistence of resting spores in the soil for up to 20 years.

RESULTS

To provide insight into the biology of the pathogen and its interaction with its primary host B. napus, we produced a draft genome of P. brassicae pathotypes 3 and 6 (Pb3 and Pb6) that differ in their host range. Pb3 is highly virulent on B. napus (but also infects other Brassica species) while Pb6 infects only vegetable Brassica crops. Both the Pb3 and Pb6 genomes are highly compact, each with a total size of 24.2 Mb, and contain less than 2 % repetitive DNA. Clustering of genome-wide single nucleotide polymorphisms (SNP) of Pb3, Pb6 and three additional re-sequenced pathotypes (Pb2, Pb5 and Pb8) shows a high degree of correlation of cluster grouping with host range. The Pb3 genome features significant reduction of intergenic space with multiple examples of overlapping untranslated regions (UTRs). Dependency on the host for essential nutrients is evident from the loss of genes for the biosynthesis of thiamine and some amino acids and the presence of a wide range of transport proteins, including some unique to P. brassicae. The annotated genes of Pb3 include those with a potential role in the regulation of the plant growth hormones cytokinin and auxin. The expression profile of Pb3 genes, including putative effectors, during infection and their potential role in manipulation of host defence is discussed.

CONCLUSION

The P. brassicae genome sequence reveals a compact genome, a dependency of the pathogen on its host for some essential nutrients and a potential role in the regulation of host plant cytokinin and auxin. Genome annotation supported by RNA sequencing reveals significant reduction in intergenic space which, in addition to low repeat content, has likely contributed to the P. brassicae compact genome.

Local diversity of heathland Cercozoa explored by in-depth sequencing

Cercozoa are abundant free-living soil protozoa and quantitatively important in soil food webs; yet, targeted high-throughput sequencing (HTS) has not yet been applied to this group. Here we describe the development of a targeted assay to explore Cercozoa using HTS, and we apply this assay to measure Cercozoan community response to drought in a Danish climate manipulation experiment (two sites exposed to artificial drought, two unexposed). Based on a comparison of the hypervariable regions of the 18S ribosomal DNA of 193 named Cercozoa, we concluded that the V4 region is the most suitable for group-specific diversity analysis. We then designed a set of highly specific primers (encompassing ~270 bp) for 454 sequencing. The primers captured all major cercozoan groups; and >95% of the obtained sequences were from Cercozoa. From 443 350 high-quality short reads (>300 bp), we recovered 1585 operational taxonomic units defined by >95% V4 sequence similarity. Taxonomic annotation by phylogeny enabled us to assign >95% of our reads to order level and ~85% to genus level despite the presence of a large, hitherto unknown diversity. Over 40% of the annotated sequences were assigned to Glissomonad genera, whereas the most common individually named genus was the euglyphid Trinema. Cercozoan diversity was largely resilient to drought, although we observed a community composition shift towards fewer testate amoebae.

Nitrate Storage and Dissimilatory Nitrate Reduction by Eukaryotic Microbes

The microbial nitrogen cycle is one of the most complex and environmentally important element cycles on Earth and has long been thought to be mediated exclusively by prokaryotic microbes. Rather recently, it was discovered that certain eukaryotic microbes are able to store nitrate intracellularly and use it for dissimilatory nitrate reduction in the absence of oxygen. The paradigm shift that this entailed is ecologically significant because the eukaryotes in question comprise global players like diatoms, foraminifers, and fungi. This review article provides an unprecedented overview of nitrate storage and dissimilatory nitrate reduction by diverse marine eukaryotes placed into an eco-physiological context. The advantage of intracellular nitrate storage for anaerobic energy conservation in oxygen-depleted habitats is explained and the life style enabled by this metabolic trait is described. A first compilation of intracellular nitrate inventories in various marine sediments is presented, indicating that intracellular nitrate pools vastly exceed porewater nitrate pools. The relative contribution by foraminifers to total sedimentary denitrification is estimated for different marine settings, suggesting that eukaryotes may rival prokaryotes in terms of dissimilatory nitrate reduction. Finally, this review article sketches some evolutionary perspectives of eukaryotic nitrate metabolism and identifies open questions that need to be addressed in future investigations.

Evolutionary Origins of Rhizarian Parasites

The SAR group (Stramenopila, Alveolata, Rhizaria) is one of the largest clades in the tree of eukaryotes and includes a great number of parasitic lineages. Rhizarian parasites are obligate and have devastating effects on commercially important plants and animals but despite this fact, our knowledge of their biology and evolution is limited. Here, we present rhizarian transcriptomes from all major parasitic lineages in order to elucidate their evolutionary relationships using a phylogenomic approach. Our results suggest that Ascetosporea, parasites of marine invertebrates, are sister to the novel clade Apofilosa. The phytomyxean plant parasites branch sister to the vampyrellid algal ectoparasites in the novel clade Phytorhiza. They also show that Ascetosporea + Apofilosa + Retaria + Filosa + Phytorhiza form a monophyletic clade, although the branching pattern within this clade is difficult to resolve and appears to be model-dependent. Our study does not support the monophyly of the rhizarian parasitic lineages (Endomyxa), suggesting independent origins for rhizarian animal and plant parasites.

Whole-genome amplification: a useful approach to characterize new genes in unculturable protozoan parasites such as Bonamia exitiosa

Bonamia exitiosa is an intracellular parasite (Haplosporidia) that has been associated with mass mortalities in oyster populations in the Southern hemisphere. This parasite was recently detected in the Northern hemisphere including Europe. Some representatives of the Bonamia genus have not been well categorized yet due to the lack of genomic information. In the present work, we have applied Whole-Genome Amplification (WGA) technique in order to characterize the actin gene in the unculturable protozoan B. exitiosa. This is the first protein coding gene described in this species. Molecular analysis revealed that B. exitiosa actin is more similar to Bonamia ostreae actin gene-1. Actin phylogeny placed the Bonamia sp. infected oysters in the same clade where the herein described B. exitiosa actin resolved, offering novel information about the classification of the genus. Our results showed that WGA methodology is a promising and valuable technique to be applied to unculturable protozoans whose genomic material is limited.

The Plasmodiophora brassicae genome reveals insights in its life cycle and ancestry of chitin synthases

Plasmodiophora brassicae causes clubroot, a major disease of Brassica oil and vegetable crops worldwide. P. brassicae is a Plasmodiophorid, obligate biotrophic protist in the eukaryotic kingdom of Rhizaria. Here we present the 25.5 Mb genome draft of P. brassicae, developmental stage-specific transcriptomes and a transcriptome of Spongospora subterranea, the Plasmodiophorid causing powdery scab on potato. Like other biotrophic pathogens both Plasmodiophorids are reduced in metabolic pathways. Phytohormones contribute to the gall phenotypes of infected roots. We report a protein (PbGH3) that can modify auxin and jasmonic acid. Plasmodiophorids contain chitin in cell walls of the resilient resting spores. If recognized, chitin can trigger defense responses in plants. Interestingly, chitin-related enzymes of Plasmodiophorids built specific families and the carbohydrate/chitin binding (CBM18) domain is enriched in the Plasmodiophorid secretome. Plasmodiophorids chitin synthases belong to two families, which were present before the split of the eukaryotic Stramenopiles/Alveolates/Rhizaria/Plantae and Metazoa/Fungi/Amoebozoa megagroups, suggesting chitin synthesis to be an ancient feature of eukaryotes. This exemplifies the importance of genomic data from unexplored eukaryotic groups, such as the Plasmodiophorids, to decipher evolutionary relationships and gene diversification of early eukaryotes.