Molecular diversity of the syndinean genus Euduboscquella based on single-cell PCR analysis

Comparative assessment of the intragenomic variations of dinoflagellate Tripos species through single-cell sequencing

Tripos is a large dinoflagellate genus widely distributed in the world's oceans. Morphology-based species identification is inconclusive due to high morphological intraspecific variability. Metabarcoding analysis has been demonstrated to be effective for species identification and tracking their spatiotemporal dynamics. However, accumulating evidence suggests high levels of intragenomic variations (IGVs) are common in many algae, leading to concerns about overinterpretation of molecular diversity in metabarcoding studies. In this project, we evaluated and compared IGVs in Tripos species by conducting the first high-throughput sequencing (HTS) of 18S rDNA V4 of Tripos single cells. High numbers of haplotypes (19-172) were identified in each of the 30 Tripos cells. Each cell contained one dominant haplotype with high relative abundance and many haplotypes with lower abundances. Thus, the presence of multiple minor haplotypes substantially overestimate the molecular diversity identified in metabarcoding analysis, which encompass not only interspecific and intraspecific diversities, but high levels of IGVs.

Compositional Graphical Lasso Resolves the Impact of Parasitic Infection on Gut Microbial Interaction Networks in a Zebrafish Model

Understanding how microbes interact with each other is key to revealing the underlying role that microorganisms play in the host or environment and to identifying microorganisms as an agent that can potentially alter the host or environment. For example, understanding how the microbial interactions associate with parasitic infection can help resolve potential drug or diagnostic test for parasitic infection. To unravel the microbial interactions, existing tools often rely on graphical models to infer the conditional dependence of microbial abundances to represent their interactions. However, current methods do not simultaneously account for the discreteness, compositionality, and heterogeneity inherent to microbiome data. Thus, we build a new approach called "compositional graphical lasso" upon existing tools by incorporating the above characteristics into the graphical model explicitly. We illustrate the advantage of compositional graphical lasso over current methods under a variety of simulation scenarios and on a benchmark study, the Tara Oceans Project. Moreover, we present our results from the analysis of a dataset from the Zebrafish Parasite Infection Study, which aims to gain insight into how the gut microbiome and parasite burden covary during infection, thus uncovering novel putative methods of disrupting parasite success. Our approach identifies changes in interaction degree between infected and uninfected individuals for three taxa, Photobacterium, Gemmobacter, and Paucibacter, which are inversely predicted by other methods. Further investigation of these method-specific taxa interaction changes reveals their biological plausibility. In particular, we speculate on the potential pathobiotic roles of Photobacterium and Gemmobacter in the zebrafish gut, and the potential probiotic role of Paucibacter. Collectively, our analyses demonstrate that compositional graphical lasso provides a powerful means of accurately resolving interactions between microbiota and can thus drive novel biological discovery.

Contrasting diversity patterns of prokaryotes and protists over time and depth at the San-Pedro Ocean Time series

Community dynamics are central in microbial ecology, yet we lack studies comparing diversity patterns among marine protists and prokaryotes over depth and multiple years. Here, we characterized microbes at the San-Pedro Ocean Time series (2005-2018), using SSU rRNA gene sequencing from two size fractions (0.2-1 and 1-80 μm), with a universal primer set that amplifies from both prokaryotes and eukaryotes, allowing direct comparisons of diversity patterns in a single set of analyses. The 16S + 18S rRNA gene composition in the small size fraction was mostly prokaryotic (>92%) as expected, but the large size fraction unexpectedly contained 46-93% prokaryotic 16S rRNA genes. Prokaryotes and protists showed opposite vertical diversity patterns; prokaryotic diversity peaked at mid-depth, protistan diversity at the surface. Temporal beta-diversity patterns indicated prokaryote communities were much more stable than protists. Although the prokaryotic communities changed monthly, the average community stayed remarkably steady over 14 years, showing high resilience. Additionally, particle-associated prokaryotes were more diverse than smaller free-living ones, especially at deeper depths, contributed unexpectedly by abundant and diverse SAR11 clade II. Eukaryotic diversity was strongly correlated with the diversity of particle-associated prokaryotes but not free-living ones, reflecting that physical associations result in the strongest interactions, including symbioses, parasitism, and decomposer relationships.

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

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

The Windblown: Possible Explanations for Dinophyte DNA in Forest Soils

Dinophytes are widely distributed in marine- and fresh-waters, but have yet to be conclusively documented in terrestrial environments. Here, we evaluated the presence of these protists from an environmental DNA metabarcoding dataset of Neotropical rainforest soils. Using a phylogenetic placement approach with a reference alignment and tree, we showed that the numerous sequencing reads that were phylogenetically placed as dinophytes did not correlate with taxonomic assignment, environmental preference, nutritional mode, or dormancy. All the dinophytes in the soils are rather windblown dispersal units of aquatic species and are not biologically active residents of terrestrial environments.

An assessment on the intrapopulational and intraindividual genetic diversity in LSU rDNA in the harmful algal blooms-forming dinoflagellate Margalefidinium (= Cochlodinium) fulvescens based on clonal cultures and bloom samples from Jiaozhou Bay, China

Large subunit ribosomal DNA (LSU rDNA) sequences have been increasingly used to infer the phylogeny and species identity of organisms, a few previous studies, however, have observed high intraspecific and even intraindividual variability in LSU rDNA in some dinoflagellate species due to, assumably, large copy numbers of rDNA in dinoflagellates. Since the copy number of LSU rDNA varies tremendously among dinoflagellate species, the intraspecific and intraindividual diversity for a species of particular interest thus needs to be investigated individually. As a toxic and HABs-forming dinoflagellate, Margalefidinium (= Cochlodinium) fulvescens has been observed to approach blooming density in Jiaozhou Bay, China since 2015 after numerous blooms having been reported from other countries. In trying to identify the source of this newly observed HABs-forming species in China by sequencing the LSU rDNA for both field samples and clonal cultures, we noticed and thus further investigated high intrapopulational and intraindividual genetic diversities of the dinoflagellate. The D1-D6 region of the LSU rDNA (1,435 bases) was amplified from 7 field samples (pooled cells) and 11 clonal cultures, cloned, sequenced, and analyzed phylogenetically for 2,341 sequences obtained. All the numbers of sequences obtained from each clonal culture were far less than the estimated rDNA copy number in M. fulvescens. In the clone library, only one unique sequence was contained in all samples as the most dominant sequence. We found high intrapopulational and intraindividual genetic diversity in M. fulvescens as reflected in the number of polymorphic sites and unique sequences in the clone library for different field samples and clonal cultures in comparison to other species. The mean number of nucleotide differences of each sequence from different field samples and clonal cultures were 6.43 and 4.42 bases, respectively, with the highest being 132 bases, nearly 10%. The sequences with highest variability may be easily annotated as different species if they were obtained from environmental genomic studies because sequence-based species identification in meta-barcoding studies often use "97% identity" threshold. Based on that the mean and overall intrapopulational genetic diversity calculated for 7 field samples was equivalent to the mean and overall intraindividual variability for 11 clonal cultures in indices of genetic diversity, together with the result of AMOVA analysis, we infer that the variability within individual cells (i.e. variability among LSU rDNA polymorphic copies) caused both the intraindividual and intrapopulational genetic diversities observed in the M. fulvescens population, and a higher interpopulational diversity may exist among different geographic populations. The results provide an insightful basis for such a comprehensive interpopulational comparison and important implications for identifying species and establishing new taxa based on the similarity comparison to reference sequences deposited in databases.

Temporal Variability and Ecological Interactions of Parasitic Marine Syndiniales in Coastal Protist Communities

Syndiniales are a ubiquitous group of protist parasites that infect and kill a wide range of hosts, including harmful bloom-forming dinoflagellates. Despite the importance of parasitism as an agent of plankton mortality, parasite-host dynamics remain poorly understood, especially over time, hindering the inclusion of parasitism in food web and ecosystem models. For a full year in the Skidaway River Estuary (Georgia), we employed weekly 18S rRNA sampling and co-occurrence network analysis to characterize temporal parasite-host infection dynamics of Syndiniales. Over the year, Syndiniales exhibited strong temporal variability, with higher relative abundance from June to October (7 to 28%) than other months in the year (0.01% to 6%). Nonmetric dimensional scaling of Syndiniales composition revealed tight clustering in June to October that coincided with elevated temperatures (23 to 31°C), though in general, abiotic factors poorly explained composition (canonical correspondence analysis [CCA] and partial least-squares [PLS]) and were less important in the network than biotic relationships. Syndiniales amplicon sequence variants (ASVs) were well represented in the co-occurrence network (20% of edges) and had significant positive associations (Spearman r > 0.7), inferred to be putative parasite-host relationships, with known dinoflagellate hosts (e.g., Akashiwo and Gymnodinium) and other protist groups (e.g., ciliates, radiolarians, and diatoms). Positive associations rarely involved a single Syndiniales and dinoflagellate species, implying flexible parasite-host infection dynamics. These findings provide insight into the temporal dynamics of Syndiniales over a full year and reinforce the importance of single-celled parasites in driving plankton population dynamics. Further empirical work is needed to confirm network interactions and to incorporate parasitism within the context of ecosystem models.IMPORTANCE Protist parasites in the marine alveolate group, Syndiniales, have been observed within infected plankton host cells for decades, and recently, global-scale efforts (Tara Ocean exploration) have confirmed their importance within microbial communities. Yet, protist parasites remain enigmatic, particularly with respect to their temporal dynamics and parasite-host interactions. We employed weekly 18S amplicon surveys over a full year in a coastal estuary, revealing strong temporal shifts in Syndiniales parasites, with highest relative abundance during warmer summer to fall months. Though influenced by temperature, Syndiniales population dynamics were also driven by a high frequency of biological interactions with other protist groups, as determined through co-occurrence network analysis. Parasitic interactions implied by the network highlighted a range of confirmed (dinoflagellates) and putative (diatoms) interactions and suggests parasites may be less selective in their preferred hosts. Understanding parasite-host dynamics over space and time will improve our ability to include parasitism as a loss term in microbial food web models.

Cryptic species in the parasitic Amoebophrya species complex revealed by a polyphasic approach

As critical primary producers and recyclers of organic matter, the diversity of marine protists has been extensively explored by high-throughput barcode sequencing. However, classification of short metabarcoding sequences into traditional taxonomic units is not trivial, especially for lineages mainly known by their genetic fingerprints. This is the case for the widespread Amoebophrya ceratii species complex, parasites of their dinoflagellate congeners. We used genetic and phenotypic characters, applied to 119 Amoebophrya individuals sampled from the same geographic area, to construct practical guidelines for species delineation that could be applied in DNA/RNA based diversity analyses. Based on the internal transcribed spacer (ITS) regions, ITS2 compensatory base changes (CBC) and genome k-mer comparisons, we unambiguously defined eight cryptic species among closely related ribotypes that differed by less than 97% sequence identity in their SSU rDNA. We then followed the genetic signatures of these parasitic species during a three-year survey of Alexandrium minutum blooms. We showed that these cryptic Amoebophrya species co-occurred and shared the same ecological niche. We also observed a maximal ecological fitness for parasites having narrow to intermediate host ranges, reflecting a high cost for infecting a broader host range. This study suggests that a complete taxonomic revision of these parasitic dinoflagellates is long overdue to understand their diversity and ecological role in the marine plankton.

Protist Interactions and Community Structure During Early Autumn in the Kerguelen Region (Southern Ocean)

This study investigated protist community composition and biotic interactions focusing on microplankton at four distinct sites around the Kerguelen Islands (Southern Ocean) after the summer phytoplankton bloom. Protist diversity in different size fractions, sampled with Niskin bottles and plankton nets, was assessed by sequencing of the V4 18S rDNA region. Combining different approaches, i.e. sequencing of different plankton size fractions, and isolation and sequencing of single cells, provided new insights into microbial interactions in protist communities. The communities displayed high variability, including short-term fluctuations in relative abundance of large protists (>35μm) highlighted by the plankton net samples. Size fractionation of protist communities showed high concentrations of free Syndiniales spores but relatively few Syndiniales associated with microplankton, suggesting low parasitic infection in early autumn. Co-variance network analyses and sequencing of individually isolated single cells highlighted the important role of Rhizaria as consumers of a wide range of different diatom taxa. The data also raised the hypothesis that different Syndiniales clades might be directly or indirectly associated with some diatom genera, thus suggesting a potentially wider host range of these parasites than has been previously reported. These associations and the potential impact on carbon fluxes are discussed.

A precedented nuclear genetic code with all three termination codons reassigned as sense codons in the syndinean Amoebophrya sp. ex Karlodinium veneficum

Amoebophrya is part of an enigmatic, diverse, and ubiquitous marine alveolate lineage known almost entirely from anonymous environmental sequencing. Two cultured Amoebophrya strains grown on core dinoflagellate hosts were used for transcriptome sequencing. BLASTx using different genetic codes suggests that Amoebophyra sp. ex Karlodinium veneficum uses the three typical stop codons (UAA, UAG, and UGA) to encode amino acids. When UAA and UAG are translated as glutamine about half of the alignments have better BLASTx scores, and when UGA is translated as tryptophan one fifth have better scores. However, the sole stop codon appears to be UGA based on conserved genes, suggesting contingent translation of UGA. Neither host sequences, nor sequences from the second strain, Amoebophrya sp. ex Akashiwo sanguinea had similar results in BLASTx searches. A genome survey of Amoebophyra sp. ex K. veneficum showed no evidence for transcript editing aside from mitochondrial transcripts. The dynein heavy chain (DHC) gene family was surveyed and of 14 transcripts only two did not use UAA, UAG, or UGA in a coding context. Overall the transcriptome displayed strong bias for A or U in third codon positions, while the tRNA genome survey showed bias against codons ending in U, particularly for amino acids with two codons ending in either C or U. Together these clues suggest contingent translation mechanisms in Amoebophyra sp. ex K. veneficum and a phylogenetically distinct instance of genetic code modification.

Further analyses of variation of ribosome DNA copy number and polymorphism in ciliates provide insights relevant to studies of both molecular ecology and phylogeny

Sequence-based approaches, such as analyses of ribosome DNA (rDNA) clone libraries and high-throughput amplicon sequencing, have been used extensively to infer evolutionary relationships and elucidate the biodiversity in microbial communities. However, recent studies demonstrate both rDNA copy number variation and intra-individual (intra-genomic) sequence variation in many organisms, which challenges the application of the rDNA-based surveys. In ciliates, an ecologically important clade of microbial eukaryotes, rDNA copy number and sequence variation are rarely studied. In the present study, we estimate the intraindividual small subunit rDNA (SSU rDNA) copy number and sequence variation in a wide range of taxa covering nine classes and 18 orders of the phylum Ciliophora. Our studies reveal that: (i) intra-individual sequence variation of SSU rDNA is ubiquitous in all groups of ciliates detected and the polymorphic level varies among taxa; (ii) there is a most common version of SSU rDNA sequence in each cell that is highly predominant and may represent the germline micronuclear template; (iii) compared with the most common version, other variant sequences differ in only 1-3 nucleotides, likely generated during macronuclear (somatic) amplification; (iv) the intra-cell sequence variation is unlikely to impact phylogenetic analyses; (v) the rDNA copy number in ciliates is highly variable, ranging from 10³ to 10⁶, with the highest record in Stentor roeselii. Overall, these analyses indicate the need for careful consideration of SSU rDNA variation in analyses of the role of ciliates in ecosystems.

Picoeukaryotic Diversity And Activity in the Northwestern Pacific Ocean Based on rDNA and rRNA High-Throughput Sequencing

Picoeukaryotes play an important role in the biogenic element cycle and energy flow in oligotrophic ecosystems. However, their biodiversity and activity are poorly studied in open ocean systems, such as the northwestern Pacific Ocean, which is characterized by a complex hydrological setting. Here, we investigated the diversity and activity of picoeukaryotes in the northwestern Pacific Ocean using high-throughput sequencing targeting the V9 region of 18S rDNA and rRNA. Our results showed that the DNA picoeukaryotic communities were mainly represented by Mamiellophyceae, MAST, MALV-II, Spirotrichea, Prymnesiophyceae, and MALV-I (69.33% of the total DNA reads), and the RNA communities were dominated by Spirotrichea, Mamiellophyceae, MAST, Pelagophyceae, and MALV-II (67.46% of the total RNA reads). The number of operational taxonomic units (OTUs) was significantly affected by temperature and salinity, and was decreased with the increasing nutrient concentration both in the DNA and RNA surveys. Significant differences were observed in the community composition between DNA-based and RNA-based molecular approaches, and these differences were mainly attributed to Mamiellophyceae, Spirotrichea, and Pelagophyceae. The RNA: DNA ratio was used as a proxy for relative metabolic activity of the individual OTUs. We found that the relative metabolic activities of Mamiellophyceae, Spirotrichea, and Pelagophyceae species in the northwestern Pacific Ocean were highly affected by the nutrient concentration, i.e., the NO₃ + NO₂ and SiO₃ concentration. Overall, our study shed light on picoeukaryotic diversity and distribution in the northwestern Pacific Ocean and revealed the correlation between the diversity, relative metabolic activities of marine picoeukaryotes, and the environmental factors.

Do All Dinoflagellates have an Extranuclear Spindle?

The syndinean dinoflagellates are a diverse assemblage of alveolate endoparasites that branch basal to the core dinoflagellates. Because of their phylogenetic position, the syndineans are considered key model microorganisms in understanding early evolution in the dinoflagellates. Closed mitosis with an extranuclear spindle that traverses the nucleus in cytoplasmic grooves or tunnels is viewed as one of the morphological features shared by syndinean and core dinoflagellates. Here we describe nuclear morphology and mitosis in the syndinean dinoflagellate Amoebophrya sp. from Akashiwo sanguinea, a member of the A. ceratii complex, as revealed by protargol silver impregnation, DNA specific fluorochromes, and transmission electron microscopy. Our observations show that not all species classified as dinoflagellates have an extranuclear spindle. In Amoebophrya sp. from A. sanguinea, an extranuclear microtubule cylinder located in a depression in the nuclear surface during interphase moves into the nucleoplasm via sequential membrane fusion events and develops into an entirely intranuclear spindle. Results suggest that the intranuclear spindle of Amoebophrya spp. may have evolved from an ancestral extranuclear spindle and indicate the need for taxonomic revision of the Amoebophryidae.

Euduboscquella costata n. sp. (Dinoflagellata, Syndinea), an Intracellular Parasite of the Ciliate Schmidingerella arcuata: Morphology, Molecular Phylogeny, Life Cycle, Prevalence, and Infection Intensity

The syndinean dinoflagellate Euduboscquella costata n. sp., an intracellular parasite of the tintinnid ciliate Schmidingerella arcuata, was discovered from Korean coastal water in November of 2013. Euduboscquella costata parasitized in about 62% of the host population, with infection intensity (= number of trophonts in a single host cell) ranging from 1 to 8. Based on morphology and nuclear 18S ribosomal RNA gene sequences, the parasite is new to science. Euduboscquella costata n. sp. had an infection cycle typical of the genus, but had morphological and developmental features that distinguished it from congeneric species. These features include: (1) episome of the trophont with 25-40 grooves converging toward the center of the shield; (2) a narrow, funnel-shaped lamina pharyngea extending from the margin of the episomal shield to the nucleus; (3) persistence of grooves during extracellular development (sporogenesis); (4) a single food vacuole during sporogenesis; (5) separation of sporocytes early in sporogenesis, regardless of type of spore formed; and (6) dinospore size (ca. 14 μm in length) and shape (bulbous episome with narrower, tapering hyposome). After sporogenesis, E. costata produced four different types of spore that showed completely identical 18S rRNA gene sequences. The gene sequence was completely identical with a previously reported population, Euduboscquella sp. ex S. arcuata, from Assawoman Bay, USA, indicating that the two populations are likely conspecific. Favella ehrenbergii, a widely recorded tintinnid known to host Euduboscquella spp., co-occurred with S. arcuata, but was not infected by E. costata in field samples or during short-term, cross-infection experiments.

Significance of plankton community structure and nutrient availability for the control of dinoflagellate blooms by parasites: a modeling approach

Dinoflagellate blooms are frequently observed under temporary eutrophication of coastal waters after heavy rains. Growth of these opportunistic microalgae is believed to be promoted by sudden input of nutrients and the absence or inefficiency of their natural enemies, such as grazers and parasites. Here, numerical simulations indicate that increasing nutrient availability not only promotes the formation of dinoflagellate blooms but can also stimulate their control by protozoan parasites. Moreover, high abundance of phytoplankton other than dinoflagellate hosts might have a significant dilution effect on the control of dinoflagellate blooms by parasites, either by resource competition with dinoflagellates (thus limiting the number of hosts available for infection) or by affecting numerical-functional responses of grazers that consume free-living parasite stages. These outcomes indicate that although both dinoflagellates and their protozoan parasites are directly affected by nutrient availability, the efficacy of the parasitic control of dinoflagellate blooms under temporary eutrophication depends strongly on the structure of the plankton community as a whole.

Phytoplankton diversity and community composition along the estuarine gradient of a temperate macrotidal ecosystem: combined morphological and molecular approaches

Microscopical and molecular analyses were used to investigate the diversity and spatial community structure of spring phytoplankton all along the estuarine gradient in a macrotidal ecosystem, the Baie des Veys (eastern English Channel). Taxa distribution at high tide in the water column appeared to be mainly driven by the tidal force which superimposed on the natural salinity gradient, resulting in a two-layer flow within the channel. Lowest taxa richness and abundance were found in the bay where Teleaulax-like cryptophytes dominated. A shift in species composition occurred towards the mouth of the river, with the diatom Asterionellopsis glacialis dramatically accumulating in the bottom waters of the upstream brackish reach. Small thalassiosiroid diatoms dominated the upper layer river community, where taxa richness was higher. Through the construction of partial 18S rDNA clone libraries, the microeukaryotic diversity was further explored for three samples selected along the surface salinity gradient (freshwater - brackish - marine). Clone libraries revealed a high diversity among heterotrophic and/or small-sized protists which were undetected by microscopy. Among them, a rich variety of Chrysophyceae and other lineages (e.g. novel marine stramenopiles) are reported here for the first time in this transition area. However, conventional microscopy remains more efficient in revealing the high diversity of phototrophic taxa, low in abundances but morphologically distinct, that is overlooked by the molecular approach. The differences between microscopical and molecular analyses and their limitations are discussed here, pointing out the complementarities of both approaches, for a thorough phytoplankton community description.

Cyttarocylis ampulla, a polymorphic tintinnid ciliate of the marine plankton

Tintinnid species are traditionally distinguished via lorica features. Recently, sequencing has revealed polymorphism, i.e., genetically identical individuals with distinct lorica morphologies. One such polymorphic species is Cyttarocylis ampulla; individuals can display lorica morphologies of formally different species of Cyttarocylis and Petalotricha, well-represented in the literature. We compiled and analysed a global database of species records to determine if there is a main form and if different morphotypes have distinct temporal or spatial distributions. The two genera show very similar widespread distributions but with some statistical evidence of spatial segregation. Examining co-occurrence among the common 'species' we found most were rarely found alone, only 6-14% of the records for all species except for 2 forms: C. eucecryphalus and P. ampulla reported alone in 34% and 43%, respectively, of their records. We identify them as the main forms and analysed data of global distributions, spatial distribution across the Mediterranean in summer and winter and temporal distributions from a site in the Adriatic. The two main forms show frequent co-occurrence, similar lack of strong seasonality and widespread geographic distributions. We tentatively conclude that the different lorica morphologies may only reflect conditions of high temporally variability such as quantities and composition of prey. Directions for further research are suggested.

Dinoflagellate phylogeny revisited: using ribosomal proteins to resolve deep branching dinoflagellate clades

The alveolates are composed of three major lineages, the ciliates, dinoflagellates, and apicomplexans. Together these 'protist' taxa play key roles in primary production and ecology, as well as in illness of humans and other animals. The interface between the dinoflagellate and apicomplexan clades has been an area of recent discovery, blurring the distinction between these two clades. Moreover, phylogenetic analysis has yet to determine the position of basal dinoflagellate clades hence the deepest branches of the dinoflagellate tree currently remain unresolved. Large-scale mRNA sequencing was applied to 11 species of dinoflagellates, including strains of the syndinean genera Hematodinium and Amoebophrya, parasites of crustaceans and dinoflagellates, respectively, to optimize and update the dinoflagellate tree. From the transcriptome-scale data a total of 73 ribosomal protein-coding genes were selected for phylogeny. After individual gene orthology assessment, the genes were concatenated into a >15,000 amino acid alignment with 76 taxa from dinoflagellates, apicomplexans, ciliates, and the outgroup heterokonts. Overall the tree was well resolved and supported, when the data was subsampled with gblocks or constraint trees were tested with the approximately unbiased test. The deepest branches of the dinoflagellate tree can now be resolved with strong support, and provides a clearer view of the evolution of the distinctive traits of dinoflagellates.

Phylogenetic investigations on ten genera of tintinnid ciliates (Ciliophora: Spirotrichea: Tintinnida), based on small subunit ribosomal RNA gene sequences

Tintinnida is a diverse taxon that accommodates over 1,500 morphospecies, which is an important component of marine planktonic food webs. However, evolutionary relationships of tintinnids are poorly known because molecular data of most groups within this order are lacking. In our study, the small subunit (SSU) rRNA genes representing 10 genera, 5 families of Tintinnida were sequenced, including the first SSU rRNA gene sequences for Coxliella, Dadayiella, Epiplocyloides, and Protorhabdonella, and phylogenetic trees were constructed to assess their intergeneric relationships. Phylogenies inferred from different methods showed that (1) Three newly sequenced Eutintinnus species fell into Eutintinnus clade forming a sister group to the clade containing Amphorides, Steenstrupiella, Amphorellopsis, and Salpingella; (2) Surprisingly, the genetic distances between Amphorides amphora and Amphorellopsis acuta population 1 was even smaller than that between the two populations of Amphorellopsis acuta, casting doubt on the validity of Amphorides and Amphorellopsis as presently defined; (3) The SSU rRNA sequences of Dadayiella ganymedes and Parundella aculeata were almost identical. Therefore, Parundella ganymedes novel combination is proposed; (4) Coxliella, which is currently assigned within Metacylididae, branched instead with some Tintinnopsis species. Furthermore, the validation of Coxliella, which was considered to be a "questionable" genus, was confirmed based on evidences from morphology, ecology, and molecular data; (5) Protorhabdonella and Rhabdonella showed rather low intergeneric distance and grouped together with strong support suggesting that Rhabdonellidae is a well-defined taxon; and (6) Epiplocyloides branched with species in Cyttarocylididae indicating their close relationship.