The utility of automated analysis of inter-simple sequence repeat (ISSR) loci for resolving relationships in the Canary Island species of Tolpis (Asteraceae)

Plants of oceanic islands, often remarkably divergent morphologically from continental relatives, are useful models for studying evolution and speciation because evolution is telescoped in time and space. Prior studies revealed little DNA sequence variation within the clade of ca. 10 Canary Island species of Tolpis, which precluded resolving species relationships. The present study assessed the utility of automated analysis of inter-simple sequence repeat (ISSR) loci for resolving relationships within the clade using 264 individuals from 36 populations of all recognized species and three undescribed morphological variants. Similarity (Dice coefficient) and Fitch parsimony were used to generate neighbor-joining (NJ) and strict consensus trees (MP), respectively. All individuals of the morphologically distinct endemic species formed clusters in both trees. There is also support for clusters of two undescribed variants in the NJ tree. Individuals from a morphologically variable complex consisting primarily of two species are not well resolved at population or species levels. The NJ and MP trees are not congruent at deeper levels, including relationships among species. Results are interpreted in terms of the biology of the species, and the utility of automated analysis of ISSR markers for interpreting patterns of evolution of Tolpis in the Canary Islands is discussed.     

Testing the potential of DNA barcoding in vertebrate radiations: the case of the littoral cichlids (Pisces,Perciformes, Cichlidae) from Lake Tanganyika

We obtained 398 cytochrome c oxidase subunit I barcodes of 96 morphospecies of Lake Tanganyika (LT) cichlids from the littoral zone. The potential of DNA barcoding in these fishes was tested using both species identification and species delineation methods. The best match (BM) and best close match (BCM) methods were used to evaluate the overall identification success. For this, three libraries were analysed in which the specimens were categorized into Operational Taxonomic Units (OTU) in three alternative ways: (A) morphologically distinct, including undescribed, species, (B) valid species and (C) complexes of morphologically similar or closely related species. For libraries A, B and C, 73, 73 and 96% (BM) and 72, 70 and 94% (BCM) of the specimens were correctly identified. Additionally, the potential of two species delineation methods was tested. The General Mixed Yule Coalescent (GMYC) analysis suggested 70 hypothetical species, while the Automatic Barcode Gap Discovery (ABGD) method revealed 115 putative species. Although the ABGD method had a tendency to oversplit, it outperformed the GMYC analysis in retrieving the species. In most cases where ABGD suggested oversplitting, this was due to intraspecific geographical variation. The failure of the GMYC method to retrieve many species could be attributed to discrepancies between mitochondrial gene trees and the evolutionary histories of LT cichlid species. Littoral LT cichlids have complex evolutionary histories that include instances of hybridization, introgression and rapid speciation. Nevertheless, although the utility of DNA barcoding in identification is restricted to the level of complexes, it has potential for species discovery in cichlid radiations.     

DNA-based species delineation in tropical beetles using mitochondrial and nuclear markers

DNA barcoding has been successfully implemented in the identification of previously described species, and in the process has revealed several cryptic species. It has been noted that such methods could also greatly assist in the discovery and delineation of undescribed species in poorly studied groups, although to date the feasibility of such an approach has not been examined explicitly. Here, we investigate the possibility of using short mitochondrial and nuclear DNA sequences to delimit putative species in groups lacking an existing taxonomic framework. We focussed on poorly known tropical water beetles (Coleoptera: Dytiscidae, Hydrophilidae) from Madagascar and dung beetles (Scarabaeidae) in the genus Canthon from the Neotropics. Mitochondrial DNA sequence variation proved to be highly structured, with >95% of the observed variation existing between discrete sets of very closely related genotypes. Sequence variation in nuclear 28S rRNA among the same individuals was lower by at least an order of magnitude, but 16 different genotypes were found in water beetles and 12 genotypes in Canthon, differing from each other by a minimum of two base pairs. The distribution of these 28S rRNA genotypes in individuals exactly matched the distribution of mtDNA clusters, suggesting that mtDNA patterns were not misleading because of introgression. Moreover, in a few cases where sequence information was available in GenBank for morphologically defined species of Canthon, these matched some of the DNA-based clusters. These findings demonstrate that clusters of close relatives can be identified readily in the sequence variation obtained in field collected samples, and that these clusters are likely to correspond to either previously described or unknown species. The results suggest that DNA-assisted taxonomy will not require more than a short fragment of mtDNA to provide a largely accurate picture of species boundaries in these groups. Applied on a large scale, this DNA-based approach could greatly improve the rate of species discovery in the large assemblages of insects that remain undescribed.     

Radiation, diversity, and host-plant interactions among island and continental legume-feeding psyllids

Island archipelagos and insect-plant associations have both independently provided many useful systems for evolutionary study. The arytainine psyllid (Sternorrhyncha: Hemiptera) radiation on broom (Fabaceae: Genisteae) in the Canary Island archipelago provides a discrete system for examining the speciation of highly host-specific phytophagous insects in an island context. Phylogenetic reconstructions based on three datasets (adult and nymph morphological characters, and two mitochondrial DNA regions: part of the small subunit rRNA, and part of cytochrome oxidase I, cytochrome oxidase II and the intervening tRNA leucine) are generally consistent. The combined molecular tree provides a well-supported estimate of psyllid relationships and shows that there have been several colonizations of the Macaronesian islands but that only one has resulted in a significant radiation. Psyllid diversification has apparently been constrained by the presence of suitable host groups within the genistoid legumes, and the diversity, distribution, and abundance of those groups. The phylogeny, by indicating pairs of sister species, allows putative mechanisms of speciation to be assessed. The most common conditions associated with psyllid speciation are geographical allopatry with a host switch to closely related hosts (six examples), or geographical allopatry on the same host (four examples). Where allopatric speciation involves a host switch, these have all been to related hosts. There is some evidence that switches between unrelated host plants may be more likely in sympatry. Only one sister pair (Aryrtainilla cytisi and A. telonicola) and the putative host races of Arytinnis modica are sympatric but on unrelated hosts, which may be a necessary condition for sympatric speciation in these insects. Where several psyllids share the same host, resources appear to be partitioned by ecological specialization and differing psyllid phenology.     

DNA-based taxonomy for associating adults and larvae in multi-species assemblages of chafers (Coleoptera: Scarabaeidae)

DNA sequences provide a universal character system in taxonomy for associating all developmental stages of organisms, but ambiguity arises due to genetic variation within species. The problem is compounded where target groups are less well studied or incompletely represented in DNA databases. Here we investigate the utility of DNA for larval-adult species associations within chafer (Coleoptera: Scarabaeidae) communities from four sites in the tropical lowlands of Nepal. We sequenced ca. 1600 bp of mitochondrial cox1 and rrnL and 700 bp of nuclear 28S rRNA from 250 larval and adult specimens. Individuals were grouped into putative species using statistical parsimony analysis and population aggregation analysis (PAA), whereby specimens from each locality were grouped according to the presence of diagnostic nucleotides. In addition, species membership was determined based on shifts in branching rates on clock-constrained trees to detect the putative transition from speciation to population coalescence patterns. Using these two methods we delineated between 48 and 56 groups, of which 16-20 were composed of larval and adult individuals. Nuclear and mtDNA-based groups were highly congruent; variation of 28S rRNA within groups was very low, while one widespread 28S rRNA genotype was universally found in a paraphyletic group of five mtDNA clusters. Linnean names could be assigned to 19 groups, and hence between 86.1% and 92.7% of larval specimens could be associated to species by their membership in clearly delineated groups that contained fully identified adults. The remaining larvae were delineated as five species, four of which could be assigned to Anomala or Adoretus based on their phylogenetic position. We conclude that the sequence variation was highly structured in this complex assemblage of chafers and that any given individual (larva or adult) can be readily associated with a particular DNA group using the criterion of diagnos ability. The association of different developmental stages therefore becomes a matter of determining the extent of the DNA-based groups, rather than matching of sequences from adult and larval individuals. This indicates the need for a purely sequence-based taxonomic system when associating different life stages via DNA.     

Mitochondrial DNA sequence variation in the spruce budworm species complex (Choristoneura: Lepidoptera)

A combination of polymerase-chain-reaction amplification and automated DNA sequencing was used to survey variation in a species complex of pest insects, the spruce budworms (Choristoneura fumiferana species group), and an outgroup species, C. rosaceana. We sequenced an mtDNA region of 1,573 bp that extends from the middle of cytochrome oxidase subunit I (COI) through tRNA leucine (UUR) to the end of cytochrome oxidase subunit II. In addition, we examined levels of intraspecific variation within a 470-bp region of the COI gene. Choristoneura fumiferana clearly represented the oldest lineage within its species group, with 2.7%-2.9% sequence divergence from the other species. In contrast, the four remaining species (C. pinus, C. biennis, C. occidentalis, and C. orae) had closely related or identical mtDNA, with < 1% divergence among most of their haplotypes. Despite its older lineage and widespread geographic distribution, C. fumiferana showed significantly lower intraspecific genetic diversity than did C. occidentalis. Choristoneura orae shared haplotypes with C. occidentalis and C. biennis, and species-level separation of these three species was not supported. Two divergent, uncommon haplotypes were also found in C. occidentalis and C. biennis. The divergent haplotype in C. biennis had an unusually high number of inferred amino acid replacements, suggesting selective differences between mitochondrial DNA haplotypes. Transition:transversion ratios in Choristoneura paralleled those found in Drosophila; transition:transversion ratios were highest in closely related sequences but decreased with increasing sequence divergence. Nucleotide composition showed an A+T bias that was near the high end of the range known for insects. This work illustrates the potential utility of direct DNA sequencing in assessing population structures, species limits, and phylogenetic relationships among organisms that have not previously been subjected to DNA analysis.      

Limited performance of DNA barcoding in a diverse community of tropical butterflies

DNA 'barcoding' relies on a short fragment of mitochondrial DNA to infer identification of specimens. The method depends on genetic diversity being markedly lower within than between species. Closely related species are most likely to share genetic variation in communities where speciation rates are rapid and effective population sizes are large, such that coalescence times are long. We assessed the applicability of DNA barcoding (here the 5' half of the cytochrome c oxidase I) to a diverse community of butterflies from the upper Amazon, using a group with a well-established morphological taxonomy to serve as a reference. Only 77% of species could be accurately identified using the barcode data, a figure that dropped to 68% in species represented in the analyses by more than one geographical race and at least one congener. The use of additional mitochondrial sequence data hardly improved species identification, while a fragment of a nuclear gene resolved issues in some of the problematic species. We acknowledge the utility of barcodes when morphological characters are ambiguous or unknown, but we also recommend the addition of nuclear sequence data, and caution that species-level identification rates might be lower in the most diverse habitats of our planet.     

Molecular phylogenetics at the population/species interface in cave spiders of the southern Appalachians (Araneae:Nesticidae:Nesticus)

This paper focuses on the relationship between population genetic structure and speciation mechanisms in a monophyletic species group of Appalachian cave spiders (Nesticus). Using mtDNA sequence data gathered from 256 individuals, I analyzed patterns of genetic variation within and between populations for three pairs of closely related sister species. Each sister-pair comparison involves taxa with differing distributional and ecological attributes; if these ecological attributes are reflected in basic demographic differences, then speciation might proceed differently across these sister taxa comparisons. Both frequency-based and gene tree analyses reveal that the genetic structure of the Nesticus species studied is characterized by similar and essentially complete population subdivision, regardless of differences in general ecology. These findings contrast with results of prior genetic studies of cave-dwelling arthropods that have typically revealed variation in population structure corresponding to differences in general ecology. Species fragmentation through both extrinsic and intrinsic evolutionary forces has resulted in discrete, perhaps independent, populations within morphologically defined species. Large sequence divergence values observed between populations suggest that this independence may extend well into the past. These patterns of mtDNA genealogical structure and divergence imply that species as morphological lineages are currently more inclusive than basal evolutionary or phylogenetic units, a suggestion that has important implications for the study of speciation mechanisms.      

Molecular evidence for ten species and Oligo-Miocene vicariance within a nominal Australian gecko species (<Emphasis Type="Italic">Crenadactylus ocellatus</Emphasis>, Diplodactylidae)

Background  

Molecular studies have revealed that many putative 'species' are actually complexes of multiple morphologically conservative, but genetically divergent 'cryptic species'. In extreme cases processes such as non-adaptive diversification (speciation without divergent selection) could mask the existence of ancient lineages as divergent as ecologically and morphologically diverse radiations recognised as genera or even families in related groups. The identification of such ancient, but cryptic, lineages has important ramifications for conservation, biogeography and evolutionary biology. Herein, we use an integrated multilocus genetic dataset (allozymes, mtDNA and nuclear DNA) to test whether disjunct populations of the widespread nominal Australian gecko species Crenadactylus ocellatus include distinct evolutionary lineages (species), and to examine the timing of diversification among these populations.     

Resolution of the Acanthopagrus black seabream complex based on mitochondrial and amplified fragment-length polymorphism analyses

In this study, DNA analyses were employed to verify the identity of six morphologically similar species that occur in the coastal waters of Taiwan: the black seabream complex (Acanthopagrus latus, Acanthopagrus schlegelii, Acanthopagrus sivicolus, Acanthopagrus taiwanensis, Acanthopagrus chinshira and Acanthopagrus pacificus). Amplified fragment-length polymorphism (AFLP) analyses clearly distinguished the same six species that are morphologically diagnosable based on subtle differences in scale counts and anal-fin colouration. In contrast, mitochondrial DNA analyses based on cytochrome b gene sequences did not distinguish individuals of A. schlegelii and A. sivicolus, reflecting either historical introgression or recent speciation and incomplete sorting of their mitochondrial lineages. Phylogenetic relationships among these six north-west Pacific Ocean species of Acanthopagrus analysed using AFLP data were consistent with scale rows above the lateral line (TRac), sperm ultrastructure and geographical distribution. The study provides molecular tools for future research relevant to improved management of these resources, and an increased understanding of the evolutionary history of this radiation.     

Deep genetic divergences among morphologically similar and parapatric Skistodiaptomus (Copepoda: Calanoida: Diaptomidae) challenge the hypothesis of Pleistocene speciation

We used mitochondrial [cytochrome c oxidase subunit I (CO I ), cytochrome b , and 16S] and nuclear [internal transcribed spacer (ITS) phylogenies of Skistodiaptomus copepods to test hypotheses of Pleistocene divergence and speciation within the genus. Mitochondrial (mt)DNA sequence divergences do not support hypotheses for Pleistocene speciation and instead suggest much more ancient speciation events in the genus. Skistodiaptomus oregonensis and Skistodiaptomus pygmaeus (i.e. two morphologically similar and parapatric species) exhibited uncorrected mtDNA sequence divergences exceeding 20%. Similarly, we identified three divergent clades of Skistodiaptomus pallidus that exhibited mtDNA sequence divergences exceeding 15%, suggesting that even intraspecific divergence within this morphospecies predates the Pleistocene. We found clear evidence of CO I pseudogenes in S. pygmaeus , but their presence did not lead to significant overestimates of sequence divergences for this gene. Substitution saturation and strong purifying selection have most likely led to underestimates of sequence divergences and divergence times among Skistodiaptomus . The widespread phenomenon of morphological stasis among genetically divergent copepod groups indicates that speciation often occurs with little or no morphological change. Instead, morphological evolution may occur idiosyncratically after speciation and create discordant patterns of morphological similarity, shared ancestry and divergence time. Cryptic species complexes are therefore common in copepods, and morphological species concepts underestimate their true species diversity.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 150–165.     

HYBRID SPECIATION AND INDEPENDENT EVOLUTION IN LINEAGES OF ALPINE BUTTERFLIES

The power of hybridization between species to generate variation and fuel adaptation is poorly understood despite long‐standing interest. There is, however, increasing evidence that hybridization often generates biodiversity, including via hybrid speciation. We tested the hypothesis of hybrid speciation in butterflies occupying extreme, high‐altitude habitats in four mountain ranges in western North America with an explicit, probabilistic model, and genome‐wide DNA sequence data. Using this approach, in concert with ecological experiments and observations and morphological data, we document three lineages of hybrid origin. These lineages have different genome admixture proportions and distinctive trait combinations that suggest unique and independent evolutionary histories.     

Coevolution of Retrovarium n. gen. (Digenea: Cryptogonimidae) in Lutjanidae and Haemulidae (Perciformes) in the Indo-West Pacific

We describe 11 new species of cryptogonimids belonging to Retrovarium n. gen., from eight species of Lutjanidae and one species of Haemulidae, from the Great Barrier Reef, French Polynesia and the Maldives. We also transfer Neoparacryptogonimus saccatus (Manter, 1963) and Neoparacryptogonimus sphericus Nahhas, Sey & Nishimoto, 1998 to Retrovarium. The morphologically based taxonomic approach was augmented with DNA sequence data from three nuclear ribosomal DNA regions (28S, ITS1 and ITS2) to explore the species integrity, biogeographic distribution and evolution of the species recognised here. Sequencing included multiple replicates and revealed 11 distinct genotypes which corroborated our morphologically based hypotheses of putative species present in the system. There was no intraspecific variation and all three rDNA regions differed between every combination of species. Two species exhibited wide geographic ranges, having identical rDNA sequences between the Great Barrier Reef and the Maldives, localities separated by over 9600 km. One host species, Symphorus nematophorus, proved to be exceptionally rich, harbouring six species. Minimum evolution analyses were conducted on each of the rDNA datasets independently; minimum evolution, maximum likelihood and Bayesian inference analyses were conducted on a combined (28S, ITS1 and ITS2) dataset for sequence comparison purposes and to explore the evolutionary history of these parasites. To examine the coevolutionary history of this complex, assessment of phylogenetic relationships between the 23 species of Lutjanidae and two species of Haemulidae collected during this survey was performed with data from 16S and cytochrome b mtDNA using Bayesian inference analysis. Despite the high host specificity observed in most of the species, mapping of the parasites on the host phylogeny revealed an absence of strict coevolution or co-descent within this complex. Overall, Retrovarium appears to have had an exceptionally patchy radiation, failing to infect many taxa, infecting species with no readily discernible pattern, and radiating dramatically within one species.     

The genetics and ecology of sympatric speciation: A case study

Mating occurs on the larval host plant in allRhagoletis species (Diptera: Tephritidae). We show how this attribute, when coupled with certain differences in other biological traits, strongly influences the mode of speciation. In species of thesuavis species group, host shifts have never occurred during speciation, and larvae feed in the husks of any walnut species(Juglans spp.), which are highly toxic. Taxa are allopatric or parapatric and exhibit deep phylogenetic nodes suggesting relatively ancient speciation events. Traits responsible for species and mate recognition, particularly in parapatric species, are morphologically distinct and strongly sexually dimorphic. All aspects of their biology, genetics and distribution are consistent with a slow rate of allopatric speciation followed by morphological divergence in secondary contact. In contrast, speciation in thepomonella species group has always involved a shift to a new, usually unrelated, non-toxic host, and all taxa within these groups are sympatric, monophagous and morphologically indistinguishable from one another. Phylogenetic nodes are very shallow, indicating recent sympatric speciation. Sympatric divergence is promoted by genetic variation which allows a portion of the original species to shift to a new habitat or host. Evidence suggests that changes in a few key loci responsible for host selection and fitness on a new host may initiate host shifts. By exploiting different habitats, competition for resources between diverging populations is reduced or avoided. We provide evidence that in phytophagous and parasitic insects sufficient intrinsic barriers to gene flow can evolve between sister populations as they adapt to different habitats or hosts to allow each population to establish independent evolutionary lineages in sympatry.     

Coexisting Cyclic Parthenogens Comprise a Holocene Species Flock in Eubosmina

Background

Mixed breeding systems with extended clonal phases and weak sexual recruitment are widespread in nature but often thought to impede the formation of discrete evolutionary clusters. Thus, cyclic parthenogens, such as cladocerans and rotifers, could be predisposed to “species problems” and a lack of discrete species. However, species flocks have been proposed for one cladoceran group, Eubosmina, where putative species are sympatric, and there is a detailed paleolimnological record indicating a Holocene age. These factors make the Eubosmina system suitable for testing the hypotheses that extended clonal phases and weak sexual recruitment inhibit speciation. Although common garden experiments have revealed a genetic component to the morphotypic variation, the evolutionary significance of the morphotypes remains controversial.

Methodology/Principal Findings

In the present study, we tested the hypothesis of a single polymorphic species (i.e., mixing occurs but selection maintains genes for morphology) in four northern European lakes where the morphotypes coexist. Our evidence is based on nuclear DNA sequence, mitochondrial DNA sequence, and morphometric analysis of coexisting morphotypes. We found significant genetic differentiation, genealogical exclusivity, and morphometric differentiation for coexisting morphotypes.

Conclusions

We conclude that the studied morphotypes represent a group of young species undergoing speciation with apparent reproductive barriers despite coexistence in the freshwater pelagic zone.     

Mating trials validate the use of DNA barcoding to reveal cryptic speciation of a marine bryozoan taxon

Despite increasing threats to the marine environment, only a fraction of the biodiversity of the oceans has been described, owing in part to the widespread occurrence of cryptic species. DNA-based barcoding through screening of an orthologous reference gene has been proposed as a powerful tool to uncover biological diversity in the face of dwindling taxonomic expertise and the limitations of traditional species identification. Although DNA barcoding should be particularly useful in the sea, given the prevalence of marine cryptic species, the link between taxa identified through DNA barcodes and reproductively isolated taxa (biological species) has rarely been explicitly tested. Here, we use an integrated framework comparing breeding compatibility, morphology and mitochondrial (cytochrome c oxidase 1) and nuclear (elongation factor-1-alpha) DNA sequence variation among globally distributed samples of the cosmopolitan marine bryozoan Celleporella hyalina (L.). Our results reveal that C. hyalina comprises numerous deep, mostly allopatric, genetic lineages that are reproductively isolated, yet share very similar morphology, indicating rampant cryptic speciation. The close correspondence between genetic lineages and reproductively isolated taxa in the context of minimal morphological change suggests that DNA barcoding will play a leading role in uncovering the hidden biodiversity of the oceans and that the sole use of morphologically based taxonomy would grossly underestimate the number of marine species.     

Natural hybridization and introgression in sympatric Ligularia species (Asteraceae, Senecioneae)

The difficulty in clarifying species of genus Ligularia Cass.has been attributed to rapid and continuous allopatric speciation in small and isolated populations,combined with interspecific diploid hybridization in the Qinghai-Tibetan Plateau and adjacent areas.However,no concrete example has been reported to prove this hypothesis.We studied a natural mixed population of six species of Ligularia in which some individuals were morphologically intermediate between L.subspicata and L.nelumbifolia.Based on DNA sequences (trnC-E trnL-rpL32,trnQ5'rpsl6,trnK-rpsl6,and internal transcribed spacer) and inter-simple sequence repeat data,we concluded that putative hybrids are primarily products of hybridization between L.nelumbifolia and L.subspicata.The other four species or additional,unknown species may also be involved in hybridization.This hybridization is bidirectional but asymmetrical.Hybrid individuals were mostly the first generation,but F2 and later-generation hybrids were also present.Moreover,the backcrossed individuals detected indicate that natural gene flow occurs among at least three Ligularia species.Hybrids may become stabilized to form new species or may function as intermediates in evolutionary diversification.     

Rampant host- and defensive phenotype-associated diversification in a goldenrod gall midge

Natural selection can play an important role in the genetic divergence of populations and their subsequent speciation. Such adaptive diversification, or ecological speciation, might underlie the enormous diversity of plant-feeding insects that frequently experience strong selection pressures associated with host plant use as well as from natural enemies. This view is supported by increasing documentation of host-associated (genetic) differentiation in populations of plant-feeding insects using alternate hosts. Here, we examine evolutionary diversification in a single nominal taxon, the gall midge Asteromyia carbonifera (O.S.), with respect to host plant use and gall phenotype. Because galls can be viewed as extended defensive phenotypes of the midges, gall morphology is likely to be a reflection of selective pressures by enemies. Using phylogenetic and comparative analyses of mtDNA and nuclear sequence data, we find evidence that A. carbonifera populations are rapidly diversifying along host plant and gall morphological lines. At a broad scale, geography explains surprisingly little genetic variation, and there is little evidence of strict co-cladogenesis with their Solidago hosts. Gall morphology is relatively labile, distinct gall morphs have evolved repeatedly and colonized multiple hosts, and multiple genetically and morphologically distinct morphs frequently coexist on a single host plant species. These results suggest that Asteromyia carbonifera is in the midst of an adaptive radiation driven by multitrophic selective pressures. Similar complex community pressures are likely to play a role in the diversification of other herbivorous insect groups.