Plant species delimitation: A comparison of morphological and molecular markers

Abstract

Species delimitation is fundamental in many areas of biology. Despite its importance, there is no agreement on criteria for species delimitation mostly due to divergence on the point of view adopted by the different biological disciplines. Two main groups of diagnostic characters are commonly used to distinguish species: the traditional morphological ones and the molecular ones. Field species recognition and sampling are generally based on morphological characters, but they can either fail to discriminate species and mask the presence of cryptic species or discriminate different species while in reality there is only one. To overcome this problem it is common to compare clusters obtained on the basis of the observed polymorphism of both characters, and to analyse their agreement. Here we compile a set of studies that have examined species delimitation with both markers. This provides a review of the different morphological and molecular markers, and of the sampling strategy and clustering methodology generally employed to delimitate species. Some conclusions are drawn with regard to species delimitation, when comparing diagnostic morphological and molecular markers.     

Failure to Replicate a Genetic Association May Provide Important Clues About Genetic Architecture

Replication has become the gold standard for assessing statistical results from genome-wide association studies. Unfortunately this replication requirement may cause real genetic effects to be missed. A real result can fail to replicate for numerous reasons including inadequate sample size or variability in phenotype definitions across independent samples. In genome-wide association studies the allele frequencies of polymorphisms may differ due to sampling error or population differences. We hypothesize that some statistically significant independent genetic effects may fail to replicate in an independent dataset when allele frequencies differ and the functional polymorphism interacts with one or more other functional polymorphisms. To test this hypothesis, we designed a simulation study in which case-control status was determined by two interacting polymorphisms with heritabilities ranging from 0.025 to 0.4 with replication sample sizes ranging from 400 to 1600 individuals. We show that the power to replicate the statistically significant independent main effect of one polymorphism can drop dramatically with a change of allele frequency of less than 0.1 at a second interacting polymorphism. We also show that differences in allele frequency can result in a reversal of allelic effects where a protective allele becomes a risk factor in replication studies. These results suggest that failure to replicate an independent genetic effect may provide important clues about the complexity of the underlying genetic architecture. We recommend that polymorphisms that fail to replicate be checked for interactions with other polymorphisms, particularly when samples are collected from groups with distinct ethnic backgrounds or different geographic regions.     

Phylogenetic analysis and intraspecific variation: performance of parsimony,likelihood, and distance methods

Intraspecific variation is abundant in all types of systematic characters but is rarely addressed in simulation studies of phylogenetic method performance. We compared the accuracy of 15 phylogenetic methods using simulations to (1) determine the most accurate method(s) for analyzing polymorphic data (under simplified conditions) and (2) test if generalizations about the performance of phylogenetic methods based on previous simulations of fixed (nonpolymorphic) characters are robust to a very different evolutionary model that explicitly includes intraspecific variation. Simulated data sets consisted of allele frequencies that evolved by genetic drift. The phylogenetic methods included eight parsimony coding methods, continuous maximum likelihood, and three distance methods (UPGMA, neighbor joining, and Fitch-Margoliash) applied to two genetic distance measures (Nei's and the modified Cavalli-Sforza and Edwards chord distance). Two sets of simulations were performed. The first examined the effects of different branch lengths, sample sizes (individuals sampled per species), numbers of characters, and numbers of alleles per locus in the eight-taxon case. The second examined more extensively the effects of branch length in the four-taxon, two-allele case. Overall, the most accurate methods were likelihood, the additive distance methods (neighbor joining and Fitch-Margoliash), and the frequency parsimony method. Despite the use of a very different evolutionary model in the present article, many of the results are similar to those from simulations of fixed characters. Similarities include the presence of the "Felsenstein zone," where methods often fail, which suggests that long-branch attraction may occur among closely related species through genetic drift. Differences between the results of fixed and polymorphic data simulations include the following: (1) UPGMA is as accurate or more accurate than nonfrequency parsimony methods across nearly all combinations of branch lengths, and (2) likelihood and the additive distance methods are not positively misled under any combination of branch lengths tested (even when the assumptions of the methods are violated and few characters are sampled). We found that sample size is an important determinant of accuracy and affects the relative success of methods (i.e., distance and likelihood methods outperform parsimony at small sample sizes). Attempts to generalize about the behavior of phylogenetic methods should consider the extreme examples offered by fixed-mutation models of DNA sequence data and genetic-drift models of allele frequencies.     

The efficacy of short tandem repeat polymorphisms versus single-nucleotide polymorphisms for resolving population structure

Accurately resolving population structure in a sample is important for both linkage and association studies. In this study we investigated the power of single-nucleotide polymorphisms (SNPs) in detecting population structure in a sample of 286 unrelated individuals. We varied the number of SNPs to determine how many are required to approach the degree of resolution obtained with the Collaborative Study on the Genetics of Alcoholism (COGA) short tandem repeat polymorphisms (STRPs). In addition, we selected SNPs with varying minor allele frequencies (MAFs) to determine whether low or high frequency SNPs are more efficient in resolving population structure. We conclude that a set of at least 100 evenly spaced SNPs with MAFs of 40-50% is required to resolve population structure in this dataset. If SNPs with lower MAFs are used, then more than 250 SNPs may be required to obtain reliable results.     

CLADES: A classification‐based machine learning method for species delimitation from population genetic data

Species are considered to be the basic unit of ecological and evolutionary studies. As multilocus genomic data are increasingly available, there have been considerable interests in the use of DNA sequence data to delimit species. In this study, we show that machine learning can be used for species delimitation. Our method treats the species delimitation problem as a classification problem for identifying the category of a new observation on the basis of training data. Extensive simulation is first conducted over a broad range of evolutionary parameters for training purposes. Each pair of known populations is combined to form training samples with a label of “same species” or “different species”. We use support vector machine (SVM) to train a classifier using a set of summary statistics computed from training samples as features. The trained classifier can classify a test sample to two outcomes: “same species” or “different species”. Given multilocus genomic data of multiple related organisms or populations, our method (called CLADES) performs species delimitation by first classifying pairs of populations. CLADES then delimits species by maximizing the likelihood of species assignment for multiple populations. CLADES is evaluated through extensive simulation and also tested on real genetic data. We show that CLADES is both accurate and efficient for species delimitation when compared with existing methods. CLADES can be useful especially when existing methods have difficulty in delimitation, for example with short species divergence time and gene flow.     

DNA and Morphology Unite Two Species and 10 Million Year Old Fossils

Species definition and delimitation is a non-trivial problem in evolutionary biology that is particularly problematic for fossil organisms. This is especially true when considering the continuity of past and present species, because species defined in the fossil record are not necessarily equivalent to species defined in the living fauna. Correctly assigned fossil species are critical for sensitive downstream analysis (e.g., diversification studies and molecular-clock calibration). The marine snail genus Alcithoe exemplifies many of the problems with species identification. The paucity of objective diagnostic characters, prevalence of morphological convergence between species and considerable variability within species that are observed in Alcithoe are typical of a broad range of fossilised organisms. Using a synthesis of molecular and morphometric approaches we show that two taxa currently recognised as distinct are morphological variants of a single species. Furthermore, we validate the fossil record for one of these morphotypes by finding a concordance between the palaeontological record and divergence time of the lineage inferred using molecular-clock analysis. This work demonstrates the utility of living species represented in the fossil record as candidates for molecular-clock calibration, as the veracity of fossil species assignment can be more rigorously tested.     

CoMuS: simulating coalescent histories and polymorphic data from multiple species

The simultaneous analysis of intra‐ and interspecies variation is challenging mainly because our knowledge about patterns of polymorphisms where both intra‐ and interspecies samples coexist is limited. In this study, we present CoMuS (Coalescent of Multiple Species), a multispecies coalescent software that can simulate intra‐ and interspecies polymorphisms. CoMuS supports a variety of speciation models and demographic scenarios related to the history of each species. In CoMuS, speciation can be accompanied by either instant or gradual isolation between sister species. Sampling may also occur in the past, and thus, we can study simultaneously extinct and extant species. Our software supports both the infinite‐ and the finite‐site model, with substitution rate heterogeneity among sites and a user‐defined proportion of invariable sites. We demonstrate the usage of CoMuS in various applications: species delimitation, software testing, model selection and parameter inference involving present‐day and ancestral samples, comparison between gradual and instantaneous isolation models, estimation of speciation time between human and chimpanzee using both intra‐ and interspecies variation. We expect that CoMuS will be particularly useful for studies where species have been separated recently from their common ancestor and phenomena such as incomplete lineage sorting or introgression still occur.     

Costs of carnivory: tooth fracture in Pleistocene and Recent carnivorans

Large, carnivorous mammals often break their teeth, probably as a result of tooth to bone contact that occurs when carcasses are consumed more fully, a behaviour likely to occur under conditions of food stress. Recent studies of Pleistocene predators revealed high numbers of teeth broken in life, suggesting that carcass utilization and, consequently, food competition was more intense in the past than at present. However, the putative association between diet and tooth fracture frequency was based on a small sample of large, highly carnivorous species. In the present study, a greater diversity of extant carnivorans is sampled, including insectivorous, omnivorous, and carnivorous forms, ranging in size from weasels to tigers. Species that habitually consume hard foods (bones, shells) had the highest fracture frequencies, followed by carnivores, and then insectivorous and/or omnivorous species. Predator and prey sizes were not associated with tooth fracture frequency, but more aggressive species did break their teeth more often. Comparison of the modern sample with five Pleistocene species confirms the previous finding of higher tooth breakage in the past, although some extant species have fracture frequencies that approach those of extinct species. Thus, the Pleistocene predator guild appears to have been characterized by relatively high levels of competition that are rarely observed today.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 68–81.     

Using geometric morphometrics for integrative taxonomy: an examination of head shapes of milksnakes (genus Lampropeltis)

Species discovery and identification has long relied on traditional morphometric analyses, although molecular methods for species delimitation are becoming increasing popular and important. Despite an increase in studies that rely solely on molecular data to differentiate between species, additional evidence that supports genealogically‐based species delimitation is desirable at least for field and museum identification of species and is part of an integrative approach to taxonomy. The present study uses geometric morphometric (GM) analyses to examine six species of milksnake (genus Lampropeltis) that have recently been delimited based on multilocus data in a coalescent framework. Landmarks are plotted onto the dorsal view of 487 specimens and canonical variate analysis (CVA) is used to determine whether the differences in head shape of these six species can be used to correctly classify specimens. For five of the six species, CVA accurately classifies individuals >70% of the time. The present study illustrates that, although GM‐based analyses may not correctly differentiate between species 100% of the time, GM methods can be useful for detecting shape differences between species and help to corroborate species delimitation. © 2015 The Linnean Society of London     

Molecular and morphological delimitation of Australian <Emphasis Type="Italic">Triops</Emphasis> species (Crustacea: Branchiopoda: Notostraca)—large diversity and little morphological differentiation

Australia has a very rich and diverse large branchiopod fauna with approximately 140 described or provisionally delimited species, but only one species of Triops, Triops australiensis (Spencer and Hall 1895), is currently recognized. Previous studies identified extensive genetic diversity within T. australiensis that suggested the presence of cryptic species. Herein, we employed an integrative approach to taxonomy to delimit putative species, integrating COI and EF1α sequence data and morphological data. Putative species were initially delimited based on COI by two computational approaches (GMYC and ABGD). The results were interpreted in the light of several species concepts, with particular emphasis on reproductive isolation. Twenty to 27 genetic lineages were delimited. Of these, up to 26 represent species following an evolutionary or phylogenetic species concept. Eighteen are biological species, though reproductive isolation could not be unambiguously established for allopatric species or species without known males. The level of co-occurrences was exceptionally high for Triops, with up to three syntopic and six sympatric species. Species delimitation was impeded by extensive overlap between intraspecific variability and interspecific variation in the genetic as well as morphological datasets. Without prior delimitation of putative species via COI, morphological delimitation would have been impossible. A potential explanation for the morphological variability is the retention of ancestral polymorphisms over long periods of time and across multiple speciation events without subsequent differentiation.     

Cytological differences and chromosomal rearrangements in four members of the Anopheles dirus complex (Diptera: Culicidae)

A reference photomap of the larval salivary gland, polytene chromosomes of the Anopheles dirus complex (species A) is presented. Samples of species A, B, C, and D from natural populations in Thailand were compared to this standard map using the larval progeny of wild-caught females. All species show differences in their chromosome banding patterns involving band size, number, and shape, particularly at the free ends of the X, 2R, and 2L. These differences provide useful diagnostic characters for separating members of the species complex. However, overall banding patterns are conservative in the group: species A, B, and C are virtually homosequential. Species D is highly polymorphic for a single paracentric inversion in each of the four autosomal arms and has a fixed inversion on the X chromosome. This same X chromosome inversion occurs at low frequency in species A.     

Is population subdivision different from speciation? From phylogeography to species delimitation

Species‐level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model‐based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.     

Species delimitation: inferring gaps in morphology across geography

Species are commonly delimited on the basis of gaps in patterns of morphological variation, but there seems to be little recent work on methods to objectively assess such gaps. Here, we introduce a statistical approach that uses measurements of continuous morphological characters and geographic variation in those characters to (i) measure the strength of the evidence for the existence of a gap in morphological variation between two hypothesized species and (ii) examine if a gap in morphological variation between two hypothesized species can be explained by an alternative hypothesis of geographic variation within a species. This approach is based on recent developments in analyses of multivariate normal mixtures, estimates of multivariate tolerance regions, and principal coordinates of neighboring matrices. We demonstrate the application of the approach by examining previously proposed hypotheses of species limits in the plant genus Escallonia. We discuss the main features of the method, including potential limitations, in relation to other approaches that use gaps in morphological variation as a criterion for species delimitation. The method we propose can help strengthen the link between the theory and practice of species delimitation by increasing the transparency and consistency of taxonomic decisions based on morphology, thus contributing to integrative approaches for species delimitation that consider morphological and geographic data on an equal footing with other kinds of information.     

Population-based resequencing of experimentally evolved populations reveals the genetic basis of body size variation in Drosophila melanogaster

Body size is a classic quantitative trait with evolutionarily significant variation within many species. Locating the alleles responsible for this variation would help understand the maintenance of variation in body size in particular, as well as quantitative traits in general. However, successful genome-wide association of genotype and phenotype may require very large sample sizes if alleles have low population frequencies or modest effects. As a complementary approach, we propose that population-based resequencing of experimentally evolved populations allows for considerable power to map functional variation. Here, we use this technique to investigate the genetic basis of natural variation in body size in Drosophila melanogaster. Significant differentiation of hundreds of loci in replicate selection populations supports the hypothesis that the genetic basis of body size variation is very polygenic in D. melanogaster. Significantly differentiated variants are limited to single genes at some loci, allowing precise hypotheses to be formed regarding causal polymorphisms, while other significant regions are large and contain many genes. By using significantly associated polymorphisms as a priori candidates in follow-up studies, these data are expected to provide considerable power to determine the genetic basis of natural variation in body size.     

Integration of conflict into integrative taxonomy: fitting hybridization in species delimitation of Mesocarabus (Coleoptera: Carabidae)

In species differentiation, characters may not diverge synchronously, and there are also processes that shuffle character states in lineages descendant from a common ancestor. Species are thus expected to show some degree of incongruence among characters; therefore, taxonomic delimitation can benefit from integrative approaches and objective strategies that account for character conflict. We illustrate the potential of exploiting conflict for species delimitation in a study case of ground beetles of the subgenus Carabus (Mesocarabus), where traditional taxonomy does not accurately delimit species. The molecular phylogenies of four mitochondrial and three nuclear genes, cladistic analysis of the aedeagus, ecological niche divergence and morphometry of pronotal shape in more than 500 specimens of Mesocarabus show that these characters are not fully congruent. For these data, a three‐step operational strategy is proposed for species delimitation by (i) delineating candidate species based on the integration of incongruence among conclusive lines of evidence, (ii) corroborating candidate species with inconclusive lines of evidence and (iii) refining a final species proposal based on an integrated characterization of candidate species based on the evolutionary analysis of incongruence. This procedure provided a general understanding of the reticulate process of hybridization and introgression acting on Mesocarabus and generated the hypothesis of seven Mesocarabus species, including two putative hybrid lineages. Our work emphasizes the importance of incorporating critical analyses of character and phylogenetic conflict to infer both the evolutionary history and species boundaries through an integrative taxonomic approach.