skip to main content


Title: Geography shapes the phylogeny of frailejones (Espeletiinae Cuatrec., Asteraceae): a remarkable example of recent rapid radiation in sky islands
Background

The páramo ecosystem, located above the timberline in the tropical Andes, has been the setting for some of the most dramatic plant radiations, and it is one of the world’s fastest evolving and most diverse high-altitude ecosystems. Today 144+ species of frailejones (subtribe Espeletiinae Cuatrec., Asteraceae) dominate the páramo. Frailejones have intrigued naturalists and botanists, not just for their appealing beauty and impressive morphological diversity, but also for their remarkable adaptations to the extremely harsh environmental conditions of the páramo. Previous attempts to reconstruct the evolutionary history of this group failed to resolve relationships among genera and species, and there is no agreement regarding the classification of the group. Thus, our goal was to reconstruct the phylogeny of the frailejones and to test the influence of the geography on it as a first step to understanding the patterns of radiation of these plants.

Methods

Field expeditions in 70 páramos of Colombia and Venezuela resulted in 555 collected samples from 110 species. Additional material was obtained from herbarium specimens. Sequence data included nrDNA (ITS and ETS) and cpDNA (rpl16), for an aligned total of 2,954 bp. Fragment analysis was performed with AFLP data using 28 primer combinations and yielding 1,665 fragments. Phylogenies based on sequence data were reconstructed under maximum parsimony, maximum likelihood and Bayesian inference. The AFLP dataset employed minimum evolution analyses. A Monte Carlo permutation test was used to infer the influence of the geography on the phylogeny.

Results

Phylogenies reconstructed suggest that most genera are paraphyletic, but the phylogenetic signal may be misled by hybridization and incomplete lineage sorting. A tree with all the available molecular data shows two large clades: one of primarily Venezuelan species that includes a few neighboring Colombian species; and a second clade of only Colombian species. Results from the Monte Carlo permutation test suggests a very strong influence of the geography on the phylogenetic relationships. Venezuelan páramos tend to hold taxa that are more distantly-related to each other than Colombian páramos, where taxa are more closely-related to each other.

Conclusions

Our data suggest the presence of two independent radiations: one in Venezuela and the other in Colombia. In addition, the current generic classification will need to be deeply revised. Analyses show a strong geographic structure in the phylogeny, with large clades grouped in hotspots of diversity at a regional scale, and in páramo localities at a local scale. Differences in the degrees of relatedness between sympatric species of Venezuelan and Colombian páramos may be explained because of the younger age of the latter páramos, and the lesser time for speciation of Espeletiinae in them.

 
more » « less
NSF-PAR ID:
10022837
Author(s) / Creator(s):
 ;  
Publisher / Repository:
PeerJ
Date Published:
Journal Name:
PeerJ
Volume:
5
ISSN:
2167-8359
Page Range / eLocation ID:
Article No. e2968
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Summary

    The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny.

    We utilize fossil data and restriction‐site associatedDNAsequencing (RAD‐seq) for 632 individuals representing nearly 250Quercusspecies to infer a time‐calibrated phylogeny of the world's oaks. We use a reversible‐jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among‐clade sampling biases. We then map the > 20 000RAD‐seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny.

    Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity.

    The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

     
    more » « less
  2. Abstract

    To examine phylogenetic heterogeneity in turtle evolution, we collected thousands of high-confidence single-copy orthologs from 19 genome assemblies representative of extant turtle diversity and estimated a phylogeny with multispecies coalescent and concatenated partitioned methods. We also collected next-generation sequences from 26 turtle species and assembled millions of biallelic markers to reconstruct phylogenies based on annotated regions from the western painted turtle (Chrysemys picta bellii) genome (coding regions, introns, untranslated regions, intergenic, and others). We then measured gene tree-species tree discordance, as well as gene and site heterogeneity at each node in the inferred trees, and tested for temporal patterns in phylogenomic conflict across turtle evolution. We found strong and consistent support for all bifurcations in the inferred turtle species phylogenies. However, a number of genes, sites, and genomic features supported alternate relationships between turtle taxa. Our results suggest that gene tree-species tree discordance in these data sets is likely driven by population-level processes such as incomplete lineage sorting. We found very little effect of substitutional saturation on species tree topologies, and no clear phylogenetic patterns in codon usage bias and compositional heterogeneity. There was no correlation between gene and site concordance, node age, and DNA substitution rate across most annotated genomic regions. Our study demonstrates that heterogeneity is to be expected even in well-resolved clades such as turtles, and that future phylogenomic studies should aim to sample as much of the genome as possible in order to obtain accurate phylogenies for assessing conservation priorities in turtles. [Discordance; genomes; phylogeny; turtles.]

     
    more » « less
  3. Abstract Background

    Core genome phylogenies are widely used to build the evolutionary history of individual prokaryote species. By using hundreds or thousands of shared genes, these approaches are the gold standard to reconstruct the relationships of large sets of strains. However, there is growing evidence that bacterial strains exchange DNA through homologous recombination at rates that vary widely across prokaryote species, indicating that core genome phylogenies might not be able to reconstruct true phylogenies when recombination rate is high. Few attempts have been made to evaluate the robustness of core genome phylogenies to recombination, but some analyses suggest that reconstructed trees are not always accurate.

    Results

    In this study, we tested the robustness of core genome phylogenies to various levels of recombination rates. By analyzing simulated and empirical data, we observed that core genome phylogenies are relatively robust to recombination rates; nevertheless, our results suggest that many reconstructed trees are not completely accurate even when bootstrap supports are high. We found that some core genome phylogenies are highly robust to recombination whereas others are strongly impacted by it, and we identified that the robustness of core genome phylogenies to recombination is highly linked to the levels of selective pressures acting on a species. Stronger selective pressures lead to less accurate tree reconstructions, presumably because selective pressures more strongly bias the routes of DNA transfers, thereby causing phylogenetic artifacts.

    Conclusions

    Overall, these results have important implications for the application of core genome phylogenies in prokaryotes.

     
    more » « less
  4. Abstract

    Rapid species radiations present difficulties for phylogenetic reconstruction due to lack of phylogenetic information and processes such as deep coalescence/incomplete lineage sorting and hybridization. Phylogenomic data can overcome some of these difficulties. In this study, we use anchored hybrid enrichment (AHE) nuclear phylogenomic data and mitochondrial genomes recovered from AHE bycatch with several concatenated and coalescent approaches to reconstruct the poorly resolved radiation of the New Zealand cicada species in the generaKikihiaDugdale andMaoricicadaDugdale. Compared with previous studies using only three to five Sanger‐sequenced genes, we find increased resolution across our phylogenies, but several branches remain unresolved due to topological conflict among genes. Some nodes that are strongly supported by traditional support measures like bootstraps and posterior probabilities still show significant gene and site concordance conflict. In addition, we find strong mito‐nuclear discordance; likely the result of interspecific hybridization events in the evolutionary history ofKikihiaandMaoricicada.

     
    more » « less
  5. Abstract Background and Aims

    Pooideae grasses contain some of the world’s most important crop and forage species. Although much work has been conducted on understanding the genetic basis of trait diversification within a few annual Pooideae, comparative studies at the subfamily level are limited by a lack of perennial models outside ‘core’ Pooideae. We argue for development of the perennial non-core genus Melica as an additional model for Pooideae, and provide foundational data regarding the group’s biogeography and history of character evolution.

    Methods

    Supplementing available ITS and ndhF sequence data, we built a preliminary Bayesian-based Melica phylogeny, and used it to understand how the genus has diversified in relation to geography, climate and trait variation surveyed from various floras. We also determine biomass accumulation under controlled conditions for Melica species collected across different latitudes and compare inflorescence development across two taxa for which whole genome data are forthcoming.

    Key Results

    Our phylogenetic analyses reveal three strongly supported geographically structured Melica clades that are distinct from previously hypothesized subtribes. Despite less geographical affinity between clades, the two sister ‘Ciliata’ and ‘Imperfecta’ clades segregate from the more phylogenetically distant ‘Nutans’ clade in thermal climate variables and precipitation seasonality, with the ‘Imperfecta’ clade showing the highest levels of trait variation. Growth rates across Melica are positively correlated with latitude of origin. Variation in inflorescence morphology appears to be explained largely through differences in secondary branch distance, phyllotaxy and number of spikelets per secondary branch.

    Conclusions

    The data presented here and in previous studies suggest that Melica possesses many of the necessary features to be developed as an additional model for Pooideae grasses, including a relatively fast generation time, perenniality, and interesting variation in physiology and morphology. The next step will be to generate a genome-based phylogeny and transformation tools for functional analyses.

     
    more » « less