Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad.

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (PST-FST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. PST-FST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

Differential gene expression reveals mechanisms related to habitat divergence between hybridizing orchids from the Neotropical coastal plains.

BACKGROUND: Closely related hybridizing species are ideal systems for identifying genomic regions underlying adaptive divergence. Although gene expression plays a central role in determining ecologically-based phenotypic differences, few studies have inferred the role of gene expression for adaptive divergence in Neotropical systems. In this study, we conduct genome-wide expression analysis alongside soil elemental analysis in sympatric and allopatric populations of Epidendrum fulgens and E. puniceoluteum (Orchidaceae), which occur in contrasting adjacent habitats in the Neotropical coastal plains. RESULTS: These species were highly differentiated by their gene expression profiles, as determined by 18-21% of transcripts. Gene ontology (GO) terms associated with reproductive processes were enriched according to comparisons between species in both allopatric and sympatric populations. Species showed differential expression in genes linked to salt and waterlogging tolerance according to comparisons between species in sympatry, and biological processes related to environmental stimulus appeared as representative among those transcripts associated with edaphic characteristics in each sympatric zone. Hybrids, in their turn, were well differentiated from E. fulgens, but exhibited a similar gene expression profile to flooding-tolerant E. puniceolutem. When compared with parental species, hybrids showed no transcripts with additive pattern of expression and increased expression for almost all transgressive transcripts. CONCLUSIONS: This study sheds light on general mechanisms promoting ecological differentiation and assortative mating, and suggests candidate genes, such as those encoding catalase and calcium-dependent protein kinase, underling adaptation to harsh edaphic conditions in the Neotropical coastal plains. Moreover, it demonstrates that differential gene expression plays a central role in determining ecologically-based phenotypic differences among co-occurring species and their hybrids.

Dispersal and local persistence shape the genetic structure of a widespread Neotropical plant species with a patchy distribution.

BACKGROUND AND AIMS: Isolated populations constitute an ideal laboratory to study the consequences of intraspecific divergence, because intrinsic incompatibilities are more likely to accumulate under reduced gene flow. Here, we use a widespread bromeliad with a patchy distribution, Pitcairnia lanuginosa, as a model to infer processes driving Neotropical diversification and, thus, to improve our understanding of the origin and evolutionary dynamics of biodiversity in this highly speciose region. METHODS: We assessed the timing of lineage divergence, genetic structural patterns and historical demography of P. lanuginosa, based on microsatellites, and plastid and nuclear sequence data sets using coalescent analyses and an Approximate Bayesian Computation framework. Additionally, we used species distribution models (SDMs) to independently estimate potential changes in habitat suitability. KEY RESULTS: Despite morphological uniformity, plastid and nuclear DNA data revealed two distinct P. lanuginosa lineages that probably diverged through dispersal from the Cerrado to the Central Andean Yungas, following the final uplift of the Andes, and passed through long-term isolation with no evidence of migration. Microsatellite data indicate low genetic diversity and high levels of inbreeding within populations, and restricted gene flow among populations, which are likely to be a consequence of bottlenecks (or founder events), and high selfing rates promoting population persistence in isolation. SDMs showed a slight expansion of the suitable range for P. lanuginosa lineages during the Last Glacial Maximum, although molecular data revealed a signature of older divergence. Pleistocene climatic oscillations thus seem to have played only a minor role in the diversification of P. lanuginosa, which probably persisted through adverse conditions in riparian forests. CONCLUSIONS: Our results imply drift as a major force shaping the evolution of P. lanuginosa, and suggest that dispersal events have a prominent role in connecting Neotropical open and forest biomes.

Temperature modulation of thermal tolerance of a CAM-tank bromeliad and the relationship with acid accumulation in different leaf regions.

Physiological changes that increase plant performance during exposure to high temperatures may play an inverse role during exposure to low temperatures. The objective of this study was to test variations in photosystem II response to heat and cold stress in the leaves of a bromeliad with crassulacean acid metabolism submitted to high or low temperatures. Leaves were maintained under constant temperatures of 10 and 35°C and used to examine possible relationships among physiological responses to high and low temperatures and organic acid accumulation. We also tested if distinct parts of bromeliad leaves show differences in photosynthetic thermotolerance. The samples from leaves maintained at 35°C showed greater heat tolerance values, while those from leaves maintained at 10°C showed lower cold tolerance values. Our results identified a strong negative relationship between the organic acid accumulation and thermal tolerance of bromeliad leaves that largely explained the differences in thermal tolerance among groups. One of these differences occurred among regions of a single leaf, with the base showing critical heat values of up to 8°C higher than the top region, suggesting a possible partitioning of leaf response among its regions. Differences in thermal tolerance were also observed between sampling times, with higher values observed in the morning.

Speciation with gene flow between two Neotropical sympatric species (Pitcairnia spp.: Bromeliaceae).

The study of mechanisms that generate new species is considered fundamental for broad areas of ecology and evolution. Speciation is a continuous process in which reproductive isolation is established, and it is of fundamental importance to understand the origins of the adaptations that contribute to this process. Hybrid zones are considered natural laboratories for the study of speciation and represent ideal systems for such studies. Here, we investigated genomic differentiation between hybridizing Neotropical species Pitcairnia staminea (G. Lodd.) and P. albiflos (Herb.). Using thousands of SNPs genotyped through RAD-seq, we estimate effective population sizes, interspecific gene flow, as well as time of divergence between these two sister species and identify candidate genomic regions for positive selection that may be related to reproductive isolation. We selected different scenarios of speciation and tested them by using approximate Bayesian computation (ABC); we found evidence of divergence with gradual reduction in gene flow between these species over time, compatible with the hypothesis of speciation with gene flow between these Pitcairnia species. The parameter estimates obtained through ABC suggested that the effective population size of P. albiflos was around three times larger than that of P. staminea. Our divergence date estimates showed that these two species diverged during the Pliocene (4.7 Mya; CI = 1.3-8.5 Mya), which has likely allowed this species to accumulate genome-wide differences. We also detected a total of 17 of 4165 loci which showed signatures of selection with high genetic differentiation (F ST > 0.85), 12 of these loci were annotated in de novo assembled transcriptomes of both species, and 4 candidate genes were identified to be putatively involved in reproductive isolation. These four candidate genes were previously associated with the function of pollen development, pollen tube germination and orientation, abiotic stress, and flower scent in plants, suggesting an interplay between pre- and postpollination barriers in the evolution of reproductive isolation between such species.

Deforestation is the turning point for the spreading of a weedy epiphyte: an IBM approach.

The rapid spread of many weeds into intensely disturbed landscapes is boosted by clonal growth and self-fertilization strategies, which conversely increases the genetic structure of populations. Here, we use empirical and modeling approaches to evaluate the spreading dynamics of Tillandsia recurvata (L.) L. populations, a common epiphytic weed with self-reproduction and clonal growth widespread in dry forests and deforested landscapes in the American continent. We introduce the TRec model, an individual-based approach to simulate the spreading of T. recurvata over time and across landscapes subjected to abrupt changes in tree density with the parameters adjusted according to the empirical genetic data based on microsatellites genotypes. Simulations with this model showed that the strong spatial genetic structure observed from empirical data in T. recurvata can be explained by a rapid increase in abundance and gene flow followed by stabilization after ca. 25 years. TRec model's results also indicate that deforestation is a turning point for the rapid increase in both individual abundance and gene flow among T. recurvata subpopulations occurring in formerly dense forests. Active reforestation can, in turn, reverse such a scenario, although with a milder intensity. The genetic-based study suggests that anthropogenic changes in landscapes may strongly affect the population dynamics of species with 'weedy' traits.