Barcoding Against a Paradox? Combined Molecular Species Delineations Reveal Multiple Cryptic Lineages in Elusive Meiofaunal
Total Page:16
File Type:pdf, Size:1020Kb
Jörger et al. BMC Evolutionary Biology 2012, 12:245 http://www.biomedcentral.com/1471-2148/12/245 RESEARCH ARTICLE Open Access Barcoding against a paradox? Combined molecular species delineations reveal multiple cryptic lineages in elusive meiofaunal sea slugs Katharina M Jörger1,2*, Jon L Norenburg3, Nerida G Wilson4 and Michael Schrödl1,2 Abstract Background: Many marine meiofaunal species are reported to have wide distributions, which creates a paradox considering their hypothesized low dispersal abilities. Correlated with this paradox is an especially high taxonomic deficit for meiofauna, partly related to a lower taxonomic effort and partly to a high number of putative cryptic species. Molecular-based species delineation and barcoding approaches have been advocated for meiofaunal biodiversity assessments to speed up description processes and uncover cryptic lineages. However, these approaches show sensitivity to sampling coverage (taxonomic and geographic) and the success rate has never been explored on mesopsammic Mollusca. Results: We collected the meiofaunal sea-slug Pontohedyle (Acochlidia, Heterobranchia) from 28 localities worldwide. With a traditional morphological approach, all specimens fall into two morphospecies. However, with a multi-marker genetic approach, we reveal multiple lineages that are reciprocally monophyletic on single and concatenated gene trees in phylogenetic analyses. These lineages are largely concordant with geographical and oceanographic parameters, leading to our primary species hypothesis (PSH). In parallel, we apply four independent methods of molecular based species delineation: General Mixed Yule Coalescent model (GMYC), statistical parsimony, Bayesian Species Delineation (BPP) and Automatic Barcode Gap Discovery (ABGD). The secondary species hypothesis (SSH) is gained by relying only on uncontradicted results of the different approaches (‘minimum consensus approach’), resulting in the discovery of a radiation of (at least) 12 mainly cryptic species, 9 of them new to science, some sympatric and some allopatric with respect to ocean boundaries. However, the meiofaunal paradox still persists in some Pontohedyle species identified here with wide coastal and trans-archipelago distributions. Conclusions: Our study confirms extensive, morphologically cryptic diversity among meiofauna and accentuates the taxonomic deficit that characterizes meiofauna research. We observe for Pontohedyle slugs a high degree of morphological simplicity and uniformity, which we expect might be a general rule for meiofauna. To tackle cryptic diversity in little explored and hard-to-sample invertebrate taxa, at present, a combined approach seems most promising, such as multi-marker-barcoding (i.e., molecular systematics using mitochondrial and nuclear markers and the criterion of reciprocal monophyly) combined with a minimum consensus approach across independent methods of molecular species delineation to define candidate species. * Correspondence: [email protected] 1Mollusca Department, Bavarian State Collection of Zoology, Münchhausenstr.21, 81247, München, Germany 2Department Biology II, Ludwig-Maximilians-Universität, BioZentrum Martinsried, Großhadernerstr. 2, 82152, Planegg-Martinsried, Germany Full list of author information is available at the end of the article © 2012 Jörger et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Jörger et al. BMC Evolutionary Biology 2012, 12:245 Page 2 of 18 http://www.biomedcentral.com/1471-2148/12/245 Background possible inadvertent inclusion of nuclear mitochondrial Sediment-associated marine meiofaunal organisms in- pseudogenes (= nonfunctional copies of mtDNA in the habit one of the largest ecosystems on earth – sediment- nucleus, or numts) [19], and other mitochondria-specific covered ocean floors and beaches – and comprise a pitfalls such as reduced effective population size or incon- major part of marine biodiversity [1]. However, only a sistent recombination [20]. The risk of incorrect species small fraction of the predicted species richness currently delineation due to incomplete lineage sorting or intro- is known to science [1-4] and recent surveys have shown gression can be reduced by analyzing independent loci a high number of new, undescribed species even in well- [21], which is generally considered superior to single- studied areas (see [4]). Minute body sizes often prohibit gene approaches [22]. We chose a barcoding approach direct visual identification in the field; instead, morpho- based on three molecular markers that have been de- logical identification generally requires time-consuming, monstrated to provide good resolution for species de- technologically sophisticated anatomical studies. Add- lineation in some Mollusca [23-25]. We included, in itionally, taxonomy frequently is complicated by mor- addition to mitochondrial Cytochrome c Oxidase subunit phological convergence and/or pronounced intraspecific I (COI) and 16S rRNA, nuclear 28S rRNA (large rib- variation (e.g., [3,5]). In Acochlidia, the most diverse osomal subunit – LSU), which has performed well for group of meiofaunal slugs, the Microhedylacea, shows species separation and was suggested as a signature ‘regressive evolution’ [6], exhibiting highly simplified gene fragment for a DNA taxonomy system for meio- organ systems and little morphological diversity even at benthos [8,26]. higher taxonomic levels [7]. Thus, it is challenging to Any method of species delineation is sensitive to use only morphology to delimit species boundaries in sampling [27], and rarity is almost universal when dealing meiofaunal slugs. In consequence of the fragmentary with invertebrates [28]. Rarity is not only a theoretical knowledge of meiofaunal taxonomy, this fauna is fre- problem in species delineation methods, but hinders as- quently neglected in conservation and biogeography, sessment of genetic variability [16,28]; if populations with and ecological analyses remain superficial despite the intermediate haplotype composition are left unsampled, undoubted importance of meiofauna; e.g., in the food barcoding and molecular-based species delineation ap- chain [8]. proaches tend to overestimate species [18,27]. With large For taxon-specific analyses, DNA-barcoding has been parts of the worldwide meiofauna still unexplored, and advocated as a fast and efficient way to reduce the taxo- patchy, discrete distributions being characteristic for nomic deficit and automate taxon determination for meiofaunal taxa [29], the present-day knowledge of this ecological research [3,5,8]. DNA-barcoding in its simple, fauna is prone to incomplete sampling. The rapidly similarity-based form of species identification [9] is not spreading biodiversity crisis with the destruction of habi- predictive; it fails if no identical sequence has yet been tats and high extinction rates calls for quick surveys and determined and deposited in a voucher database, or if realistic data for efficient conservation strategies (e.g., [16], no limit in species boundaries has been established and references therein). Currently, most molecular species [10,11]. In well-known taxa with good sampling cover- delineation approaches have been conducted on large age, identification rates via DNA-barcoding can be quite datasets with dense sampling coverage or on local scales successful (e.g.,[12,13]), but in case of meiofauna finding (e.g., [30-35]), with few exceptions using small datasets in identical sequences in public databases for a newly integrative approaches (e.g., [36]). Barcoding approaches collected mollusk or other under-investigated taxon is using COI trees for defining species clusters and revealing not expected to become the rule for decades to come. gaps between intra- and intercluster distance; multi-locus The application of the typical barcoding approach for tree-based methods with or without using diagnostic char- species delineation - COI in conjunction with a com- acters; and a couple of newly developed, tree-independent parison of pairwise distances - in Mollusca has resulted methods for species delineation, all serve as methods for in mixed reports: although the identification success DNA taxonomy. However, as a solution to address the with known taxa was generally high (e.g., [12,14,15]), the challenges of the taxonomic impediment in problematic determination of a ‘barcoding gap’ (i.e., significant differ- taxa the power of these methods is still largely untested. ence between inter- and intraspecific variation) and thus Here, we performed thorough phylogenetic analyses of a delimiting threshold has been problematic, especially all three molecular markers and integrate available add- above local-scale approaches and in undersampled phy- itional data from morphology and geography. In parallel, logenies [12]. Doubts have also arisen concerning spe- we applied four different methods of species delineation: cies identification and delimitation based on single-locus 1) The General Mixed Yule Coalescent model (GMYC) DNA sequences, which frequently result in problematic [32,33] is a maximum likelihood approach, able to dis- under- or overestimation of species [16-18]. Mitochon- criminate between population and speciation patterns drial markers, especially, came under criticism