Molecular characterization of the myxosporean associated with parasitic encephalitis of farmed Atlantic salmon Salmo salar in Ireland

During seasonal epizootics of neurologic disease and mass mortality in the summers of 1992, 1993 and 1994 on a sea-farm in Ireland, Atlantic salmon Salmo salar smolts suffered from encephalitis associated with infection by a neurotropic parasite. Based on ultrastructural studies, this neurotropic parasite was identified as an intercellular presporogonic multicellular developmental stage of a histozoic myxosporean, possibly a Myxobolus species. In order to generate sequence data for phylogenetic comparisons to substantiate the present morphological identification of this myxosporean in the absence of detectable sporogony, polymerase chain reaction (PCR), Southern blot hybridization, dideoxynucleotide chain-termination DNA sequencing, and in situ hybridization (ISH) were used in concert to characterize segments of the small subunit ribosomal RNA (SSU rRNA) gene. Oligonucleotide primers were created from sequences of the SSU rRNA gene of M. cerebralis and were employed in PCR experiments using DNA extracted from formalin-fixed paraffin-embedded tissue sections of brains from Atlantic salmon smolts in which the myxosporean had been detected by light microscopy. Five segments of the SSU rRNA gene of the myxosporean, ranging in length from 187 to 287 base pairs, were amplified, detected by hybridization with sequence-specific probes, and sequenced. Consensus sequences from these segments were aligned to create a partial sequence of the SSU rRNA gene of the myxosporean. Assessments of sequence identity were made between this partial sequence and sequences of SSU rRNA genes from 7 myxosporeans, including Ceratomyxa shasta, Henneguya doori, M. arcticus, M. cerebralis, M. insidiosus, M. neurobius, and M. squamalis. The partial SSU rRNA gene sequence from the myxosporean had more sequence identity with SSU rRNA gene sequences from neurotropic and myotropic species of Myxobolus than to those from epitheliotropic species of Myxobolus or Henneguya, or the enterotropic species of Ceratomyxa, and was identical to regions of the SSU rRNA gene of M. cerebralis. Digoxigenin-labeled oligonucleotide DNA probes complementary to multiple segments of the SSU rRNA gene of M. cerebralis hybridized with DNA of the parasite in histologic sections of brain in ISH experiments, demonstrating definitively that the segments of genome amplified were from the organisms identified by histology and ultrastructural analysis. Based on sequence data derived entirely from genetic material of extrasporogonic stages, the SSU rDNA sequence identity discovered in this study supports the hypothesis that the myxosporean associated with encephalitis of farmed Atlantic salmon smolts is a neurotropic species of the genus Myxobolus, with sequences identical to those of M. cerebralis.     

Recognition of chimeric small-subunit ribosomal DNAs composed of genes from uncultivated microorganisms.

When PCR was used to recover small-subunit (SSU) rRNA genes from a hot spring cyanobacterial mat community, chimeric SSU rRNA sequences which exhibited little or no secondary structural abnormality were recovered. They were revealed as chimeras of SSU rRNA genes of uncultivated species through separate phylogenetic analysis of short sequence domains.     

Phylogenetic diversity and whole-cell hybridization of oxymonad flagellates from the hindgut of the wood-feeding lower termite Reticulitermes flavipes

SSU rRNA genes of oxymonad protists from the hindgut of the wood-feeding termite Reticulitermes flavipes were PCR-amplified using a newly designed oxymonad-specific forward primer and a newly designed reverse primer specific for termite gut flagellates. After cloning, the clone library was sorted into four groups by RFLP analysis and nearly full-length SSU rRNA gene sequences were obtained for representative clones from each group. Phylogenetic analysis revealed that sequences of all four groups formed a monophyletic cluster with the only other existing SSU rRNA gene sequence of oxymonads. Using whole-cell hybridization with clone-specific fluorescently labeled probes, each of the four clone groups could be assigned to a specific morphotype, which were identified as Dinenympha gracilis, Dinenympha fimbriata, and so-far undescribed species of Pyrsonympha and Dinenympha. Our results demonstrate that the morphological variety of oxymonads is not caused by the presence of different developmental stages of the same organism, but that the various morphotypes represent different species.     

Protistan microbial observatory in the Cariaco Basin,Caribbean. I. Pyrosequencing vs Sanger insights into species richness

Microbial diversity and distribution are topics of intensive research. In two companion papers in this issue, we describe the results of the Cariaco Microbial Observatory (Caribbean Sea, Venezuela). The Basin contains the largest body of marine anoxic water, and presents an opportunity to study protistan communities across biogeochemical gradients. In the first paper, we survey 18S ribosomal RNA (rRNA) gene sequence diversity using both Sanger- and pyrosequencing-based approaches, employing multiple PCR primers, and state-of-the-art statistical analyses to estimate microbial richness missed by the survey. Sampling the Basin at three stations, in two seasons, and at four depths with distinct biogeochemical regimes, we obtained the largest, and arguably the least biased collection of over 6000 nearly full-length protistan rRNA gene sequences from a given oceanographic regime to date, and over 80 000 pyrosequencing tags. These represent all major and many minor protistan taxa, at frequencies globally similar between the two sequence collections. This large data set provided, via the recently developed parametric modeling, the first statistically sound prediction of the total size of protistan richness in a large and varied environment, such as the Cariaco Basin: over 36 000 species, defined as almost full-length 18S rRNA gene sequence clusters sharing over 99% sequence homology. This richness is a small fraction of the grand total of known protists (over 100 000–500 000 species), suggesting a degree of protistan endemism.     

Methodology of protistan discovery: from rRNA detection to quality scanning electron microscope images

Each year, thousands of new protistan 18S rRNA sequences are detected in environmental samples. Many of these sequences are molecular signatures of new protistan species, classes, and/or kingdoms that have never been seen before. The main goal of this study was to enable visualization of these novel organisms and to conduct quality ultrastructural examination. We achieved this goal by modifying standard procedures for cell fixation, fluorescence in situ hybridization, and scanning electron microscopy (SEM) and by making these methodologies work in concert. As a result, the same individual cell can now be detected by 18S rRNA-targeted fluorochrome-labeled probes and then viewed by SEM to reveal its diagnostic morphological characteristics. The method was successfully tested on a wide range of protists (alveolates, stramenopiles, kinetoplastids, and cryptomonads). The new methodology thus opens a way for fine microscopy studies of many organisms previously known exclusively by their 18S rRNA sequences.     

Can Abundance of Protists Be Inferred from Sequence Data: A Case Study of Foraminifera

Protists are key players in microbial communities, yet our understanding of their role in ecosystem functioning is seriously impeded by difficulties in identification of protistan species and their quantification. Current microscopy-based methods used for determining the abundance of protists are tedious and often show a low taxonomic resolution. Recent development of next-generation sequencing technologies offered a very powerful tool for studying the richness of protistan communities. Still, the relationship between abundance of species and number of sequences remains subjected to various technical and biological biases. Here, we test the impact of some of these biological biases on sequence abundance of SSU rRNA gene in foraminifera. First, we quantified the rDNA copy number and rRNA expression level of three species of foraminifera by qPCR. Then, we prepared five mock communities with these species, two in equal proportions and three with one species ten times more abundant. The libraries of rDNA and cDNA of the mock communities were constructed, Sanger sequenced and the sequence abundance was calculated. The initial species proportions were compared to the raw sequence proportions as well as to the sequence abundance normalized by rDNA copy number and rRNA expression level per species. Our results showed that without normalization, all sequence data differed significantly from the initial proportions. After normalization, the congruence between the number of sequences and number of specimens was much better. We conclude that without normalization, species abundance determination based on sequence data was not possible because of the effect of biological biases. Nevertheless, by taking into account the variation of rDNA copy number and rRNA expression level we were able to infer species abundance, suggesting that our approach can be successful in controlled conditions.     

Phylogenetic analysis of ribosomal RNA operons from uncultivated coastal marine bacterioplankton

Analyses of small subunit ribosomal RNA genes (SSU rDNAs) have significantly influenced our understanding of the composition of aquatic microbial assemblages. Unfortunately, SSU rDNA sequences often do not have sufficient resolving power to differentiate closely related species. To address this general problem for uncultivated bacterioplankton taxa, we analysed and compared sequences of polymerase chain reaction (PCR)-generated and bacterial artificial chromosome (BAC)-derived clones that contained most of the SSU rDNAs, the internal transcribed spacer (ITS) and the large subunit ribosomal RNA gene (LSU rDNA). The phylogenetic representation in the rRNA operon PCR library was similar to that reported previously in coastal bacterioplankton SSU rDNA libraries. We observed good concordance between the phylogenetic relationships among coastal bacterioplankton inferred from SSU or LSU rDNA sequences. ITS sequences confirmed the close intragroup relationships among members of the SAR11, SAR116 and SAR86 clades that were predicted by SSU and LSU rDNA sequence analyses. We also found strong support for homologous recombination between the ITS regions of operons from the SAR11 clade.     

A discontinuous small subunit ribosomal RNA in Tetrahymena pyriformis mitochondria

We show here that in the mitochondria of Tetrahymena pyriformis, the small subunit (SSU) rRNA is discontinuous, being comprised of two separate components which we term "alpha" (a novel low molecular weight RNA, approximately equal to 200 nucleotides long) and "beta" (a previously described 14 S RNA). The SSU alpha rRNA has been sequenced in its entirety; it represents the immediate 5'-terminal domain of conventional SSU rRNA. The sequences at the ends of the SSU beta rRNA have also been determined; they show that this molecule corresponds to the 3'-terminal 7/8 of conventional SSU rRNA. A 2.5-kilobase pair XbaI restriction fragment of T. pyriformis mitochondrial DNA which contains the SSU alpha and SSU beta rRNA genes was cloned and its complete nucleotide sequence was determined. This revealed that the genes encoding the two segments of SSU rRNA are separated by a 54-base pair (A + T)-rich spacer. The alpha and beta sequences can be fitted to a generalized secondary structure model for eubacterial 16 S rRNA, with the two RNA species associating through long range interactions to form base-paired regions characteristic of SSU rRNA. In this model, the spacer is situated in a region of pronounced primary and secondary structural variation among SSU rRNAs. The significance of these findings with respect to rRNA biosynthesis and processing and the possible evolutionary relationship between spacers and variable regions in rRNA genes is discussed.     

Comparative microbial diversity in the gastrointestinal tracts of food animal species

Molecular tools based on small subunit (SSU) rDNA gene sequencesoffer a powerful and rapid tool for the analysis of complexmicrobial communities found in the gastrointestinal tracts (GIT)of food animal species. Extensive comparative sequence analysisof SSU rRNA molecules representing a wide diversity of organismsshows that different regions of the molecule vary in sequenceconservation. Oligonucleotides complementing regions of universallyconversed SSU rRNA sequences are used as universal probes, whilethose complementing more variable regions of sequence are usefulas selective probes targeting species, genus, or phylogeneticgroups. Different approaches derive different information andthis is highly dependent on the type of target nucleic acidemployed and the conceptual and technical basis used for nucleicacid probe design. Generally these approaches can be dividedinto DNA-based methods employing empirically characterized probesand rRNA-based methods based on comparative sequence analysisfor design and interpretation of "rational" probes. Polymerasechain reaction (PCR) based techniques can also be applied tothe analysis of microbial communities in the GIT. Direct cloningof SSU rDNA genes amplified from these complex communities canbe used to determine the extent of diversity in these GIT communities.Denaturing gradient gel electrophoresis (DGGE) is another powerfultool for profiling microbial diversity of microbial communitiesin GI tracts. Sequence analysis of the excised DGGE ampliconscan then be used to presumptively identify predominant bacterialspecies. Examples of how these molecular approaches are beingused to study the microbial diversity of communities from steersfed different diets, swine fed probiotics, and Atlantic salmonfed aquaculture diets are presented.     

Complete sequence and structure of ribosomal RNA gene of Heterosporis anguillarum

The ribosomal RNA (rRNA) gene region of the microsporidium Heterosporis anguillarum has been examined. Complete DNA sequence data (4060 bp, GenBank Accession No. AF402839) of the rRNA gene of H. anguillarum are presented for the small subunit gene (SSU rRNA: 1359 bp), the internal transcribed spacer (ITS: 37 bp), and the large subunit gene (LSU rRNA: 2664 bp). The secondary structures of the H. anguillarum SSU and LSU rRNA genes are constructed and described. This is the first complete sequence of an rRNA gene published for a fish-infecting microsporidian species. In the phylogenetic analysis, the sequences, including partial SSU rRNA, ITS, and partial LSU rRNA sequences of the fish-infecting microsporidia, were aligned and analysed. The taxonomic position of H. anguillarum as suggested by Lom et al. (2000; Dis Aquat Org 43:225-231) is confirmed in this paper.     

Confidence intervals of similarity values determined for cloned SSU rRNA genes from environmental samples

The goal of this research was to investigate the influence of the error rate of sequence determination on the differentiation of cloned SSU rRNA gene sequences for assessment of community structure. SSU rRNA cloned sequences from groundwater samples that represent different bacterial divisions were sequenced multiple times with the same sequencing primer. From comparison of sequence alignments with unedited data, confidence intervals were obtained from both a 'double binomial' model of sequence comparison and by non-parametric methods. The results indicated that similarity values below 0.9946 are likely derived from dissimilar sequences at a confidence level of 0.95, and not sequencing errors. The results confirmed that screening by direct sequence determination could be reliably used to differentiate at the species level. However, given sequencing errors comparable to those seen in this study, sequences with similarities above 0.9946 should be treated as the same sequence if a 95% confidence is desired.     

Extremely acidophilic protists from acid mine drainage host Rickettsiales-lineage endosymbionts that have intervening sequences in their 16S rRNA genes

During a molecular phylogenetic survey of extremely acidic (pH < 1), metal-rich acid mine drainage habitats in the Richmond Mine at Iron Mountain, Calif., we detected 16S rRNA gene sequences of a novel bacterial group belonging to the order Rickettsiales in the Alphaproteobacteria. The closest known relatives of this group (92% 16S rRNA gene sequence identity) are endosymbionts of the protist Acanthamoeba. Oligonucleotide 16S rRNA probes were designed and used to observe members of this group within acidophilic protists. To improve visualization of eukaryotic populations in the acid mine drainage samples, broad-specificity probes for eukaryotes were redesigned and combined to highlight this component of the acid mine drainage community. Approximately 4% of protists in the acid mine drainage samples contained endosymbionts. Measurements of internal pH of the protists showed that their cytosol is close to neutral, indicating that the endosymbionts may be neutrophilic. The endosymbionts had a conserved 273-nucleotide intervening sequence (IVS) in variable region V1 of their 16S rRNA genes. The IVS does not match any sequence in current databases, but the predicted secondary structure forms well-defined stem loops. IVSs are uncommon in rRNA genes and appear to be confined to bacteria living in close association with eukaryotes. Based on the phylogenetic novelty of the endosymbiont sequences and initial culture-independent characterization, we propose the name "Candidatus Captivus acidiprotistae." To our knowledge, this is the first report of an endosymbiotic relationship in an extremely acidic habitat.     

Evolution of the RNA polymerase B' subunit gene (rpoB') in Halobacteriales: a complementary molecular marker to the SSU rRNA gene

Many prokaryotes have multiple ribosomal RNA operons. Generally, sequence differences between small subunit (SSU) rRNA genes are minor (<1%) and cause little concern for phylogenetic inference or environmental diversity studies. For Halobacteriales, an order of extremely halophilic, aerobic Archaea, within-genome SSU rRNA sequence divergence can exceed 5%, rendering phylogenetic assignment problematic. The RNA polymerase B' subunit gene (rpoB') is a single-copy conserved gene that may be an appropriate alternative phylogenetic marker for Halobacteriales. We sequenced a fragment of the rpoB' gene from 21 species, encompassing 15 genera of Halobacteriales. To examine the utility of rpoB' as a phylogenetic marker in Halobacteriales, we investigated three properties of rpoB' trees: the variation in resolution between trees inferred from the rpoB' DNA and RpoB' protein alignment, the degree of mutational saturation between taxa, and congruence with the SSU rRNA tree. The rpoB' DNA and protein trees were for the most part congruent and consistently recovered two well-supported monophyletic groups, the clade I and clade II haloarchaea, within a collection of less well resolved Halobacteriales lineages. A comparison of observed versus inferred numbers of substitution revealed mutational saturation in the rpoB' DNA data set, particularly between more distant species. Thus, the RpoB' protein sequence may be more reliable than the rpoB' DNA sequence for inferring Halobacteriales phylogeny. AU tests of tree selection indicated the trees inferred from rpoB' DNA and protein alignments were significantly incongruent with the SSU rRNA tree. We discuss possible explanations for this incongruence, including tree reconstruction artifact, differential paralog sampling, and lateral gene transfer. This is the first study of Halobacteriales evolution based on a marker other than the SSU rRNA gene. In addition, we present a valuable phylogenetic framework encompassing a broad diversity of Halobacteriales, in which novel sequences can be inserted for evolutionary, ecological, or taxonomic investigations.     

Identification and isolation of a Castellaniella species important during biostimulation of an acidic nitrate- and uranium-contaminated aquifer

Immobilization of uranium in groundwater can be achieved through microbial reduction of U(VI) to U(IV) upon electron donor addition. Microbial community structure was analyzed in ethanol-biostimulated and control sediments from a high-nitrate (>130 mM), low-pH, uranium-contaminated site in Oak Ridge, TN. Analysis of small subunit (SSU) rRNA gene clone libraries and polar lipid fatty acids from sediments revealed that biostimulation resulted in a general decrease in bacterial diversity. Specifically, biostimulation resulted in an increase in the proportion of Betaproteobacteria (10% of total clones in the control sediment versus 50 and 79% in biostimulated sediments) and a decrease in the proportion of Gammaproteobacteria and Acidobacteria. Clone libraries derived from dissimilatory nitrite reductase genes (nirK and nirS) were also dominated by clones related to Betaproteobacteria (98% and 85% of total nirK and nirS clones, respectively). Within the nirK libraries, one clone sequence made up 59 and 76% of sequences from biostimulated sediments but only made up 10% of the control nirK library. Phylogenetic analysis of SSU rRNA and nirK gene sequences from denitrifying pure cultures isolated from the site indicate that all belong to a Castellaniella species; nearly identical sequences also constituted the majority of biostimulated SSU rRNA and nirK clone libraries. Thus, by combining culture-independent with culture-dependent techniques, we were able to link SSU rRNA clone library information with nirK sequence data and conclude that a potentially novel Castellaniella species is important for in situ nitrate removal at this site.     

Synthetic Statistical Approach Reveals a High Degree of Richness of Microbial Eukaryotes in an Anoxic Water Column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.     

Evolution according to large ribosomal subunit RNA

Evolutionary trees were constructed, by distance methods, from an alignment of 225 complete large subunit (LSU) rRNA sequences, representing Eucarya, Archaea, Bacteria, plastids, and mitochondria. A comparison was made with trees based on sets of small subunit (SSU) rRNA sequences. Trees constructed on the set of 172 species and organelles for which the sequences of both molecules are known had a very similar topology, at least with respect to the divergence order of large taxa such as the eukaryotic kingdoms and the bacterial divisions. However, since there are more than ten times as many SSU as LSU rRNA sequences, it is possible to select many SSU rRNA sequence sets of equivalent size but different species composition. The topologies of these trees showed considerable differences according to the particular species set selected.The effect of the dataset and of different distance correction methods on tree topology was tested for both LSU and SSU rRNA by repetitive random sampling of a single species from each large taxon. The impact of the species set on the topology of the resulting consensus trees is much lower using LSU than using SSU rRNA. This might imply that LSU rRNA is a better molecule for studying wide-range relationships. The mitochondria behave clearly as a monophyletic group, clustering with the Proteobacteria. Gram-positive bacteria appear as two distinct groups, which are found clustered together in very few cases. Archaea behave as if monophyletic in most cases, but with a low confidence.Abbreviations LSU rRNA large subunit ribosomal RNA - SSU rRNA small subunit ribosomal RNA - JC Jukes and Cantor - JN Jin and Nei Correspondence to: R. De Wachter     

Exploring Microbial Diversity and Taxonomy Using SSU rRNA Hypervariable Tag Sequencing

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA. As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe. However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process. The resulting taxonomic census provides information on both composition and diversity of the microbial community. To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent. The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match. The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the “rare biosphere” than does capillary sequencing of the full-length gene. In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation. This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.     

Investigating Microbial Eukaryotic Diversity from a Global Census: Insights from a Comparison of Pyrotag and Full-Length Sequences of 18S rRNA Genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.     

Capturing diversity of marine heterotrophic protists: one cell at a time

Recent applications of culture-independent, molecular methods have revealed unexpectedly high diversity in a variety of functional and phylogenetic groups of microorganisms in the ocean. However, none of the existing research tools are free from significant limitations, such as PCR and cloning biases, low phylogenetic resolution and others. Here, we employed novel, single-cell sequencing techniques to assess the composition of small (<10 μm diameter), heterotrophic protists from the Gulf of Maine. Single cells were isolated by flow cytometry, their genomes amplified, and 18S rRNA marker genes were amplified and sequenced. We compared the results to traditional environmental PCR cloning of sorted cells. The diversity of heterotrophic protists was significantly higher in the library of single amplified genomes (SAGs) than in environmental PCR clone libraries of the 18S rRNA gene, obtained from the same coastal sample. Libraries of SAGs, but not clones contained several recently discovered, uncultured groups, including picobiliphytes and novel marine stramenopiles. Clone, but not SAG, libraries contained several large clusters of identical and nearly identical sequences of Dinophyceae, Cercozoa and Stramenopiles. Similar results were obtained using two alternative primer sets, suggesting that PCR biases may not be the only explanation for the observed patterns. Instead, differences in the number of 18S rRNA gene copies among the various protist taxa probably had a significant role in determining the PCR clone composition. These results show that single-cell sequencing has the potential to more accurately assess protistan community composition than previously established methods. In addition, the creation of SAG libraries opens opportunities for the analysis of multiple genes or entire genomes of the uncultured protist groups.     

Amplification of mitochondrial small subunit ribosomal DNA of polypores and its potential for phylogenetic analysis

There has been a systematic need to seek adequate phylogenetic markers that can be applied in phylogenetic analyses of fungal taxa at various levels. The mitochondrial small subunit ribosomal DNA (mt SSU rDNA) is generally considered to be one of the molecules that are appropriate for phylogenetic analyses at a family level. In order to obtain universal primers for polypores of Hymenomycetes, mt SSU rRNA genes were cloned from Bjerkandera adusta, Ganoderma lucidum, Phlebiopsis gigantea, and Phellinus laevigatus and their sequences were determined. Based on the conserved sequences of cloned genes from polypores and Agrocybe aegerita, PCR primers were designed for amplification and sequencing of mt SSU rDNAs. New primers allowed effective amplification and sequencing of almost full-sized genes from representative species of polypores and related species. Phylogenetic relationships were resolved quite efficiently by mt SSU rDNA sequences, and they proved to be more useful in phylogenetic reconstruction of Ganoderma than nuclear internal transcribed spacer (ITS) rDNA sequences.