A Clustering Optimization Strategy for Molecular Taxonomy Applied to Planktonic Foraminifera SSU rDNA

Species composition of arbuscular mycorrhizal communities changes with elevation in the Andes of South Ecuador

Arbuscular mycorrhizal fungi (AMF) are the most prominent mycobionts of plants in the tropics, yet little is known about their diversity, species compositions and factors driving AMF distribution patterns. To investigate whether elevation and associated vegetation type affect species composition, we sampled 646 mycorrhizal samples in locations between 1000 and 4000 m above sea level (masl) in the South of Ecuador. We estimated diversity, distribution and species compositions of AMF by cloning and Sanger sequencing the 18S rDNA (the section between AML1 and AML2) and subsequent derivation of fungal OTUs based on 99% sequence similarity. In addition, we analyzed the phylogenetic structure of the sites by computing the mean pairwise distance (MPD) and the mean nearest taxon difference (MNTD) for each elevation level. It revealed that AMF species compositions at 1000 and 2000 masl differ from 3000 and 4000 masl. Lower elevations (1000 and 2000 masl) were dominated by members of Glomeraceae, whereas Acaulosporaceae were more abundant in higher elevations (3000 and 4000 masl). Ordination of OTUs with respect to study sites revealed a correlation to elevation with a continuous turnover of species from lower to higher elevations. Most of the abundant OTUs are not endemic to South Ecuador. We also found a high proportion of rare OTUs at all elevations: 79-85% of OTUs occurred in less than 5% of the samples. Phylogenetic community analysis indicated clustering and evenness for most elevation levels indicating that both, stochastic processes and habitat filtering are driving factors of AMF community compositions.

TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy

Microbial taxonomy is increasingly influenced by genome-based computational methods. Yet such analyses can be complex and require expert knowledge. Here we introduce TYGS, the Type (Strain) Genome Server, a user-friendly high-throughput web server for genome-based prokaryote taxonomy, connected to a large, continuously growing database of genomic, taxonomic and nomenclatural information. It infers genome-scale phylogenies and state-of-the-art estimates for species and subspecies boundaries from user-defined and automatically determined closest type genome sequences. TYGS also provides comprehensive access to nomenclature, synonymy and associated taxonomic literature. Clinically important examples demonstrate how TYGS can yield new insights into microbial classification, such as evidence for a species-level separation of previously proposed subspecies of Salmonella enterica. TYGS is an integrated approach for the classification of microbes that unlocks novel scientific approaches to microbiologists worldwide and is particularly helpful for the rapidly expanding field of genome-based taxonomic descriptions of new genera, species or subspecies.

Surface ocean metabarcoding confirms limited diversity in planktonic foraminifera but reveals unknown hyper-abundant lineages

Since the advent of DNA metabarcoding surveys, the planktonic realm is considered a treasure trove of diversity, inhabited by a small number of abundant taxa, and a hugely diverse and taxonomically uncharacterized consortium of rare species. Here we assess if the apparent underestimation of plankton diversity applies universally. We target planktonic foraminifera, a group of protists whose known morphological diversity is limited, taxonomically resolved and linked to ribosomal DNA barcodes. We generated a pyrosequencing dataset of ~100,000 partial 18S rRNA foraminiferal sequences from 32 size fractioned photic-zone plankton samples collected at 8 stations in the Indian and Atlantic Oceans during the Tara Oceans expedition (2009–2012). We identified 69 genetic types belonging to 41 morphotaxa in our metabarcoding dataset. The diversity saturated at local and regional scale as well as in the three size fractions and the two depths sampled indicating that the diversity of foraminifera is modest and finite. The large majority of the newly discovered lineages occur in the small size fraction, neglected by classical taxonomy. These unknown lineages dominate the bulk [>0.8 µm] size fraction, implying that a considerable part of the planktonic foraminifera community biomass has its origin in unknown lineages.

Alignment-free sequence comparison: benefits, applications, and tools

Alignment-free sequence analyses have been applied to problems ranging from whole-genome phylogeny to the classification of protein families, identification of horizontally transferred genes, and detection of recombined sequences. The strength of these methods makes them particularly useful for next-generation sequencing data processing and analysis. However, many researchers are unclear about how these methods work, how they compare to alignment-based methods, and what their potential is for use for their research. We address these questions and provide a guide to the currently available alignment-free sequence analysis tools.

Nomenclature for the Nameless: A Proposal for an Integrative Molecular Taxonomy of Cryptic Diversity Exemplified by Planktonic Foraminifera

Investigations of biodiversity, biogeography, and ecological processes rely on the identification of "species" as biologically significant, natural units of evolution. In this context, morphotaxonomy only provides an adequate level of resolution if reproductive isolation matches morphological divergence. In many groups of organisms, morphologically defined species often disguise considerable genetic diversity, which may be indicative of the existence of cryptic species. The diversity hidden by morphological species can be disentangled through genetic surveys, which also provide access to data on the ecological distribution of genetically circumscribed units. These units can be identified by unique DNA sequence motifs and allow studies of evolutionary and ecological processes at different levels of divergence. However, the nomenclature of genetically circumscribed units within morphological species is not regulated and lacks stability. This represents a major obstacle to efforts to synthesize and communicate data on genetic diversity for multiple stakeholders. We have been confronted with such an obstacle in our work on planktonic foraminifera, where the stakeholder community is particularly diverse, involving geochemists, paleoceanographers, paleontologists, and biologists, and the lack of stable nomenclature beyond the level of formal morphospecies prevents effective transfer of knowledge. To circumvent this problem, we have designed a stable, reproducible, and flexible nomenclature system for genetically circumscribed units, analogous to the principles of a formal nomenclature system. Our system is based on the definition of unique DNA sequence motifs collocated within an individual, their typification (in analogy with holotypes), utilization of their hierarchical phylogenetic structure to define levels of divergence below that of the morphospecies, and a set of nomenclature rules assuring stability. The resulting molecular operational taxonomic units remain outside the domain of current nomenclature codes, but are linked to formal morphospecies as regulated by the codes. Subsequently, we show how this system can be applied to classify genetically defined units using the SSU rDNA marker in planktonic foraminifera and we highlight its potential use for other groups of organisms where similarly high levels of connectivity between molecular and formal taxonomies can be achieved.

Divergence Times and Phylogenetic Patterns of Sebacinales, a Highly Diverse and Widespread Fungal Lineage

Patterns of geographic distribution and composition of fungal communities are still poorly understood. Widespread occurrence in terrestrial ecosystems and the unique richness of interactions of Sebacinales with plants make them a target group to study evolutionary events in the light of nutritional lifestyle. We inferred diversity patterns, phylogenetic structures and divergence times of Sebacinales with respect to their nutritional lifestyles by integrating data from fossil-calibrated phylogenetic analyses. Relaxed molecular clock analyses indicated that Sebacinales originated late Permian within Basidiomycota, and their split into Sebacinaceae and Serendipitaceae nom. prov. likely occurred during the late Jurassic and the early Cretaceous, coinciding with major diversifications of land plants. In Sebacinaceae, diversification of species with ectomycorrhizal lifestyle presumably started during the Paleocene. Lineage radiations of the core group of ericoid and cavendishioid mycorrhizal Sebacinales started probably in the Eocene, coinciding with diversification events of their hosts. The diversification of Sebacinales with jungermannioid interactions started during the Oligocene, and occurred much later than the diversification of their hosts. Sebacinales communities associated either with ectomycorrhizal plants, achlorophyllous orchids, ericoid and cavendishioid Ericaceae or liverworts were phylogenetically clustered and globally distributed. Major Sebacinales lineage diversifications started after the continents had drifted apart. We also briefly discuss dispersal patterns of extant Sebacinales.

Phylogenetic classification of yeasts and related taxa within Pucciniomycotina

Most small genera containing yeast species in the Pucciniomycotina (Basidiomycota, Fungi) are monophyletic, whereas larger genera including Bensingtonia, Rhodosporidium, Rhodotorula, Sporidiobolus and Sporobolomyces are polyphyletic. With the implementation of the “One Fungus = One Name” nomenclatural principle these polyphyletic genera were revised. Nine genera, namely Bannoa, Cystobasidiopsis, Colacogloea, Kondoa, Erythrobasidium, Rhodotorula, Sporobolomyces, Sakaguchia and Sterigmatomyces, were emended to include anamorphic and teleomorphic species based on the results obtained by a multi-gene phylogenetic analysis, phylogenetic network analyses, branch length-based methods, as well as morphological, physiological and biochemical comparisons. A new class Spiculogloeomycetes is proposed to accommodate the order Spiculogloeales. The new families Buckleyzymaceae with Buckleyzyma gen. nov., Chrysozymaceae with Chrysozyma gen. nov., Microsporomycetaceae with Microsporomyces gen. nov., Ruineniaceae with Ruinenia gen. nov., Symmetrosporaceae with Symmetrospora gen. nov., Colacogloeaceae and Sakaguchiaceae are proposed. The new genera Bannozyma, Buckleyzyma, Fellozyma, Hamamotoa, Hasegawazyma, Jianyunia, Rhodosporidiobolus, Oberwinklerozyma, Phenoliferia, Pseudobensingtonia, Pseudohyphozyma, Sampaiozyma, Slooffia, Spencerozyma, Trigonosporomyces, Udeniozyma, Vonarxula, Yamadamyces and Yunzhangia are proposed to accommodate species segregated from the genera Bensingtonia, Rhodosporidium, Rhodotorula, Sporidiobolus and Sporobolomyces. Ballistosporomyces is emended and reintroduced to include three Sporobolomyces species of the sasicola clade. A total of 111 new combinations are proposed in this study.

Towards an integrated phylogenetic classification of the Tremellomycetes

Families and genera assigned to Tremellomycetes have been mainly circumscribed by morphology and for the yeasts also by biochemical and physiological characteristics. This phenotype-based classification is largely in conflict with molecular phylogenetic analyses. Here a phylogenetic classification framework for the Tremellomycetes is proposed based on the results of phylogenetic analyses from a seven-genes dataset covering the majority of tremellomycetous yeasts and closely related filamentous taxa. Circumscriptions of the taxonomic units at the order, family and genus levels recognised were quantitatively assessed using the phylogenetic rank boundary optimisation (PRBO) and modified general mixed Yule coalescent (GMYC) tests. In addition, a comprehensive phylogenetic analysis on an expanded LSU rRNA (D1/D2 domains) gene sequence dataset covering as many as available teleomorphic and filamentous taxa within Tremellomycetes was performed to investigate the relationships between yeasts and filamentous taxa and to examine the stability of undersampled clades. Based on the results inferred from molecular data and morphological and physiochemical features, we propose an updated classification for the Tremellomycetes. We accept five orders, 17 families and 54 genera, including seven new families and 18 new genera. In addition, seven families and 17 genera are emended and one new species name and 185 new combinations are proposed. We propose to use the term pro tempore or pro tem. in abbreviation to indicate the species names that are temporarily maintained.

Toward accurate molecular identification of species in complex environmental samples: testing the performance of sequence filtering and clustering methods

Metabarcoding has the potential to become a rapid, sensitive, and effective approach for identifying species in complex environmental samples. Accurate molecular identification of species depends on the ability to generate operational taxonomic units (OTUs) that correspond to biological species. Due to the sometimes enormous estimates of biodiversity using this method, there is a great need to test the efficacy of data analysis methods used to derive OTUs. Here, we evaluate the performance of various methods for clustering length variable 18S amplicons from complex samples into OTUs using a mock community and a natural community of zooplankton species. We compare analytic procedures consisting of a combination of (1) stringent and relaxed data filtering, (2) singleton sequences included and removed, (3) three commonly used clustering algorithms (mothur, UCLUST, and UPARSE), and (4) three methods of treating alignment gaps when calculating sequence divergence. Depending on the combination of methods used, the number of OTUs varied by nearly two orders of magnitude for the mock community (60–5068 OTUs) and three orders of magnitude for the natural community (22–22191 OTUs). The use of relaxed filtering and the inclusion of singletons greatly inflated OTU numbers without increasing the ability to recover species. Our results also suggest that the method used to treat gaps when calculating sequence divergence can have a great impact on the number of OTUs. Our findings are particularly relevant to studies that cover taxonomically diverse species and employ markers such as rRNA genes in which length variation is extensive.

Metabarcoding vs. morphological identification to assess diatom diversity in environmental studies

Diatoms are frequently used for water quality assessments; however, identification to species level is difficult, time-consuming and needs in-depth knowledge of the organisms under investigation, as nonhomoplastic species-specific morphological characters are scarce. We here investigate how identification methods based on DNA (metabarcoding using NGS platforms) perform in comparison to morphological diatom identification and propose a workflow to optimize diatom fresh water quality assessments. Diatom diversity at seven different sites along the course of the river system Odra and Lusatian Neisse from the source to the mouth is analysed with DNA and morphological methods, which are compared. The NGS technology almost always leads to a higher number of identified taxa (270 via NGS vs. 103 by light microscopy LM), whose presence could subsequently be verified by LM. The sequence-based approach allows for a much more graduated insight into the taxonomic diversity of the environmental samples. Taxa retrieval varies considerably throughout the river system, depending on species occurrences and the taxonomic depth of the reference databases. Mostly rare taxa from oligotrophic parts of the river systems are less well represented in the reference database used. A workflow for DNA-based NGS diatom identification is presented. 28 000 diatom sequences were evaluated. Our findings provide evidence that metabarcoding of diatoms via NGS sequencing of the V4 region (18S) has a great potential for water quality assessments and could complement and maybe even improve the identification via light microscopy.

SSU rDNA divergence in planktonic foraminifera: molecular taxonomy and biogeographic implications

The use of planktonic foraminifera in paleoceanography requires taxonomic consistency and precise assessment of the species biogeography. Yet, ribosomal small subunit (SSUr) DNA analyses have revealed that most of the modern morpho-species of planktonic foraminifera are composed of a complex of several distinct genetic types that may correspond to cryptic or pseudo-cryptic species. These genetic types are usually delimitated using partial sequences located at the 3'end of the SSUrDNA, but typically based on empirical delimitation. Here, we first use patristic genetic distances calculated within and among genetic types of the most common morpho-species to show that intra-type and inter-type genetic distances within morpho-species may significantly overlap, suggesting that genetic types have been sometimes inconsistently defined. We further apply two quantitative and independent methods, ABGD (Automatic Barcode Gap Detection) and GMYC (General Mixed Yule Coalescent) to a dataset of published and newly obtained partial SSU rDNA for a more objective assessment of the species status of these genetic types. Results of these complementary approaches are highly congruent and lead to a molecular taxonomy that ranks 49 genetic types of planktonic foraminifera as genuine (pseudo)cryptic species. Our results advocate for a standardized sequencing procedure allowing homogenous delimitations of (pseudo)cryptic species. On the ground of this revised taxonomic framework, we finally provide an integrative taxonomy synthesizing geographic, ecological and morphological differentiations that can occur among the genuine (pseudo)cryptic species. Due to molecular, environmental or morphological data scarcities, many aspects of our proposed integrative taxonomy are not yet fully resolved. On the other hand, our study opens up the potential for a correct interpretation of environmental sequence datasets.

Phylogeography of the tropical planktonic foraminifera lineage globigerinella reveals isolation inconsistent with passive dispersal by ocean currents

Morphologically defined species of marine plankton often harbor a considerable level of cryptic diversity. Since many morphospecies show cosmopolitan distribution, an understanding of biogeographic and evolutionary processes at the level of genetic diversity requires global sampling. We use a database of 387 single-specimen sequences of the SSU rDNA of the planktonic foraminifera Globigerinella as a model to assess the biogeographic and phylogenetic distributions of cryptic diversity in marine microplankton on a global scale. Our data confirm the existence of multiple, well isolated genetic lineages. An analysis of their abundance and distribution indicates that our sampling is likely to approximate the actual total diversity. Unexpectedly, we observe an uneven allocation of cryptic diversity among the phylogenetic lineages. We show that this pattern is neither an artifact of sampling intensity nor a function of lineage age. Instead, we argue that it reflects an ongoing speciation process in one of the three major lineages. Surprisingly, four of the six genetic types in the hyperdiverse lineage are biogeographically restricted to the Indopacific. Their mutual co-occurrence and their hierarchical phylogenetic structure provide no evidence for an origin through sudden habitat fragmentation and their limitation to the Indopacific challenges the view of a global gene flow within the warm-water provinces. This phenomenon shows that passive dispersal is not sufficient to describe the distribution of plankton diversity. Rather, these organisms show differentiated distribution patterns shaped by species interactions and reflecting phylogenetic contingency with unique histories of diversification rates.

New metrics for comparison of taxonomies reveal striking discrepancies among species delimitation methods in Madascincus lizards

Delimiting and describing species is fundamental to numerous biological disciplines such as evolution, macroecology, and conservation. Delimiting species as independent evolutionary lineages may and often does yield different outcomes depending on the species criteria applied, but methods should be chosen that minimize the inference of objectively erroneous species limits. Several protocols exploit single-gene or multi-gene coalescence statistics, assignment tests or other rationales related to nuclear DNA (nDNA) allele sharing to automatically delimit species. We apply seven different species delimitation protocols to a taxonomically confusing group of Malagasy lizards (Madascincus), and compare the resulting taxonomies with two newly developed metrics: the Taxonomic index of congruence C tax which quantifies the congruence between two taxonomies, and the Relative taxonomic resolving power index R tax which quantifies the potential of an approach to capture a high number of species boundaries. The protocols differed in the total number of species proposed, between 9 and 34, and were also highly incongruent in placing species boundaries. The Generalized Mixed Yule-Coalescent approach captured the highest number of potential species boundaries but many of these were clearly contradicted by extensive nDNA admixture between sympatric mitochondrial DNA (mtDNA) haplotype lineages. Delimiting species as phenotypically diagnosable mtDNA clades failed to detect two cryptic species that are unambiguous due to a lack of nDNA gene flow despite sympatry. We also consider the high number of species boundaries and their placement by multi-gene Bayesian species delimitation as poorly reliable whereas the Bayesian assignment test approach provided a species delimitation highly congruent with integrative taxonomic practice. The present study illustrates the trade-off in taxonomy between reliability (favored by conservative approaches) and resolving power (favored by inflationist approaches). Quantifying excessive splitting is more difficult than quantifying excessive lumping, suggesting a priority for conservative taxonomies in which errors are more liable to be detected and corrected by subsequent studies.

Genome sequence-based species delimitation with confidence intervals and improved distance functions

BACKGROUND

For the last 25 years species delimitation in prokaryotes (Archaea and Bacteria) was to a large extent based on DNA-DNA hybridization (DDH), a tedious lab procedure designed in the early 1970s that served its purpose astonishingly well in the absence of deciphered genome sequences. With the rapid progress in genome sequencing time has come to directly use the now available and easy to generate genome sequences for delimitation of species. GBDP (Genome Blast Distance Phylogeny) infers genome-to-genome distances between pairs of entirely or partially sequenced genomes, a digital, highly reliable estimator for the relatedness of genomes. Its application as an in-silico replacement for DDH was recently introduced. The main challenge in the implementation of such an application is to produce digital DDH values that must mimic the wet-lab DDH values as close as possible to ensure consistency in the Prokaryotic species concept.

RESULTS

Correlation and regression analyses were used to determine the best-performing methods and the most influential parameters. GBDP was further enriched with a set of new features such as confidence intervals for intergenomic distances obtained via resampling or via the statistical models for DDH prediction and an additional family of distance functions. As in previous analyses, GBDP obtained the highest agreement with wet-lab DDH among all tested methods, but improved models led to a further increase in the accuracy of DDH prediction. Confidence intervals yielded stable results when inferred from the statistical models, whereas those obtained via resampling showed marked differences between the underlying distance functions.

CONCLUSIONS

Despite the high accuracy of GBDP-based DDH prediction, inferences from limited empirical data are always associated with a certain degree of uncertainty. It is thus crucial to enrich in-silico DDH replacements with confidence-interval estimation, enabling the user to statistically evaluate the outcomes. Such methodological advancements, easily accessible through the web service at http://ggdc.dsmz.de, are crucial steps towards a consistent and truly genome sequence-based classification of microorganisms.

Vertical niche partitioning between cryptic sibling species of a cosmopolitan marine planktonic protist

A large portion of the surface-ocean biomass is represented by microscopic unicellular plankton. These organisms are functionally and morphologically diverse, but it remains unclear how their diversity is generated. Species of marine microplankton are widely distributed because of passive transport and lack of barriers in the ocean. How does speciation occur in a system with a seemingly unlimited dispersal potential? Recent studies using planktonic foraminifera as a model showed that even among the cryptic genetic diversity within morphological species, many genetic types are cosmopolitan, lending limited support for speciation by geographical isolation. Here we show that the current two-dimensional view on the biogeography and potential speciation mechanisms in the microplankton may be misleading. By depth-stratified sampling, we present evidence that sibling genetic types in a cosmopolitan species of marine microplankton, the planktonic foraminifer Hastigerina pelagica, are consistently separated by depth throughout their global range. Such strong separation between genetically closely related and morphologically inseparable genetic types indicates that niche partitioning in marine heterotrophic microplankton can be maintained in the vertical dimension on a global scale. These observations indicate that speciation along depth (depth-parapatric speciation) can occur in vertically structured microplankton populations, facilitating diversification without the need for spatial isolation.

Species delimitation and global biosecurity

Species delimitation directly impacts on global biosecurity. It is a critical element in the decisions made by national governments in regard to the flow of trade and to the biosecurity measures imposed to protect countries from the threat of invasive species. Here we outline a novel approach to species delimitation, “tip to root”, for two highly invasive insect pests, Bemisia tabaci (sweetpotato whitefly) and Lymantria dispar (Asian gypsy moth). Both species are of concern to biosecurity, but illustrate the extremes of phylogenetic resolution that present the most complex delimitation issues for biosecurity; B. tabaci having extremely high intra-specific genetic variability and L. dispar composed of relatively indistinct subspecies. This study tests a series of analytical options to determine their applicability as tools to provide more rigorous species delimitation measures and consequently more defensible species assignments and identification of unknowns for biosecurity. Data from established DNA barcode datasets (COI), which are becoming increasingly considered for adoption in biosecurity, were used here as an example. The analytical approaches included the commonly used Kimura two-parameter (K2P) inter-species distance plus four more stringent measures of taxon distinctiveness, (1) Rosenberg’s reciprocal monophyly, (P(AB)),1 (2) Rodrigo’s (P(randomly distinct)),2 (3) genealogical sorting index, (gsi),3 and (4) General mixed Yule-coalescent (GMYC).4,5 For both insect datasets, a comparative analysis of the methods revealed that the K2P distance method does not capture the same level of species distinctiveness revealed by the other three measures; in B. tabaci there are more distinct groups than previously identified using the K2P distances and for L. dipsar far less variation is apparent within the predefined subspecies. A consensus for the results from P(AB), P(randomly distinct) and gsi offers greater statistical confidence as to where genetic limits might be drawn. In the species cases here, the results clearly indicate that there is a need for more gene sampling to substantiate either the new cohort of species indicated for B. tabaci or to detect the established subspecies taxonomy of L. dispar. Given the ease of use through the Geneious species delimitation plugins, similar analysis of such multi-gene datasets would be easily accommodated. Overall, the tip to root approach described here is recommended where careful consideration of species delimitation is required to support crucial biosecurity decisions based on accurate species identification.

Species delimitation in taxonomically difficult fungi: the case of Hymenogaster

BACKGROUND

False truffles are ecologically important as mycorrhizal partners of trees and evolutionarily highly interesting as the result of a shift from epigeous mushroom-like to underground fruiting bodies. Since its first description by Vittadini in 1831, inappropriate species concepts in the highly diverse false truffle genus Hymenogaster has led to continued confusion, caused by a large variety of prevailing taxonomical opinions.

METHODOLOGY

In this study, we reconsidered the species delimitations in Hymenogaster based on a comprehensive collection of Central European taxa comprising more than 140 fruiting bodies from 20 years of field work. The ITS rDNA sequence dataset was subjected to phylogenetic analysis as well as clustering optimization using OPTSIL software.

CONCLUSIONS

Among distinct species concepts from the literature used to create reference partitions for clustering optimization, the broadest concept resulted in the highest agreement with the ITS data. Our results indicate a highly variable morphology of H. citrinus and H. griseus, most likely linked to environmental influences on the phenology (maturity, habitat, soil type and growing season). In particular, taxa described in the 19(th) century frequently appear as conspecific. Conversely, H. niveus appears as species complex comprising seven cryptic species with almost identical macro- and micromorphology. H. intermedius and H. huthii are described as novel species, each of which with a distinct morphology intermediate between two species complexes. A revised taxonomy for one of the most taxonomically difficult genera of Basidiomycetes is proposed, including an updated identification key. The (semi-)automated selection among species concepts used here is of importance for the revision of taxonomically problematic organism groups in general.