Bacteriome dataset from the rhizosphere of trees in a Pinus pinaster and Pinus halepensis dominated forest subjected to drought conditions.

The Mediterranean basin is drastically affected by intense and frequent droughts, which jeopardize the diversity and survival of its forest, for example, Pinus pinaster forests. The dynamics of the bacterial communities inhabiting the rhizosphere of Pinus pinaster and other plants from a pine dominated forest under contrasting hydric conditions was monitored. The forest was located in Sierra de Oria (southeast Spain), and it was mainly composed by P. pinaster, P. halepensis, woody shrub species and herbaceous plants. 18 trees visually belonging to P. pinaster located along the perimeter and across the forest were selected for the analysis. All the trees were separated at least 50 m each other. Although all of them belonged to P. pinaster morphologically according to visual identification, the genotyping of the roots confirmed that they corresponded to P. pinaster, P. halepensis, and other plant species different from genus Pinus, although in the last case it was not possible to identify the plant species. At a distance less than 50 cm from the trunk, the litter and topsoil were removed, and the soil closely attached to non-suberified roots (rhizosphere soil) was collected (depth of 5-25 cm). Sampling was carried out in two seasons with contrasting temperature and rainfall patterns: on July 18, 2017 (summer) and April 24, 2018 (spring). After rhizosphere soil DNA and RNA extraction (and cDNA synthesis), a metabarcoding approach was followed by sequencing the V3-V4 hypervariable regions of the 16S rRNA gene and its derived transcripts by Illumina MiSeq platform. Sequencing reads were bioinformatically processed; specifically, they were filtered, trimmed, clustered into ASV (Amplicon Sequence Variants), and taxonomically identified. As a result, a total of 1,123,209 and 1,089,359 quality sequences were obtained from DNA and RNA-derived libraries, which resulted in 5,241 and 5,231 ASVs, respectively. Total communities (DNA) were mainly dominated by phyla Proteobacteria, Acidobacteria, Actinobacteria, Verrucomicrobia and Bacteroidetes in summer and spring, while potentially active populations (RNA libraries) were rich in Proteobacteria, Acidobacteria, Candidate division WPS-1, Actinobacteria and Verrucomicrobia both in summer and spring. On the other hand, DNA libraries were mainly dominated by genera Sphingomonas and acidobacterial groups Gp4 and Gp6, while potentially active bacteria (RNA) were rich in acidobacterial Gp3, Gp4, Gp6 and Phenylobacterium, although their relative abundance depended on the considered season. This dataset can provide valuable information about bacterial candidates which could be used as bioindicators of drought conditions. In addition to shifts in the bacterial relative abundance due to seasonal changes, the ratio RNA-based cDNA:DNA could be calculated as proxy of the potential activity of bacterial taxa. Moreover, these data could aid in developing bioformulations based on microorganisms which could be resistant to desiccation and involved in the drought resistance mechanisms of the host plant.

Correlating the above- and belowground genotype of Pinus pinaster trees and rhizosphere bacterial communities under drought conditions.

Increasing temperatures along with severe droughts are factors that may jeopardize the survival of the forests in the Mediterranean basin. In this region, Pinus pinaster is a common conifer species, that has been used as a model species in evolutionary studies due to its adaptive response to changing environments. Although its drought tolerance mechanisms are already known, knowledge about the dynamics of its root microbiota is still scarce. We aimed to decipher the structural (bacterial abundance), compositional, functional and associative changes of the P. pinaster rhizosphere bacterial communities in spring and summer, at DNA and RNA level (environmental DNA, live and dead cells, and those synthesizing proteins). A fundamental aspect of root microbiome-based approaches is to guarantee the correct origin of the samples. Thus, we assessed the genotype of host needles and roots from which rhizosphere samples were obtained. For more than 50% of the selected trees, genotype discrepancies were found and in three cases the plant species could not be determined. Rhizosphere bacterial communities were homogeneous with respect to diversity and structural levels regardless of the host genotype in both seasons. Nonetheless, significant changes were seen in the taxonomic profiles depending on the season. Seasonal changes were also evident in the bacterial co-occurrence patterns, both in DNA and RNA libraries. While spring communities switched to more complex networks, summer populations resulted in more compartmentalized networks, suggesting that these communities were facing a disturbance. These results may mirror the future status of bacterial communities in a context of climate change. A keystone hub was ascribed to the genus Phenylobacterium in the functional network calculated for summer. Overall, it is important to validate the origin and identity of plant samples in any plant-microbiota study so that more reliable ecological analyses are performed.

Organ-specific metabolic responses to drought in Pinus pinaster Ait.

Drought is an important driver of plant survival, growth, and distribution. Water deficit affects different pathways of metabolism, depending on plant organ. While previous studies have mainly focused on the metabolic drought response of a single organ, analysis of metabolic differences between organs is essential to achieve an integrated understanding of the whole plant response. In this work, untargeted metabolic profiling was used to examine the response of roots, stems, adult and juvenile needles from Pinus pinaster Ait. full-sib individuals, subjected to a moderate and long lasting drought period. Cyclitols content showed a significant alteration, in response to drought in all organs examined, but other metabolites increased or decreased differentially depending on the analyzed organ. While a high number of flavonoids were only detected in aerial organs, an induction of the glutathione pathway was mainly detected in roots. This result may reflect different antioxidant mechanisms activated in aerial organs and roots. Metabolic changes were more remarkable in roots than in the other organs, highlighting its prominent role in the response to water stress. Significant changes in flavonoids and ascorbate metabolism were also observed between adult and juvenile needles, consistent with previously proven differential functional responses between the two developmental stages. Genetic polymorphisms in candidate genes coding for a Myb1 transcription factor and a malate dehydrogenase (EC 1.1.1.37) were associated with different concentration of phenylalanine, phenylpropanoids and malate, respectively. The results obtained will support further research on metabolites and genes potentially involved in functional mechanisms related to drought tolerance in trees.

Epigenetic regulation of adaptive responses of forest tree species to the environment.

Epigenetic variation is likely to contribute to the phenotypic plasticity and adaptative capacity of plant species, and may be especially important for long-lived organisms with complex life cycles, including forest trees. Diverse environmental stresses and hybridization/polyploidization events can create reversible heritable epigenetic marks that can be transmitted to subsequent generations as a form of molecular "memory". Epigenetic changes might also contribute to the ability of plants to colonize or persist in variable environments. In this review, we provide an overview of recent data on epigenetic mechanisms involved in developmental processes and responses to environmental cues in plant, with a focus on forest tree species. We consider the possible role of forest tree epigenetics as a new source of adaptive traits in plant breeding, biotechnology, and ecosystem conservation under rapid climate change.

Annotated genetic linkage maps of Pinus pinaster Ait. from a Central Spain population using microsatellite and gene based markers.

BACKGROUND: Pinus pinaster Ait. is a major resin producing species in Spain. Genetic linkage mapping can facilitate marker-assisted selection (MAS) through the identification of Quantitative Trait Loci and selection of allelic variants of interest in breeding populations. In this study, we report annotated genetic linkage maps for two individuals (C14 and C15) belonging to a breeding program aiming to increase resin production. We use different types of DNA markers, including last-generation molecular markers. RESULTS: We obtained 13 and 14 linkage groups for C14 and C15 maps, respectively. A total of 211 and 215 markers were positioned on each map and estimated genome length was between 1,870 and 2,166 cM respectively, which represents near 65% of genome coverage. Comparative mapping with previously developed genetic linkage maps for P. pinaster based on about 60 common markers enabled aligning linkage groups to this reference map. The comparison of our annotated linkage maps and linkage maps reporting QTL information revealed 11 annotated SNPs in candidate genes that co-localized with previously reported QTLs for wood properties and water use efficiency. CONCLUSIONS: This study provides genetic linkage maps from a Spanish population that shows high levels of genetic divergence with French populations from which segregating progenies have been previously mapped. These genetic maps will be of interest to construct a reliable consensus linkage map for the species. The importance of developing functional genetic linkage maps is highlighted, especially when working with breeding populations for its future application in MAS for traits of interest.

Towards decoding the conifer giga-genome.

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world's conifer forests.

EuroPineDB: a high-coverage web database for maritime pine transcriptome.

BACKGROUND: Pinus pinaster is an economically and ecologically important species that is becoming a woody gymnosperm model. Its enormous genome size makes whole-genome sequencing approaches are hard to apply. Therefore, the expressed portion of the genome has to be characterised and the results and annotations have to be stored in dedicated databases. DESCRIPTION: EuroPineDB is the largest sequence collection available for a single pine species, Pinus pinaster (maritime pine), since it comprises 951 641 raw sequence reads obtained from non-normalised cDNA libraries and high-throughput sequencing from adult (xylem, phloem, roots, stem, needles, cones, strobili) and embryonic (germinated embryos, buds, callus) maritime pine tissues. Using open-source tools, sequences were optimally pre-processed, assembled, and extensively annotated (GO, EC and KEGG terms, descriptions, SNPs, SSRs, ORFs and InterPro codes). As a result, a 10.5× P. pinaster genome was covered and assembled in 55 322 UniGenes. A total of 32 919 (59.5%) of P. pinaster UniGenes were annotated with at least one description, revealing at least 18 466 different genes. The complete database, which is designed to be scalable, maintainable, and expandable, is freely available at: http://www.scbi.uma.es/pindb/. It can be retrieved by gene libraries, pine species, annotations, UniGenes and microarrays (i.e., the sequences are distributed in two-colour microarrays; this is the only conifer database that provides this information) and will be periodically updated. Small assemblies can be viewed using a dedicated visualisation tool that connects them with SNPs. Any sequence or annotation set shown on-screen can be downloaded. Retrieval mechanisms for sequences and gene annotations are provided. CONCLUSIONS: The EuroPineDB with its integrated information can be used to reveal new knowledge, offers an easy-to-use collection of information to directly support experimental work (including microarray hybridisation), and provides deeper knowledge on the maritime pine transcriptome.

"Contrasting patterns of selection at Pinus pinaster Ait. Drought stress candidate genes as revealed by genetic differentiation analyses".

The importance of natural selection for shaping adaptive trait differentiation among natural populations of allogamous tree species has long been recognized. Determining the molecular basis of local adaptation remains largely unresolved, and the respective roles of selection and demography in shaping population structure are actively debated. Using a multilocus scan that aims to detect outliers from simulated neutral expectations, we analyzed patterns of nucleotide diversity and genetic differentiation at 11 polymorphic candidate genes for drought stress tolerance in phenotypically contrasted Pinus pinaster Ait. populations across its geographical range. We compared 3 coalescent-based methods: 2 frequentist-like, including 1 approach specifically developed for biallelic single nucleotide polymorphisms (SNPs) here and 1 Bayesian. Five genes showed outlier patterns that were robust across methods at the haplotype level for 2 of them. Two genes presented higher F(ST) values than expected (PR-AGP4 and erd3), suggesting that they could have been affected by the action of diversifying selection among populations. In contrast, 3 genes presented lower F(ST) values than expected (dhn-1, dhn2, and lp3-1), which could represent signatures of homogenizing selection among populations. A smaller proportion of outliers were detected at the SNP level suggesting the potential functional significance of particular combinations of sites in drought-response candidate genes. The Bayesian method appeared robust to low sample sizes, flexible to assumptions regarding migration rates, and powerful for detecting selection at the haplotype level, but the frequentist-like method adapted to SNPs was more efficient for the identification of outlier SNPs showing low differentiation. Population-specific effects estimated in the Bayesian method also revealed populations with lower immigration rates, which could have led to favorable situations for local adaptation. Outlier patterns are discussed in relation to the different genes' putative involvement in drought tolerance responses, from published results in transcriptomics and association mapping in P. pinaster and other related species. These genes clearly constitute relevant candidates for future association studies in P. pinaster.