Inferring phylogeny despite incomplete lineage sorting

Cryptic speciation in a biodiversity hotspot: multilocus molecular data reveal new velvet worm species from Western Australia (Onychophora : Peripatopsidae : Kumbadjena)

There is a yet uncovered multitude of species to be found among Western Australian Onychophora. Kumbadjena, one of the two genera that reside in this region, has been previously suggested to house an extensive species complex. Morphology alone has not been able to elucidate the diversity in this genus and has instead muddled species delineations. Topologies and species delimitation analyses resulting from the sequences of two mitochondrial ribosomal markers (12S rRNA and 16S rRNA), one nuclear ribosomal marker (18S rRNA), and one mitochondrial protein-coding gene (cytochrome c oxidase subunit I) are indicative of several undescribed species. Fixed diagnostic nucleotide changes in the highly conserved sequences of 18S rRNA warrant distinction of three new species of Kumbadjena: K. toolbrunupensis, sp. nov., K. karricola, sp. nov., and K. extrema, sp. nov. The geographic distributions of the proposed species suggest that Kumbadjena is another example of short-range endemism, a common occurrence in the flora and fauna of the region. The extensive biodiversity and endemism in the region necessitates conservation to preserve the species and processes that promote speciation harboured by Western Australia.

Rooting phylogenetic trees under the coalescent model using site pattern probabilities

BACKGROUND

Phylogenetic tree inference is a fundamental tool to estimate ancestor-descendant relationships among different species. In phylogenetic studies, identification of the root - the most recent common ancestor of all sampled organisms - is essential for complete understanding of the evolutionary relationships. Rooted trees benefit most downstream application of phylogenies such as species classification or study of adaptation. Often, trees can be rooted by using outgroups, which are species that are known to be more distantly related to the sampled organisms than any other species in the phylogeny. However, outgroups are not always available in evolutionary research.

METHODS

In this study, we develop a new method for rooting species tree under the coalescent model, by developing a series of hypothesis tests for rooting quartet phylogenies using site pattern probabilities. The power of this method is examined by simulation studies and by application to an empirical North American rattlesnake data set.

RESULTS

The method shows high accuracy across the simulation conditions considered, and performs well for the rattlesnake data. Thus, it provides a computationally efficient way to accurately root species-level phylogenies that incorporates the coalescent process. The method is robust to variation in substitution model, but is sensitive to the assumption of a molecular clock.

CONCLUSIONS

Our study establishes a computationally practical method for rooting species trees that is more efficient than traditional methods. The method will benefit numerous evolutionary studies that require rooting a phylogenetic tree without having to specify outgroups.

Genomics clarifies taxonomic boundaries in a difficult species complex

Efforts to taxonomically delineate species are often confounded with conflicting information and subjective interpretation. Advances in genomic methods have resulted in a new approach to taxonomic identification that stands to greatly reduce much of this conflict. This approach is ideal for species complexes, where divergence times are recent (evolutionarily) and lineages less well defined. The California Roach/Hitch fish species complex is an excellent example, experiencing a convoluted geologic history, diverse habitats, conflicting species designations and potential admixture between species. Here we use this fish complex to illustrate how genomics can be used to better clarify and assign taxonomic categories. We performed restriction-site associated DNA (RAD) sequencing on 255 Roach and Hitch samples collected throughout California to discover and genotype thousands of single nucleotide polymorphism (SNPs). Data were then used in hierarchical principal component, admixture, and FST analyses to provide results that consistently resolved a number of ambiguities and provided novel insights across a range of taxonomic levels. At the highest level, our results show that the CA Roach/Hitch complex should be considered five species split into two genera (4 + 1) as opposed to two species from distinct genera (1 +1). Subsequent levels revealed multiple subspecies and distinct population segments within identified species. At the lowest level, our results indicate Roach from a large coastal river are not native but instead introduced from a nearby river. Overall, this study provides a clear demonstration of the power of genomic methods for informing taxonomy and serves as a model for future studies wishing to decipher difficult species questions. By allowing for systematic identification across multiple scales, taxonomic structure can then be tied to historical and contemporary ecological, geographic or anthropogenic factors.

The Kalanchoë genome provides insights into convergent evolution and building blocks of crassulacean acid metabolism

Crassulacean acid metabolism (CAM) is a water-use efficient adaptation of photosynthesis that has evolved independently many times in diverse lineages of flowering plants. We hypothesize that convergent evolution of protein sequence and temporal gene expression underpins the independent emergences of CAM from C3 photosynthesis. To test this hypothesis, we generate a de novo genome assembly and genome-wide transcript expression data for Kalanchoë fedtschenkoi, an obligate CAM species within the core eudicots with a relatively small genome (~260 Mb). Our comparative analyses identify signatures of convergence in protein sequence and re-scheduling of diel transcript expression of genes involved in nocturnal CO2 fixation, stomatal movement, heat tolerance, circadian clock, and carbohydrate metabolism in K. fedtschenkoi and other CAM species in comparison with non-CAM species. These findings provide new insights into molecular convergence and building blocks of CAM and will facilitate CAM-into-C3 photosynthesis engineering to enhance water-use efficiency in crops.

Taxonomic placement of Paphiopedilum rungsuriyanum (Cypripedioideae; Orchidaceae) based on morphological, cytological and molecular analyses

BACKGROUND

Paphiopedilum rungsuriyanum from Northern Laos was discovered and described in 2014. It is characterized by having miniature tessellated leaves, a flower having a helmet shaped lip with a V-shaped neckline, and a semi-lunate, 3-dentate staminode with an umbo. These morphological features distinguish P. rungsuriyanum from the other known sections/subgenera of Paphiopedilum, making it difficult to group with existing infrageneric units.

RESULTS

Paphiopedilum rungsuriyanum has chromosome number of 2n = 26. Fluorescence in situ hybridization study demonstrates that there are two 45S rDNA signals in the telomeric region of chromosomes, and more than 20 5S rDNA signals dispersed signals in the pericentromeric and centromeric regions. Phylogenetic analyses based on four nuclear (i.e. ITS, ACO, DEF4 and RAD51) and four plastid (i.e. atpI-atpH, matK, trnS-trnfM and ycf1) gene regions indicate that P. rungsuriyanum is nested in subgenus Paphiopedilum and is a sister to section Paphiopedilum.

CONCLUSIONS

The results in combination with karyomorphological, rDNA FISH patterns, morphological and phylogenetic analyses suggest a new section Laosianum to accommodate this species in the current sectional circumscription of subgenus Paphiopedilum.

Multilocus inference of species trees and DNA barcoding

The unprecedented amount of data resulting from next-generation sequencing has opened a new era in phylogenetic estimation. Although large datasets should, in theory, increase phylogenetic resolution, massive, multilocus datasets have uncovered a great deal of phylogenetic incongruence among different genomic regions, due both to stochastic error and to the action of different evolutionary process such as incomplete lineage sorting, gene duplication and loss and horizontal gene transfer. This incongruence violates one of the fundamental assumptions of the DNA barcoding approach, which assumes that gene history and species history are identical. In this review, we explain some of the most important challenges we will have to face to reconstruct the history of species, and the advantages and disadvantages of different strategies for the phylogenetic analysis of multilocus data. In particular, we describe the evolutionary events that can generate species tree—gene tree discordance, compare the most popular methods for species tree reconstruction, highlight the challenges we need to face when using them and discuss their potential utility in barcoding. Current barcoding methods sacrifice a great amount of statistical power by only considering one locus, and a transition to multilocus barcodes would not only improve current barcoding methods, but also facilitate an eventual transition to species-tree-based barcoding strategies, which could better accommodate scenarios where the barcode gap is too small or inexistent.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Species and hybrids in the genus Diaphanosoma Fischer, 1850 (Crustacea: Branchiopoda: Cladocera)

Cladocerans are well-studied planktonic crustaceans, especially those of the genus Daphnia in which interesting evolutionary questions have been addressed on speciation processes. The aim of the present study is to demonstrate that other genera of cladocerans show similar levels of cryptic diversity, intraspecific gene flow, and thus become useful model systems for comparison. In order to do so, we chose the genus Diaphanosoma, widespread in tropical and temperate areas. We started with a survey of species diversity in the genus Diaphanosoma in Asia using a morphological approach, then obtained sequences from a mitochondrial and a nuclear marker from multiple individuals of different species, performed tests on DNA taxonomy and molecular phylogenies, and assessed the role of hybridization in explaining the cases of mitonuclear discordance. The results are that cryptic diversity occurs in Diaphanosoma, and mitonuclear discordance was found in about 6% of the sequenced animals. Past hybridization is supported as the most likely explanation for the discordance: no evidence was found of first generation hybrids with heterozygous sequences. Our analysis on patterns of genetic diversity in Diaphanosoma supports similarities and differences with what is known in Daphnia.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Unconstrained Diameters for Deep Coalescence

The minimizing-deep-coalescence (MDC) approach infers a median (species) tree for a given set of gene trees under the deep coalescence cost. This cost accounts for the minimum number of deep coalescences needed to reconcile a gene tree with a species tree where the leaf-genes are mapped to the leaf-species through a function called leaf labeling. In order to better understand the MDC approach we investigate here the diameter of a gene tree, which is an important property of the deep coalescence cost. This diameter is the maximal deep coalescence costs for a given gene tree under all leaf labelings for each possible species tree topology. While we prove that this diameter is generally infinite, this result relies on the diameter's unrealistic assumption that species trees can be of infinite size. Providing a more practical definition, we introduce a natural extension of the gene tree diameter that constrains the species tree size by a given constant. For this new diameter, we describe an exact formula, present a complete classification of the trees yielding this diameter, derive formulas for its mean and variance, and demonstrate its ability using comparative studies.

IDXL: Species Tree Inference Using Internode Distance and Excess Gene Leaf Count

We propose an extension of the distance matrix methods NJst and ASTRID to infer species trees from incongruent gene trees having Incomplete Lineage Sorting. Both approaches consider the average internode distance (ID) between individual taxa pairs as the distance measure. The measure ID does not use the root of a tree, and thus may not always infer the relative position of a taxon with respect to the root. We define a novel distance measure excess gene leaf count (XL) between individual couplets. The XL measure is computed using the root of a tree. It is proved to be additive, and is shown to infer the relative order of divergence among individual couplets better. We propose a novel method IDXL which uses both the XL and ID measures for species tree construction. IDXL is shown to perform better than NJst and other distance matrix approaches for most of the biological and simulated datasets. Having the same computational complexity as NJst, IDXL can be applied for species tree inference on large-scale biological datasets.

Inferring rooted species trees from unrooted gene trees using approximate Bayesian computation

Methods for inferring species trees from gene trees motivated by incomplete lineage sorting typically use either rooted gene trees to infer a rooted species tree, or use unrooted gene trees to infer an unrooted species tree, which is then typically rooted using one or more outgroups. Theoretically, however, it has been known since 2011 that it is possible to consistently infer the root of the species tree directly from unrooted gene trees without assuming an outgroup. Here, we use approximate Bayesian computation to infer the root of the species tree from unrooted gene trees assuming the multispecies coalescent model. It is hoped that this approach will be useful in cases where an appropriate outgroup is difficult to find and gene trees do not follow a molecular clock. We use approximate Bayesian computation to infer the root of the species tree from unrooted gene trees. This approach could also be useful when there is prior information that makes a small number of root locations plausible in an unrooted species tree.

Isolation and molecular characterization of Acanthamoeba from patients with keratitis in Spain

In order to improve our knowledge on the epidemiology of amoebic keratitis, as well as the identification of Acanthamoeba isolates, we have isolated Acanthamoeba spp. from five symptomatic patients in Spain in the present study. All isolates were grown in axenic liquid medium, with only one exception. The morphology of these isolates were characterized by optical and scanning electron microscopy. Their structural features corresponded to those of amphizoic amoebae (namely Acanthamoeba spp.). The molecular characterization of the five Acanthamoeba isolates yielded six sequences. Almost complete 18S rRNA gene sequences (>2000bp) were obtained from three isolates and partial sequences (∼1500bp) from the other two. A robust phylogenetic analysis based on the almost complete 18S rRNA sequence showed that four isolates belonged to the T4 genotype and the other one to the T3 genotype. However, all isolates were identified as T4 genotype using the ASA.S1 fragment. As previously suggested by other researchers, only a robust phylogenetic approach may be helpful in identifying Acanthamoeba genotypes. In addition, new data on the phylogenetic relationships among the Acanthamoeba genotypes is provided and discussed.

Molecular markers provide evidence for a broad-fronted recolonisation of the widespread calcareous grassland species Sanguisorba minor from southern and cryptic northern refugia

Calcareous grasslands belong to the most species-rich and endangered habitats in Europe. However, little is known about the origin of the species typically occurring in these grasslands. In this study we analysed the glacial and post-glacial history of Sanguisorba minor, a typical plant species frequently occurring in calcareous grasslands. The study comprised 38 populations throughout the whole distribution range of the species across Europe. We used molecular markers (AFLP) and applied Bayesian cluster analysis as well as spatial principal components analysis (sPCA) to identify glacial refugia and post-glacial migration routes to Central Europe. Our study revealed significant differences in the level of genetic variation and the occurrence of rare fragments within populations of S. minor and a distinct separation of eastern and western lineages. The analyses uncovered traditional southern but also cryptic northern refugia and point towards a broad fronted post-glacial recolonisation. Based on these results we postulate that incomplete lineage sorting may have contributed to the detected pattern of genetic variation and that S. minor recolonised Central Europe post-glacially from Iberia and northern glacial refugia in France, Belgium or Germany. Our results highlight the importance of refugial areas for the conservation of intraspecific variation in calcareous grassland species.

The befuddling nature of mouse lemur hands and feet at Bezà Mahafaly, SW Madagascar

The reddish-gray mouse lemur (Microcebus griseorufus) possesses striking phenotypic and behavioral variation. This project investigates differences in autopod proportions in neighboring populations of M. griseorufus from the Special Reserve at Bezà Mahafaly in southwest Madagascar. One population resides in an environment generally preferred by M. griseorufus-a spiny forest with large-trunked trees, vertically-oriented supports, and more open ground, while the other resides in a gallery forest with abundant small, often horizontal peripheral branches in high canopy. We demonstrate significant interpopulation differences in autopod morphophology despite no evidence of divergence in mitochondrial cytochrome b. We test two hypotheses regarding ultimate causation. The first, based on the Fine Branch Arborealism Hypothesis (FBAH), holds that autopod differences are related to different locomotor practices in the two environments, and the second, based on the Narrow Niche Hypothesis (NNH), holds that the observed differences reflect a relaxation (from ancestral to descendant conditions) of selective pressure for terrestrial locomotion and/or use of large, vertical supports combined with positive selection for locomoting in peripheral branch settings. Our data conform well to FBAH expectations and show some support for the NNH. Individuals from the gallery forest possess disproportionally long posterior digits that facilitate locomotion on small, flexible canopy supports while individuals from the spiny forest possess shorter posterior digits and a longer pollex/hallux that increase functional grasping diameter for large vertical supports and facilitate efficient ground locomotion. Focal individual data confirm differences in how often individuals descend to the ground and use vertical supports. We further show that predispersal juveniles, like adults, possess autopod morphologies suited to their natal forest. We explore two proximate mechanisms that could generate these cheiridial differences. The first posits an in vivo plastic response to different locomotor behaviors, the second posits differences that manifest in early development.

Challenges in Species Tree Estimation Under the Multispecies Coalescent Model

The multispecies coalescent (MSC) model has emerged as a powerful framework for inferring species phylogenies while accounting for ancestral polymorphism and gene tree-species tree conflict. A number of methods have been developed in the past few years to estimate the species tree under the MSC. The full likelihood methods (including maximum likelihood and Bayesian inference) average over the unknown gene trees and accommodate their uncertainties properly but involve intensive computation. The approximate or summary coalescent methods are computationally fast and are applicable to genomic datasets with thousands of loci, but do not make an efficient use of information in the multilocus data. Most of them take the two-step approach of reconstructing the gene trees for multiple loci by phylogenetic methods and then treating the estimated gene trees as observed data, without accounting for their uncertainties appropriately. In this article we review the statistical nature of the species tree estimation problem under the MSC, and explore the conceptual issues and challenges of species tree estimation by focusing mainly on simple cases of three or four closely related species. We use mathematical analysis and computer simulation to demonstrate that large differences in statistical performance may exist between the two classes of methods. We illustrate that several counterintuitive behaviors may occur with the summary methods but they are due to inefficient use of information in the data by summary methods and vanish when the data are analyzed using full-likelihood methods. These include (i) unidentifiability of parameters in the model, (ii) inconsistency in the so-called anomaly zone, (iii) singularity on the likelihood surface, and (iv) deterioration of performance upon addition of more data. We discuss the challenges and strategies of species tree inference for distantly related species when the molecular clock is violated, and highlight the need for improving the computational efficiency and model realism of the likelihood methods as well as the statistical efficiency of the summary methods.

Epigean gammarids survived millions of years of severe climatic fluctuations in high latitude refugia throughout the Western Carpathians

Isolated glacial refugia have been documented in Central Europe for a number of taxa, but conclusive evidence for epigean aquatic species has remained elusive. Using molecular data (mitochondrial and nuclear markers), we compared the spatial patterns of lineage diversity of the widely distributed Gammarus fossarum species complex between two adjacent biogeographically and geomorphologically distinct Central European regions: the Bohemian Massif and the Western Carpathians. We investigated if the observed patterns of spatial diversity are more likely to stem from historical or present-day factors. Phylogenetic and phylogeographic analyses revealed eight phylogenetically diverse lineages: two exhibiting local signatures of recent demographic expansion inhabit both regions, while the other six display a relict distributional pattern and are found only in the Western Carpathians. Molecular dating indicates that these lineages are old and probably diverged throughout the Miocene (7-18Ma). Furthermore, their distribution does not seem to be constrained by the present boundaries of river catchments or topography. The contrasting spatial patterns of diversity observed between the two regions thus more likely result from historical rather than contemporaneous or recent factors. Our results indicate that despite the high latitude and proximity to the Pleistocene ice sheets, the Western Carpathians functioned as long-term glacial refugia for permanent freshwater fauna, allowing the uninterrupted survival of ancient lineages through millions of years of drastic climatic fluctuations.

Simultaneous speciation in the European high mountain flowering plant genus Facchinia (Minuartia s.l., Caryophyllaceae) revealed by genotyping-by-sequencing

Understanding the relative importance of different mechanisms of speciation in a given lineage requires fully resolved interspecific relationships. Using Facchinia, a genus of seven species centred in the European Alps, we explore whether the polytomy found by Sanger sequencing analyses of standard nuclear (ITS) and plastid markers (trnQ-rps16) is a hard or soft polytomy by substantially increasing the amount of DNA sequence data, generated by genotyping-by-sequencing. In comparison to 142 phylogenetically informative sites in the Sanger sequences the GBS sequences yielded 3363 phylogenetically informative sites after exclusion of apparently oversaturated SNPs. Maximum parsimony, maximum likelihood, NeighborNet, SVDquartets and Astral-II analyses all resulted in phylogenetic trees (and networks) in which interspecific relationships were largely unresolved. After excluding incomplete lineage sorting, hybridisation and oversaturation of characters as possible causes for lack of phylogenetic resolution, we conclude that the polytomy obtained most likely represents a hard polytomy. We hypothesize that diversification of Facchinia is best interpreted as the result of multiple simultaneous vicariance in response to climatic changes during the Early Quaternary.

A genomic assessment of species boundaries and hybridization in a group of highly polymorphic anoles (distichus species complex)

Delimiting young species is one of the great challenges of systematic biology, particularly when the species in question exhibit little morphological divergence. Anolis distichus, a trunk anole with more than a dozen subspecies that are defined primarily by dewlap color, may actually represent several independent evolutionary lineages. To test this, we utilized amplified fragment length polymorphisms (AFLP) genome scans and genetic clustering analyses in conjunction with a coalescent‐based species delimitation method. We examined a geographically widespread set of samples and two heavily sampled hybrid zones. We find that genetic divergence is associated with a major biogeographic barrier, the Hispaniolan paleo‐island boundary, but not with dewlap color. Additionally, we find support for hypotheses regarding colonization of two Hispaniolan satellite islands and the Bahamas from mainland Hispaniola. Our results show that A. distichus is composed of seven distinct evolutionary lineages still experiencing a limited degree of gene flow. We suggest that A. distichus merits taxonomic revision, but that dewlap color cannot be relied upon as the primary diagnostic character.

Deep reticulation and incomplete lineage sorting obscure the diploid phylogeny of rain-lilies and allies (Amaryllidaceae tribe Hippeastreae)

Hybridization is a frequent and important force in plant evolution. Next-generation sequencing (NGS) methods offer new possibilities for clade resolution and ambitious sampling of gene genealogies, yet difficulty remains in detecting deep reticulation events using currently available methods. We reconstructed the phylogeny of diploid representatives of Amaryllidaceae tribe Hippeastreae to test the hypothesis of ancient hybridizations preceding the radiation of its major subclade, Hippeastrinae. Through hybrid enrichment of DNA libraries and NGS, we obtained data for 18 nuclear loci through a curated assembly approach and nearly complete plastid genomes for 35 ingroup taxa plus 5 outgroups. Additionally, we obtained alignments for 39 loci through an automated assembly algorithm. These data were analyzed with diverse phylogenetic methods, including concatenation, coalescence-based species tree estimation, Bayesian concordance analysis, and network reconstructions, to provide insights into the evolutionary relationships of Hippeastreae. Causes for gene tree heterogeneity and cytonuclear discordance were examined through a Bayesian posterior predictive approach (JML) and coalescent simulations. Two major clades were found, Hippeastrinae and Traubiinae, as previously reported. Our results suggest the presence of two major nuclear lineages in Hippeastrinae characterized by different chromosome numbers: (1) Tocantinia and Hippeastrum with 2n=22, and (2) Eithea, Habranthus, Rhodophiala, and Zephyranthes mostly with 2n=12, 14, and 18. Strong cytonuclear discordance was confirmed in Hippeastrinae, and a network scenario with at least six hybridization events is proposed to reconcile nuclear and plastid signals, along a backbone that may also have been affected by incomplete lineage sorting at the base of each major subclade.

Evolutionary and domestication history of Cucurbita (pumpkin and squash) species inferred from 44 nuclear loci

Phylogenetics can facilitate the study of plant domestication by resolving sister relationships between crops and their wild relatives, thereby identifying the ancestors of cultivated plants. Previous phylogenetic studies of the six Cucurbita crop lineages (pumpkins and squashes) and their wild relatives suggest histories of deep coalescence that complicate uncovering the genetic origins of the six crop taxa. We investigated the evolution of wild and domesticated Cucurbita using the most comprehensive and robust molecular-based phylogeny for Cucurbita to date based on 44 loci derived from introns of single-copy nuclear genes. We discovered novel relationships among Cucurbita species and recovered the first Cucurbita tree with well-supported resolution within species. Cucurbita comprises a clade of mesophytic annual species that includes all six crop taxa and a grade of xerophytic perennial species that represent the ancestral xerophytic habit of the genus. Based on phylogenetic resolution within-species we hypothesize that the magnitude of domestication bottlenecks varies among Cucurbita crop lineages. Our phylogeny clarifies how wild Cucurbita species are related to the domesticated taxa. We find close relationships between two wild species and crop lineages not previously identified. Expanded geographic sampling of key wild species is needed for improved understanding of the evolution of domesticated Cucurbita.

Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content

A comprehensive and meaningful phylogenetic hypothesis for the commercially important coffee genus (Coffea) has long been a key objective for coffee researchers. For molecular studies, progress has been limited by low levels of sequence divergence, leading to insufficient topological resolution and statistical support in phylogenetic trees, particularly for the major lineages and for the numerous species occurring in Madagascar. We report here the first almost fully resolved, broadly sampled phylogenetic hypothesis for coffee, the result of combining genotyping-by-sequencing (GBS) technology with a newly developed, lab-based workflow to integrate short read next-generation sequencing for low numbers of additional samples. Biogeographic patterns indicate either Africa or Asia (or possibly the Arabian Peninsula) as the most likely ancestral locality for the origin of the coffee genus, with independent radiations across Africa, Asia, and the Western Indian Ocean Islands (including Madagascar and Mauritius). The evolution of caffeine, an important trait for commerce and society, was evaluated in light of our phylogeny. High and consistent caffeine content is found only in species from the equatorial, fully humid environments of West and Central Africa, possibly as an adaptive response to increased levels of pest predation. Moderate caffeine production, however, evolved at least one additional time recently (between 2 and 4Mya) in a Madagascan lineage, which suggests that either the biosynthetic pathway was already in place during the early evolutionary history of coffee, or that caffeine synthesis within the genus is subject to convergent evolution, as is also the case for caffeine synthesis in coffee versus tea and chocolate.

Genomewide SNP data reveal cryptic phylogeographic structure and microallopatric divergence in a rapids-adapted clade of cichlids from the Congo River

The lower Congo River is a freshwater biodiversity hot spot in Africa characterized by some of the world's largest rapids. However, little is known about the evolutionary forces shaping this diversity, which include numerous endemic fishes. We investigated phylogeographic relationships in Teleogramma, a small clade of rheophilic cichlids, in the context of regional geography and hydrology. Previous studies have been unable to resolve phylogenetic relationships within Teleogramma due to lack of variation in nuclear genes and discrete morphological characters among putative species. To sample more broadly across the genome, we analysed double-digest restriction-associated sequencing (ddRAD) data from 53 individuals across all described species in the genus. We also assessed body shape and mitochondrial variation within and between taxa. Phylogenetic analyses reveal previously unrecognized lineages and instances of microallopatric divergence across as little as ~1.5 km. Species ranges appear to correspond to geographic regions broadly separated by major hydrological and topographic barriers, indicating these features are likely important drivers of diversification. Mitonuclear discordance indicates one or more introgressive hybridization events, but no clear evidence of admixture is present in nuclear genomes, suggesting these events were likely ancient. A survey of female fin patterns hints that previously undetected lineage-specific patterning may be acting to reinforce species cohesion. These analyses highlight the importance of hydrological complexity in generating diversity in certain freshwater systems, as well as the utility of ddRAD-Seq data in understanding diversification processes operating both below and above the species level.

Phylogenetic Relationships of Japanese Auritibicen Species (Hemiptera: Cicadidae: Cryptotympanini) Inferred from Mitochondrial and Nuclear Gene Sequences

We investigated the phylogenetic relationships and divergence times within the genus Auritibicen(Cicadidae: Cicadinae: Cryptotympanini), analyzing five Japanese species (A. japonicus, A. bihamatus,A. kyushyuensis, A. esakii and A. flammatus) and three species from East Asian mainland and Taiwan (A. atrofasciatus, A. intermedius and A. chujoi) using mitochondrial cytochrome oxidase subunit I (COI) and nuclear elongation factor 1-alpha (EF-1a) gene sequences. Although the EF-1a gene tree did not resolve the relationships among these Auritibicen species, the trees based on COI gene and the combined data set showed that Japanese taxa comprised three distinct lineages: the individual species A. flammatus and A. bihamatus, and the A. japonicus group, comprising A. japonicus, A. esakii and A. kyushyuensis from Japan and A. intermedius from Korea. In A. kyushyuensis, which comprises three populations in Kyushu, western Honshu and Shikoku, the specimens from western Honshu and Shikoku were closely related to each other, but not to the specimen from Kyushu; instead, they were sister to the Korean A. intermedius. The incongruence between the gene tree and species tree necessitates further population genetic and morphological studies to confirm the classification and species status of the western Honshu and Shikoku populations of A. kyushyuensis, which were originally described as two independent species. Divergence time estimation suggested that the most recent common ancestor of Auritibicen species studied dated back to the late Pliocene and that the species of the A. japonicus group diverged during the mid Pleistocene. Thus, the Pleistocene climatic fluctuation may have promoted the divergence of the Auritibicen species.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.

Species Tree Estimation from Genome-Wide Data with guenomu

The history of particular genes and that of the species that carry them can be different for a variety of reasons. In particular, gene trees and species trees can differ due to well-known evolutionary processes such as gene duplication and loss, lateral gene transfer, or incomplete lineage sorting. Species tree reconstruction methods have been developed to take this incongruence into account; these can be divided grossly into supertree and supermatrix approaches. Here we introduce a new Bayesian hierarchical model that we have recently developed and implemented in the program guenomu. The new model considers multiple sources of gene tree/species tree disagreement. Guenomu takes as input posterior distributions of unrooted gene tree topologies for multiple gene families, in order to estimate the posterior distribution of rooted species tree topologies.

The D3-D5 region of large subunit ribosomal DNA provides good resolution of German limnic and terrestrial nematode communities

Reliable and well-developed DNA barcode databases are indispensable for the identification of microscopic life. However, effectiveness of molecular barcoding in identifying terrestrial specimens, and nematodes in particular, has received little attention. In this study, ca 600 ribosomal large subunit DNA fragments (D3-D5 region) were successfully amplified for 79 limnic and terrestrial nematode species sampled at 147 locations across Germany. Distinctive DNA motifs in the LSU region were identified in 80% of all species examined. For 13 supposedly single morphospecies, 2-7 LSU barcode groups were detected with a wide range of intraspecific variations (0.09-7.9%). This region seems to be more suitable for the assessment of limno-terrestrial nematode diversity than the frequently used mitochondrial gene COI, as amplification success of the latter fragment is low for several nematode species. Our reference database for nematodes may serve as a starting point for applied and fundamental studies for these ubiquitous, ecologically highly relevant, organisms.

Comprehensive Phylogenetic Analysis of Bovine Non-aureus Staphylococci Species Based on Whole-Genome Sequencing

Non-aureus staphylococci (NAS), a heterogeneous group of a large number of species and subspecies, are the most frequently isolated pathogens from intramammary infections in dairy cattle. Phylogenetic relationships among bovine NAS species are controversial and have mostly been determined based on single-gene trees. Herein, we analyzed phylogeny of bovine NAS species using whole-genome sequencing (WGS) of 441 distinct isolates. In addition, evolutionary relationships among bovine NAS were estimated from multilocus data of 16S rRNA, hsp60, rpoB, sodA, and tuf genes and sequences from these and numerous other single genes/proteins. All phylogenies were created with FastTree, Maximum-Likelihood, Maximum-Parsimony, and Neighbor-Joining methods. Regardless of methodology, WGS-trees clearly separated bovine NAS species into five monophyletic coherent clades. Furthermore, there were consistent interspecies relationships within clades in all WGS phylogenetic reconstructions. Except for the Maximum-Parsimony tree, multilocus data analysis similarly produced five clades. There were large variations in determining clades and interspecies relationships in single gene/protein trees, under different methods of tree constructions, highlighting limitations of using single genes for determining bovine NAS phylogeny. However, based on WGS data, we established a robust phylogeny of bovine NAS species, unaffected by method or model of evolutionary reconstructions. Therefore, it is now possible to determine associations between phylogeny and many biological traits, such as virulence, antimicrobial resistance, environmental niche, geographical distribution, and host specificity.

Evolutionary scalpels for dissecting tumor ecosystems

Amidst the growing literature on cancer genomics and intratumor heterogeneity, essential principles in evolutionary biology recur time and time again. Here we use these principles to guide the reader through major advances in cancer research, highlighting issues of "hit hard, hit early" treatment strategies, drug resistance, and metastasis. We distinguish between two frameworks for understanding heterogeneous tumors, both of which can inform treatment strategies: (1) The tumor as diverse ecosystem, a Darwinian population of sometimes-competing, sometimes-cooperating cells; (2) The tumor as tightly integrated, self-regulating organ, which may hijack developmental signals to restore functional heterogeneity after treatment. While the first framework dominates literature on cancer evolution, the second framework enjoys support as well. Throughout this review, we illustrate how mathematical models inform understanding of tumor progression and treatment outcomes. Connecting models to genomic data faces computational and technical hurdles, but high-throughput single-cell technologies show promise to clear these hurdles. This article is part of a Special Issue entitled: Evolutionary principles - heterogeneity in cancer?, edited by Dr. Robert A. Gatenby.

Evolutionary history of the sequestrate genus Rossbeevera (Boletaceae) reveals a new genus Turmalinea and highlights the utility of ITS minisatellite-like insertions for molecular identification

The sequestrate (truffle-like) basidiomycete genera Rossbeevera, Chamonixia, and Octaviania are closely related to the epigeous mushroom genera Leccinum and Leccinellum. In order to elucidate the properties and placement of several undescribed sequestrate taxa in the group and to reveal the evolutionary history of Rossbeevera and its allies, we conducted phylogenetic analyses based on three nuclear (ITS, nLSU, EF-1α) and two mitochondrial DNA loci (ATP6 and mtSSU) as well as precise morphological observations. Phylogenetic analyses of three nuclear loci suggest a complex evolutionary history with sequestrate fruiting bodies present in several clades, including a previously unrecognized sister clade to Rossbeevera. Here we propose a new sequestrate genus, Turmalinea, with four new species and one new subspecies as well as two new species of Rossbeevera. The three-locus nuclear phylogeny resolves species-level divergence within the Rossbeevera-Turmalinea lineage, whereas a separate phylogeny based on two mitochondrial genes corresponds to geographic distance within each species-level lineage and suggests incomplete lineage sorting (ILS) and gene introgression within several intraspecific lineages of Rossbeevera. Furthermore, topological incongruence among the three nuclear single-locus phylogenies suggests that ancient speciation within Rossbeevera probably involved considerable ILS. We also found an unusually long, minisatellite-like insertion within the ITS2 in all Rossbeevera and Turmalinea species. A barcode gap analysis demonstrates that the insertion is more informative for discrimination at various taxonomic levels than the rest of the ITS region and could therefore serve as a unique molecular barcode for these genera.

On the origins of Balkan endemics: the complex evolutionary history of the Cyanus napulifer group (Asteraceae)

BACKGROUND AND AIMS

The Balkan Peninsula is one of the most important centres of plant diversity in Europe. Here we aim to fill the gap in the current knowledge of the evolutionary processes and factors modelling this astonishing biological richness by applying multiple approaches to the Cyanus napulifer group.

METHODS

To reconstruct the mode of diversification within the C. napulifer group and to uncover its relationships with potential relatives with x = 10 from Europe and Northern Africa, we examined variation in genetic markers (amplified fragment length polymorphisms [AFLPs]; 460 individuals), relative DNA content (4',6-diamidino-2-phenylindole [DAPI] flow cytometry, 330 individuals) and morphology (multivariate morphometrics, 40 morphological characters, 710 individuals). To elucidate its evolutionary history, we analysed chloroplast DNA (cpDNA) sequences of the genus Cyanus deposited in the GenBank database.

KEY RESULTS

The AFLPs revealed a suite of closely related entities with variable levels of differentiation. The C. napulifer group formed a genetically well-defined unit. Samples outside the group formed strongly diversified and mostly species-specific genetic lineages with no further geographical patterns, often characterized also by a different DNA content. AFLP analysis of the C. napulifer group revealed extensive radiation and split it into nine allopatric (sub)lineages with varying degrees of congruence among genetic, DNA-content and morphological patterns. Genetic admixture was usually detected in contact zones between genetic lineages. Plastid data indicated extensive maintenance of ancestral variation across Cyanus perennials.

CONCLUSION

The C. napulifer group is an example of a rapidly and recently diversified plant group whose genetic lineages have evolved in spatio-temporal isolation on the topographically complex Balkan Peninsula. Adaptive radiation, accompanied in some cases by long-term isolation and hybridization, has contributed to the formation of this species complex and its mosaic pattern.

Genotyping-by-sequencing provides the discriminating power to investigate the subspecies of Daucus carota (Apiaceae)

BACKGROUND

The majority of the subspecies of Daucus carota have not yet been discriminated clearly by various molecular or morphological methods and hence their phylogeny and classification remains unresolved. Recent studies using 94 nuclear orthologs and morphological characters, and studies employing other molecular approaches were unable to distinguish clearly many of the subspecies. Fertile intercrosses among traditionally recognized subspecies are well documented. We here explore the utility of single nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing (GBS) to serve as an effective molecular method to discriminate the subspecies of the D. carota complex.

RESULTS

We used GBS to obtain SNPs covering all nine Daucus carota chromosomes from 162 accessions of Daucus and two related genera. To study Daucus phylogeny, we scored a total of 10,814 or 38,920 SNPs with a maximum of 10 or 30 % missing data, respectively. To investigate the subspecies of D. carota, we employed two data sets including 150 accessions: (i) rate of missing data 10 % with a total of 18,565 SNPs, and (ii) rate of missing data 30 %, totaling 43,713 SNPs. Consistent with prior results, the topology of both data sets separated species with 2n = 18 chromosome from all other species. Our results place all cultivated carrots (D. carota subsp. sativus) in a single clade. The wild members of D. carota from central Asia were on a clade with eastern members of subsp. sativus. The other subspecies of D. carota were in four clades associated with geographic groups: (1) the Balkan Peninsula and the Middle East, (2) North America and Europe, (3) North Africa exclusive of Morocco, and (4) the Iberian Peninsula and Morocco. Daucus carota subsp. maximus was discriminated, but neither it, nor subsp. gummifer (defined in a broad sense) are monophyletic.

CONCLUSIONS

Our study suggests that (1) the morphotypes identified as D. carota subspecies gummifer (as currently broadly circumscribed), all confined to areas near the Atlantic Ocean and the western Mediterranean Sea, have separate origins from sympatric members of other subspecies of D. carota, (2) D. carota subsp. maximus, on two clades with some accessions of subsp. carota, can be distinguished from each other but only with poor morphological support, (3) D. carota subsp. capillifolius, well distinguished morphologically, is an apospecies relative to North African populations of D. carota subsp. carota, (4) the eastern cultivated carrots have origins closer to wild carrots from central Asia than to western cultivated carrots, and (5) large SNP data sets are suitable for species-level phylogenetic studies in Daucus.