Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability.

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

Genetic control of functional traits related to photosynthesis and water use efficiency in Pinus pinaster Ait. drought response: integration of genome annotation, allele association and QTL detection for candidate gene identification.

BACKGROUND: Understanding molecular mechanisms that control photosynthesis and water use efficiency in response to drought is crucial for plant species from dry areas. This study aimed to identify QTL for these traits in a Mediterranean conifer and tested their stability under drought. RESULTS: High density linkage maps for Pinus pinaster were used in the detection of QTL for photosynthesis and water use efficiency at three water irrigation regimes. A total of 28 significant and 27 suggestive QTL were found. QTL detected for photochemical traits accounted for the higher percentage of phenotypic variance. Functional annotation of genes within the QTL suggested 58 candidate genes for the analyzed traits. Allele association analysis in selected candidate genes showed three SNPs located in a MYB transcription factor that were significantly associated with efficiency of energy capture by open PSII reaction centers and specific leaf area. CONCLUSIONS: The integration of QTL mapping of functional traits, genome annotation and allele association yielded several candidate genes involved with molecular control of photosynthesis and water use efficiency in response to drought in a conifer species. The results obtained highlight the importance of maintaining the integrity of the photochemical machinery in P. pinaster drought response.

Maritime Pine Rootstock Genotype Modulates Gene Expression Associated with Stress Tolerance in Grafted Stems.

Climate change-induced hazards, such as drought, threaten forest resilience, particularly in vulnerable regions such as the Mediterranean Basin. Maritime pine (Pinus pinaster Aiton), a model species in Western Europe, plays a crucial role in the Mediterranean forest due to its genetic diversity and ecological plasticity. This study characterizes transcriptional profiles of scion and rootstock stems of four P. pinaster graft combinations grown under well-watered conditions. Our grafting scheme combined drought-sensitive and drought-tolerant genotypes for scions (GAL1056: drought-sensitive scion; and Oria6: drought-tolerant scion) and rootstocks (R1S: drought-sensitive rootstock; and R18T: drought-tolerant rootstock). Transcriptomic analysis revealed expression patterns shaped by genotype provenance and graft combination. The accumulation of differentially expressed genes (DEGs) encoding proteins, involved in defense mechanisms and pathogen recognition, was higher in drought-sensitive scion stems and also increased when grafted onto drought-sensitive rootstocks. DEGs involved in drought tolerance mechanisms were identified in drought-tolerant genotypes as well as in drought-sensitive scions grafted onto drought-tolerant rootstocks, suggesting their establishment prior to drought. These mechanisms were associated with ABA metabolism and signaling. They were also involved in the activation of the ROS-scavenging pathways, which included the regulation of flavonoid and terpenoid metabolisms. Our results reveal DEGs potentially associated with the conifer response to drought and point out differences in drought tolerance strategies. These findings suggest genetic trade-offs between pine growth and defense, which could be relevant in selecting more drought-tolerant Pinus pinaster trees.

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.

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.

Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms.

Adaptation of long-living forest trees to respond to environmental changes is essential to secure their performance under adverse conditions. Water deficit is one of the most significant stress factors determining tree growth and survival. Maritime pine (Pinus pinaster Ait.), the main source of softwood in southwestern Europe, is subjected to recurrent drought periods which, according to climate change predictions for the years to come, will progressively increase in the Mediterranean region. The mechanisms regulating pine adaptive responses to environment are still largely unknown. The aim of this work was to go a step further in understanding the molecular mechanisms underlying maritime pine response to water stress and drought tolerance at the whole plant level. A global transcriptomic profiling of roots, stems, and needles was conducted to analyze the performance of siblings showing contrasted responses to water deficit from an ad hoc designed full-sib family. Although P. pinaster is considered a recalcitrant species for vegetative propagation in adult phase, the analysis was conducted using vegetatively propagated trees exposed to two treatments: well-watered and moderate water stress. The comparative analyses led us to identify organ-specific genes, constitutively expressed as well as differentially expressed when comparing control versus water stress conditions, in drought-sensitive and drought-tolerant genotypes. Different response strategies can point out, with tolerant individuals being pre-adapted for coping with drought by constitutively expressing stress-related genes that are detected only in latter stages on sensitive individuals subjected to drought.

Analysis of DNA Cytosine Methylation Patterns Using Methylation-Sensitive Amplification Polymorphism (MSAP).

Different molecular techniques have been developed to study either the global level of methylated cytosines or methylation at specific gene sequences. One of them is the methylation-sensitive amplified polymorphism technique (MSAP) which is a modification of amplified fragment length polymorphism (AFLP). It has been used to study methylation of anonymous CCGG sequences in different fungi, plants, and animal species. The main variation of this technique resides on the use of isoschizomers with different methylation sensitivity (such as HpaII and MspI) as a frequent-cutter restriction enzyme. For each sample, MSAP analysis is performed using both EcoRI/HpaII- and EcoRI/MspI-digested samples. A comparative analysis between EcoRI/HpaII and EcoRI/MspI fragment patterns allows the identification of two types of polymorphisms: (1) methylation-insensitive polymorphisms that show common EcoRI/HpaII and EcoRI/MspI patterns but are detected as polymorphic amplified fragments among samples and (2) methylation-sensitive polymorphisms which are associated with the amplified fragments that differ in their presence or absence or in their intensity between EcoRI/HpaII and EcoRI/MspI patterns. This chapter describes a detailed protocol of this technique and discusses the modifications that can be applied to adjust the technology to different species of interest.

New insights on glyphosate mode of action in nodular metabolism: Role of shikimate accumulation.

The short-term effects of the herbicide glyphosate (1.25-10 mM) on the growth, nitrogen fixation, carbohydrate metabolism, and shikimate pathway were investigated in leaves and nodules of nodulated lupine plants. All glyphosate treatments decreased nitrogenase activity rapidly (24 h) after application, even at the lowest and sublethal dose used (1.25 mM). This early effect on nitrogenase could not be related to either damage to nitrogenase components (I and II) or limitation of carbohydrates supplied by the host plant. In fact, further exposure to increasing glyphosate concentrations (5 mM) and greater time after exposure (5 days) decreased nodule starch content and sucrose synthase (SS; EC 2.4.1.13) activity but increased sucrose content within the nodule. These effects were accompanied by a great inhibition of the activity of phosphoenolpyruvate carboxylase (PEPC; EC 4.1.1.31). There were remarkable and rapid effects on the increase of shikimic and protocatechuic (PCA) acids in nodules and leaves after herbicide application. On the basis of the role of shikimic acid and PCA in the regulation of PEPC, as potent competitive inhibitors, this additional effect provoked by glyphosate on 5-enolpyruvylshikimic-3-phosphate synthase enzyme (EPSPS; EC 2.5.1.19) inhibition would divert most PEP into the shikimate pathway, depriving energy substrates to bacteroids to maintain nitrogen fixation. These findings provide a new explanation for the effectiveness of glyphosate as a herbicide in other plant tissues, for the observed differences in tolerance among species or cultivars, and for the transitory effects on glyphosate-resistant transgenic crops under several environmental conditions.

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.

Alterations induced by glyphosate on lupin photosynthetic apparatus and nodule ultrastructure and some oxygen diffusion related proteins.

The effects of glyphosate on protein metabolism, mesophyll cell ultrastructure and nodule ultrastructure and functioning of Lupinus albus cv. Multolupa inoculated with Bradyrhizobium sp. (Lupinus) were investigated. Young leaves and nodules were especially affected because these organs act as sinks of the herbicide. The alterations on nodular and chloroplast ultrastructure varied depending on herbicide concentration and time of exposure. After 3 days of 2.5 mM glyphosate application some toxic effects were detected. The most important alterations on nodules were the progressive cellular degradation of plant and bacteroidal cytosol and the rupture of bacteroidal membrane, whilst the peribacteroid membrane of the symbiosomes was preserved. This is the first report on the effect of glyphosate on legume-nodule ultrastructure. Glyphosate inhibited B. sp. (Lupinus) growth at concentrations higher than 62.5 microM. In the mesophyll cells, gradual disorganization of grana and intergrana was observed, loosing the parallel alignment with the chloroplast axis. As in nodules, degradation of membrane systems was observed, with the deformation, and even the rupture, of the tonoplast. These progressive effects were similar to those described in senescence processes. The adverse effects produced on infected zone can be due both to a direct effect of the herbicide on microsymbiont and to an indirect effect of glyphosate action on photosynthetic apparatus. Glyphosate produced changes in nodule cytosol and bacteroid proteins content and polypeptide pattern of leaves and nodules. With respect to proteins related to the oxygen diffusion mechanism, a large decrease in leghemoglobin and glycoproteins (recognized by antibodies MAC236 and MAC265) content was detected, which suggests that the oxygen diffusion mechanisms were also affected by glyphosate.

Nucleotide polymorphisms in a pine ortholog of the Arabidopsis degrading enzyme cellulase KORRIGAN are associated with early growth performance in Pinus pinaster.

We have carried out a candidate-gene-based association genetic study in Pinus pinaster Aiton and evaluated the predictive performance for genetic merit gain of the most significantly associated genes and single nucleotide polymorphisms (SNPs). We used a second generation 384-SNP array enriched with candidate genes for growth and wood properties to genotype mother trees collected in 20 natural populations covering most of the European distribution of the species. Phenotypic data for total height, polycyclism, root-collar diameter and biomass were obtained from a replicated provenance-progeny trial located in two sites with contrasting environments (Atlantic vs Mediterranean climate). General linear models identified strong associations between growth traits (total height and polycyclism) and four SNPs from the korrigan candidate gene, after multiple testing corrections using false discovery rate. The combined genomic breeding value predictions assessed for the four associated korrigan SNPs by ridge regression-best linear unbiased prediction (RR-BLUP) and cross-validation accounted for up to 8 and 15% of the phenotypic variance for height and polycyclic growth, respectively, and did not improve adding SNPs from other growth-related candidate genes. For root-collar diameter and total biomass, they accounted for 1.6 and 1.1% of the phenotypic variance, respectively, but increased to 15 and 4.1% when other SNPs from lp3.1, lp3.3 and cad were included in RR-BLUP models. These results point towards a desirable integration of candidate-gene studies as a means to pre-select relevant markers, and aid genomic selection in maritime pine breeding programs.

Evidence of intense chromosomal shuffling during conifer evolution.

While recent advances have been gained on genome evolution in angiosperm lineages, virtually nothing is known about karyotype evolution in the other group of seed plants, the gymnosperms. Here we used high density gene-based linkage mapping to compare the karyotype structure of two families of conifers (the most abundant group of gymnosperms) separated around 290 million years ago: Pinaceae and Cupressaceae. We propose for the first time a model based on the fusion of 20 ancestral chromosomal blocks that may have shaped the modern karyotpes of Pinaceae (with n=12) and Cupressaceae (with n=11). The considerable difference in modern genome organization between these two lineages contrasts strongly with the remarkable level of synteny already reported within the Pinaceae. It also suggests a convergent evolutionary mechanism of chromosomal block shuffling that has shaped the genomes of the spermatophytes.

Isoschizomers and amplified fragment length polymorphism for the detection of specific cytosine methylation changes.

Different molecular techniques have been developed to study either the global level of methylated cytosines or methylation at specific gene sequences. One of them is a modification of the Amplified Fragment Length Polymorphism (AFLP) technique that has been used to study methylation of anonymous CCGG sequences in different fungi, plant and animal species. The main variation of this technique is based on the use of isoschizomers with different methylation sensitivity (such as HpaII and MspI) as a frequent cutter restriction enzyme. For each sample, AFLP analysis is performed using both EcoRI/HpaII and EcoRI/MspI digested samples. Comparative analysis between EcoRI/HpaII and EcoRI/MspI fragment patterns allows the identification of two types of polymorphisms: (1) "Methylation-insensitive polymorphisms" that show common EcoRI/HpaII and EcoRI/MspI patterns but are detected as polymorphic amplified fragments among samples; and (2) "Methylation-sensitive polymorphisms" that are associated with amplified fragments differing in their presence or absence or in their intensity between EcoRI/HpaII and EcoRI/MspI patterns. This chapter describes a detailed protocol of this technique and discusses modifications that can be applied to adjust the technology to different species of interest.

Putative porin of Bradyrhizobium sp. (Lupinus) bacteroids induced by glyphosate.

Application of glyphosate (N-[phosphonomethyl] glycine) to Bradyrhizobium sp. (Lupinus)-nodulated lupin plants caused modifications in the protein pattern of bacteroids. The most significant change was the presence of a 44-kDa polypeptide in bacteroids from plants treated with the higher doses of glyphosate employed (5 and 10 mM). The polypeptide has been characterized by the amino acid sequencing of its N terminus and the isolation and nucleic acid sequencing of its encoding gene. It is putatively encoded by a single gene, and the protein has been identified as a putative porin. Protein modeling revealed the existence of several domains sharing similarity to different porins, such as a transmembrane beta-barrel. The protein has been designated BLpp, for Bradyrhizobium sp. (Lupinus) putative porin, and would be the first porin described in Bradyrhizobium sp. (Lupinus). In addition, a putative conserved domain of porins has been identified which consists of 87 amino acids, located in the BLpp sequence 30 amino acids downstream of the N-terminal region. In bacteroids, mRNA of the BLpp gene shows a basal constitutive expression that increases under glyphosate treatment, and the expression of the gene is seemingly regulated at the transcriptional level. By contrast, in free-living bacteria glyphosate treatment leads to an inhibition of BLpp mRNA accumulation, indicating a different effect of glyphosate on BLpp gene expression in bacteroids and free-living bacteria. The possible role of BLpp in a metabolite interchange between Bradyrhizobium and lupin is discussed.