Phylogeny and species reassessment of Hyalopterus (Aphididae,Aphidinae)

The broadly distributed genus Hyalopterus currently comprises three formally recognized species that are highly similar morphologically and hence difficult to be identified with certainty. This group has undergone multiple revisions in the past century, but none of these has assessed species from Asia, which has hampered our understanding of the species diversity within this genus. Based on a comprehensive data set from morphological data and host-associated data, and by coalescent-based delimitation approaches, the Hyalopterus species boundaries, distribution and diversity were clarified here to further reveal the composition of the species. Two single-locus (ML-GMYC and mPTP) and two multilocus (BPP and STACEY) delimitation methods were conducted based on extensive sampling. Then, the phylogenetic relationships and morphological divergence were assessed. Our data strongly supported that the number of recognized species in Hyalopterus had likely been underestimated. The phylogenetic analyses recovered four major clades, which corresponded to distinct host-plant preferences. Also, the morphological analyses showed significant differentiation for only one of the newly recognized candidate species uncovered by the delimitation approaches, suggesting the existence of at least two independent evolutionary lineages within Hyalopterus arundiniformis, which showed different patterns of host association. Moreover, based on our data, the taxonomic misidentification of H. arundiniformis in China was corrected here. This study lays the groundwork for the thorough taxonomic revision of Hyalopterus and for future evolutionary studies and underlines the importance of an integrated framework for species determination.     

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.     

Genomic data reveal deep genetic structure but no support for current taxonomic designation in a grasshopper species complex

Taxonomy has traditionally relied on morphological and ecological traits to interpret and classify biological diversity. Over the last decade, technological advances and conceptual developments in the field of molecular ecology and systematics have eased the generation of genomic data and changed the paradigm of biodiversity analysis. Here we illustrate how traditional taxonomy has led to species designations that are supported neither by high throughput sequencing data nor by the quantitative integration of genomic information with other sources of evidence. Specifically, we focus on Omocestus antigai and Omocestus navasi, two montane grasshoppers from the Pyrenean region that were originally described based on quantitative phenotypic differences and distinct habitat associations (alpine vs. Mediterranean‐montane habitats). To validate current taxonomic designations, test species boundaries, and understand the factors that have contributed to genetic divergence, we obtained phenotypic (geometric morphometrics) and genome‐wide SNP data (ddRADSeq) from populations covering the entire known distribution of the two taxa. Coalescent‐based phylogenetic reconstructions, integrative Bayesian model‐based species delimitation, and landscape genetic analyses revealed that populations assigned to the two taxa show a spatial distribution of genetic variation that do not match with current taxonomic designations and is incompatible with ecological/environmental speciation. Our results support little phenotypic variation among populations and a marked genetic structure that is mostly explained by geographic distances and limited population connectivity across the abrupt landscapes characterizing the study region. Overall, this study highlights the importance of integrative approaches to identify taxonomic units and elucidate the evolutionary history of species.     

Integration of molecular, ecological, morphological and endosymbiont data for species delimitation within the Pnigalio soemius complex (Hymenoptera: Eulophidae)

Integrative taxonomy is a recently developed approach that uses multiple lines of evidence such as molecular, morphological, ecological and geographical data to test species limits, and it stands as one of the most promising approaches to species delimitation in taxonomically difficult groups. The Pnigalio soemius complex (Hymenoptera: Eulophidae) represents an interesting taxonomical and ecological study case, as it is characterized by a lack of informative morphological characters, deep mitochondrial divergence, and is susceptible to infection by parthenogenesis‐inducing Rickettsia. We tested the effectiveness of an integrative taxonomy approach in delimiting species within the P. soemius complex. We analysed two molecular markers (COI and ITS2) using different methods, performed multivariate analysis on morphometric data and exploited ecological data such as host–plant system associations, geographical separation, and the prevalence, type and effects of endosymbiont infection. The challenge of resolving different levels of resolution in the data was met by setting up a formal procedure of data integration within and between conflicting independent lines of evidence. An iterative corroboration process of multiple sources of data eventually indicated the existence of several cryptic species that can be treated as stable taxonomic hypotheses. Furthermore, the integrative approach confirmed a trend towards host specificity within the presumed polyphagous P. soemius and suggested that Rickettsia could have played a major role in the reproductive isolation and genetic diversification of at least two species.     

Multivariate species boundaries and conservation of harlequin poison frogs

In this study, we present an iterative method for delimiting species under the general lineage concept (GLC) based on the multivariate clustering of morphological, ecological and genetic data. Our rationale is that distinct multivariate groups correspond to evolutionarily independent metapopulation lineages because they reflect the common signal of different secondary defining properties (environmental and genetic distinctiveness, phenotypic diagnosability, etc.) that imply the existence of barriers preventing or limiting gene exchange. We applied this method to study a group of endangered poison frogs, the Oophaga histrionica complex. In our study case, we used next‐generation targeted amplicon sequencing to obtain a robust genetic data set that we combined with patterns of morphological and ecological features. Our analyses revealed the existence of at least five different species in the histrionica complex (three, new to science), some of them, occurring in small isolated populations outside any protected areas. The lineage delimitation proposed here has important conservation implications as it revealed that some of the Oophaga species should be considered among the most vulnerable of the Neotropical frogs. More broadly, our study exemplifies how multiple‐amplicon and multivariate statistical techniques can be integrated to successfully identify species and their boundaries.     

Cryptic species in the Neotropical fish genus Curimatopsis (Teleostei,Characiformes)

Detritivores of the fish family Curimatidae are assigned to eight genera, one of which, the Curimatopsis, with only five species, is the least speciose genus and sister to other seven genera in the family. Ongoing morphological investigations reveal, however, the likely existence of additional species. In this study, fifty‐one specimens of Curimatopsis from multiple rivers of the Amazon, Paraguay and Suriname drainages were identified morphologically according to the present species concepts and then barcoded using the universal cytochrome c oxidase subunit I (COI) mitochondrial marker. Species delimitation analyses were conducted using Bayesian methods through the general mixed Yule‐coalescent analysis combined with conventional likelihood, genetic distance and haplotypic diversity approaches. We found eleven well‐supported clusters that represent four of the named species and seven cryptic, undescribed species of Curimatopsis. Our results show a clear delimitation of species boundaries constrained by distinct Amazonian river ecotones that may have promoted intrageneric lineage diversification. This is the first of a series of genetic studies applicable to future taxonomic, phylogenetic and evolutionary studies across the Curimatidae.     

Integrative insect taxonomy based on morphology,mitochondrial DNA,and hyperspectral reflectance profiling

Integrative taxonomy is considered a reliable taxonomic approach of closely related and cryptic species by integrating different sources of taxonomic data (genetic, ecological, and morphological characters). In order to infer the boundaries of seven species of the evacanthine leafhopper genus Bundera Distant, 1908 (Hemiptera: Cicadellidae), an integrated analysis based on morphology, mitochondrial DNA, and hyperspectral reflectance profiling (37 spectral bands from 411–870 nm) was conducted. Despite their morphological similarities, the genetic distances of the cytochrome c oxidase subunit I (COI) gene among the tested species are relatively large (5.8–17.3%). The species‐specific divergence of five morphologically similar species (Bundera pellucida and Bundera spp. 1–4) was revealed in mitochondrial DNA data and reflectance profiling. A key to identifying males is provided, and their morphological characters are described. Average reflectance profiles from the dorsal side of specimens were classified based on linear discriminant analysis. Cross‐validation of reflectance‐based classification revealed that the seven species could be distinguished with 91.3% classification accuracy. This study verified the feasibility of using hyperspectral imaging data in insect classification, and our work provides a good example of using integrative taxonomy in studies of closely related and cryptic species. © 2015 The Linnean Society of London     

Species detection and individual assignment in species delimitation: can integrative data increase efficacy?

Statistical species delimitation usually relies on singular data, primarily genetic, for detecting putative species and individual assignment to putative species. Given the variety of speciation mechanisms, singular data may not adequately represent the genetic, morphological and ecological diversity relevant to species delimitation. We describe a methodological framework combining multivariate and clustering techniques that uses genetic, morphological and ecological data to detect and assign individuals to putative species. Our approach recovers a similar number of species recognized using traditional, qualitative taxonomic approaches that are not detected when using purely genetic methods. Furthermore, our approach detects groupings that traditional, qualitative taxonomic approaches do not. This empirical test suggests that our approach to detecting and assigning individuals to putative species could be useful in species delimitation despite varying levels of differentiation across genetic, phenotypic and ecological axes. This work highlights a critical, and often overlooked, aspect of the process of statistical species delimitation—species detection and individual assignment. Irrespective of the species delimitation approach used, all downstream processing relies on how individuals are initially assigned, and the practices and statistical issues surrounding individual assignment warrant careful consideration.     

Geographic variation in hybridization and ecological differentiation between three syntopic,morphologically similar species of montane lizards

To understand factors shaping species boundaries in closely related taxa, a powerful approach is to compare levels of genetic admixture at multiple points of contact and determine how this relates to intrinsic and extrinsic factors, such as genetic, morphological and ecological differentiation. In the Australian Alps, the threatened alpine bog skink Pseudemoia cryodroma co‐occurs with two morphologically and ecologically similar congeners, P. entrecasteauxii and P. pagenstecheri, and all three species are suspected to hybridize. We predicted that the frequency of hybridization should be negatively correlated with genetic divergence, morphological differentiation and microhabitat separation. We tested this hypothesis using a mitochondrial locus, 13 microsatellite loci, morphological and microhabitat data and compared results across three geographically isolated sites. Despite strong genetic structure between species, we detected hybridization between all species pairs, including evidence of backcrossed individuals at the two sites where all three species are syntopic. Hybridization frequencies were not consistently associated with genetic, morphological or ecological differentiation. Furthermore, P. entrecasteauxii and P. pagenstecheri only hybridized at the two sites where they are syntopic with P. cryodroma, but not at the largest site where P. cryodroma was not recorded, suggesting that P. cryodroma may serve as a bridging species. This study reveals the complex dynamics within a three species hybrid zone and provides a baseline for assessing the impact of climate change and anthropogenic habitat modification on future hybridization frequencies.     

Species boundaries in the messy middle—A genome‐scale validation of species delimitation in a recently diverged lineage of coastal fog desert lichen fungi

Species delimitation among closely related species is challenging because traditional phenotype‐based approaches, for example, using morphology, ecological, or chemical characteristics, may not coincide with natural groupings. With the advent of high‐throughput sequencing, it has become increasingly cost‐effective to acquire genome‐scale data which can resolve previously ambiguous species boundaries. As the availability of genome‐scale data has increased, numerous species delimitation analyses, such as BPP and SNAPP+Bayes factor delimitation (BFD*), have been developed to delimit species boundaries. However, even empirical molecular species delimitation approaches can be biased by confounding evolutionary factors, for example, hybridization/introgression and incomplete lineage sorting, and computational limitations. Here, we investigate species boundaries and the potential for micro‐endemism in a lineage of lichen‐forming fungi, Niebla Rundel & Bowler, in the family Ramalinaceae by analyzing single‐locus and genome‐scale data consisting of (a) single‐locus species delimitation analysis using ASAP, (b) maximum likelihood‐based phylogenetic tree inference, (c) genome‐scale species delimitation models, e.g., BPP and SNAPP+BFD, and (d) species validation using the genealogical divergence index (gdi). We specifically use these methods to cross‐validate results between genome‐scale and single‐locus datasets, differently sampled subsets of genomic data and to control for population‐level genetic divergence. Our species delimitation models tend to support more speciose groupings that were inconsistent with traditional taxonomy, supporting a hypothesis of micro‐endemism, which may include morphologically cryptic species. However, the models did not converge on robust, consistent species delimitations. While the results of our analysis are somewhat ambiguous in terms of species boundaries, they provide a valuable perspective on how to use these empirical species delimitation methods in a nonmodel system. This study thus highlights the challenges inherent in delimiting species, particularly in groups such as Niebla, with complex, relatively recent phylogeographic histories.     

Integrated species delimitation and conservation implications of an endangered weevil Pachyrhynchus sonani (Coleoptera: Curculionidae) in Green and Orchid Islands of Taiwan

Oceanic islands are productive habitats for generating new species and high endemism, which is primarily due to their geographical isolation, smaller population sizes and local adaptation. However, the short divergence times and subtle morphological or ecological divergence of insular organisms may obscure species identity, so the cryptic endemism on islands may be underestimated. The endangered weevil Pachyrhynchus sonani Kôno (Coleoptera: Curculionidae: Entiminae: Pachyrhynchini) is endemic to Green Island and Orchid Island of the Taiwan‐Luzon Archipelago and displays widespread variation in coloration and host range, thus raising questions regarding its species boundaries and degree of cryptic diversity. We tested the species boundaries of P. sonani using an integrated approach that combined morphological (body size and shape, genital shape, coloration and cuticular scale), genetic (four genes and restriction site‐associated DNA sequencing, RAD‐seq) and ecological (host range and distribution) diversity. The results indicated that all the morphological datasets for male P. sonani, except for the colour spectrum, reveal overlapping but statistically significant differences between islands. In contrast, the morphology of the female P. sonani showed minimum divergence between island populations. The populations of P. sonani on the two islands were significantly different in their host ranges, and the genetic clustering and phylogenies of P. sonani established two valid evolutionary species. Integrated species delimitation combining morphological, molecular and ecological characters supported two distinct species of P. sonani from Green Island and Orchid Island. The Green Island population was described as P. jitanasaius sp.n. Chen & Lin, and it is recommended that its threatened conservation status be recognized. Our findings suggest that the inter‐island speciation of endemic organisms inhabiting both islands may be more common than previously thought, and they highlight the possibility that the cryptic diversity of small oceanic islands may still be largely underestimated.     

Assessing the diversity and distribution of Cephalothrix species (Nemertea: Palaeonemertea) in European waters by comparing different species delimitation methods

Soft‐bodied marine taxa, like ribbon worms (Nemertea), often lack clear diagnostic morphological characters impeding traditional species delimitation. Therefore, recent studies concentrated on molecular genetic methods to solve taxonomic issues. Different delimitation methods were employed to explore species boundaries and the presence of cryptic species. However, the performance of the different delimitation methods needs to be tested. A particularly promising nemertean genus in this regard is the palaeonemertean genus Cephalothrix that is commonly found in European waters. In order to gain information on the number and distribution of European cephalotrichids and to test different tree‐based and non‐tree‐based delimitation methods, we analyzed a dataset comprising the barcoding region of the mitochondrial cytochrome c oxidase subunit I (COI) of 215 European Cephalothrix specimens, of which 78 were collected for this study. Our results show the presence of 12–13 European lineages of which several can be assigned to known European species. Analyzing a second dataset comprising 74 additional sequences from the Pacific and the Atlantic Oceans helped identify some of the unassigned European specimens. One resulting clade seems to represent a non‐native introduced Cephalothrix species, while another has never been recorded from Europe before. In our analysis, especially the tree‐based methods and the phylogenetic analysis proved to be a useful tool when delimiting species. It remains unclear whether the different identified clades result from cryptic speciation or from a high genetic variability of the COI gene.     

A genomic evaluation of taxonomic trends through time in coast horned lizards (genus Phrynosoma)

Determining the boundaries between species and deciding when to describe new species are challenging practices that are particularly difficult in groups with high levels of geographic variation. The coast horned lizards (Phrynosoma blainvillii, Phrynosoma cerroense and P. coronatum) have an extensive geographic distribution spanning many distinctive ecological regions ranging from northern California to the Cape Region of Baja California, Mexico, and populations differ substantially with respect to external morphology across much of this range. The number of taxa recognized in the group has been reevaluated by herpetologists over 20 times during the last 180 years, and typically without the aid of explicit species delimitation methods, resulting in a turbulent taxonomy containing anywhere from one to seven taxa. In this study, we evaluate taxonomic trends through time by ranking 15 of these species delimitation models (SDMs) using coalescent analyses of nuclear loci and SNPs in a Bayesian model comparison framework. Species delimitation models containing more species were generally favoured by Bayesian model selection; however, several three‐species models outperformed some four‐ and five‐species SDMs, and the top‐ranked model, which contained five species, outperformed all SDMs containing six species. Model performance peaked in the 1950s based on marginal likelihoods estimated from nuclear loci and SNPs. Not surprisingly, SDMs based on genetic data outperformed morphological taxonomies when using genetic data alone to evaluate models. The de novo estimation of population structure favours a three‐population model that matches the currently recognized integrative taxonomy containing three species. We discuss why Bayesian model selection might favour models containing more species, and why recognizing more than three species might be warranted.     

DNA‐based species delimitation separates highly divergent populations within morphologically coherent clades of poorly dispersing beetles

We compared species‐level entities recovered using distance, tree‐based, and DNA‐character based methods with morphologically defined species in poorly dispersing lyponiine net‐winged beetles. The phylogenetic relationships were investigated using the cytochrome c oxidase 1–2 mtDNA fragment. We identified 31 species using a morphology‐based concept and additional candidate species were delimited by the general mixed Yule‐coalescent method and barcoding threshold within the morphologically coherent lineages owing to high genetic divergence (up to 10.97% within morphologically defined species, highest density 2.5–6.5%). Genetic divergence was positively correlated with geographical distance: lower in continental China where no apparent dispersal barriers are present (r² = 0.172, P ≤ 0.001, average increase of genetic distance of 0.32% per 100 km) and much higher in Japan (r² = 0.490, P ≤ 0.001, 1.81% per 100 km). We hypothesize that low dispersal propensity contributed to the high level of intraspecific, geographically structured divergence. DNA‐based methods suggested a high number of morphologically undistinguishable species. The observed patterns agree with the model of neutral evolution. The poor dispersers produce gradually divergent populations across the range of the morphologically coherent lineages. We pose the question of what size of range and level of genetic difference justify formal acceptance of a species without morphological divergence from both the taxonomical and conservation management view points. © 2015 The Linnean Society of London     

Delimiting cryptic pathogen species causing apple Valsa canker with multilocus data

Fungal diseases are posing tremendous threats to global economy and food safety. Among them, Valsa canker, caused by fungi of Valsa and their Cytospora anamorphs, has been a serious threat to fruit and forest trees and is one of the most destructive diseases of apple in East Asia, particularly. Accurate and robust delimitation of pathogen species is not only essential for the development of effective disease control programs, but also will advance our understanding of the emergence of plant diseases. However, species delimitation is especially difficult in Valsa because of the high variability of morphological traits and in many cases the lack of the teleomorph. In this study, we delimitated species boundary for pathogens causing apple Valsa canker with a multifaceted approach. Based on three independent loci, the internal transcribed spacer (ITS), β‐tubulin (Btu), and translation elongation factor‐1 alpha (EF1α), we inferred gene trees with both maximum likelihood and Bayesian methods, estimated species tree with Bayesian multispecies coalescent approaches, and validated species tree with Bayesian species delimitation. Through divergence time estimation and ancestral host reconstruction, we tested the possible underlying mechanisms for fungal speciation and host‐range change. Our results proved that two varieties of the former morphological species V. mali represented two distinct species, V. mali and V. pyri, which diverged about 5 million years ago, much later than the divergence of their preferred hosts, excluding a scenario of fungi–host co‐speciation. The marked different thermal preferences and contrasting pathogenicity in cross‐inoculation suggest ecological divergences between the two species. Apple was the most likely ancestral host for both V. mali and V. pyri. Host‐range expansion led to the occurrence of V. pyri on both pear and apple. Our results also represent an example in which ITS data might underestimate species diversity.     

Integrative analyses of Nervilia (Orchidaceae) section Linervia reveal further undescribed cryptic diversity in Thailand

The delimitation of cryptic species is necessary to accurately classify and appropriately conserve biodiversity. Integrative analyses can be incisive in detecting and circumscribing cryptic diversity, especially in species complexes whose members are delineated by minor or overlapping morphological variation. We adopt an integrative approach to assess species relationships and resolve species boundaries in the taxonomically difficult Nervilia adolphi/punctata species alliance of N. sect. Linervia, an Old World complex of reduced, one-flowered terrestrial orchids that is both species-rich and poorly known in tropical and warm temperate Asia. We sampled 12 of the 27 known species of the alliance in Asia, including all four species reported from Thailand and a further 20 plants collected in that country that could not be satisfactorily identified using morphology alone. Phylogenetic analyses using one nuclear (ITS) and two plastid (matK and trnL-F) markers confirmed both N. sect. Linervia and the alliance itself as monophyletic, and corroborated 11 of the 12 sampled species; N. punctata proved polyphyletic, with the Thai samples referred to this Indonesian species falling sister to the Himalayan N. mackinnonii. The 20 unidentified Thai samples formed three distinct, strongly supported clades. STACEY, a Bayesian coalescence approach to species delimitation, resolved the same three clusters, but provided evidence suggesting that one comprised two distinct sub-clades. Building on this genetic evidence, we identify subtle morphological differences and invoke a diagnosable species concept to circumscribe three previously unrecognized cryptic species from Thailand. This objective approach to species delimitation validates ostensibly minor morphological differences as a basis for differentiating species within the alliance, paving the way for a global analysis of species boundaries throughout the genus as a whole.     

Gene flow and species delimitation in fishes of Western North America: Flannelmouth (Catostomus latipinnis) and Bluehead sucker (C. Pantosteus discobolus)

The delimitation of species boundaries, particularly those obscured by reticulation, is a critical step in contemporary biodiversity assessment. It is especially relevant for conservation and management of indigenous fishes in western North America, represented herein by two species with dissimilar life histories codistributed in the highly modified Colorado River (i.e., flannelmouth sucker, Catostomus latipinnis; bluehead sucker, C. (Pantosteus) discobolus). To quantify phylogenomic patterns and examine proposed taxonomic revisions, we first employed double‐digest restriction site‐associated DNA sequencing (ddRAD), yielding 39,755 unlinked SNPs across 139 samples. These were subsequently evaluated with multiple analytical approaches and by contrasting life history data. Three phylogenetic methods and a Bayesian assignment test highlighted similar phylogenomic patterns in each, but with considerable difference in presumed times of divergence. Three lineages were detected in bluehead sucker, supporting elevation of C. (P.) virescens to species status and recognizing C. (P.) discobolus yarrowi (Zuni bluehead sucker) as a discrete entity. Admixture in the latter necessitated a reevaluation of its contemporary and historic distributions, underscoring how biodiversity identification can be confounded by complex evolutionary histories. In addition, we defined three separate flannelmouth sucker lineages as ESUs (evolutionarily significant units), given limited phenotypic and genetic differentiation, contemporary isolation, and lack of concordance (per the genealogical concordance component of the phylogenetic species concept). Introgression was diagnosed in both species, with the Little Colorado and Virgin rivers in particular. Our diagnostic methods, and the agreement of our SNPs with previous morphological, enzymatic, and mitochondrial work, allowed us to partition complex evolutionary histories into requisite components, such as isolation versus secondary contact.