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

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.

Towards high-throughput parallel imaging and single-cell transcriptomics of microbial eukaryotic plankton

Single-cell transcriptomics has the potential to provide novel insights into poorly studied microbial eukaryotes. Although several such technologies are available and benchmarked on mammalian cells, few have been tested on protists. Here, we applied a microarray single-cell sequencing (MASC-seq) technology, that generates microscope images of cells in parallel with capturing their transcriptomes, on three species representing important plankton groups with different cell structures; the ciliate Tetrahymena thermophila, the diatom Phaeodactylum tricornutum, and the dinoflagellate Heterocapsa sp. Both the cell fixation and permeabilization steps were adjusted. For the ciliate and dinoflagellate, the number of transcripts of microarray spots with single cells were significantly higher than for background spots, and the overall expression patterns were correlated with that of bulk RNA, while for the much smaller diatom cells, it was not possible to separate single-cell transcripts from background. The MASC-seq method holds promise for investigating "microbial dark matter", although further optimizations are necessary to increase the signal-to-noise ratio.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

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.

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.

Multigene phylogenetics of euglenids based on single-cell transcriptomics of diverse phagotrophs

Euglenids are a well-known group of single-celled eukaryotes, with phototrophic, osmotrophic and phagotrophic members. Phagotrophs represent most of the phylogenetic diversity of euglenids, and gave rise to the phototrophs and osmotrophs, but their evolutionary relationships are poorly understood. Symbiontids, in contrast, are anaerobes that are alternatively inferred to be derived euglenids, or a separate euglenozoan group. Most phylogenetic studies of euglenids have examined the SSU rDNA only, which is often highly divergent. Also, many phagotrophic euglenids (and symbiontids) are uncultured, restricting collection of other molecular data. We generated transcriptome data for 28 taxa, mostly using a single-cell approach, and conducted the first multigene phylogenetic analyses of euglenids to include phagotrophs and symbiontids. Euglenids are recovered as monophyletic, with symbiontids forming an independent branch within Euglenozoa. Spirocuta, the clade of flexible euglenids that contains both the phototrophs (Euglenophyceae) and osmotrophs (Aphagea), is robustly resolved, with the ploeotid Olkasia as its sister group, forming the new taxon Olkaspira. Ploeotids are paraphyletic, although Ploeotiidae (represented by Ploeotia spp.), Lentomonas, and Keelungia form a robust clade (new taxon Alistosa). Petalomonadida branches robustly as sister to other euglenids in outgroup-rooted analyses. Within Spirocuta, Euglenophyceae is a robust clade that includes Rapaza, and Anisonemia is a well-supported monophyletic group containing Anisonemidae (Anisonema and Dinema spp.), 'Heteronema II' (represented by H. vittatum), and a clade of Neometanema plus Aphagea. Among 'peranemid' phagotrophs, Chasmostoma branches with included Urceolus, and Peranema with the undescribed 'Jenningsia II', while other relationships are weakly supported and consequently the closest sister group to Euglenophyceae remains unresolved. Our results are inconsistent with recent inferences that Entosiphon is the evolutionarily pivotal sister either to other euglenids, or to Spirocuta. At least three transitions between posterior and anterior flagellar gliding occurred in euglenids, with the phylogenetic positions and directions of those transitions remaining ambiguous.

The architecture of the Plasmodiophora brassicae nuclear and mitochondrial genomes

Plasmodiophora brassicae is a soil-borne pathogen that attacks roots of cruciferous plants causing clubroot disease. The pathogen belongs to the Plasmodiophorida order in Phytomyxea. Here we used long-read SMRT technology to clarify the P. brassicae e3 genomic constituents along with comparative and phylogenetic analyses. Twenty contigs representing the nuclear genome and one mitochondrial (mt) contig were generated, together comprising 25.1 Mbp. Thirteen of the 20 nuclear contigs represented chromosomes from telomere to telomere characterized by [TTTTAGGG] sequences. Seven active gene candidates encoding synaptonemal complex-associated and meiotic-related protein homologs were identified, a finding that argues for possible genetic recombination events. The circular mt genome is large (114,663 bp), gene dense and intron rich. It shares high synteny with the mt genome of Spongospora subterranea, except in a unique 12 kb region delimited by shifts in GC content and containing tandem minisatellite- and microsatellite repeats with partially palindromic sequences. De novo annotation identified 32 protein-coding genes, 28 structural RNA genes and 19 ORFs. ORFs predicted in the repeat-rich region showed similarities to diverse organisms suggesting possible evolutionary connections. The data generated here form a refined platform for the next step involving functional analysis, all to clarify the complex biology of P. brassicae.

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.

A single-cell genome perspective on studying intracellular associations in unicellular eukaryotes

Single-cell genomics (SCG) methods provide a unique opportunity to analyse whole genome information at the resolution of an individual cell. While SCG has been extensively used to investigate bacterial and archaeal genomes, the technique has been rarely used to access the genetic makeup of uncultivated microbial eukaryotes. In this regard, the use of SCG can provide a wealth of information; not only do the methods allow exploration of the genome, they can also help elucidate the relationship between the cell and intracellular entities extant in nearly all eukaryotes. SCG enables the study of total eukaryotic cellular DNA, which in turn allows us to better understand the evolutionary history and diversity of life, and the physiological interactions that define complex organisms.

This article is part of a discussion meeting issue ‘Single cell ecology’.

Quartet-Based Computations of Internode Certainty Provide Robust Measures of Phylogenetic Incongruence

Incongruence, or topological conflict, is prevalent in genome-scale data sets. Internode certainty (IC) and related measures were recently introduced to explicitly quantify the level of incongruence of a given internal branch among a set of phylogenetic trees and complement regular branch support measures (e.g., bootstrap, posterior probability) that instead assess the statistical confidence of inference. Since most phylogenomic studies contain data partitions (e.g., genes) with missing taxa and IC scores stem from the frequencies of bipartitions (or splits) on a set of trees, IC score calculation typically requires adjusting the frequencies of bipartitions from these partial gene trees. However, when the proportion of missing taxa is high, the scores yielded by current approaches that adjust bipartition frequencies in partial gene trees differ substantially from each other and tend to be overestimates. To overcome these issues, we developed three new IC measures based on the frequencies of quartets, which naturally apply to both complete and partial trees. Comparison of our new quartet-based measures to previous bipartition-based measures on simulated data shows that: (1) on complete data sets, both quartet-based and bipartition-based measures yield very similar IC scores; (2) IC scores of quartet-based measures on a given data set with and without missing taxa are more similar than the scores of bipartition-based measures; and (3) quartet-based measures are more robust to the absence of phylogenetic signal and errors in phylogenetic inference than bipartition-based measures. Additionally, the analysis of an empirical mammalian phylogenomic data set using our quartet-based measures reveals the presence of substantial levels of incongruence for numerous internal branches. An efficient open-source implementation of these quartet-based measures is freely available in the program QuartetScores (https://github.com/lutteropp/QuartetScores).

Comparative and Functional Algal Genomics

Over 100 whole-genome sequences from algae are published or soon to be published. The rapidly increasing availability of these fundamental resources is changing how we understand one of the most diverse, complex, and understudied groups of photosynthetic eukaryotes. Genome sequences provide a window into the functional potential of individual algae, with phylogenomics and functional genomics as tools for contextualizing and transferring knowledge from reference organisms into less well-characterized systems. Remarkably, over half of the proteins encoded by algal genomes are of unknown function, highlighting the volume of functional capabilities yet to be discovered. In this review, we provide an overview of publicly available algal genomes, their associated protein inventories, and their quality, with a summary of the statuses of protein function understanding and predictions.

Inferring Ancient Relationships with Genomic Data: A Commentary on Current Practices

Contemporary phylogeneticists enjoy an embarrassment of riches, not only in the volumes of data now available, but also in the diversity of bioinformatic tools for handling these data. Here, I discuss a subset of these tools I consider well-suited to the task of inferring ancient relationships with coding sequence data in particular, encompassing data generation, orthology assignment, alignment and gene tree inference, supermatrix construction, and analysis under the best-fitting models applicable to large-scale datasets. Throughout, I compare and critique methods, considering both their theoretical principles and the details of their implementation, and offering practical tips on usage where appropriate. I also entertain different motivations for analyzing what are almost always originally DNA sequence data as codons, amino acids, and higher-order recodings. Although presented in a linear order, I see value in using the diversity of tools available to us to assess the sensitivity of clades of biological interest to different gene and taxon sets and analytical modes, which can be an indication of the presence of systematic error, of which a few forms remain poorly controlled by even the best available inference methods.

Phylogenomics and Morphological Reconstruction of Arcellinida Testate Amoebae Highlight Diversity of Microbial Eukaryotes in the Neoproterozoic

Life was microbial for the majority of Earth's history, but as very few microbial lineages leave a fossil record, the Precambrian evolution of life remains shrouded in mystery. Shelled (testate) amoebae stand out as an exception with rich documented diversity in the Neoproterozoic as vase-shaped microfossils (VSMs). While there is general consensus that most of these can be attributed to the Arcellinida lineage in Amoebozoa, it is still unclear whether they can be used as key fossils for interpretation of early eukaryotic evolution. Here, we present a well-resolved phylogenomic reconstruction based on 250 genes, obtained using single-cell transcriptomic techniques from a representative selection of 19 Arcellinid testate amoeba taxa. The robust phylogenetic framework enables deeper interpretations of evolution in this lineage and demanded an updated classification of the group. Additionally, we performed reconstruction of ancestral morphologies, yielding hypothetical ancestors remarkably similar to existing Neoproterozoic VSMs. We demonstrate that major lineages of testate amoebae were already diversified before the Sturtian glaciation (720 mya), supporting the hypothesis that massive eukaryotic diversification took place in the early Neoproterozoic and congruent with the interpretation that VSM are arcellinid testate amoebae.

Hemimastigophora is a novel supra-kingdom-level lineage of eukaryotes

Almost all eukaryote life forms have now been placed within one of five to eight supra-kingdom-level groups using molecular phylogenetics^1-4. The 'phylum' Hemimastigophora is probably the most distinctive morphologically defined lineage that still awaits such a phylogenetic assignment. First observed in the nineteenth century, hemimastigotes are free-living predatory protists with two rows of flagella and a unique cell architecture^5-7; to our knowledge, no molecular sequence data or cultures are currently available for this group. Here we report phylogenomic analyses based on high-coverage, cultivation-independent transcriptomics that place Hemimastigophora outside of all established eukaryote supergroups. They instead comprise an independent supra-kingdom-level lineage that most likely forms a sister clade to the 'Diaphoretickes' half of eukaryote diversity (that is, the 'stramenopiles, alveolates and Rhizaria' supergroup (Sar), Archaeplastida and Cryptista, as well as other major groups). The previous ranking of Hemimastigophora as a phylum understates the evolutionary distinctiveness of this group, which has considerable importance for investigations into the deep-level evolutionary history of eukaryotic life-ranging from understanding the origins of fundamental cell systems to placing the root of the tree. We have also established the first culture of a hemimastigote (Hemimastix kukwesjijk sp. nov.), which will facilitate future genomic and cell-biological investigations into eukaryote evolution and the last eukaryotic common ancestor.

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.

Rhizarian 'Novel Clade 10' Revealed as Abundant and Diverse Planktonic and Terrestrial Flagellates, including Aquavolon n. gen

Rhizarian ‘Novel Clade 10’ (NC10) is frequently detected by 18S rRNA gene sequencing studies in freshwater planktonic samples. We describe a new genus and two species of eukaryovorous biflagellate protists, Aquavolon hoantrani n. gen. n. sp. and A. dientrani n. gen. n. sp., which represent the first morphologically characterized members of NC10, here named Aquavolonida ord. nov. The slightly metabolic cells possess naked heterodynamic flagella, whose kinetosomes lie at a right angle to each other and are connected by at least one fibril. Unlike their closest known relative Tremula longifila, they rotate around their longitudinal axis when swimming and only very rarely glide on surfaces. Screening of a wide range of environmental DNA extractions with lineage‐specific PCR primers reveals that Aquavolonida consists of a large radiation of protists, which are most diversified in freshwater planktonic habitats and as yet undetected in marine environments. Earlier‐branching lineages in Aquavolonida include less frequently detected organisms from soils and freshwater sediments. The 18S rRNA gene phylogeny suggests that Aquavolonida forms a common evolutionary lineage with tremulids and uncharacterized ‘Novel Clade 12’, which likely represents one of the deepest lineages in the Rhizaria, separate from Cercozoa (Filosa), Endomyxa, and Retaria.

Encystation: the most prevalent and underinvestigated differentiation pathway of eukaryotes

Not long ago, protists were considered one of four eukaryote kingdoms, but recent gene-based phylogenies show that they contribute to all nine eukaryote subdomains. The former kingdoms of animals, plants and fungi are now relegated to lower ranks within subdomains. Most unicellular protists respond to adverse conditions by differentiating into dormant walled cysts. As cysts, they survive long periods of starvation, drought and other environmental threats, only to re-emerge when conditions improve. For protists pathogens, the resilience of their cysts can prevent successful treatment or eradication of the disease. In this context, effort has been directed towards understanding the molecular mechanisms that control encystation. We here firstly summarize the prevalence of encystation across protists and next focus on Amoebozoa, where most of the health-related issues occur. We review current data on processes and genes involved in encystation of the obligate parasite Entamoeba histolytica and the opportunistic pathogen Acanthamoeba. We show how the cAMP-mediated signalling pathway that controls spore and stalk cell encapsulation in Dictyostelium fruiting bodies could be retraced to a stress-induced pathway controlling encystation in solitary Amoebozoa. We highlight the conservation and prevalence of cAMP signalling genes in Amoebozoan genomes and the suprisingly large and varied repertoire of proteins for sensing and processing environmental signals in individual species.

Alpha-actinin of the chlorarchiniophyte Bigelowiella natans

The genome of the chlorarchiniophyte Bigelowiella natans codes for a protein annotated as an α-actinin-like protein. Analysis of the primary sequence indicate that this protein has the same domain structure as other α-actinins, a N-terminal actin-binding domain and a C-terminal calmodulin-like domain. These two domains are connected by a short rod domain, albeit long enough to form a single spectrin repeat. To analyse the functional properties of this protein, the full-length protein as well as the separate domains were cloned and isolated. Characerisation showed that the protein is capable of cross-linking actin filaments into dense bundles, probably due to dimer formation. Similar to human α-actinin, calcium-binding occurs to the most N-terminal EF-hand motif in the calmodulin-like C-terminal domain. The results indicate that this Bigelowiella protein is a proper α-actinin, with all common characteristics of a typical α-actinin.

The immunophilin repertoire of Plasmodiophora brassicae and functional analysis of PbCYP3 cyclophilin

Plasmodiophora brassicae is a soil-borne pathogen that belongs to Rhizaria, an almost unexplored eukaryotic organism group. This pathogen requires a living host for growth and multiplication, which makes molecular analysis further complicated. To broaden our understanding of a plasmodiophorid such as P. brassicae, we here chose to study immunophilins, a group of proteins known to have various cellular functions, including involvement in plant defense and pathogen virulence. Searches in the P. brassicae genome resulted in 20 putative immunophilins comprising of 11 cyclophilins (CYPs), 7 FK506-binding proteins (FKBPs) and 2 parvulin-like proteins. RNAseq data showed that immunophilins were differentially regulated in enriched life stages such as germinating spores, maturing spores, and plasmodia, and infected Brassica hosts (B. rapa, B. napus and B. oleracea). PbCYP3 was highly induced in all studied life stages and during infection of all three Brassica hosts, and hence was selected for further analysis. PbCYP3 was heterologously expressed in Magnaporthe oryzae gene-inactivated ΔCyp1 strain. The new strain ΔCyp1+ overexpressing PbCYP3 showed increased virulence on rice compared to the ΔCyp1 strain. These results suggest that the predicted immunophilins and particularly PbCYP3 are activated during plant infection. M. oryzae is a well-studied fungal pathogen and could be a valuable tool for future functional studies of P. brassicae genes, particularly elucidating their role during various infection phases.