Understanding the spectacular failure of DNA barcoding in willows (Salix): Does this result from a trans-specific selective sweep?

Multi-locus phylogenetics, lineage sorting, and reticulation in Pinus subsection Australes

PREMISE OF THE STUDY

Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long-lived, wind-pollinated plants with large population sizes and weak reproductive barriers.

METHODS

We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low-copy-number nuclear genes and high-copy-number plastomes (Hyb-Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation.

KEY RESULTS

Concatenation- and coalescent-based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed-cone pines) and Oocarpae (the egg-cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non-monophyly of the Attenuatae in the plastid tree was equivocal.

CONCLUSIONS

Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.

Complete chloroplast genome of seven Fritillaria species, variable DNA markers identification and phylogenetic relationships within the genus

Fritillaria spp. constitute important traditional Chinese medicinal plants. Xinjiang is one of two diversity hotspots in China in which eight Fritillaria species occur, two of which are endemic to the region. Furthermore, the phylogenetic relationships of Xinjiang Fritillaria species (including F. yuminensis) within the genus are unclear. In the present study, we sequenced the chloroplast (cp) genomes of seven Fritillaria species in Xinjiang using the Illumina HiSeq platform, with the aim of assessing the global structural patterns of the seven cp genomes and identifying highly variable cp DNA sequences. These were compared to previously sequenced Fritillaria cp genomes. Phylogenetic analysis was then used to evaluate the relationships of the Xinjiang species and assess the evolution of an undivided stigma. The seven cp genomes ranged from 151,764 to 152,112 bp, presenting a traditional quadripartite structure. The gene order and gene content of the seven cp genomes were identical. A comparison of the 13 cp genomes indicated that the structure is highly conserved. Ten highly divergent regions were identified that could be valuable in phylogenetic and population genetic studies. The phylogenetic relationships of the 13 Fritillaria species inferred from the protein-coding genes, large single-copy, small single-copy, and inverted repeat regions were identical and highly resolved. The phylogenetic relationships of the species corresponded with their geographic distribution patterns, in that the north group (consisting of eight species from Xinjiang and Heilongjiang in North China) and the south group (including six species from South China) were basically divided at 40°N. Species with an undivided stigma were not monophyletic, suggesting that this trait might have evolved several times in the genus.

A time and a place for everything: phylogenetic history and geography as joint predictors of oak plastome phylogeny

Owing to high rates of introgressive hybridization, the plastid genome is poorly suited to fine-scale DNA barcoding and phylogenetic studies of the oak genus (Quercus, Fagaceae). At the tips of the oak plastome phylogeny, recent gene migration and reticulation generally cause topology to reflect geographic structure, while deeper branches reflect lineage divergence. In this study, we quantify the simple and partial effects of geographic proximity and nucleome-inferred phylogenetic history on oak plastome phylogeny at different evolutionary scales. Our study compares pairwise phylogenetic distances based on complete plastome sequences, pairwise phylogenetic distances from nuclear restriction site-associated DNA sequences (RADseq), and pairwise geographic distances for 34 individuals of the white oak clade representing 24 North American and Eurasian species. Within the North American white oak clade alone, phylogenetic history has essentially no effect on plastome variation, while geography explains 11%–21% of plastome phylogenetic variance. However, across multiple continents and clades, phylogeny predicts 30%–41% of plastome variation, geography 3%–41%. Tipwise attenuation of phylogenetic informativeness in the plastome means that in practical terms, plastome data has little use in solving phylogenetic questions, but can still be a useful barcoding or phylogenetic marker for resolving questions among major clades.

What can we do with 1000 plastid genomes?

The plastid genome of plants is the smallest and most gene-rich of the three genomes in each cell and the one generally present in the highest copy number. As a result, obtaining plastid DNA sequence is a particularly cost-effective way of discovering genetic information about a plant. Until recently, the sequence information gathered in this way was generally limited to small portions of the genome amplified by polymerase chain reaction, but recent advances in sequencing technology have stimulated a substantial rate of increase in the sequencing of complete plastid genomes. Within the last year, the number of complete plastid genomes accessible in public sequence repositories has exceeded 1000. This sudden flood of data raises numerous challenges in data analysis and interpretation, but also offers the keys to potential insights across large swathes of plant biology. We examine what has been learnt so far, what more could be learnt if we look at the data in the right way, and what we might gain from the tens of thousands more genome sequences that will surely arrive in the next few years. The most exciting new discoveries are likely to be made at the interdisciplinary interfaces between molecular biology and ecology.

DNA barcoding evaluation and implications for phylogenetic relationships in Lauraceae from China

Lauraceae are an important component of tropical and subtropical forests and have major ecological and economic significance. Owing to lack of clear-cut morphological differences between genera and species, this family is an ideal case for testing the efficacy of DNA barcoding in the identification and discrimination of species and genera. In this study, we evaluated five widely recommended plant DNA barcode loci matK, rbcL, trnH–psbA, ITS2 and the entire ITS region for 409 individuals representing 133 species, 12 genera from China. We tested the ability of DNA barcoding to distinguish species and as an alternative tool for correcting species misidentification. We also used the rbcL+matK+trnH–psbA+ITS loci to investigate the phylogenetic relationships of the species examined. Among the gene regions and their combinations, ITS was the most efficient for identifying species (57.5%) and genera (70%). DNA barcoding also had a positive role for correcting species misidentification (10.8%). Furthermore, based on the results of the phylogenetic analyses, Chinese Lauraceae species formed three supported monophyletic clades, with the Cryptocarya group strongly supported (PP = 1.00, BS = 100%) and the clade including the Persea group, Laureae and Cinnamomum also receiving strong support (PP = 1.00, BS = 98%), whereas the Caryodaphnopsis–Neocinnamomum received only moderate support (PP = 1.00 and BS = 85%). This study indicates that molecular barcoding can assist in screening difficult to identify families like Lauraceae, detecting errors of species identification, as well as helping to reconstruct phylogenetic relationships. DNA barcoding can thus help with large-scale biodiversity inventories and rare species conservation by improving accuracy, as well as reducing time and costs associated with species identification.

DNA barcoding and species delimitation of Chaitophorinae (Hemiptera, Aphididae)

Chaitophorinae aphids are widespread across Eurasia and North America, and include some important agricultural and horticultural pests. So, accurate rapid species identification is very important. Here, we used three mitochondrial genes and one endosymbiont gene to calculate and analyze the genetic distances within different datasets. For species delimitation, two distance-based methods were employed, threshold with NJ (neighbor-joining) and ABGD (Automatic Barcode Gap Discovery), and two tree-based approaches, GMYC (General Mixed Yule Coalescent) and PTP (Poisson Tree Process). The genetic interspecific divergence was clearly larger than the intraspecific divergence for four molecular markers. COI and COII genes were found to be more suitable for Chaitophorinae DNA barcoding. For species delimitation, at least one distance-based method combined with one tree-based method would be preferable. Based on the data for Chaitophorussaliniger and Laingiapsammae, DNA barcoding may also reveal geographical variation.

Testing the Efficacy of DNA Barcodes for Identifying the Vascular Plants of Canada

Their relatively slow rates of molecular evolution, as well as frequent exposure to hybridization and introgression, often make it difficult to discriminate species of vascular plants with the standard barcode markers (rbcL, matK, ITS2). Previous studies have examined these constraints in narrow geographic or taxonomic contexts, but the present investigation expands analysis to consider the performance of these gene regions in discriminating the species in local floras at sites across Canada. To test identification success, we employed a DNA barcode reference library with sequence records for 96% of the 5108 vascular plant species known from Canada, but coverage varied from 94% for rbcL to 60% for ITS2 and 39% for matK. Using plant lists from 27 national parks and one scientific reserve, we tested the efficacy of DNA barcodes in identifying the plants in simulated species assemblages from six biogeographic regions of Canada using BLAST and mothur. Mean pairwise distance (MPD) and mean nearest taxon distance (MNTD) were strong predictors of barcode performance for different plant families and genera, and both metrics supported ITS2 as possessing the highest genetic diversity. All three genes performed strongly in assigning the taxa present in local floras to the correct genus with values ranging from 91% for rbcL to 97% for ITS2 and 98% for matK. However, matK delivered the highest species discrimination (~81%) followed by ITS2 (~72%) and rbcL (~44%). Despite the low number of plant taxa in the Canadian Arctic, DNA barcodes had the least success in discriminating species from this biogeographic region with resolution ranging from 36% with rbcL to 69% with matK. Species resolution was higher in the other settings, peaking in the Woodland region at 52% for rbcL and 87% for matK. Our results indicate that DNA barcoding is very effective in identifying Canadian plants to a genus, and that it performs well in discriminating species in regions where floristic diversity is highest.

From writing to reading the encyclopedia of life

Prologue 'As the study of natural science advances, the language of scientific description may be greatly simplified and abridged. This has already been done by Linneaus and may be carried still further by other invention. The descriptions of natural orders and genera may be reduced to short definitions, and employment of signs, somewhat in the manner of algebra, instead of long descriptions. It is more easy to conceive this, than it is to conceive with what facility, and in how short a time, a knowledge of all the objects of natural history may ultimately be acquired; and that which is now considered learning and science, and confined to a few specially devoted to it, may at length be universally possessed in every civilized country and in every rank of life'. J. C. Louden 1829. Magazine of natural history, vol. 1: This article is part of the themed issue 'From DNA barcodes to biomes'.

Telling plant species apart with DNA: from barcodes to genomes

Land plants underpin a multitude of ecosystem functions, support human livelihoods and represent a critically important component of terrestrial biodiversity-yet many tens of thousands of species await discovery, and plant identification remains a substantial challenge, especially where material is juvenile, fragmented or processed. In this opinion article, we tackle two main topics. Firstly, we provide a short summary of the strengths and limitations of plant DNA barcoding for addressing these issues. Secondly, we discuss options for enhancing current plant barcodes, focusing on increasing discriminatory power via either gene capture of nuclear markers or genome skimming. The former has the advantage of establishing a defined set of target loci maximizing efficiency of sequencing effort, data storage and analysis. The challenge is developing a probe set for large numbers of nuclear markers that works over sufficient phylogenetic breadth. Genome skimming has the advantage of using existing protocols and being backward compatible with existing barcodes; and the depth of sequence coverage can be increased as sequencing costs fall. Its non-targeted nature does, however, present a major informatics challenge for upscaling to large sample sets.This article is part of the themed issue 'From DNA barcodes to biomes'.

Understanding invasion history and predicting invasive niches using genetic sequencing technology in Australia: case studies from Cucurbitaceae and Boraginaceae

Part of the challenge in dealing with invasive plant species is that they seldom represent a uniform, static entity. Often, an accurate understanding of the history of plant introduction and knowledge of the real levels of genetic diversity present in species and populations of importance is lacking. Currently, the role of genetic diversity in promoting the successful establishment of invasive plants is not well defined. Genetic profiling of invasive plants should enhance our understanding of the dynamics of colonization in the invaded range. Recent advances in DNA sequencing technology have greatly facilitated the rapid and complete assessment of plant population genetics. Here, we apply our current understanding of the genetics and ecophysiology of plant invasions to recent work on Australian plant invaders from the Cucurbitaceae and Boraginaceae. The Cucurbitaceae study showed that both prickly paddy melon (Cucumis myriocarpus) and camel melon (Citrullus lanatus) were represented by only a single genotype in Australia, implying that each was probably introduced as a single introduction event. In contrast, a third invasive melon, Citrullus colocynthis, possessed a moderate level of genetic diversity in Australia and was potentially introduced to the continent at least twice. The Boraginaceae study demonstrated the value of comparing two similar congeneric species; one, Echium plantagineum, is highly invasive and genetically diverse, whereas the other, Echium vulgare, exhibits less genetic diversity and occupies a more limited ecological niche. Sequence analysis provided precise identification of invasive plant species, as well as information on genetic diversity and phylogeographic history. Improved sequencing technologies will continue to allow greater resolution of genetic relationships among invasive plant populations, thereby potentially improving our ability to predict the impact of these relationships upon future spread and better manage invaders possessing potentially diverse biotypes and exhibiting diverse breeding systems, life histories and invasion histories.

Evolutionary histories determine DNA barcoding success in vascular plants: seven case studies using intraspecific broad sampling of closely related species

Background

Four plastid regions, rpoB, rpoC1, matK, and trnH-psbA, have been recommended as DNA barcodes for plants. Their success in delimiting species boundaries depends on the existence of a clear-cut difference between inter- and intraspecific variability. We tested the ability of these regions to discriminate among closely related species in seven genera of flowering plants with different generation times (trees, perennials, and annuals). To ensure a maximum coverage of intraspecific diversity, and therefore to better evaluate the resolution power of each barcode, we applied a population genetics approach by sampling three to 45 individuals per species over a wide geographical range.

Results

All possible combinations between loci were analysed, which showed that using more than one locus does not always improve the resolution power. The trnH-psbA locus was most effective at discriminating among closely related species (Acer, Lonicera, Geranium, and Veronica), singly or in combination. For Salix, Adenostyles, and Gentiana, the best results were obtained with the combination of matK, rpoB, and trnH-psbA. No barcoding gap was found within six genera analysed, excepting Lonicera. This is due to shared polymorphisms among species, combined with very divergent sequences within species. These genetic patterns reflect incomplete lineage sorting and hybridization events followed by chloroplast capture.

Conclusions

Our results strongly suggest that adding trnH-psbA to the two obligate DNA barcodes proposed by the CBOL plant-working group (matK and rbcL) should be mandatory for closely related species. In our sampling, generation time had no influence on DNA barcoding success, as the best and worst identification successes were found for the two tree genera (Acer, 64 % success and Salix, 86 % failure). Evolutionary histories are the main factor influencing DNA barcoding success in the studied genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0678-0) contains supplementary material, which is available to authorized users.

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation-interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.

Assessing the performance of DNA barcoding using posterior predictive simulations

Accurate estimates of biodiversity are required for research in a broad array of biological subdisciplines including ecology, evolution, systematics, conservation and biodiversity science. The use of statistical models and genetic data, particularly DNA barcoding, has been suggested as an important tool for remedying the large gaps in our current understanding of biodiversity. However, the reliability of biodiversity estimates obtained using these approaches depends on how well the statistical models that are used describe the evolutionary process underlying the genetic data. In this study, we utilize data from the Barcode of Life Database and posterior predictive simulations to assess the performance of DNA barcoding under commonly used substitution models. We demonstrate that the success of DNA barcoding varies widely across DNA substitution models and that model choice has a substantial impact on the number of operational taxonomic units identified (changing results by ~4-31%). Additionally, we demonstrate that the widely followed practice of a priori assuming the Kimura 2-parameter model for DNA barcoding is statistically unjustified and should be avoided. Using both data-based and inference-based test statistics, we detect variation in model performance across taxonomic groups, clustering algorithms, genetic divergence thresholds and substitution models. Taken together, these results illustrate the importance of considering both model selection and model adequacy in studies quantifying biodiversity.

Promise and Challenge of DNA Barcoding in Venus Slipper (Paphiopedilum)

Orchidaceae are one of the largest families of flowering plants, with over 27,000 species described and all orchids are listed in CITES. Moreover, the seedlings of orchid species from the same genus are similar. The objective of DNA barcoding is rapid, accurate, and automated species identification, which may be used to identify illegally traded endangered species from vegetative specimens of Paphiopedilum (Venus slipper), a flagship group for plant conservation with high ornamental and commercial values. Here, we selected eight chloroplast barcodes and nrITS to evaluate their suitability in Venus slippers. The results indicate that all tested barcodes had no barcoding gap and the core plant barcodes showed low resolution for the identification of Venus slippers (18.86%). Of the single-locus barcodes, nrITS is the most efficient for the species identification of the genus (52.27%), whereas matK + atpF-atpH is the most efficient multi-locus combination (28.97%). Therefore, we recommend the combination of matK + atpF-atpH + ITS as a barcode for Venus slippers. Furthermore, there is an upper limit of resolution of the candidate barcodes, and only half of the taxa with multiple samples were identified successfully. The low efficiency of these candidate barcodes in Venus slippers may be caused by relatively recent speciation, the upper limit of the barcodes, and/or the sampling density. Although the discriminatory power is relatively low, DNA barcoding may be a promising tool to identify species involved in illegal trade, which has broad applications and is valuable for orchid conservation.

Salix transect of Europe: latitudinal patterns in willow diversity from Greece to arctic Norway

Background

Willows (Salix spp.) are ecosystem "foundation species" that are hosts to large numbers of associated insects. Determining their patterns of distribution across Europe is therefore of interest for understanding the spatial distribution of associated fauna. The aim of this study was to record species composition at multiple sites on a long latitudinal gradient (megatransect) across Europe as a baseline for the future detailed analysis of insect fauna at these sites. In this way we used willow stands as comparable mesocosms in which to study floristic and faunistic changes with latitude across Europe.

New information

To determine spatial patterning of  an ecologically important group on a latitudinal gradient across Europe, we sampled willows at the stand level in 42 sites, approximately 100 km apart, from the Aegean (38.8°N) to the Arctic Ocean (70.6°N), but at a similar longitude (21.2 to 26.1°E). The sites were predominantly lowland (elevations 1 to 556 metres amsl, median = 95 m) and wet (associated with rivers, lakes, drainage ditches or wet meadows). The median number of willow taxa (species and hybrids) per stand was four, and varied from one to nine. There is a progressive increase in willow diversity from south to north with the median number of taxa per stand in southern Europe being three, and in northern Europe six. A total of 20 willow species were recorded, along with 12 hybrids. The most widespread willow in the transect was Salixalba L. (occurring in 20 sites out of 42) followed by S.triandra L. (15 sites), S.caprea L., S.phylicifolia L. (14 sites) and S.myrsinifolia Salisb., Salix×fragilis L. (13 sites). Voucher specimens from this study are deposited in the herbaria of the Natural History Museum (BM) and the Royal Botanic Gardens Kew (K). These samples provide a "snapshot" of willow diversity along a latitudinal gradient and an indication of the geographically changing taxonomic diversity that is presented to willow-feeding herbivores across Europe. It is anticipated that further papers will examine the insect fauna collected from these sites as part of this study.

Genetic and morphological evidence for introgression between three species of willows

Background

Hybridization and introgression are said to occur relatively frequently in plants, and in particular among different species of willows. However, data on the actual frequency of natural hybridization and introgression is rare. Here, we report the first fine-scale genetic analysis of a contact zone shared between the three basket willow species, Salix dasyclados, S. schwerinii and S. viminalis in the vicinity of the Lake Baikal in Southern Siberia. Individuals were sampled in fourteen populations and classified as pure species or hybrids based on a set of morphological characters. They were then genotyped at 384 nuclear SNP and four chloroplast SSR loci. The STRUCTURE and NewHybrids softwares were used to estimate the frequency and direction of hybridization using genotypic data at the nuclear SNP loci.

Results

As many as 19 % of the genotyped individuals were classified as introgressed individuals and these were mainly encountered in the centre of the contact zone. All introgressed individuals were backcrosses to S. viminalis or S. schwerinii and no F1 or F2 hybrids were found. The rest of the genotyped individuals were classified as pure species and formed two clusters, one with S. schwerinii individuals and the other with S. viminalis and S. dasyclados individuals. The two clusters were significantly genetically differentiated, with F_ST = 0.333 (0.282–0.382, p < 0.001). In contrast, for the chloroplast haplotypes, no genetic differentiation was observed as they were completely shared between the species. Based on morphological classification only 5 % of the individuals were classified as introgressed individuals, which was much less than what was detected using genotypic data.

Conclusions

We have discovered a new willow hybrid zone with relatively high frequency of introgressed individuals. The low frequency of F1 hybrids indicates that ongoing hybridization is limited, which could be because of the presence of reproductive barriers or simply because the conditions are not favorable for hybridization. We further conclude that in order to get a complete picture of the species composition of a hybrid zone it is necessary to use a combination of morphological characters and genetic data from both nuclear and chloroplast markers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0461-7) contains supplementary material, which is available to authorized users.

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high-throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole-genome sequencing project, while the others were obtained by means of long-range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference-based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes ('Brindisino' artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty-nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short 'variable' cp regions. The evaluation of the molecular resources obtained from this study led us to support the 'super-barcode' theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.

The Miocene to Pleistocene colonization of the Philippine archipelago by Begonia sect. Baryandra (Begoniaceae)

PREMISE OF THE STUDY

One third of the species-rich Philippine flora is endemic, and most of the islands in the archipelago have never been connected to a continental region. We currently lack any well-sampled angiosperm phylogenies that span the archipelago, prohibiting the formation of informed hypotheses as to the evolution of this rich and highly endemic flora.

METHODS

We produced time-calibrated phylogenetic trees from both nuclear (ITS) and chloroplast (ndhA intron, ndhF-rpl32 spacer, rpl32-trnL spacer, trnC-trnD spacer) regions of 41 species of Begonia sect. Baryandra, all except one endemic to the Philippines. Historical biogeography was reconstructed across the chloroplast phylogeny using a Bayesian binary method of character optimization. Comparison of phylogenies from the two genomes permitted insight into the prevalence of hybridization in the group.

KEY RESULTS

The Philippine archipelago was colonized by Begonia sect. Baryandra in the late Miocene, via long-distance dispersal from western Malesia and a point of entry likely to be in the northwestern region of the archipelago. Palawan, Luzon, and Panay all bear early-branching lineages from this initial colonization. There have been Plio-Pleistocene dispersals from these islands into Borneo and Mindanao. Hybridization was common between species as evidenced by haplotype sharing and phylogenetic incongruence.

CONCLUSIONS

The phylogenies show a high degree of geographic structure, which millions of years of exposure to typhoons have not blurred, showing long-term species and population stability. The recent dispersals to Mindanao are congruent with the geologically recent arrival of the island at its current latitude in the southern Philippines.

Phylogenetic relationships of American willows (Salix L., Salicaceae)

Salix L. is the largest genus in the family Salicaceae (450 species). Several classifications have been published, but taxonomic subdivision has been under continuous revision. Our goal is to establish the phylogenetic structure of the genus using molecular data on all American willows, using three DNA markers. This complete phylogeny of American willows allows us to propose a biogeographic framework for the evolution of the genus. Material was obtained for the 122 native and introduced willow species of America. Sequences were obtained from the ITS (ribosomal nuclear DNA) and two plastid regions, matK and rbcL. Phylogenetic analyses (parsimony, maximum likelihood, Bayesian inference) were performed on the data. Geographic distribution was mapped onto the tree. The species tree provides strong support for a division of the genus into two subgenera, Salix and Vetrix. Subgenus Salix comprises temperate species from the Americas and Asia, and their disjunction may result from Tertiary events. Subgenus Vetrix is composed of boreo-arctic species of the Northern Hemisphere and their radiation may coincide with the Quaternary glaciations. Sixteen species have ambiguous positions; genetic diversity is lower in subg. Vetrix. A molecular phylogeny of all species of American willows has been inferred. It needs to be tested and further resolved using other molecular data. Nonetheless, the genus clearly has two clades that have distinct biogeographic patterns.

DNA barcoding evaluation and its taxonomic implications in the species-rich genus Primula L. in China

The genus Primula is extremely diverse in the east Himalaya-Hengduan Mountains (HHM) in China as a result of rapid radiation. In order to overcome the difficulty of morphological classification of this genus, we surveyed three plastid regions (rbcL, matK, and trnH-psbA) and two nuclear markers (ITS and ITS2) from 227 accessions representing 66 Primula species across 18 sections, to assess their discriminatory power as barcodes. We found that ITS alone or combined with plastid regions showed the best discrimination across different infrageneric ranks and at species level. We suggest rbcL + matK + ITS as the first choice at present to barcode Primula plants. Although the present barcoding combination performed poorly in many closely related species of Primula, it still provided many new insights into current Primula taxonomy, such as the underlying presence of cryptic species, and several potential improper taxonomic treatments. DNA barcoding is one useful technique in the integrative taxonomy of the genus Primula, but it still requires further efforts to improve its effectiveness in some taxonomically challenging groups.

Efficiency of ITS sequences for DNA barcoding in Passiflora (Passifloraceae)

DNA barcoding is a technique for discriminating and identifying species using short, variable, and standardized DNA regions. Here, we tested for the first time the performance of plastid and nuclear regions as DNA barcodes in Passiflora. This genus is a largely variable, with more than 900 species of high ecological, commercial, and ornamental importance. We analyzed 1034 accessions of 222 species representing the four subgenera of Passiflora and evaluated the effectiveness of five plastid regions and three nuclear datasets currently employed as DNA barcodes in plants using barcoding gap, applied similarity-, and tree-based methods. The plastid regions were able to identify less than 45% of species, whereas the nuclear datasets were efficient for more than 50% using “best match” and “best close match” methods of TaxonDNA software. All subgenera presented higher interspecific pairwise distances and did not fully overlap with the intraspecific distance, and similarity-based methods showed better results than tree-based methods. The nuclear ribosomal internal transcribed spacer 1 (ITS1) region presented a higher discrimination power than the other datasets and also showed other desirable characteristics as a DNA barcode for this genus. Therefore, we suggest that this region should be used as a starting point to identify Passiflora species.

Phylogeny of Salix subgenus Salix s.l. (Salicaceae): delimitation, biogeography, and reticulate evolution

Background

The taxonomy and systematics of Salix subgenus Salix s.l. is difficult. The reliability and evolutionary implications of two important morphological characters (number of stamens, and morphology of bud scales) used in subgeneric classification within Salix remain untested, and a disjunct Old–New World distribution pattern of a main clade of subgenus Salix s.l., revealed by a previous study, lacks a reasonable explanation. To study these questions, we conducted phylogenetic analyses based on 4,688 bp of sequence data from four plastid (rbcL, trnD–T, matK, and atpB–rbcL) and two nuclear markers (ETS and ITS) covering all subgenera of Salix, and all sections of subgenus Salix s.l.

Results

Subgenus Salix came out as para- or polyphyletic in both nrDNA and plastid trees. The plastid phylogeny successfully resolved relationships among the major clades of Salix, but resolution within subgenus Salix s.l. remained low. Nevertheless, three monophyletic groups were identifiable in subgenus Salix s.l.: the ‘main clade’ of subgenus Salix s.l., with New and Old World species being reciprocally monophyletic; the section Triandroides clade; and the subgenus Pleuradenia clade. While nrDNA regions showed higher resolution within subgenus Salix s.l., they failed to resolve subgeneric relationships. Extensive, statistically significant gene-tree incongruence was detected across nrDNA–plastid as well as nrDNA ETS–ITS phylogenies, suggesting reticulate evolution or hybridization within the group. The results were supported by network analyses. Ancestral-state reconstructions indicated that multiple stamens and free bud scales represent the plesiomorphic states within Salix, and that several significant shifts in stamen number and bud scale morphology have occurred.

Conclusions

Subgenus Salix s.l. is not monophyletic, and the evolutionary history of the subgenus has involved multiple reticulation events that may mainly be due to hybridization. The delimitation of subgenus Salix s.l. should be redefined by excluding section Triandrae and subgenus Pleuradenia from it. The evolutionary lability of bud-scale morphology and stamen number means that these characters are unreliable bases for classification. The disjunct Old–New World distribution of subgenus Salix s.l. appears to be linked to the profound climatic cooling during the Tertiary, which cut off gene exchange between New and Old World lineages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0311-7) contains supplementary material, which is available to authorized users.

A comparative analysis of genetic differentiation across six shared willow host species in leaf- and bud-galling sawflies

Genetic divergence and speciation in plant-feeding insects could be driven by contrasting selection pressures imposed by different plant species and taxa. While numerous examples of host-associated differentiation (HAD) have been found, the overall importance of HAD in insect diversification remains unclear, as few studies have investigated its frequency in relation to all speciation events. One promising way to infer the prevalence and repeatability of HAD is to estimate genetic differentiation in multiple insect taxa that use the same set of hosts. To this end, we measured and compared variation in mitochondrial COI and nuclear ITS2 sequences in population samples of leaf-galling Pontania and bud-galling Euura sawflies (Hymenoptera: Tenthredinidae) collected from six Salix species in two replicate locations in northern Fennoscandia. We found evidence of frequent HAD in both species complexes, as individuals from the same willow species tended to cluster together on both mitochondrial and nuclear phylogenetic trees. Although few fixed differences among the putative species were found, hierarchical AMOVAs showed that most of the genetic variation in the samples was explained by host species rather than by sampling location. Nevertheless, the levels of HAD measured across specific pairs of host species were not correlated in the two focal galler groups. Hence, our results support the hypothesis of HAD as a central force in herbivore speciation, but also indicate that evolutionary trajectories are only weakly repeatable even in temporally overlapping radiations of related insect taxa.

Testing evolutionary hypotheses for DNA barcoding failure in willows

The goal of DNA barcoding is to enable the rapid identification of taxa from short diagnostic DNA sequence profiles. But how feasible is this objective when many evolutionary processes, such as hybridization and selective sweeps, cause alleles to be shared among related taxa? In this issue of Molecular Ecology, Percy et al. (2014) test the full suite of seven candidate plant barcoding loci in a broad geographic sample of willow species. They show exceptional plastid haplotype sharing between species across continents, with most taxa not possessing a unique barcode sequence. Using population genetic and molecular dating analyses, they implicate hybridization and selective sweeps, but not incomplete lineage sorting, as the historical processes causing widespread haplotype sharing among willow taxa. This study represents an exceptional case of how poorly barcoding can perform, and highlights methodological issues using universal organellar regions for species identification.