Foliar mycobiome remains unaltered under urban air-pollution but differentially express stress-related genes.

Air pollution caused by tropospheric ozone contributes to the decline of forest ecosystems; for instance, sacred fir, Abies religiosa (Kunth) Schltdl. & Cham. forests in the peri-urban region of Mexico City. Individual trees within these forests exhibit variation in their response to ozone exposure, including the severity of visible symptoms in needles. Using RNA-Seq metatranscriptomic data and ITS2 metabarcoding, we investigated whether symptom variation correlates with the taxonomic and functional composition of fungal mycobiomes from needles collected in this highly polluted area in the surroundings of Mexico City. Our findings indicate that ozone-related symptoms do not significantly correlate with changes in the taxonomic composition of fungal mycobiomes. However, genes coding for 30 putative proteins were differentially expressed in the mycobiome of asymptomatic needles, including eight genes previously associated with resistance to oxidative stress. These results suggest that fungal communities likely play a role in mitigating the oxidative burst caused by tropospheric ozone in sacred fir. Our study illustrates the feasibility of using RNA-Seq data, accessible from global sequence repositories, for the characterization of fungal communities associated with plant tissues, including their gene expression.

Stochastic processes and changes in evolutionary rate are associated with diversification in a lineage of tropical hard pines (Pinus).

The study of the patterns of polymorphism and molecular evolution among closely related species is key to understanding the evolutionary forces involved in the diversification of lineages. This point is a big challenge in species with slow evolutionary rates, long life cycles, and ancient, shared polymorphisms such as conifers. Under the premise of divergence in a stepwise migration process, we expect clinal geographical patterns of purifying selection efficiency, and genetic structure related to latitude or longitude. If migration is accompanied by changes in the environment, we could further expect a role of positive selection in driving species divergence. Here, we infer patterns of polymorphism, efficiency of purifying selection, and molecular evolution using a dataset of 161 nuclear genes (∼71 Kb) in a lineage of hard pines from North America, the Caribbean, Mexico, and Central America presumed to have migrated from North America toward lower latitudes with tropical conditions. Under the premise of differences in selective pressures, we also look for possible signals of positive selection. To test our hypothesis, first we estimated different indices to infer patterns of polymorphism and efficiency of purifying selection (Ka, Ks, Ka/Ks, dN, dS, dN/dS, and dxy) and compared these metrics across five clades. Also, we investigated possible clinal patterns in these indices and morphological traits (needle length and cone length). Then we inferred genetic structure and environmental differences among species to test for possible signals of positive selection using phylogenetic methods in specific clades. We found differences among clades using Ka, Ks, and Ka/Ks with a relaxation of purifying selection, especially in the Elliotti and Patula clades. We also found environmental differences related to geographic distance, and among clades suggesting differences in selective pressures. The indices Ks, dxy, and needle length had relationships with geography but not ovulate cone length. Finally, we found that most analyzed genes are under purifying selection, but there was an exception of faster evolutionary rate in some pine species, suggesting the possible action of positive selection in divergence. Our study indicated that stochastic processes have played a key role in the diversification of the group, with a possible input of positive selection in pines from Mexico and Central America.

Evaluating pollution-related damage and restoration success in urban forests with participatory monitoring and digital tools.

Peri-urban forest monitoring requires indicators of vegetation damage. An example is the sacred fir (Abies religiosa) forests surrounding Mexico City, which have been heavily exposed to tropospheric ozone, a harmful pollutant, for over 4 decades. We developed a participatory monitoring system with which local community members and scientists generated data on ozone tree damage. Santa Rosa Xochiac rangers (13) used the digital tool KoboToolBox to record ozone damage to trees, tree height, tree ages, tree condition, tree position, and whether the tree had been planted. Thirty-five percent of the trees (n = 1765) had ozone damage. Younger trees had a lower percentage of foliage damaged by ozone than older trees (p < 0.0001), and asymptomatic trees tended to be younger (p < 0.0001). Symptomatic trees were taller than asymptomatic trees of the same age (R2 c  = 0.43, R2 m  = 0.27). Involving local communities facilitated forest monitoring and using digital technology improved data quality. This participatory system can be used to monitor forest condition change over time and thus aids restoration efforts driven by government or local communities' interests, facilitating local decision-making.

Evaluación del daño relacionado a la contaminación y del éxito de la restauración de los bosques urbanos con un monitoreo participativo y herramientas digitales Resumen El monitoreo de los bosques periurbanos requiere indicadores de daños en la vegetación. Un ejemplo son los bosques de abeto (Abies religiosa) que rodean la Ciudad de México, pues desde hace más de 4 décadas han estado expuestos al ozono troposférico, un contaminante nocivo. Desarrollamos un sistema de monitoreo participativo con el que miembros de la comunidad local y científicos generaron datos sobre los daños causados por el ozono en los árboles. Los guardabosques de Santa Rosa Xochiac (133) utilizaron la herramienta digital KoboToolBox para registrar los daños causados por el ozono en los árboles, su altura, edad, estado, posición y si eran árboles plantados. El 35% de los árboles (n = 1,765) presentó daños por ozono. Los árboles más jóvenes tenían un menor porcentaje de follaje dañado por el ozono que los árboles más viejos (p<0.0001), y los árboles asintomáticos tendían a ser más jóvenes (p<0.0001). Los árboles sintomáticos eran más altos que los asintomáticos de la misma edad (R2 c = 0.43, R2 m = 0.27). La participación de las comunidades locales facilitó el monitoreo forestal y el uso de tecnología digital mejoró la calidad de los datos. Este sistema participativo puede utilizarse para monitorear los cambios en el estado de los bosques a lo largo del tiempo y contribuir a los esfuerzos de restauración impulsados por el gobierno o las comunidades locales, facilitando la toma de decisiones a nivel local.

Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico.

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Comparative analysis of differential gene expression indicates divergence in ontogenetic strategies of leaves in two conifer genera.

In land plants, heteroblasty broadly refers to a drastic change in morphology during growth through ontogeny. Juniperus flaccida and Pinus cembroides are conifers of independent lineages known to exhibit leaf heteroblasty between the juvenile and adult life stage of development. Juvenile leaves of P. cembroides develop spirally on the main stem and appear decurrent, flattened, and needle-like; whereas adult photosynthetic leaves are triangular or semi-circular needle-like, and grow in whorls on secondary or tertiary compact dwarf shoots. By comparison, J. flaccida juvenile leaves are decurrent and needle-like, and adult leaves are compact, short, and scale-like. Comparative analyses were performed to evaluate differences in anatomy and gene expression patterns between developmental phases in both species. RNA from 12 samples was sequenced and analyzed with available software. They were assembled de novo from the RNA-Seq reads. Following assembly, 63,741 high-quality transcripts were functionally annotated in P. cembroides and 69,448 in J. flaccida. Evaluation of the orthologous groups yielded 4140 shared gene families among the four references (adult and juvenile from each species). Activities related to cell division and development were more abundant in juveniles than adults in P. cembroides, and more abundant in adults than juveniles in J. flaccida. Overall, there were 509 up-regulated and 81 down-regulated genes in the juvenile condition of P. cembroides and 14 up-regulated and 22 down-regulated genes in J. flaccida. Gene interaction network analysis showed evidence of co-expression and co-localization of up-regulated genes involved in cell wall and cuticle formation, development, and phenylpropanoid pathway, in juvenile P. cembroides leaves. Whereas in J. flaccida, differential expression and gene interaction patterns were detected in genes involved in photosynthesis and chloroplast biogenesis. Although J. flaccida and P. cembroides both exhibit leaf heteroblastic development, little overlap was detected, and unique genes and pathways were highlighted in this study.

Polygenic adaptation and negative selection across traits, years and environments in a long-lived plant species (Pinus pinaster Ait., Pinaceae).

A decade of genetic association studies in multiple organisms suggests that most complex traits are polygenic; that is, they have a genetic architecture determined by numerous loci, each with small effect-size. Thus, determining the degree of polygenicity and its variation across traits, environments and time is crucial to understand the genetic basis of phenotypic variation. We applied multilocus approaches to estimate the degree of polygenicity of fitness-related traits in a long-lived plant (Pinus pinaster Ait., maritime pine) and to analyse this variation across environments and years. We evaluated five categories of fitness-related traits (survival, height, phenology, functional, and biotic-stress response) in a clonal common-garden network planted in contrasted environments (over 20,500 trees). Most of the analysed traits showed evidence of local adaptation based on Qst -Fst comparisons. We further observed a remarkably stable degree of polygenicity, averaging 6% (range of 0%-27%), across traits, environments and years. We detected evidence of negative selection, which could explain, at least partially, the high degree of polygenicity. Because polygenic adaptation can occur rapidly, our results suggest that current predictions on the capacity of natural forest tree populations to adapt to new environments should be revised, especially in the current context of climate change.

Divergence with gene flow is driven by local adaptation to temperature and soil phosphorus concentration in teosinte subspecies (Zea mays parviglumis and Zea mays mexicana).

Patterns of genomic divergence between hybridizing taxa can be heterogeneous along the genome. Both differential introgression and local adaptation may contribute to this pattern. Here, we analysed two teosinte subspecies, Zea mays ssp. parviglumis and ssp. mexicana, to test whether their divergence has occurred in the face of gene flow and to infer which environmental variables have been important drivers of their ecological differentiation. We generated 9,780 DArTseqTM SNPs for 47 populations, and used an additional data set containing 33,454 MaizeSNP50 SNPs for 49 populations. With these data, we inferred features of demographic history and performed genome wide scans to determine the number of outlier SNPs associated with climate and soil variables. The two data sets indicate that divergence has occurred or been maintained despite continuous gene flow and/or secondary contact. Most of the significant SNP associations were to temperature and to phosphorus concentration in the soil. A large proportion of these candidate SNPs were located in regions of high differentiation that had been identified previously as putative inversions. We therefore propose that genomic differentiation in teosintes has occurred by a process of adaptive divergence, with putative inversions contributing to reduced gene flow between locally adapted populations.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding.

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.

Local effects drive heterozygosity-fitness correlations in an outcrossing long-lived tree.

Heterozygosity-fitness correlations (HFCs) have been used to understand the complex interactions between inbreeding, genetic diversity and evolution. Although frequently reported for decades, evidence for HFCs was often based on underpowered studies or inappropriate methods, and hence their underlying mechanisms are still under debate. Here, we used 6100 genome-wide single nucleotide polymorphisms (SNPs) to test for general and local effect HFCs in maritime pine (Pinus pinaster Ait.), an iconic Mediterranean forest tree. Survival was used as a fitness proxy, and HFCs were assessed at a four-site common garden under contrasting environmental conditions (total of 16 288 trees). We found no significant correlations between genome-wide heterozygosity and fitness at any location, despite variation in inbreeding explaining a substantial proportion of the total variance for survival. However, four SNPs (including two non-synonymous mutations) were involved in significant associations with survival, in particular in the common gardens with higher environmental stress, as shown by a novel heterozygosity-fitness association test at the species-wide level. Fitness effects of SNPs involved in significant HFCs were stable across maritime pine gene pools naturally growing in distinct environments. These results led us to dismiss the general effect hypothesis and suggested a significant role of heterozygosity in specific candidate genes for increasing fitness in maritime pine. Our study highlights the importance of considering the species evolutionary and demographic history and different spatial scales and testing environments when assessing and interpreting HFCs.

Molecular proxies for climate maladaptation in a long-lived tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Development and implementation of a highly-multiplexed SNP array for genetic mapping in maritime pine and comparative mapping with loblolly pine.

BACKGROUND: Single nucleotide polymorphisms (SNPs) are the most abundant source of genetic variation among individuals of a species. New genotyping technologies allow examining hundreds to thousands of SNPs in a single reaction for a wide range of applications such as genetic diversity analysis, linkage mapping, fine QTL mapping, association studies, marker-assisted or genome-wide selection. In this paper, we evaluated the potential of highly-multiplexed SNP genotyping for genetic mapping in maritime pine (Pinus pinaster Ait.), the main conifer used for commercial plantation in southwestern Europe. RESULTS: We designed a custom GoldenGate assay for 1,536 SNPs detected through the resequencing of gene fragments (707 in vitro SNPs/Indels) and from Sanger-derived Expressed Sequenced Tags assembled into a unigene set (829 in silico SNPs/Indels). Offspring from three-generation outbred (G2) and inbred (F2) pedigrees were genotyped. The success rate of the assay was 63.6% and 74.8% for in silico and in vitro SNPs, respectively. A genotyping error rate of 0.4% was further estimated from segregating data of SNPs belonging to the same gene. Overall, 394 SNPs were available for mapping. A total of 287 SNPs were integrated with previously mapped markers in the G2 parental maps, while 179 SNPs were localized on the map generated from the analysis of the F2 progeny. Based on 98 markers segregating in both pedigrees, we were able to generate a consensus map comprising 357 SNPs from 292 different loci. Finally, the analysis of sequence homology between mapped markers and their orthologs in a Pinus taeda linkage map, made it possible to align the 12 linkage groups of both species. CONCLUSIONS: Our results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in maritime pine, a conifer species that has a genome seven times the size of the human genome. This SNP-array will be extended thanks to recent sequencing effort using new generation sequencing technologies and will include SNPs from comparative orthologous sequences that were identified in the present study, providing a wider collection of anchor points for comparative genomics among the conifers.

Diverging patterns of mitochondrial and nuclear DNA diversity in subarctic black spruce: imprint of a founder effect associated with postglacial colonization.

High-latitude ecotonal populations at the species margins may exhibit altered patterns of genetic diversity, resulting from more or less recent founder events and from bottleneck effects in response to climate oscillations. Patterns of genetic diversity were investigated in nine populations of the conifer black spruce (Picea mariana [Mill.] BSP.) in northwestern Québec, Canada, using seed-dispersed mitochondrial (mt) DNA and nuclear (nc) DNA. mtDNA diversity (mitotypes) was assessed at three loci, and ncDNA diversity was estimated for nine expressed sequence tag polymorphism (ESTP) loci. Sampling included populations from the boreal forest and the southern and northern subzones of the subarctic forest-tundra, a fire-born ecotone. For ncDNA, populations from all three vegetation zones were highly diverse with little population differentiation (thetaN = 0.014); even the northernmost populations showed no loss of rare alleles. Patterns of mitotype diversity were strikingly different: within-population diversity and population differentiation were high for boreal forest populations [expected heterozygosity per locus (HE) = 0.58 and thetaM = 0.529], but all subarctic populations were fixed for a single mitotype (HE = 0). This lack of variation suggests a founder event caused by long-distance seed establishment during postglacial colonization, consistent with palaeoecological data. The estimated movement of seeds alone (effective number of migrants per generation, NmM < 2) was much restricted compared to that estimated from nuclear variants, which including pollen movement (NmN > 17). This could account for the conservation of a founder imprint in the mtDNA of subarctic black spruce. After reduction, presumably in the early Holocene, the diversity in ncDNA would have been replenished rapidly by pollen-mediated gene flow, and maintained subsequently through vegetative layering during the current cooler period covering the last 3000 years.