Species delimitation and phylogeography of the Abies chensiensis complex inferred from morphological and molecular data

Species delimitation in the Abies chensiensis complex, consisting of A. chensiensis, A. recurvata, A. ernestii and A. ernestii var. salouenensis, has been the subject of a long‐term dispute, as various taxonomic combinations have been proposed. We combined different lines of evidence, including morphological characters and mitochondrial and plastid DNA sequence data, to assess species delimitation. Plastid DNA sequence variation was highly consistent with the morphological result and separated A. recurvata from the other three taxa, but the mitochondrial DNA genealogy was more strongly related to geographical distribution rather than species delimitation, as proposed in previous studies. On the basis of morphological characters and plastid DNA sequence variation, the current species status of A. recurvata was well supported, and we also accepted A. ernestii and A. ernestii var. salouenensis as varieties of A. chensiensis. Furthermore, the phylogeographical history of A. chensiensis was surveyed by molecular data and ecological niche modelling. On the basis of molecular data, the current phylogeographical pattern of the A. chensiensis complex was probably shaped by recent rapid expansions from multiple refugia in the Hengduan Mountains and habitat contraction after the last glaciation. This hypothesis was supported by the results of ecological niche modelling and previous fossil records. Our findings suggest that these fir taxa might have continued to expand their range after the largest glaciation until the end of the last glaciation, but have contracted since then with the increasing temperature and habitat changes. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 175–188.     

Integrative taxonomy of the primitively segmented spider genus Ganthela (Araneae: Mesothelae: Liphistiidae): DNA barcoding gap agrees with morphology

Species delimitation is difficult for taxa in which the morphological characters are poorly known because of the rarity of adult morphs or sexes, and in cryptic species. In primitively segmented spiders, family Liphistiidae, males are often unknown, and female genital morphology – usually species‐specific in spiders – exhibits considerable intraspecific variation. Here, we report on an integrative taxonomic study of the liphistiid genus Ganthela Xu & Kuntner, 2015, endemic to south‐east China, where males are only available for two of the seven morphological species (two known and five undescribed). We obtained DNA barcodes (cytochrome c oxidase subunit I gene, COI) for 51 newly collected specimens of six morphological species and analysed them using five species‐delimitation methods: DNA barcoding gap, species delimitation plugin [P ID(Liberal)], automatic barcode gap discovery (ABGD), generalized mixed Yule‐coalescent model (GMYC), and statistical parsimony (SP). Whereas the first three agreed with the morphology, GMYC and SP indicate several additional species. We used the consensus results to delimit and diagnose six Ganthela species, which in addition to the type species Ganthela yundingensis Xu, 2015, completes the revision of the genus. Although multi‐locus phylogenetic approaches may be needed for complex taxonomic delimitations, our results indicate that even single‐locus analyses based on the COI barcodes, if integrated with morphological and geographical data, may provide sufficiently reliable species delimitation. © 2015 The Linnean Society of London     

Discordance between morphological and molecular species boundaries among Caribbean species of the reef sponge Callyspongia

Sponges are among the most species‐rich and ecologically important taxa on coral reefs, yet documenting their diversity is difficult due to the simplicity and plasticity of their morphological characters. Genetic attempts to identify species are hampered by the slow rate of mitochondrial sequence evolution characteristic of sponges and some other basal metazoans. Here we determine species boundaries of the Caribbean coral reef sponge genus Callyspongia using a multilocus, model‐based approach. Based on sequence data from one mitochondrial (COI), one ribosomal (28S), and two single‐copy nuclear protein‐coding genes, we found evolutionarily distinct lineages were not concordant with current species designations in Callyspongia. While C. fallax, C. tenerrima, and C. plicifera were reciprocally monophyletic, four taxa with different morphologies (C. armigera, C. longissima, C. eschrichtii, and C. vaginalis) formed a monophyletic group and genetic distances among these taxa overlapped distances within them. A model‐based method of species delimitation supported collapsing these four into a single evolutionary lineage. Variation in spicule size among these four taxa was partitioned geographically, not by current species designations, indicating that in Callyspongia, these key taxonomic characters are poor indicators of genetic differentiation. Taken together, our results suggest a complex relationship between morphology and species boundaries in sponges.     

Range expansion during the Pleistocene drove morphological radiation of the fir genus (Abies,Pinaceae) in the Qinghai‐Tibet Plateau and Himalayas

Range expansion caused by climate oscillations in the past probably promoted morphological radiation in a few plant groups. In this study, we aim to test this hypothesis through phylogeographical analysis of the cold‐tolerant fir genus (Abies) in the Qinghai‐Tibet Plateau (QTP) and Himalayas, where it comprises 12 described species. We examined sequence variation in two maternally inherited mitochondrial (mt) DNA fragments (nad5‐4 and nad7‐1) and two paternally inherited plastid DNA fragments (trnS‐G and trnL‐F) for 733 individuals from 75 populations of the species in a monophyletic group. Only six mtDNA haplotypes were recovered, but five were shared between multiple species and one occurred at a high frequency, providing strong evidence of range expansion. Forty‐three plastid DNA haplotypes were detected, 19 of which were shared between species and three occurred at high frequency. Network, mismatch and Bayesian skyline plot analyses of all plastid DNA haplotypes from this clade clearly suggested range expansion. This expansion was dated as having occurred during the longest and most extensive glaciation in the Pleistocene. Our results therefore supported the range expansion hypothesis for this clade of Abies during the Pleistocene; expansion probably drove the morphological radiation of the clade in the QTP and Himalayas, although it remains unclear whether the different morphotypes should be acknowledged as independent, reproductively isolated species. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 444–453.     

Allozyme variation and population differentiation of the Aconitum delavayi complex (Ranunculaceae) in the Hengduan Mountains of China

The Aconitum delavayi complex is a group of four climbing species with trisect-leaves occurring in the Hengduan Mountains. The species of this complex are highly localized on very narrow regions with quite small population sizes. Because of rapid environmental changes recently in the Hengduan Mountains, this complex shows complicated morphological variability, which makes it difficult to delimit species. In the present study, 10 enzyme systems coding for 14 putative loci were employed to detect the interspecific and intraspecific genetic variation of the complex. In addition to low genetic diversity within all eight populations surveyed, the results indicate that A. episcopale is a distinct species because of high genetic identities among its three populations. Very low genetic divergence among populations of A. stapfianum and A. delavayi suggests that the two species should be treated as a single one.     

Integrative species delimitation reveals cryptic diversity in the southern Appalachian Antrodiaetus unicolor (Araneae: Antrodiaetidae) species complex

Although species delimitation can be highly contentious, the development of reliable methods to accurately ascertain species boundaries is an imperative step in cataloguing and describing Earth's quickly disappearing biodiversity. Spider species delimitation remains largely based on morphological characters; however, many mygalomorph spider populations are morphologically indistinguishable from each other yet have considerable molecular divergence. The focus of our study, the Antrodiaetus unicolor species complex containing two sympatric species, exhibits this pattern of relative morphological stasis with considerable genetic divergence across its distribution. A past study using two molecular markers, COI and 28S, revealed that A. unicolor is paraphyletic with respect to A. microunicolor. To better investigate species boundaries in the complex, we implement the cohesion species concept and use multiple lines of evidence for testing genetic exchangeability and ecological interchangeability. Our integrative approach includes extensively sampling homologous loci across the genome using a RADseq approach (3RAD), assessing population structure across their geographic range using multiple genetic clustering analyses that include structure , principal components analysis and a recently developed unsupervised machine learning approach (Variational Autoencoder). We evaluate ecological similarity by using large‐scale ecological data for niche‐based distribution modelling. Based on our analyses, we conclude that this complex has at least one additional species as well as confirm species delimitations based on previous less comprehensive approaches. Our study demonstrates the efficacy of genomic‐scale data for recognizing cryptic species, suggesting that species delimitation with one data type, whether one mitochondrial gene or morphology, may underestimate true species diversity in morphologically homogenous taxa with low vagility.     

DEFINING TAXON BOUNDARIES IN MEMBERS OF THE MORPHOLOGICALLY AND GENETICALLY PLASTIC GENES CAULERPA (CAULERPALES,CHLOROPHYTA)1

Numerous attempts to capture the morphological variability of the genus Caulerpa have resulted in an unstable classification of numerous varieties and formae. In the present study we attempted to test taxon boundaries by investigating morphological and genetic variation within and between seven taxa of Caulerpa, supposedly belonging to four species, sampled at different sites in a Philippine reef system. Using both field and culture observations, we described the relation between the variability of a set of morphological characters and ecological parameters, such as wave exposure, light intensity, and substrate type. Statistical analyses showed that the limits between two (out of three) ecads of the C. racemosa (Forsskål) J. Agardh complex were obscured by the presence of morphological plasticity. Other studied taxa of Caulerpa (i.e. C. cupressoides [Vahl] C. Agardh, C. serrulata [Forsskål] J. Agardh, and two formae of C. sertularioides [S. Gmelin] Howe) could be grouped based on morphology despite the presence of morphological plasticity. Our results indicated a strong association between light intensity and several quantitative morphological variables. Genetic diversity of these taxa was assessed by partial sequencing chloroplast rbcL and tufA genes and the ycf10‐chlB chloroplast spacer. In all phylogenetic analyses, C. serrulata, C. cupressoides, C. sertularioides, and the three ecads of C. racemosa emerged as distinct genetic units. Despite the presence of morphological plasticity and morphological convergence, a subset of morphological characters traditionally used in taxonomic delimitation still had sufficient discriminative power to recognize the terminal phylogenetic clades.     

Phylogeny,phylogeography and geographical variation in the Crocidura monax (Soricidae) species complex from the montane islands of Tanzania,with descriptions of three new species

We assess morphological and multilocus genetic variation among 11 isolated montane populations of white‐toothed shrews from Tanzania that have been referred to either Crocidura monax Thomas or C. montis Thomas. The montane sites we sampled represent ‘sky‐islands’ from two geologically distinct archipelagos (Northern Highlands and the Eastern Arc Mountains) and are a significant component of the Eastern Afromontane Biodiversity Hotspot. We used multivariate analyses of morphometric traits and phylogenetic and species‐delimitation analyses of multilocus DNA sequence data to assess species‐level diversity. Our species delimitation analyses included a novel, pairwise validation approach that avoids potential biases associated with specifying a guide tree. These analyses reveal several distinct lineages, which we treat as six allopatric species. Each species is restricted to one, two or four mountains. We use available names to recognize C. monax, C. tansaniana Hutterer and C. usambarae Dippenaar, while naming and describing three new species. Our results demonstrate the effectiveness of combining morphological and genetic data to uncover and describe hidden diversity in a cryptic mammalian system. © 2015 The Linnean Society of London     

Generic delimitation and biogeography of Maianthemum and Smilacina (Ruscaceae sensu lato): preliminary results based on partial 3′ matK gene and trnK 3′ intron sequences of cpDNA

The controversy over generic delimitation between Maianthemum and Smilacina has been unresolved for almost two centuries. Distributions of the two genera in the Northern Hemisphere also provide an excellent opportunity to further understand the disjunct distribution patterns of the North Temperate Flora. To test the generic delimitation and to investigate biogeographic patterns, we sequenced the partial 3′ matK gene and trnK 3′ intron of chloroplast DNA for 38 accessions, representing three species of Maianthemum, seven species of Smilacina, and four outgroup taxa. Maximum parsimony and neighbor-joining trees showed reciprocal monophyly of the two genera with very weak bootstrap support for each genus. Within each genus, relationships among species were poorly resolved. Despite its low resolution, this study shows that eastern Asian species of Smilacina and Maianthemum are generally more closely related to eastern North American taxa than to western ones. More detailed sampling of Smilacina from different geographic regions, especially from the two centers of diversity (southeastern Asia and Mexico/Central America), and additional sequences from cpDNA, as well as from nuclear DNA, are needed to test the reciprocal monophyly of the two genera and also to understand current distributions of disjunct taxa.     

New insights into DNA barcoding of seagrasses

Taxonomists find some plant genera challenging because of the few morphological differences or unclear characters among closely related species, which leads to the misidentification of taxa. DNA barcoding is an approach to identify species by using short orthologous DNA sequences, known as ‘DNA barcodes’. Concatenated rbcL and matK sequences are considered DNA barcodes for seagrasses. However, these markers are not applicable to all members of seagrasses at the species level, especially within the genus Halophila. Our previous studies indicated that the internal transcribed spacer (ITS) showed higher species resolution than the concatenated rbcL and matK sequences in the case of Halophila ovalis and closely related species. In this study, 26 ITS, two rbcL and two matK consensus sequences from 18 seagrass taxa belonging to four families collected in India, Vietnam, Germany, Croatia and Egypt were processed. Molecular ITS analysis resolved five clades. The results also indicate that the Cymodoceaceae family might be a non-monophyletic group. In conclusion, ITS could be applied as a DNA barcode for seagrasses instead of the rbcL/matK system previously proposed.     

Integrative taxonomy and preliminary assessment of species limits in the Liolaemus walkeri complex (Squamata,Liolaemidae) with descriptions of?three?new species from Peru

Species delimitation studies based on integrative taxonomic approaches have received considerable attention in the last few years, and have provided the strongest hypotheses of species boundaries. We used three lines of evidence (molecular, morphological, and niche envelopes) to test for species boundaries in Peruvian populations of the Liolaemus walkeri complex. Our results show that different lines of evidence and analyses are congruent in different combinations, for unambiguous delimitation of three lineages that were “hidden” within known species, and now deserve species status. Our phylogenetic analysis shows that L. walkeri, L. tacnae and the three new species are strongly separated from other species assigned to the alticolor-bibronii group. Few conventional morphological characters distinguish the new species from closely related taxa and this highlights the need to integrate other sources of data to erect strong hypothesis of species limits. A taxonomic key for known Peruvian species of the subgenus Lioalemus is provided.     

Multilocus species delimitation in Mesoamerican Scaptotrigona stingless bees (Apidae: Meliponini) supports the existence of cryptic species

Accelerating taxonomic knowledge and making accurate species identifications are critically important given the current biodiversity crisis, particularly in biodiversity hotspots such as Mesoamerica. Objective species delimitation that reduces investigator‐driven bias is fundamental to the establishment of appropriate conservation strategies, above all in managed species. Previous morphological and molecular studies on three managed stingless bee species of the genus Scaptotrigona distributed in Mexico (S. mexicana, S. pectoralis and S. hellwegeri) suggested that both S. mexicana and S. hellwegeri are cryptic species complexes. Herein we tested species delimitation by analysing sequence information of five markers (two mitochondrial: cox1 and 16S, and three nuclear: ITS1, EF1‐α, ArgK) within a Bayesian coalescent framework to test the putative species. We obtained two different hypotheses using a Generalized Mixed Yule Coalescent (GMYC) model: four (cox1) and six (16S) species. After the species validation step with the Bayesian species‐delimitation analysis (BPP), we suggest that only S. mexicana is a complex of two species with different distribution (along the Pacific and the Atlantic coasts, respectively). We highly recommend avoiding colony exchange between geographical regions in order to conserve the genetic integrity of both taxa.     

Exploring species boundaries with multiple genetic loci using empirical data from non‐biting midges

Over the past decade, molecular approaches to species delimitation have seen rapid development. However, species delimitation based on a single locus, for example, DNA barcodes, can lead to inaccurate results in cases of recent speciation and incomplete lineage sorting. Here, we compare the performance of Automatic Barcode Gap Discovery (ABGD), Bayesian Poisson tree processes (PTP), networks, generalized mixed Yule coalescent (GMYC) and Bayesian phylogenetics and phylogeography (BPP) models to delineate cryptic species previously detected by DNA barcodes within Tanytarsus (Diptera: Chironomidae) non‐biting midges. We compare the results from analyses of one mitochondrial (cytochrome c oxidase subunit I [COI]) and three nuclear (alanyl‐tRNA synthetase 1 [AATS1], carbamoyl phosphate synthetase 1 [CAD1] and 6‐phosphogluconate dehydrogenase [PGD]) protein‐coding genes. Our results show that species delimitation based on multiple nuclear DNA markers is largely concordant with morphological variation and delimitations using a single locus, for example, the COI barcode. However, ABGD, GMYC, PTP and network models led to conflicting results based on a single locus and delineate species differently than morphology. Results from BPP analyses on multiple loci correspond best with current morphological species concept. In total, 10 lineages of the Tanytarsus curticornis species complex were uncovered. Excluding a Norwegian population of Tanytarsus brundini which might have undergone recent hybridization, this suggests six hitherto unrecognized species new to science. Five distinct species are well supported in the Tanytarsus heusdensis species complex, including two species new to science.     

Phylogeny and species delimitation of North European Lumbricillus (Clitellata,Enchytraeidae)

The enchytraeid genus Lumbricillus comprises about 80 described species of clitellate worms, which are up to a few centimetres long, and they mostly inhabit the littoral zone of non‐tropical marine and brackish waters world‐wide. The phylogeny of this genus is poorly studied, but previous work has suggested that Lumbricillus is a non‐monophyletic group. In this study, species boundaries and the phylogeny of this genus is re‐assessed using more than 300 DNA‐barcoded specimens (using COI mtDNA), part of which was also sequenced for two additional mitochondrial and four nuclear molecular markers. Statistical and coalescent based applications were used for the delimitation of a total of 24 species, of which 20 were identified as belonging to 17 described morphospecies; one morphospecies was found to be a complex of four delimited species, and another four delimited species could not be matched with any described species. Furthermore, gene trees, concatenation and multispecies coalescent based species trees were estimated using Bayesian inference. The estimated phylogenies confirm a non‐monophyletic Lumbricillus as L. semifuscus is clearly excluded from the genus. Furthermore, the placement of a monophyletic clade consisting of L. arenarius, L. dubius, and an unidentified species varies between analyses; they are either found as the sister‐group to the genus Grania or as sister‐group to the remaining Lumbricillus, where the latter relationship is supported by the multispecies coalescent, which we consider as the most reliable method.     

A polyphasic approach to the delimitation of diatom species: a case study for the genus Pinnularia (Bacillariophyta)

Diatoms are one of the most abundant and arguably the most species‐rich group of protists. Diatom species delimitation has often been based exclusively on the recognition of morphological discontinuities without investigation of other lines of evidence. Even though DNA sequences and reproductive experiments have revealed several examples of (pseudo)cryptic diversity, our understanding of diatom species boundaries and diversity remains limited. The cosmopolitan pennate raphid diatom genus Pinnularia represents one of the most taxon‐rich diatom genera. In this study, we focused on the delimitation of species in one of the major clades of the genus, the Pinnularia subgibba group, based on 105 strains from a worldwide origin. We compared genetic distances between the sequences of seven molecular markers and selected the most variable pair, the mitochondrial cox1 and nuclear encoded LSU rDNA, to formulate a primary hypothesis on the species limits using three single‐locus automated species delimitation methods. We compared the DNA‐based primary hypotheses with morphology and with other available lines of evidence. The results indicate that our data set comprised 15 species of the P. subgibba group. The vast majority of these taxa have an uncertain taxonomic identity, suggesting that several may be unknown to science and/or members of (pseudo)cryptic species complexes within the P. subgibba group.     

The systematics of Lobophora (Dictyotales,Phaeophyceae) in the western Atlantic and eastern Pacific oceans: eight new species

Lobophora is a common tropical to temperate genus of brown algae found in a plethora of habitats including shallow and deep‐water coral reefs, rocky shores, mangroves, seagrass beds, and rhodoliths beds. Recent molecular studies have revealed that Lobophora species diversity has been severely underestimated. Current estimates of the species numbers range from 100 to 140 species with a suggested center of diversity in the Central Indo‐Pacific. This study used three molecular markers (cox3, rbcL, psbA), different single‐marker species delimitation methods (GMYC, ABGD, PTP), and morphological evidence to evaluate Lobophora species diversity in the Western Atlantic and the Eastern Pacific oceans. Cox3 provided the greatest number of primary species hypotheses(PSH), followed by rbcL and then psbA. GMYC species delimitation analysis was the most conservative across all three markers, followed by PTP, and then ABGD. The most informative diagnostic morphological characters were thallus thickness and number of cell layers in both the medulla and the dorsal/ventral cortices. Following a consensus approach, 14 distinct Lobophora species were identified in the Western Atlantic and five in the Eastern Pacific. Eight new species from these two oceans were herein described: L. adpressa sp. nov., L. cocoensis sp. nov., L. colombiana sp. nov., L. crispata sp. nov., L. delicata sp. nov., L. dispersa sp. nov., L. panamensis sp. nov., and L. tortugensis sp. nov. This study showed that the best approach to confidently identify Lobophora species is to analyze DNA sequences (preferably cox3 and rbcL) followed by comparative morphological and geographical assessment.     

Phylogeography of the Oenanthe hispanica–pleschanka–cypriaca complex (Aves,Muscicapidae: Saxicolinae): Diversification history of open‐habitat specialists based on climate niche models,genetic data,and morphometric data

The succession of glacials and interglacials during the Pleistocene strongly influenced the diversification and distribution patterns in birds. In contrast to species of temperate regions, open‐habitat specialists should have experienced range expansion during the longer glacial periods, while range contractions occurred during the shorter interglacials. However, only few studies have tested this prediction so far. We studied the Oenanthe hispanica–pleschanka–cypriaca (Aves, Muscicapidae: Saxicolinae) complex characteristic of open habitats in the Palearctic. Based on three mitochondrial and one Z‐linked nuclear marker, we inferred its phylogeny, historical diversification, and demography. Ecological niche modeling was used to reconstruct potential distributions during the last glacial maximum and the last interglacial. Using 19 morphological traits, we tested for morphometric differences among the different taxa. Mitochondrial markers revealed strong genetic differences between O. h. hispanica and the other taxa with a divergence event at around 1.7 million years ago. No consistent genetic differences were revealed between O. cypriaca, O. h. melanoleuca, and O. pleschanka. The latter two hybridize in contact zones, which might explain partly the lack of genetic differentiation; yet, further analyses using genomic data are needed to infer the true divergence history of the complex. Signs of population expansions in the clade comprising O. h. melanoleuca, O. pleschanka, and O. cypriaca at 90,000 years ago coincided with the last glacial as predicted. Population expansion then was also supported by ecological climate niche models. O. h. hispanica was not consistently separated from the other taxa in morphometrics. It might nonetheless warrant species status, pending further analyses.