A global blueberry phylogeny: Evolution, diversification, and biogeography of Vaccinieae (Ericaceae).

Vaccinieae is a morphologically diverse and species-rich (∼1430 species) tribe in Ericaceae. Although the majority of diversity is tropical, Vaccinieae are best known for temperate crops (i.e., blueberries, cranberries, and lingonberries) in Vaccinium. Vaccinium itself (∼500 species) has been previously suggested as highly polyphyletic and taxonomic boundaries among many of the other genera in the tribe remain uncertain. We assessed the evolutionary history of Vaccinieae with phylogenomic analyses based on a target-enrichment dataset containing 256 low-copy nuclear loci and 210 species representing 30 of the 35 genera in the tribe and 25 of the 29 sections of Vaccinium. We conducted time-calibrated biogeographic analyses and diversification analyses to explore the area of origin and global dispersal history of the tribe. The analysis recovered a temperate North American origin for Vaccinieae approximately 30 million years ago. Tropical diversity of Vaccinieae was inferred to result from multiple, independent movements into the tropics from north-temperate ancestors. Diversification rate increases corresponded to radiation into the Andes and SE Asia. The pseudo-10-locular ovary evolved once in the tribe from the five-locular state, coinciding with the diversification of a major clade that includes most Asian Vaccinium and the group from which commercial blueberries are derived (V. sect. Cyanococcus). A reconstruction from available chromosome counts suggests that a major polyploid event predated the evolution of nearly half the diversity of Vaccinieae. The extent of polyphyly in Vaccinium documented here supports the need for a generic reclassification of the tribe.

The phylogeny and global biogeography of Primulaceae based on high-throughput DNA sequence data.

The angiosperm family Primulaceae is morphologically diverse and distributed nearly worldwide. However, phylogenetic uncertainty has obstructed the identification of major morphological and biogeographic transitions within the clade. We used target capture sequencing with the Angiosperms353 probes, taxon-sampling encompassing nearly all genera of the family, tree-based sequence curation, and multiple phylogenetic approaches to investigate the major clades of Primulaceae and their relationship to other Ericales. We generated dated phylogenetic trees and conducted broad-scale biogeographic analyses as well as stochastic character mapping of growth habit. We show that Ardisia, a pantropical genus and the largest in the family, is not monophyletic, with at least 19 smaller genera nested within it. Neotropical members of Ardisia and several smaller genera form a clade, an ancestor of which arrived in the Neotropics and began diversifying about 20 Ma. This Neotropical clade is most closely related to Elingamita and Tapeinosperma, which are most diverse on islands of the Pacific. Both Androsace and Primula are non-monophyletic by the inclusion of smaller genera. Ancestral state reconstructions revealed that there have either been parallel transitions to an herbaceous habit in Primuloideae, Samolus, and at least three lineages of Myrsinoideae, or a common ancestor of nearly all Primulaceae was herbaceous. Our results provide a robust estimate of phylogenetic relationships across Primulaceae and show that a revised classification of Myrsinoideae and several other clades within the family is necessary to render all genera monophyletic.

Genomic evidence of survival near ice sheet margins for some, but not all, North American trees.

Temperate species experienced dramatic range reductions during the Last Glacial Maximum, yet refugial populations from which modern populations are descended have never been precisely located. Climate-based models identify only broad areas of potential habitat, traditional phylogeographic studies provide poor spatial resolution, and pollen records for temperate forest communities are difficult to interpret and do not provide species-level taxonomic resolution. Here we harness signals of range expansion from large genomic datasets, using a simulation-based framework to infer the precise latitude and longitude of glacial refugia in two widespread, codistributed hickories (Carya spp.) and to quantify uncertainty in these estimates. We show that one species likely expanded from close to ice sheet margins near the site of a previously described macrofossil for the genus, highlighting support for the controversial notion of northern microrefugia. In contrast, the expansion origin inferred for the second species is compatible with classic hypotheses of distant displacement into southern refugia. Our statistically rigorous, powerful approach demonstrates how refugia can be located from genomic data with high precision and accuracy, addressing fundamental questions about long-term responses to changing climates and providing statistical insight into longstanding questions that have previously been addressed primarily qualitatively.

Admixture may be extensive among hyperdominant Amazon rainforest tree species.

Admixture is a mechanism by which species of long-lived plants may acquire novel alleles. However, the potential role of admixture in the origin and maintenance of tropical plant diversity is unclear. We ask whether admixture occurs in an ecologically important clade of Eschweilera (Parvifolia clade, Lecythidaceae), which includes some of the most widespread and abundant tree species in Amazonian forests. Using target capture sequencing, we conducted a detailed phylogenomic investigation of 33 species in the Parvifolia clade and investigated specific hypotheses of admixture within a robust phylogenetic framework. We found strong evidence of admixture among three ecologically dominant species, E. coriacea, E. wachenheimii and E. parviflora, but a lack of evidence for admixture among other lineages. Accepted species were largely distinguishable from one another, as was geographic structure within species. We show that hybridization may play a role in the evolution of the most widespread and ecologically variable Amazonian tree species. While admixture occurs among some species of Eschweilera, it has not led to widespread erosion of most species' genetic or morphological identities. Therefore, current morphological based species circumscriptions appear to provide a useful characterization of the clade's lineage diversity.

Tangled banks: A landscape genomic evaluation of Wallace's Riverine barrier hypothesis for three Amazon plant species.

Wallace's Riverine Barrier hypothesis is one of the earliest biogeographic explanations for Amazon speciation, but it has rarely been tested in plants. In this study, we used three woody Amazonian plant species to evaluate Wallace's Hypothesis using tools of landscape genomics. We generated unlinked single-nucleotide polymorphism (SNP) data from the nuclear genomes of 234 individuals (78 for each plant species) across 13 sampling sites along the Rio Branco, Brazil, for Amphirrhox longifolia (8,075 SNPs), Psychotria lupulina (9,501 SNPs) and Passiflora spinosa (14,536 SNPs). Although significantly different migration rates were estimated between species, the population structure data do not support the hypothesis that the Rio Branco-an allopatric barrier for primates and birds-is a significant genetic barrier for Amphirrhox longifolia, Passiflora spinosa or Psychotria lupulina. Overall, we demonstrated that medium-sized rivers in the Amazon Basin, such as the Rio Branco, are permeable barriers to gene flow for animal-dispersed and animal-pollinated plant species.

Target sequence capture in the Brazil nut family (Lecythidaceae): Marker selection and in silico capture from genome skimming data.

Reconstructing species trees from multi-loci datasets is becoming a standard practice in phylogenetics. Nevertheless, access to high-throughput sequencing may be costly, especially with studies of many samples. The potential high cost makes a priori assessments desirable in order to make informed decisions about sequencing. We generated twelve transcriptomes for ten species of the Brazil nut family (Lecythidaceae), identified a set of putatively orthologous nuclear loci and evaluated, in silico, their phylogenetic utility using genome skimming data of 24 species. We designed the markers using MarkerMiner, and developed a script, GoldFinder, to efficiently sub-select the best makers for sequencing. We captured, in silico, all designed 354 nuclear loci and performed a maximum likelihood phylogenetic analysis on the concatenated sequence matrix. We also calculated individual gene trees with maximum likelihood and used them for a coalescent-based species tree inference. Both analyses resulted in almost identical topologies. However, our nuclear-loci phylogenies were strongly incongruent with a published plastome phylogeny, suggesting that plastome data alone is not sufficient for species tree estimation. Our results suggest that using hundreds of nuclear markers (i.e. 354) will significantly improve the Lecythidaceae species tree. The framework described here will be useful, generally, for developing markers for species tree inference.

Wide but not impermeable: Testing the riverine barrier hypothesis for an Amazonian plant species.

Wallace's riverine barrier hypothesis postulates that large rivers, such as the Amazon and its tributaries, reduce or prevent gene flow between populations on opposite banks, leading to allopatry and areas of species endemism occupying interfluvial regions. Several studies have shown that two major tributaries, Rio Branco and Rio Negro, are important barriers to gene flow for birds, amphibians and primates. No botanical studies have considered the potential role of the Rio Branco as a barrier, while a single botanical study has evaluated the Rio Negro as a barrier. We studied an Amazon shrub, Amphirrhox longifolia (A. St.-Hil.) Spreng (Violaceae), as a model to test the riverine barrier hypothesis. Twenty-six populations of A. longifolia were sampled on both banks of the Rio Branco and Rio Negro in the core Amazon Basin. Double-digest RADseq was used to identify 8,010 unlinked SNP markers from the nuclear genome of 156 individuals. Data relating to population structure support the hypothesis that the Rio Negro acted as a significant genetic barrier for A. longifolia. On the other hand, no genetic differentiation was detected among populations spanning the narrower Rio Branco, which is a tributary of the Rio Negro. This study shows that the strength of riverine barriers for Amazon plants is dependent on the width of the river separating populations and species-specific dispersal traits. Future studies of plants with contrasting life history traits will further improve our understanding of the landscape genetics and allopatric speciation history of Amazon plant diversity.

Native bees mediate long-distance pollen dispersal in a shade coffee landscape mosaic.

Coffee farms are often embedded within a mosaic of agriculture and forest fragments in the world's most biologically diverse tropical regions. Although shade coffee farms can potentially support native pollinator communities, the degree to which these pollinators facilitate gene flow for native trees is unknown. We examined the role of native bees as vectors of gene flow for a reproductively specialized native tree, Miconia affinis, in a shade coffee and remnant forest landscape mosaic. We demonstrate extensive cross-habitat gene flow by native bees, with pollination events spanning more than 1,800 m. Pollen was carried twice as far within shade coffee habitat as in nearby forest, and trees growing within shade coffee farms received pollen from a far greater number of sires than trees within remnant forest. The study shows that shade coffee habitats support specialized native pollinators that enhance the fecundity and genetic diversity of remnant native trees.

Anonymous and EST-based microsatellite DNA markers that transfer broadly across the fig tree genus (Ficus, Moraceae).

PREMISE OF THE STUDY: We developed a set of microsatellite markers for broad utility across the species-rich pantropical tree genus Ficus (fig trees). The markers were developed to study population structure, hybridization, and gene flow in neotropical species. METHODS AND RESULTS: We developed seven novel primer sets from expressed sequence tag (EST) libraries of F. citrifolia and F. popenoei (subgen. Urostigma sect. Americana) and optimized five previously developed anonymous loci for cross-species amplification. The markers were successfully tested on four species from the basal subgenus Pharmacosycea sect. Pharmacosycea (F. insipida, F. maxima, F. tonduzii, and F. yoponensis) and seven species of the derived subgenus Urostigma (F. citrifolia, F. colubrinae, F. costaricana, F. nymphaeifolia, F. obtusifolia, F. pertusa, and F. popenoei). The 12 markers amplified consistently and displayed polymorphism in all the species. CONCLUSIONS: This set of microsatellite markers is transferable across the phylogenetic breadth of Ficus, and should therefore be useful for studies of population structure and gene flow in approximately 750 fig species worldwide.

By Animal, Water, or Wind: Can Dispersal Mode Predict Genetic Connectivity in Riverine Plant Species?

Seed dispersal is crucial to gene flow among plant populations. Although the effects of geographic distance and barriers to gene flow are well studied in many systems, it is unclear how seed dispersal mediates gene flow in conjunction with interacting effects of geographic distance and barriers. To test whether distinct seed dispersal modes (i.e., hydrochory, anemochory, and zoochory) have a consistent effect on the level of genetic connectivity (i.e., gene flow) among populations of riverine plant species, we used unlinked single-nucleotide polymorphisms (SNPs) for eight co-distributed plant species sampled across the Rio Branco, a putative biogeographic barrier in the Amazon basin. We found that animal-dispersed plant species exhibited higher levels of genetic diversity and lack of inbreeding as a result of the stronger genetic connectivity than plant species whose seeds are dispersed by water or wind. Interestingly, our results also indicated that the Rio Branco facilitates gene dispersal for all plant species analyzed, irrespective of their mode of dispersal. Even at a small spatial scale, our findings suggest that ecology rather than geography play a key role in shaping the evolutionary history of plants in the Amazon basin. These results may help improve conservation and management policies in Amazonian riparian forests, where degradation and deforestation rates are high.

Pervasive alteration of tree communities in undisturbed Amazonian forests.

Amazonian rainforests are some of the most species-rich tree communities on earth. Here we show that, over the past two decades, forests in a central Amazonian landscape have experienced highly nonrandom changes in dynamics and composition. Our analyses are based on a network of 18 permanent plots unaffected by any detectable disturbance. Within these plots, rates of tree mortality, recruitment and growth have increased over time. Of 115 relatively abundant tree genera, 27 changed significantly in population density or basal area--a value nearly 14 times greater than that expected by chance. An independent, eight-year study in nearby forests corroborates these shifts in composition. Contrary to recent predictions, we observed no increase in pioneer trees. However, genera of faster-growing trees, including many canopy and emergent species, are increasing in dominance or density, whereas genera of slower-growing trees, including many subcanopy species, are declining. Rising atmospheric CO2 concentrations may explain these changes, although the effects of this and other large-scale environmental alterations remain uncertain. These compositional changes could have important impacts on the carbon storage, dynamics and biota of Amazonian forests.

Vulnerability and resilience of tropical forest species to land-use change.

We provide a cross-taxon and historical analysis of what makes tropical forest species vulnerable to extinction. Several traits have been important for species survival in the recent and distant geological past, including seed dormancy and vegetative growth in plants, small body size in mammals, and vagility in insects. For major past catastrophes, such as the five mass extinction events, large range size and vagility or dispersal were key to species survival. Traits that make some species more vulnerable to extinction are consistent across time scales. Terrestrial organisms, particularly animals, are more extinction prone than marine organisms. Plants that persist through dramatic changes often reproduce vegetatively and possess mechanisms of die back. Synergistic interactions between current anthropogenic threats, such as logging, fire, hunting, pests and diseases, and climate change are frequent. Rising temperatures threaten all organisms, perhaps particularly tropical organisms adapted to small temperature ranges and isolated by distance from suitable future climates. Mutualist species and trophic specialists may also be more threatened because of such range-shift gaps. Phylogenetically specialized groups may be collectively more prone to extinction than generalists. Characterization of tropical forest species' vulnerability to anthropogenic change is constrained by complex interactions among threats and by both taxonomic and ecological impediments, including gross undersampling of biotas and poor understanding of the spatial patterns of taxa at all scales.

A Biogeographic Barrier Test Reveals a Strong Genetic Structure for a Canopy-Emergent Amazon Tree Species.

Wallace's (1854) Riverine Barrier hypothesis is one of the earliest explanations for Amazon biotic diversification. Despite the importance of this hypothesis for explaining speciation in some animal groups, it has not been studied extensively for plant species. In this study we use a prominent Amazon tree, Buchenavia oxycarpa (Mart.) Eichler (Combretaceae), to evaluate Wallace's hypothesis along the Rio Negro, a major Amazon tributary that has driven allopatric speciation for several animal taxa. We sampled six individuals from sixteen localities along both river banks, and used a modified ddRADseq protocol to identify SNP markers. Our population genomic data revealed strong genetic structure for B. oxycarpa sampled across banks of the Rio Negro (ϕCT = 0.576, P < 0.001), supporting the hypothesis that the Rio Negro acted as a significant genetic barrier for B. oxycarpa. Our study shows that gene flow for this large and well-dispersed Amazon tree is impeded by riverine barriers, though this has not yet resulted in speciation. Future studies focused on species with different life histories, including species restricted to non-flooded forests, are needed to further advance our understanding of Amazon rivers as drivers of biotic diversification.

The complex biogeographic history of a widespread tropical tree species.

Many tropical forest tree species have broad geographic ranges, and fossil records indicate that population disjunctions in some species were established millions of years ago. Here we relate biogeographic history to patterns of population differentiation, mutational and demographic processes in the widespread rainforest tree Symphonia globulifera using ribosomal (ITS) and chloroplast DNA sequences and nuclear microsatellite (nSSR) loci. Fossil records document sweepstakes dispersal origins of Neotropical S. globulifera populations from Africa during the Miocene. Despite historical long-distance gene flow, nSSR differentiation across 13 populations from Costa Rica, Panama, Ecuador (east and west of Andes) and French Guiana was pronounced (F(ST)= 0.14, R(ST)= 0.39, P < 0.001) and allele-size mutations contributed significantly (R(ST) > F(ST)) to the divergences between cis- and trans-Andean populations. Both DNA sequence and nSSR data reflect contrasting demographic histories in lower Mesoamerica and Amazonia. Amazon populations show weak phylogeographic structure and deviation from drift-mutation equilibrium indicating recent population expansion. In Mesoamerica, genetic drift was strong and contributed to marked differentiation among populations. The genetic structure of S. globulifera contains fingerprints of drift-dispersal processes and phylogeographic footprints of geological uplifts and sweepstakes dispersal.

New interpretations of fine-scale spatial genetic structure.

Recent methodological advances permit refined inferences of evolutionary processes from the fine-scale spatial genetic structure of plant populations. In this issue of Molecular Ecology, Born et al. (2008) exploit the full power of these methods by examining effects of ancient and recent landscape histories in an African rainforest tree species. The authors first detected admixture of distinct gene pools that may have formed in Pleistocene forest refuges. Then, comparing across six study populations in Gabon, the authors found similar patterns of fine-scale spatial genetic structure despite natural and anthropogenic variation in population density. The latter results suggest that enhanced gene dispersal may compensate for low population densities in fragmented landscapes.

Complete plastome sequences from Bertholletia excelsa and 23 related species yield informative markers for Lecythidaceae.

PREMISE OF THE STUDY: The tropical tree family Lecythidaceae has enormous ecological and economic importance in the Amazon basin. Lecythidaceae species can be difficult to identify without molecular data, however, and phylogenetic relationships within and among the most diverse genera are poorly resolved. METHODS: To develop informative genetic markers for Lecythidaceae, we used genome skimming to de novo assemble the full plastome of the Brazil nut tree (Bertholletia excelsa) and 23 other Lecythidaceae species. Indices of nucleotide diversity and phylogenetic signal were used to identify regions suitable for genetic marker development. RESULTS: The B. excelsa plastome contained 160,472 bp and was arranged in a quadripartite structure. Using the 24 plastome alignments, we developed primers for 10 coding and non-coding DNA regions containing exceptional nucleotide diversity and phylogenetic signal. We also developed 19 chloroplast simple sequence repeats for population-level studies. DISCUSSION: The coding region ycf1 and the spacer rpl16-rps3 outperformed plastid DNA markers previously used for barcoding and phylogenetics. Used in a phylogenetic analysis, the matrix of 24 plastomes showed with 100% bootstrap support that Lecythis and Eschweilera are polyphyletic. The plastomes and primers presented in this study will facilitate a broad array of ecological and evolutionary studies in Lecythidaceae.

Scale and strength oak-mesophyte interactions in a transitional oak-hickory forest.

Forests in eastern North America are undergoing rapid compositional changes as they experience novel climate, disturbance, and pest conditions. One striking pattern is the replacement of canopy oaks (Quercus spp.) by mesic, fire-sensitive, shade-tolerant species like red maple (Acer rubrum). To gain insight into the successional patterns driving stand-level canopy oak replacement we ask two questions: (i) What is the spatial association of oak and mesophyte recruitment compared to oak and mesophyte overstory individuals, and (ii) How do oaks and mesophytes differentially respond to canopy openings. We analyzed census data from a 23 ha forest plot surveyed in 2003, 2008 and 2014. We show that oak recruits are negatively associated with overstory red maples and black cherries (Prunus serotina), while mesophytic recruits were positively associated with overstory oaks. Second, we found that proximity to a dead overstory tree increased growth and survival for black cherries, increased growth for red maples, but had no effect on oaks. Black cherries and red maples are therefore better suited than oaks to take advantage of canopy openings and the moderate light available under adult oaks. These same fine scale competitive processes are contributing to canopy oak replacement across eastern North America.

Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species.

High-throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double-digest restriction-associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra- and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of FST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high-throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.