Protistan diversity in the Arctic: a case of paleoclimate shaping modern biodiversity?

18S rDNA gene metabarcoding of microeukaryotes and epi-endophytes in the holobiome of seven species of large brown algae

Brown algae (Phaeophyceae) are habitat-forming species in coastal ecosystems and include kelp forests and seaweed beds that support a wide diversity of marine life. Host-associated microbial communities are an integral part of phaeophyte biology, and whereas the bacterial microbial partners have received considerable attention, the microbial eukaryotes associated with brown algae have hardly been studied. Here, we used broadly targeted "pan-eukaryotic" primers (metabarcoding) to investigate brown algal-associated eukaryotes (the eukaryome). Using this approach, we aimed to investigate the eukaryome of seven large brown algae that are important and common species in coastal ecosystems. We also aimed to assess whether these macroalgae harbor novel eukaryotic diversity and to ascribe putative functional roles to the host-associated eukaryome based on taxonomic affiliation and phylogenetic placement. We detected a significant diversity of microeukaryotic and algal lineages associated with the brown algal species investigated. The operational taxonomic units (OTUs) were taxonomically assigned to 10 of the eukaryotic major supergroups, including taxonomic groups known to be associated with seaweeds as epibionts, endobionts, parasites, and commensals. Additionally, we revealed previously unrecorded sequence types, including novel phaeophyte OTUs, particularly in the Fucus spp. samples, that may represent fucoid genomic variants, sequencing artifacts, or undescribed epi-/endophytes. Our results provide baseline data and technical insights that will be useful for more comprehensive seaweed eukaryome studies investigating the evidently lineage-rich and functionally diverse symbionts of brown algae.

Community dynamics of microbial eukaryotes in intertidal mudflats in the hypertidal Bay of Fundy

Protists (microbial eukaryotes) are a critically important but understudied group of microorganisms. They are ubiquitous, represent most of the genetic and functional diversity among eukaryotes, and play essential roles in nutrient and energy cycling. Yet, protists remain a black box in marine sedimentary ecosystems like the intertidal mudflats in the Bay of Fundy. The harsh conditions of the intertidal zone and high energy nature of tides in the Bay of Fundy provide an ideal system for gaining insights into the major food web players, diversity patterns and potential structuring influences of protist communities. Our 18S rDNA metabarcoding study quantified seasonal variations and vertical stratification of protist communities in Bay of Fundy mudflat sediments. Three 'SAR' lineages were consistently dominant (in terms of abundance, richness, and prevalence), drove overall community dynamics and formed the core microbiome in sediments. They are Cercozoa (specifically thecate, benthic gliding forms), Bacillariophyta (mainly cosmopolitan, typically planktonic diatoms), and Dinophyceae (dominated by a toxigenic, bloom-forming species). Consumers were the dominant trophic functional group and were comprised mostly of eukaryvorous and bacterivorous Cercozoa, and omnivorous Ciliophora, while phototrophs were dominated by Bacillariophyta. The codominance of Apicomplexa (invertebrate parasites) and Syndiniales (protist parasites) in parasite assemblages, coupled with broader diversity patterns, highlighted the combined marine and terrestrial influences on microbial communities inhabiting intertidal sediments. Our findings, the most comprehensive in a hypertidal benthic system, suggest that synergistic interactions of both local and regional processes (notably benthic-pelagic coupling) may drive heterogenous microbial distribution in high-energy coastal systems.

Evidence for an Independent Hydrogenosome-to-Mitosome Transition in the CL3 Lineage of Fornicates

Fornicata, a lineage of a broader and ancient anaerobic eukaryotic clade Metamonada, contains diverse taxa that are ideally suited for evolutionary studies addressing various fundamental biological questions, such as the evolutionary trajectory of mitochondrion-related organelles (MROs), the transition between free-living and endobiotic lifestyles, and the derivation of alternative genetic codes. To this end, we conducted detailed microscopic and transcriptome analyses in a poorly documented strain of an anaerobic free-living marine flagellate, PCS, in the so-called CL3 fornicate lineage. Fortuitously, we discovered that the original culture contained two morphologically similar and closely related CL3 representatives, which doubles the taxon representation within this lineage. We obtained a monoeukaryotic culture of one of them and formally describe it as a new member of the family Caviomonadidae, Euthynema mutabile gen. et sp. nov. In contrast to previously studied caviomonads, the endobiotic Caviomonas mobilis and Iotanema spirale, E. mutabile possesses an ultrastructurally discernible MRO. We sequenced and assembled the transcriptome of E. mutabile, and by sequence subtraction, obtained transcriptome data from the other CL3 clade representative present in the original PCS culture, denoted PCS-ghost. Transcriptome analyses showed that the reassignment of only one of the UAR stop codons to encode Gln previously reported from I. spirale does not extend to its free-living relatives and is likely due to a unique amino acid substitution in I. spirale’s eRF1 protein domain responsible for termination codon recognition. The backbone fornicate phylogeny was robustly resolved in a phylogenomic analysis, with the CL3 clade amongst the earliest branching lineages. Metabolic and MRO functional reconstructions of CL3 clade members revealed that all three, including I. spirale, encode homologs of key components of the mitochondrial protein import apparatus and the ISC pathway, indicating the presence of a MRO in all of them. In silico evidence indicates that the organelles of E. mutabile and PCS-ghost host ATP and H₂ production, unlike the cryptic MRO of I. spirale. These data suggest that the CL3 clade has experienced a hydrogenosome-to-mitosome transition independent from that previously documented for the lineage leading to Giardia.

Seasonal Variability of Photosynthetic Microbial Eukaryotes (<3 µm) in the Kara Sea Revealed by 18S rDNA Metabarcoding of Sediment Trap Fluxes

This survey is the first to explore the seasonal cycle of microbial eukaryote diversity (<3 µm) using the NGS method and a 10-month sediment trap (2018–2019). The long-term trap was deployed from September to June in the northwestern part of the Kara Sea. A water sample collected before the sediment trap was deployed and also analyzed. The taxonomic composition of microbial eukaryotes in the water sample significantly differed from sediment trap samples, characterized by a high abundance of Ciliophora reads and low abundance of Fungi while trap samples contained an order of magnitude less Ciliophora sequences and high contribution of Fungi. Photosynthetic eukaryotes (PEs) accounting for about 34% of total protists reads were assigned to five major divisions: Chlorophyta, Cryptophyta, Dinoflagellata, Haptophyta, and Ochrophyta. The domination of phototrophic algae was revealed in late autumn. Mamiellophyceae and Trebouxiophyceae were the predominant PEs in mostly all of the studied seasons. Micromonas polaris was constantly present throughout the September–June period in the PE community. The obtained results determine the seasonal dynamics of picoplankton in order to improve our understanding of their role in polar ecosystems.

Evidence from the resurrected family Polyrhabdinidae Kamm, 1922 (Apicomplexa: Gregarinomorpha) supports the epimerite, an attachment organelle, as a major eugregarine innovation

Background

Gregarines are a major group of apicomplexan parasites of invertebrates. The gregarine classification is largely incomplete because it relies primarily on light microscopy, while electron microscopy and molecular data in the group are fragmentary and often do not overlap. A key characteristic in gregarine taxonomy is the structure and function of their attachment organelles (AOs). AOs have been commonly classified as “mucrons” or “epimerites” based on their association with other cellular traits such as septation. An alternative proposal focused on the AOs structure, functional role, and developmental fate has recently restricted the terms “mucron” to archigregarines and “epimerite” to eugregarines.

Methods

Light microscopy and scanning and transmission electron microscopy, molecular phylogenetic analyses of ribosomal RNA genes.

Results

We obtained the first data on fine morphology of aseptate eugregarines Polyrhabdina pygospionis and Polyrhabdina cf. spionis, the type species. We demonstrate that their AOs differ from the mucron in archigregarines and represent an epimerite structurally resembling that in other eugregarines examined using electron microscopy. We then used the concatenated ribosomal operon DNA sequences (SSU, 5.8S, and LSU rDNA) of P. pygospionis to explore the phylogeny of eugregarines with a resolution superior to SSU rDNA alone. The obtained phylogenies show that the Polyrhabdina clade represents an independent, deep-branching family in the Ancoroidea clade within eugregarines. Combined, these results lend strong support to the hypothesis that the epimerite is a synapomorphic innovation of eugregarines. Based on these findings, we resurrect the family Polyrhabdinidae Kamm, 1922 and erect and diagnose the family Trollidiidae fam. n. within the superfamily Ancoroidea Simdyanov et al., 2017. Additionally, we re-describe the characteristics of P. pygospionis, emend the diagnoses of the genus Polyrhabdina, the family Polyrhabdinidae, and the superfamily Ancoroidea.

Microeukaryotic Communities Associated With the Seagrass Zostera marina Are Spatially Structured

Epibiotic microorganisms link seagrass productivity to higher trophic levels, but little is known about the processes structuring these communities, and which taxa consistently associate with seagrass. We investigated epibiotic microeukaryotes on seagrass (Zostera marina) leaves, substrates, and planktonic microeukaryotes in ten meadows in the Northeast Pacific. Seagrass epibiotic communities are distinct from planktonic and substrate communities. We found sixteen core microeukaryotes, including dinoflagellates, diatoms, and saprotrophic stramenopiles. Some likely use seagrass leaves as a substrate, others for grazing, or they may be saprotrophic organisms involved in seagrass decomposition or parasites; their relatives have been previously reported from marine sediments and in association with other hosts such as seaweeds. Core microeukaryotes were spatially structured, and none were ubiquitous across meadows. Seagrass epibiota were more spatially structured than planktonic communities, mostly due to spatial distance and changes in abiotic conditions across space. Seawater communities were relatively more similar in composition across sites and more influenced by the environmental component, but more variable over time. Core and transient taxa were both mostly structured by spatial distance and the abiotic environment, with little effect of host attributes, further indicating that those core taxa would not show a strong specific association with Z. marina.

Nutritional intake of Aplanochytrium (Labyrinthulea, Stramenopiles) from living diatoms revealed by culture experiments suggesting the new prey-predator interactions in the grazing food web of the marine ecosystem

Labyrinthuleans (Labyrinthulea, Stramenopiles) are recognized as decomposers in marine ecosystems but their nutrient sources are not fully understood. We conducted two-membered culture experiments with labyrinthuleans and diatoms to discover where labyrinthuleans obtain their nutrients from. The results showed that Aplanochytrium strains obtained nutrients by consuming living diatoms. Aplanochytrium cells did not release digestive enzymes into the medium, but adhered to diatom cells via the tip of their characteristic ectoplasmic net system to obtain nutrients from them. The chloroplast and cell contents of the diatoms shrank and were absorbed, and then the number of Aplanochytrium cells rapidly increased as multiple aplanospores were released. To estimate the effect of labyrinthulean organisms including Aplanochytrium on marine ecosystem, we explored the dataset generated by the Tara Oceans Project from a wide range of oceanic regions. The average proportion of all labyrinthulean sequences to diatom sequences at each station was about 10%, and labyrinthulids, oblongichytrids, and aplanochytrids were the major constituent genera, accounting for more than 80% of labyrinthuleans. Therefore, these groups are suggested to greatly affect the marine ecosystem. There were positive correlations between aplanochytrids and phototrophs, green algae, and diatoms. At many stations, relatively large proportions of aplanochytrid sequences were detected in the size fraction larger than their cell size. This implied that Aplanochytrium cells increased their particle size by adhering to each other and forming aggregates with diatoms that are captured by larger zooplankton in the environment, thereby bypassing the food web pathway via aplanochytrids to higher predators. The intake of nutrients from diatoms by aplanochytrids represents a newly recognized pathway in the grazing food chain in the marine ecosystem.

A new view on the morphology and phylogeny of eugregarines suggested by the evidence from the gregarine Ancora sagittata (Leuckart, 1860) Labbé, 1899 (Apicomplexa: Eugregarinida)

BACKGROUND

Gregarines are a group of early branching Apicomplexa parasitizing invertebrate animals. Despite their wide distribution and relevance to the understanding the phylogenesis of apicomplexans, gregarines remain understudied: light microscopy data are insufficient for classification, and electron microscopy and molecular data are fragmentary and overlap only partially.

METHODS

Scanning and transmission electron microscopy, PCR, DNA cloning and sequencing (Sanger and NGS), molecular phylogenetic analyses using ribosomal RNA genes (18S (SSU), 5.8S, and 28S (LSU) ribosomal DNAs (rDNAs)).

RESULTS AND DISCUSSION

We present the results of an ultrastructural and molecular phylogenetic study on the marine gregarine Ancora sagittata from the polychaete Capitella capitata followed by evolutionary and taxonomic synthesis of the morphological and molecular phylogenetic evidence on eugregarines. The ultrastructure of Ancora sagittata generally corresponds to that of other eugregarines, but reveals some differences in epicytic folds (crests) and attachment apparatus to gregarines in the family Lecudinidae, where Ancora sagittata has been classified. Molecular phylogenetic trees based on SSU (18S) rDNA reveal several robust clades (superfamilies) of eugregarines, including Ancoroidea superfam. nov., which comprises two families (Ancoridae fam. nov. and Polyplicariidae) and branches separately from the Lecudinidae; thus, all representatives of Ancoroidea are here officially removed from the Lecudinidae. Analysis of sequence data also points to possible cryptic species within Ancora sagittata and the inclusion of numerous environmental sequences from anoxic habitats within the Ancoroidea. LSU (28S) rDNA phylogenies, unlike the analysis of SSU rDNA alone, recover a well-supported monophyly of the gregarines involved (eugregarines), although this conclusion is currently limited by sparse taxon sampling and the presence of fast-evolving sequences in some species. Comparative morphological analyses of gregarine teguments and attachment organelles lead us to revise their terminology. The terms "longitudinal folds" and "mucron" are restricted to archigregarines, whereas the terms "epicystic crests" and "epimerite" are proposed to describe the candidate synapomorphies of eugregarines, which, consequently, are considered as a monophyletic group. Abolishing the suborders Aseptata and Septata, incorporating neogregarines into the Eugregarinida, and treating the major molecular phylogenetic lineages of eugregarines as superfamilies appear as the best way of reconciling recent morphological and molecular evidence. Accordingly, the diagnosis of the order Eugregarinida Léger, 1900 is updated.

Primary Structure of 28S rRNA Gene Confirms Monophyly of Free-Living Heterotrophic and Phototrophic Apicomplexans (Alveolata)

Phylogenetic analysis of large subunit ribosomal RNA (LSU rRNA or 28S rRNA) gene sequences from free-living predatory flagellates Colpodella angusta, Voromonas pontica, and Alphamonas edax (Apicomplexa) confirms their close relationship with chromerids Chromera velia and Vitrella brassicaformis, which possess a functional photosynthetic plastid. Together these organisms form a sister group to parasitic apicomplexans (coccidians and gregarines, or sporozoans sensu lato). This result agrees with the previous conclusion on monophyly of colpodellids and chromerids (chrompodellids) based on phylogenomic data. The revealed relationships demonstrate a complex pattern of acquisition, loss, or modification of plastids and transition to parasitism during alveolate evolution.

Metagenomic analysis of lacustrine viral diversity along a latitudinal transect of the Antarctic Peninsula

Environmental viruses constitute the most abundant biological entities on earth, and harbor an enormous genetic diversity. While their strong influence on the ecosystem is widely acknowledged, current knowledge about their diversity and distribution remains limited. Here we present the metagenomic study of viral communities from freshwater bodies located along a transect of the Antarctic Peninsula. These ecosystems were chosen on the basis of environmental and biogeographical variation. The results obtained indicate that the virus assemblages were diverse, and that the larger fraction represented viruses with no close relatives in the databases. Comparisons to existing metaviromes showed that the communities studied were dissimilar to other freshwater viromes including those from the Arctic. Finally, we observed no indication of there being a reduction in either viral richness or diversity estimates with increasing latitude along the studied transect, further adding to the controversy regarding the possible existence of latitudinal gradients of diversity in the microbial world.

New record and phylogenetic affinities of the oomycete Olpidiopsis feldmanni infecting Asparagopsis sp. (Rhodophyta)

A new geographic record of the oomycete Olpidiopsis feldmanni infecting the tetrasporophytic stage of the red alga Asparagopsis sp. from the Adriatic Sea, confirmed through morphological identification, allowed us to expand previous observations of this organism. Ultrastructural investigations of environmental material showed a large central vacuole and a cell wall thicker than previously reported from other basal oomycete pathogens of algae. Phylogenetic analysis closely associates O. feldmanni to O. bostrychiae concurrent with structural observations. This constitutes the first genetic characterisation of an Olpidiopsis species that was initially described before 1960, adding to the genetic data of 3 other marine Olpidiopsis species established and genetically characterised in the last 2 decades. The paper discusses concurrences of the ultrastructural observations made here and in previous studies of the marine Olpidiopsis species with those made on the freshwater species.

Is Planktonic Diversity Well Recorded in Sedimentary DNA? Toward the Reconstruction of Past Protistan Diversity

Studies based on the coupling of a paleolimnological approach and molecular tools (e.g., sequencing of sedimentary DNA) present a promising opportunity to obtain long-term data on past lacustrine biodiversity. However, certain validations are still required, such as the evaluation of DNA preservation in sediments for various planktonic taxa that do not leave any morphological diagnostic features. In this study, we focused on the diversity of planktonic unicellular eukaryotes and verified the presence of their DNA in sediment archives. We compared the molecular inventories (high-throughput sequencing of 18S ribosomal DNA) obtained from monitoring the water column with those obtained for DNA archived in the first 30 cm of sediment. Seventy-one percent of taxonomic units found in the water samples were detected in sediment samples, including pigmented taxa, such as Chlorophyta, Dinophyceae, and Chrysophyceae, phagotrophic taxa, such as Ciliophora, parasitic taxa, such as Apicomplexa and Chytridiomycota, and saprotrophs, such as Cryptomycota. Parallel analysis of 18S ribosomal RNA (rRNA) transcripts revealed the presence of living eukaryotic taxa only in the top 2 cm of sediment; although some limits exist in using RNA/DNA ratio as indicator of microbial activity, these results suggested that the sedimentary DNA mostly represented DNA from past and inactive communities. Only the diversity of a few groups, such as Cryptophyta and Haptophyta, seemed to be poorly preserved in sediments. Our overall results showed that the application of sequencing techniques to sedimentary DNA could be used to reconstruct past diversity for numerous planktonic eukaryotic groups.

Vermistella arctica n. sp. Nominates the Genus Vermistella as a Candidate for Taxon with Bipolar Distribution

A new amoebozoan species, Vermistella arctica n. sp., is described from marine habitats in the central part of Svalbard archipelago. This is the first report on Arctic amoebae belonging to the genus Vermistella Moran and Anderson, 2007, the type species of which was described from the opposite pole of the planet. Psychrophily proved in the new strains qualifies the genus Vermistella as a bipolar taxon. Molecular phylogenetic analyses based on 18S rDNA and actin sequences did not show any affinity of the genus Vermistella to Stygamoeba regulata ATCC(®) 50892(™) strain. A close phylogenetic relationship was found between Vermistella spp. and a sequence originating from an environmental sample from Cariaco basin, the largest marine permanently anoxic system in the world. Possible mechanisms of bipolar distribution are discussed.

Assessing the global phylum level diversity within the bacterial domain: A review

Microbial ecology is the study of microbes in the natural environment and their interactions with each other. Investigating the nature of microorganisms residing within a specific habitat is an extremely important component of microbial ecology. Such microbial diversity surveys aim to determine the identity, physiological preferences, metabolic capabilities, and genomic features of microbial taxa within a specific ecosystem. A comprehensive review of various aspects of microbial diversity (phylogenetic, functional, and genomic diversities) in the microbial (bacterial, archaeal, and microeukaryotic) world is clearly a daunting task that could not be aptly summarized in a single review. Here, we focus on one aspect of diversity (phylogenetic diversity) in one microbial domain (the Bacteria). We restrict our analysis to the highest taxonomic rank (phylum) and attempt to investigate the extent of global phylum level diversity within the Bacteria. We present a brief historical perspective on the subject and highlight how the adaptation of molecular biological and phylogenetic approaches has greatly expanded our view of global bacterial diversity. We also summarize recent progress toward the discovery of novel bacterial phyla, present evidences that the scope of phylum level diversity in nature has hardly been exhausted, and propose novel approaches that could greatly facilitate the discovery process of novel bacterial phyla within various ecosystems.

Intracellular diversity of the V4 and V9 regions of the 18S rRNA in marine protists (radiolarians) assessed by high-throughput sequencing

Metabarcoding is a powerful tool for exploring microbial diversity in the environment, but its accurate interpretation is impeded by diverse technical (e.g. PCR and sequencing errors) and biological biases (e.g. intra-individual polymorphism) that remain poorly understood. To help interpret environmental metabarcoding datasets, we investigated the intracellular diversity of the V4 and V9 regions of the 18S rRNA gene from Acantharia and Nassellaria (radiolarians) using 454 pyrosequencing. Individual cells of radiolarians were isolated, and PCRs were performed with generalist primers to amplify the V4 and V9 regions. Different denoising procedures were employed to filter the pyrosequenced raw amplicons (Acacia, AmpliconNoise, Linkage method). For each of the six isolated cells, an average of 541 V4 and 562 V9 amplicons assigned to radiolarians were obtained, from which one numerically dominant sequence and several minor variants were found. At the 97% identity, a diversity metrics commonly used in environmental surveys, up to 5 distinct OTUs were detected in a single cell. However, most amplicons grouped within a single OTU whereas other OTUs contained very few amplicons. Different analytical methods provided evidence that most minor variants forming different OTUs correspond to PCR and sequencing artifacts. Duplicate PCR and sequencing from the same DNA extract of a single cell had only 9 to 16% of unique amplicons in common, and alignment visualization of V4 and V9 amplicons showed that most minor variants contained substitutions in highly-conserved regions. We conclude that intracellular variability of the 18S rRNA in radiolarians is very limited despite its multi-copy nature and the existence of multiple nuclei in these protists. Our study recommends some technical guidelines to conservatively discard artificial amplicons from metabarcoding datasets, and thus properly assess the diversity and richness of protists in the environment.

Phylogeny and biogeography of an uncultured clade of snow chytrids

Numerous studies have shown that snow can contain a diverse array of algae known as 'snow algae'. Some reports also indicate that parasites of algae (e.g. chytrids) are also found in snow, but efforts to phylogenetically identify 'snow chytrids' have not been successful. We used culture-independent molecular approaches to phylogenetically identify chytrids that are common in long-lived snowpacks of Colorado and Europe. The most remarkable finding of the present study was the discovery of a new clade of chytrids that has representatives in snowpacks of Colorado and Switzerland and cold sites in Nepal and France, but no representatives from warmer ecosystems. This new clade ('Snow Clade 1' or SC1) is as deeply divergent as its sister clade, the Lobulomycetales, and phylotypes of SC1 show significant (P < 0.003) genetic-isolation by geographic distance patterns, perhaps indicating a long evolutionary history in the cryosphere. In addition to SC1, other snow chytrids were phylogenetically shown to be in the order Rhizophydiales, a group with known algal parasites and saprotrophs. We suggest that these newly discovered snow chytrids are important components of snow ecosystems where they contribute to snow food-web dynamics and the release of nutrients due to their parasitic and saprotrophic activities.

General patterns of diversity in major marine microeukaryote lineages

Microeukaryotes have vital roles for the functioning of marine ecosystems, but still some general characteristics of their current diversity and phylogeny remain unclear. Here we investigated both aspects in major oceanic microeukaryote lineages using 18S rDNA (V4–V5 hypervariable regions) sequences from public databases that derive from various marine environmental surveys. A very carefully and manually curated dataset of 8291 Sanger sequences was generated and subsequently split into 65 taxonomic groups (roughly to Class level based on KeyDNATools) prior to downstream analyses. First, we calculated genetic distances and clustered sequences into Operational Taxonomic Units (OTUs) using different distance cut-off levels. We found that most taxonomic groups had a maximum pairwise genetic distance of 0.25. Second, we used phylogenetic trees to study general evolutionary patterns. These trees confirmed our taxonomic classification and served to run Lineage Through Time (LTT) plots. LTT results indicated different cladogenesis dynamics across groups, with some displaying an early diversification and others a more recent one. Overall, our study provides an improved description of the microeukaryote diversity in the oceans in terms of genetic differentiation within groups as well as in the general phylogenetic structure. These results will be important to interpret the large amount of sequence data that is currently generated by High Throughput Sequencing technologies.

Discovery of a diverse clade of gregarine apicomplexans (Apicomplexa: Eugregarinorida) from Pacific eunicid and onuphid polychaetes, including descriptions of Paralecudina n. gen., Trichotokara japonica n. sp., and T. eunicae n. sp

Marine gregarines are poorly understood apicomplexan parasites with large trophozoites that inhabit the body cavities of marine invertebrates. Two novel species of gregarines were discovered in polychaete hosts collected in Canada and Japan. The trophozoites of Trichotokara japonica n. sp. were oval to rhomboidal shaped, and covered with longitudinal epicytic folds with a density of six to eight folds/micron. The nucleus was situated in the middle of the cell, and the mucron was elongated and covered with hair-like projections; antler-like projections also extended from the anterior tip of the mucron. The distinctively large trophozoites of Trichotokara eunicae n. sp. lacked an elongated mucron and had a tadpole-like cell shape consisting of a bulbous anterior region and a tapered tail-like posterior region. The cell surface was covered with longitudinal epicytic folds with a density of three to five folds/micron. Small subunit (SSU) rDNA sequences of both species were very divergent and formed a strongly supported clade with the recently described species Trichotokara nothriae and an environmental sequence (AB275074). This phylogenetic context combined with the morphological features of T. eunicae n. sp. required us to amend the description for Trichotokara. The sister clade to the Trichotokara clade consisted of environmental sequences and Lecudina polymorpha, which also possesses densely packed epicyctic folds (3-5 folds/micron) and a prominently elongated mucron. This improved morphological and molecular phylogenetic context justified the establishment of Paralecudina (ex. Lecudina) polymorpha n. gen. et comb.

Molecular Diversity of Fungi from Marine Oxygen-Deficient Environments (ODEs)

Molecular diversity surveys of marine fungi have demonstrated that the species richness known to date is just the tip of the iceberg and that there is a large extent of unknown fungal diversity in marine habitats. Reports of novel fungal lineages at higher taxonomic levels are documented from a large number of marine habitats, including the various marine oxygen-deficient environments (ODEs). In the past few years, a strong focus of eukaryote diversity research has been on a variety of ODEs, as these environments are considered to harbor a large number of organisms, which are highly divergent to known diversity and could provide insights into the early eukaryotic evolution. ODEs that have been targeted so far include shallow water sediments, hydrothermal vent systems, deep-sea basins, intertidal habitats, and fjords. Most, if not all, molecular diversity studies in marine ODEs have shown, that contrary to previous assumptions, fungi contribute significantly to the micro-eukaryotic community in such habitats. In this chapter, we have reanalyzed the environmental fungal sequences obtained from the molecular diversity survey in 14 different sites to obtain a comprehensive picture of fungal diversity in these marine habitats. The phylogenetic analysis of the fungal environmental sequences from various ODEs have grouped these sequences into seven distinct clades (Clade 1-7) clustering with well-known fungal taxa. Apart from this, four environmental clades (EnvClade A, B, C, and D) with exclusive environmental sequences were also identified. This has provided information on the positioning of the environmental sequences at different taxonomic levels within the major fungal phylums. The taxonomic distribution of these environmental fungal sequences into clusters and clades has also shown that they are not restricted by geographical boundaries. The distribution pattern together with the reports on the respiratory abilities of fungi under reduced oxygen conditions shows that they are highly adaptive and may have a huge ecological role in these oxygen deficient habitats.

Subulatomonas tetraspora nov. gen. nov. sp. is a member of a previously unrecognized major clade of eukaryotes

While a large number of aerobic free-living protists have been described within the last decade, the number of new anaerobic or microaerophilic microbial eukaryotic taxa has lagged behind. Here we describe a microaerophilic genus and species of amoeboflagellate isolated from a near-shore marine site off the coast at Plymouth, Massachusetts: Subulatomonas tetraspora nov. gen. nov. sp. This taxon is closely related to Breviata anathema based on both microscopical features and phylogenetic analyses of sequences of three genes: SSU-rDNA, actin, and alpha-tubulin. However, Subulatomonas tetraspora nov. gen. nov. sp. and B. anathema are morphologically distinctive, differ by 14.9% at their SSU-rDNA locus, and were isolated from marine and 'slightly brackish' environments, respectively. Phylogenetic analyses of these two taxa plus closely related sequences from environmental surveys provide support for a novel clade of eukaryotes that is distinct from the major clades including the Opisthokonta, Excavata, Amoebozoa and 'SAR' (Stramenopile, Alveolate, Rhizaria).

Diversity and vertical distribution of microbial eukaryotes in the snow, sea ice and seawater near the north pole at the end of the polar night

Our knowledge about the microorganisms living in the high Arctic Ocean is still rudimentary compared to other oceans mostly because of logistical challenges imposed by its inhospitable climate and the presence of a multi-year ice cap. We have used 18S rRNA gene libraries to study the diversity of microbial eukaryotes in the upper part of the water column (0–170 m depth), the sea ice (0–1.5 m depth) and the overlying snow from samples collected in the vicinity of the North Pole (N88°35′, E015°59) at the very end of the long polar night. We detected very diverse eukaryotes belonging to Alveolata, Fungi, Amoebozoa, Viridiplantae, Metazoa, Rhizaria, Heterokonta, and Telonemia. Different alveolates (dinoflagellates and Marine Alveolate Groups I and II species) were the most abundant and diverse in gene libraries from water and sea ice, representing 80% of the total number of clones and operational taxonomic units. Only contaminants and/or species from continental ecosystems were detected in snow, suggesting wind- and animal- or human-mediated cosmopolitan dispersal of some taxa. By contrast, sea ice and seawater samples harbored a larger and more similar inter-sample protist diversity as compared with snow. The North Pole was found to harbor distinctive eukaryotic communities along the vertical gradient with an unparalleled diversity of core dinoflagellates, largely dominant in libraries from the water column, as compared to other oceanic locations. In contrast, phototrophic organisms typical of Arctic sea ice and plankton, such as diatoms and prasinophytes, were very rare in our samples. This was most likely due to a decrease of their populations after several months of polar night darkness and to the presence of rich populations of diverse grazers. Whereas strict phototrophs were scarce, we identified a variety of likely mixotrophic taxa, which supports the idea that mixotrophy may be important for the survival of diverse protists through the long polar night.

Delimiting operational taxonomic units for assessing ciliate environmental diversity using small-subunit rRNA gene sequences

Delineating operational taxonomic units (OTUs) is a central element in any culture-independent analysis of environmental microbial eukaryotic diversity. Previous studies either have not justified their choice in sequence distance used to bin small-subunit ribosomal RNA (SSU rRNA) gene sequences amplified from environmental samples into OTUs, or have used a value based on the average across a broad sampling of microbial eukaryotes. Here, we analyse distances (320 922 pairwise comparisons) among sequences just from identified ciliates, and compare these with their taxonomic hierarchy. Our results show that no single sequence similarity value can always and unambiguously delineate species boundaries and higher taxa. Nevertheless, we suggest the use of 98% similarity to delineate ciliate OTUs because this threshold at least accounts for intra-specific polymorphism among multiple rRNA cistron copies. However, we suggest refraining from reconciling SSU rRNA gene-based OTUs and ciliate morphotypes; these OTUs should be used to analyse ciliate phylotype diversity, not ciliate species diversity.

Identification of a divergent environmental DNA sequence clade using the phylogeny of gregarine parasites (Apicomplexa) from crustacean hosts

BACKGROUND

Environmental SSU rDNA surveys have significantly improved our understanding of microeukaryotic diversity. Many of the sequences acquired using this approach are closely related to lineages previously characterized at both morphological and molecular levels, making interpretation of these data relatively straightforward. Some sequences, by contrast, appear to be phylogenetic orphans and are sometimes inferred to represent "novel lineages" of unknown cellular identity. Consequently, interpretation of environmental DNA surveys of cellular diversity rely on an adequately comprehensive database of DNA sequences derived from identified species. Several major taxa of microeukaryotes, however, are still very poorly represented in these databases, and this is especially true for diverse groups of single-celled parasites, such as gregarine apicomplexans.

METHODOLOGY/PRINCIPAL FINDINGS

This study attempts to address this paucity of dna sequence data by characterizing four different gregarine species, isolated from the intestines of crustaceans, at both morphological and molecular levels: Thiriotia pugettiae sp. n. from the graceful kelp crab (Pugettia gracilis), Cephaloidophora cf. communis from two different species of barnacles (Balanus glandula and B. balanus), Heliospora cf. longissima from two different species of freshwater amphipods (Eulimnogammarus verrucosus and E. vittatus), and Heliospora caprellae comb. n. from a skeleton shrimp (Caprella alaskana). SSU rDNA sequences were acquired from isolates of these gregarine species and added to a global apicomplexan alignment containing all major groups of gregarines characterized so far. Molecular phylogenetic analyses of these data demonstrated that all of the gregarines collected from crustacean hosts formed a very strongly supported clade with 48 previously unidentified environmental DNA sequences.

CONCLUSIONS/SIGNIFICANCE

This expanded molecular phylogenetic context enabled us to establish a major clade of intestinal gregarine parasites and infer the cellular identities of several previously unidentified environmental SSU rDNA sequences, including several sequences that have formerly been discussed broadly in the literature as a suspected "novel" lineage of eukaryotes.

Protistan microbial observatory in the Cariaco Basin, Caribbean. II. Habitat specialization

This is the second paper in a series of three that investigates eukaryotic microbial diversity and taxon distribution in the Cariaco Basin, Venezuela, the ocean's largest anoxic marine basin. Here, we use phylogenetic information, multivariate community analyses and statistical richness predictions to test whether protists exhibit habitat specialization within defined geochemical layers of the water column. We also analyze spatio-temporal distributions of protists across two seasons and two geographic sites within the basin. Non-metric multidimensional scaling indicates that these two basin sites are inhabited by distinct protistan assemblages, an observation that is supported by the minimal overlap in observed and predicted richness of sampled sites. A comparison of parametric richness estimations indicates that protistan communities in closely spaced-but geochemically different-habitats are very dissimilar, and may share as few as 5% of total operational taxonomic units (OTUs). This is supported by a canonical correspondence analysis, indicating that the empirically observed OTUs are organized along opposing gradients in oxidants and reductants. Our phylogenetic analyses identify many new clades at species to class levels, some of which appear restricted to specific layers of the water column and have a significantly nonrandom distribution. These findings suggest many pelagic protists are restricted to specific habitats, and likely diversify, at least in part due to separation by geochemical barriers.

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.

A wide diversity of previously undetected free-living relatives of diplomonads isolated from marine/saline habitats

Over the last 15 years classical culturing and environmental PCR techniques have revealed a modest number of genuinely new major lineages of protists; however, some new groups have greatly influenced our understanding of eukaryote evolution. We used culturing techniques to examine the diversity of free-living protists that are relatives of diplomonads and retortamonads, a group of evolutionary and parasitological importance. Until recently, a single organism, Carpediemonas membranifera, was the only representative of this region of the tree. We report 18 new isolates of Carpediemonas-like organisms (CLOs) from anoxic marine sediments. Only one is a previously cultured species. Eleven isolates are conspecific and were classified within a new genus, Kipferlia n. gen. The remaining isolates include representatives of three other lineages that likely represent additional undescribed genera (at least). Small-subunit ribosomal RNA gene phylogenies show that CLOs form a cloud of six major clades basal to the diplomonad-retortamonad grouping (i.e. each of the six CLO clades is potentially as phylogenetically distinct as diplomonads and retortamonads). CLOs will be valuable for tracing the evolution of diplomonad cellular features, for example, their extremely reduced mitochondrial organelles. It is striking that the majority of CLO diversity was undetected by previous light microscopy surveys and environmental PCR studies, even though they inhabit a commonly sampled environment. There is no reason to assume this is a unique situation - it is likely that undersampling at the level of major lineages is still widespread for protists.

Microbial diversity of a Brazilian coastal region influenced by an upwelling system and anthropogenic activity

Background

Upwelling systems are characterised by an intense primary biomass production in the surface (warmest) water after the outcrop of the bottom (coldest) water, which is rich in nutrients. Although it is known that the microbial assemblage plays an important role in the food chain of marine systems and that the upwelling systems that occur in southwest Brazil drive the complex dynamics of the food chain, little is known about the microbial composition present in this region.

Methodology/Principal Findings

We carried out a molecular survey based on SSU rRNA gene from the three domains of the phylogenetic tree of life present in a tropical upwelling region (Arraial do Cabo, Rio de Janeiro, Brazil). The aim was to analyse the horizontal and vertical variations of the microbial composition in two geographically close areas influenced by anthropogenic activity (sewage disposal/port activity) and upwelling phenomena, respectively. A lower estimated diversity of microorganisms of the three domains of the phylogenetic tree of life was found in the water of the area influenced by anthropogenic activity compared to the area influenced by upwelling phenomena. We observed a heterogenic distribution of the relative abundance of taxonomic groups, especially in the Archaea and Eukarya domains. The bacterial community was dominated by Proteobacteria, Cyanobacteria and Bacteroidetes phyla, whereas the microeukaryotic community was dominated by Metazoa, Fungi, Alveolata and Stramenopile. The estimated archaeal diversity was the lowest of the three domains and was dominated by uncharacterised marine Crenarchaeota that were most closely related to Marine Group I.

Conclusions/Significance

The variety of conditions and the presence of different microbial assemblages indicated that the area of Arraial do Cabo can be used as a model for detailed studies that contemplate the correlation between pollution-indicating parameters and the depletion of microbial diversity in areas close to anthropogenic activity; functional roles and geochemical processes; phylogeny of the uncharacterised diversity; and seasonal variations of the microbial assemblages.

Are low temperature habitats hot spots of microbial evolution driven by viruses?

There is an increasing body of evidence to show that viruses are important drivers of microbial evolution and that they can store a great deal of the Earth's microbial diversity in their genomes. Examination of microbial diversity in polar regions has revealed a higher than expected diversity of viruses, bacteria and eukaryotic microbes. Further, the few available studies in polar regions reveal that viral control of microbial mortality is important in these habitats. In this opinion article, we argue that strong relationships between viruses and their hosts in a range of polar habitats could be key in explaining why polar regions are in fact hot spots of microbial diversity and evolution. Further, we argue that periodic glaciations, and particularly the Neoproterozoic low-latitude glaciation, known as 'snowball Earth', could have been periods of intense diversification in aquatic refuges.

Sequence diversity and novelty of natural assemblages of picoeukaryotes from the Indian Ocean

Despite the ecological importance of marine pico-size eukaryotes, the study of their in situ diversity using molecular tools started just a few years ago. These studies have revealed that marine picoeukaryotes are very diverse and include many novel taxa. However, the amount and structure of their phylogenetic diversity and the extent of their sequence novelty still remains poorly known, as a systematic analysis has been seldom attempted. In this study, we use a coherent and carefully curated data set of 500 published 18S ribosomal DNA sequences to quantify the diversity and novelty patterns of picoeukaryotes in the Indian Ocean. Our phylogenetic tree showed many distant lineages. We grouped sequences in OTUs (operational taxonomic units) at discrete values delineated by pair-wise Jukes-Cantor (JC) distances and tree patristic distances. At a distance of 0.01, the number of OTUs observed (237/242; using JC or patristic distances, respectively) was half the number of sequences analyzed, indicating the existence of microdiverse clusters of highly related sequences. At this distance level, we estimated 600-800 OTUs using several statistical methods. The number of OTUs observed was still substantial at higher distances (39/82 at 0.20 distance) suggesting a large diversity at high-taxonomic ranks. Most sequences were related to marine clones from other sites and many were distant to cultured organisms, highlighting the huge culturing gap within protists. The novelty analysis indicated the putative presence of pseudogenes and of truly novel high-rank phylogenetic lineages. The identified diversity and novelty patterns among marine picoeukaryotes are of great importance for understanding and interpreting their ecology and evolution.

Essential genes from Arctic bacteria used to construct stable, temperature-sensitive bacterial vaccines

All bacteria share a set of evolutionarily conserved essential genes that encode products that are required for viability. The great diversity of environments that bacteria inhabit, including environments at extreme temperatures, place adaptive pressure on essential genes. We sought to use this evolutionary diversity of essential genes to engineer bacterial pathogens to be stably temperature-sensitive, and thus useful as live vaccines. We isolated essential genes from bacteria found in the Arctic and substituted them for their counterparts into pathogens of mammals. We found that substitution of nine different essential genes from psychrophilic (cold-loving) bacteria into mammalian pathogenic bacteria resulted in strains that died below their normal-temperature growth limits. Substitution of three different psychrophilic gene orthologs of ligA, which encode NAD-dependent DNA ligase, resulted in bacterial strains that died at 33, 35, and 37 degrees C. One ligA gene was shown to render Francisella tularensis, Salmonella enterica, and Mycobacterium smegmatis temperature-sensitive, demonstrating that this gene functions in both Gram-negative and Gram-positive lineage bacteria. Three temperature-sensitive F. tularensis strains were shown to induce protective immunity after vaccination at a cool body site. About half of the genes that could be tested were unable to mutate to temperature-resistant forms at detectable levels. These results show that psychrophilic essential genes can be used to create a unique class of bacterial temperature-sensitive vaccines for important human pathogens, such as S. enterica and Mycobacterium tuberculosis.

Spatio-temporal variations in protistan communities along an O/HS gradient in the anoxic Framvaren Fjord (Norway)

Despite its relevance for ecology and biodiversity, the stability of spatial microeukaryote diversity patterns in time has received only little attention using gene-based strategies, and there is little knowledge about the relation of spatial vs. temporal variation. We addressed this subject by investigating seasonal fluctuations in protistan communities in three ecologically distinct marine habitats. We analyzed 3360 eukaryote small subunit rRNA gene sequences collected along an O(2)/H(2)S gradient in a Norwegian fjord in order to reveal shifts in protistan community composition and structure in three different seasons. In all nine clone libraries, ciliates and stramenopiles accounted for the largest proportion. Yet, as expected, at the phylotype level, the protistan communities from distinct habitats differed significantly, with the number of shared phylotypes between two habitats being as low as 18%. This confirmed previous notions that environmental factors along the stratification gradient shape biodiversity patterns. Surprisingly, the intrahabitat community composition and structure varied at a comparable order of magnitude over time, with only 18-28% phylotypes shared within the same habitat. Our study demonstrates that the consideration of local fluctuations in microeukaryote diversity over time offers additional information for diversity surveys and can significantly contribute to the revelation of spatial protistan community patterns.

Massively parallel tag sequencing reveals the complexity of anaerobic marine protistan communities

Background

Recent advances in sequencing strategies make possible unprecedented depth and scale of sampling for molecular detection of microbial diversity. Two major paradigm-shifting discoveries include the detection of bacterial diversity that is one to two orders of magnitude greater than previous estimates, and the discovery of an exciting 'rare biosphere' of molecular signatures ('species') of poorly understood ecological significance. We applied a high-throughput parallel tag sequencing (454 sequencing) protocol adopted for eukaryotes to investigate protistan community complexity in two contrasting anoxic marine ecosystems (Framvaren Fjord, Norway; Cariaco deep-sea basin, Venezuela). Both sampling sites have previously been scrutinized for protistan diversity by traditional clone library construction and Sanger sequencing. By comparing these clone library data with 454 amplicon library data, we assess the efficiency of high-throughput tag sequencing strategies. We here present a novel, highly conservative bioinformatic analysis pipeline for the processing of large tag sequence data sets.

Results

The analyses of ca. 250,000 sequence reads revealed that the number of detected Operational Taxonomic Units (OTUs) far exceeded previous richness estimates from the same sites based on clone libraries and Sanger sequencing. More than 90% of this diversity was represented by OTUs with less than 10 sequence tags. We detected a substantial number of taxonomic groups like Apusozoa, Chrysomerophytes, Centroheliozoa, Eustigmatophytes, hyphochytriomycetes, Ichthyosporea, Oikomonads, Phaeothamniophytes, and rhodophytes which remained undetected by previous clone library-based diversity surveys of the sampling sites. The most important innovations in our newly developed bioinformatics pipeline employ (i) BLASTN with query parameters adjusted for highly variable domains and a complete database of public ribosomal RNA (rRNA) gene sequences for taxonomic assignments of tags; (ii) a clustering of tags at k differences (Levenshtein distance) with a newly developed algorithm enabling very fast OTU clustering for large tag sequence data sets; and (iii) a novel parsing procedure to combine the data from individual analyses.

Conclusion

Our data highlight the magnitude of the under-sampled 'protistan gap' in the eukaryotic tree of life. This study illustrates that our current understanding of the ecological complexity of protist communities, and of the global species richness and genome diversity of protists, is severely limited. Even though 454 pyrosequencing is not a panacea, it allows for more comprehensive insights into the diversity of protistan communities, and combined with appropriate statistical tools, enables improved ecological interpretations of the data and projections of global diversity.

Fungal diversity in oxygen-depleted regions of the Arabian Sea revealed by targeted environmental sequencing combined with cultivation

In order to study fungal diversity in oxygen minimum zones of the Arabian Sea, we analyzed 1440 cloned small subunit rRNA gene (18S rRNA gene) sequences obtained from environmental samples using three different PCR primer sets. Restriction fragment length polymorphism (RFLP) analyses yielded 549 distinct RFLP patterns, 268 of which could be assigned to fungi (Dikarya and zygomycetes) after sequence analyses. The remaining 281 RFLP patterns represented a variety of nonfungal taxa, even when using putatively fungal-specific primers. A substantial number of fungal sequences were closely related to environmental sequences from a range of other anoxic marine habitats, but distantly related to known sequences of described fungi. Community similarity analyses suggested distinctively different structures of fungal communities from normoxic sites, seasonally anoxic sites and permanently anoxic sites, suggesting different adaptation strategies of fungal communities to prevailing oxygen conditions. Additionally, we obtained 26 fungal cultures from the study sites, most of which were closely related (>97% sequence similarity) to well-described Dikarya. This indicates that standard cultivation mainly produces more of what is already known. However, two of these cultures were highly divergent to known sequences and seem to represent novel fungal groups on high taxonomic levels. Interestingly, none of the cultured isolates is identical to any of the environmental sequences obtained. Our study demonstrates the importance of a multiple-primer approach combined with cultivation to obtain deeper insights into the true fungal diversity in environmental samples and to enable adequate intersample comparisons of fungal communities.

Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencing-generated fragments in 16S rRNA gene-based environmental surveys

Pyrosequencing-based 16S rRNA gene surveys are increasingly utilized to study highly diverse bacterial communities, with special emphasis on utilizing the large number of sequences obtained (tens to hundreds of thousands) for species richness estimation. However, it is not yet clear how the number of operational taxonomic units (OTUs) and, hence, species richness estimates determined using shorter fragments at different taxonomic cutoffs correlates with the number of OTUs assigned using longer, nearly complete 16S rRNA gene fragments. We constructed a 16S rRNA clone library from an undisturbed tallgrass prairie soil (1,132 clones) and used it to compare species richness estimates obtained using eight pyrosequencing candidate fragments (99 to 361 bp in length) and the nearly full-length fragment. Fragments encompassing the V1 and V2 (V1+V2) region and the V6 region (generated using primer pairs 8F-338R and 967F-1046R) overestimated species richness; fragments encompassing the V3, V7, and V7+V8 hypervariable regions (generated using primer pairs 338F-530R, 1046F-1220R, and 1046F-1392R) underestimated species richness; and fragments encompassing the V4, V5+V6, and V6+V7 regions (generated using primer pairs 530F-805R, 805F-1046R, and 967F-1220R) provided estimates comparable to those obtained with the nearly full-length fragment. These patterns were observed regardless of the alignment method utilized or the parameter used to gauge comparative levels of species richness (number of OTUs observed, slope of scatter plots of pairwise distance values for short and nearly complete fragments, and nonparametric and parametric species richness estimates). Similar results were obtained when analyzing three other datasets derived from soil, adult Zebrafish gut, and basaltic formations in the East Pacific Rise. Regression analysis indicated that these observed discrepancies in species richness estimates within various regions could readily be explained by the proportions of hypervariable, variable, and conserved base pairs within an examined fragment.

Parametric models for estimating the number of classes

We consider parametric distributions intended to model heterogeneity in population size estimation, especially parametric stochastic abundance models for species richness estimation. We briefly review (conditional) maximum likelihood estimation of the number of species, and summarize the results of fitting 7 candidate models to frequency-count data, from a database of >40000 such instances, mostly arising from microbial ecology. We consider error estimation, goodness-of-fit assessment, data subsetting, and other practical matters. We find that, although the array of candidate models can be improved, finite mixtures of a small number of components (point masses or simple diffuse distributions) represent a promising direction. Finally we consider the connections between parametric models for abundance and incidence data, again noting the usefulness of finite mixture models.