Distinct expression and methylation patterns for genes with different fates following a single whole-genome duplication in flowering plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K–pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein–protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein–protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant

Cyanobacterial diversity in Salar de Huasco, a high altitude saline wetland in Northern Chile, are highly similar to Antarctic cyanobacteria

The diversity of Cyanobacteria in water and sediment samples from four representative sites of the Salar de Huasco was examined using denaturing gradient gel electrophoresis and analysis of clone libraries of 16S rRNA gene PCR products. Salar de Huasco is a high altitude (3800 m altitude) saline wetland located in the Chilean Altiplano. We analyzed samples from a tributary stream (H0) and three shallow lagoons (H1, H4, H6) that contrasted in their physicochemical conditions and associated biota. Seventy-eight phylotypes were identified in a total of 268 clonal sequences deriving from seven clone libraries of water and sediment samples. Oscillatoriales were frequently found in water samples from sites H0, H1 and H4 and in sediment samples from sites H1 and H4. Pleurocapsales were found only at site H0, while Chroococcales were recovered from sediment samples of sites H0 and H1, and from water samples of site H1. Nostocales were found in sediment samples from sites H1 and H4, and water samples from site H1 and were largely represented by sequences highly similar to Nodularia spumigena. We suggest that cyanobacterial communities from Salar de Huasco are unique - they include sequences related to others previously described from the Antarctic, along with others from diverse, but less extreme environments

Climatically stable landscapes predict patterns of genetic structure and admixture in the Californian canyon live oak

[Aim]: We studied which factors shape contemporary patterns of genetic structure, diversity and admixture in the canyon live oak (Quercus chrysolepis). Specifically, we tested two alternative hypotheses: (1) that areas with high habitat suitability and stability since the Last Glacial Maximum (LGM) sustain higher effective population sizes, resulting in increased levels of genetic diversity; and (2) that populations from areas with lower habitat stability show higher levels of genetic admixture due to their recurrent colonization by individuals originating from genetically differentiated populations. Furthermore, we analysed the relative importance of past and current habitat suitability and their additive effects on contemporary patterns of genetic structure. [Location]: California, USA. [Methods]: We sampled 160 individuals from 33 localities across the distribution range of the canyon live oak in California and then combined information from 13 nuclear microsatellite DNA markers and climate niche modelling to study patterns of genetic variation in this species. We used Bayesian clustering analyses to analyse geographical patterns of genetic structure and admixture, and circuit theory to generate isolation-by-resistance (IBR) distance matrices. [Results]: We found that the degree of genetic admixture was higher in localities with lower inferred population stability, but that genetic diversity was not associated with habitat suitability or stability. Landscape genetic analyses identified habitat stability as the primary driver of population genetic differentiation. [Main conclusions]: This study shows that habitat stability can be a major factor shaping genetic variation in wind-pollinated trees and supports the idea that stable regions contribute to genetic connectivity across different climatic periods. To our knowledge, this study is the first to report an association between patterns of genetic admixture and stability of local habitat.During this work J.O. was supported by Juan de la Cierva (MICINN), José Castillejo (ME) and Severo Ochoa (EBD) research fellowships. P.F.G. received post-doctoral support from a UCLA research award to V.L.S. This workreceived financial support from grants CGL2011-25053(MICINN) and UNCM08-1E-018 (FEDER).Peer reviewe

Data from: Influence of late Quaternary climate change on present patterns of genetic variation in valley oak, Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate-based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large-scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions

Data from: Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analyzed inter- and intra-specific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models, and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with species distribution models and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process

Data from: Genome-wide signature of local adaptation linked to variable CpG methylation in oak populations

It has long been known that adaptive evolution can occur through genetic mutations in DNA sequence, but it is unclear whether adaptive evolution can occur through analogous epigenetic mechanisms, such as through DNA methylation. If epigenetic variation contributes directly to evolution, species under threat of disease, invasive competition, climate change or other stresses would have greater stores of variation from which to draw. We looked for evidence of natural selection acting on variably methylated DNA sites using population genomic analysis across three climatologically distinct populations of valley oaks. We found patterns of genetic and epigenetic differentiations that indicate local adaptation is operating on large portions of the oak genome. While CHG methyl polymorphisms are not playing a significant role and would make poor targets for natural selection, our findings suggest that CpG methyl polymorphisms as a whole are involved in local adaptation, either directly or through linkage to regions under selection

batch_1

VCF output from Stacks

Input files used for Mantel tests

Input files for Mantel tests performed using the program ZT. These files include two indexes of niche overlap (D and I) between each pair of species and different estimates of interspecific genetic differentiation (Fst, Reynolds distance, G`st and Jost`s D) calculated for 16 nuclear (nuSSR) and 5 chloroplast (cpSSR) simple sequence repeat loci

Evolutionary lessons from California plant phylogeography

Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change