The Chinese pine genome and methylome unveil key features of conifer evolution.

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

A transcriptome-based association study of growth, wood quality, and oleoresin traits in a slash pine breeding population.

Slash pine (Pinus elliottii Engelm.) is an important timber and resin species in the United States, China, Brazil and other countries. Understanding the genetic basis of these traits will accelerate its breeding progress. We carried out a genome-wide association study (GWAS), transcriptome-wide association study (TWAS) and weighted gene co-expression network analysis (WGCNA) for growth, wood quality, and oleoresin traits using 240 unrelated individuals from a Chinese slash pine breeding population. We developed high quality 53,229 single nucleotide polymorphisms (SNPs). Our analysis reveals three main results: (1) the Chinese breeding population can be divided into three genetic groups with a mean inbreeding coefficient of 0.137; (2) 32 SNPs significantly were associated with growth and oleoresin traits, accounting for the phenotypic variance ranging from 12.3% to 21.8% and from 10.6% to 16.7%, respectively; and (3) six genes encoding PeTLP, PeAP2/ERF, PePUP9, PeSLP, PeHSP, and PeOCT1 proteins were identified and validated by quantitative real time polymerase chain reaction for their association with growth and oleoresin traits. These results could be useful for tree breeding and functional studies in advanced slash pine breeding program.

Comparative metabolomics and transcriptomics analyses provide insights into branching changes induced by heterografting in Pinus massoniana seedlings.

The investigation of dwarfing rootstocks for the establishment of high-generation seed orchards is a prospective avenue of research. In this investigation, Pinus massoniana, Pinus yunnanensis var. pygmaea (P. pygmaea), and P. elliottii seedlings were used as rootstocks for grafting with P. massoniana scions. Grafting P. massoniana onto P. pygmaea rootstock resulted in observable phenotypic alterations in lateral branches, apical buds, and needle length. Certain characteristic metabolites of rootstocks, such as fatty acyls, pregnenolones, steroids, and steroid derivatives, were found to be highly expressed in scions after grafting. RNA-seq analysis revealed MYB-related, SBP, and bHLH demonstrating a significant positive correlation, while C2H2 and Orphans exhibited negative correlations with the differential intensity of metabolites related to lipids and lipid-like molecules. This study offers valuable insights for the establishment of rootstock breeding programs.

Pine has two glutamine synthetase paralogs, GS1b.1 and GS1b.2, exhibiting distinct biochemical properties.

The enzyme glutamine synthetase (EC 6.3.1.2) is mainly responsible for the incorporation of inorganic nitrogen into organic molecules in plants. In the present work, a pine (Pinus pinaster) GS1 (PpGS1b.2) gene was identified, showing a high sequence identity with the GS1b.1 gene previously characterized in conifers. Phylogenetic analysis revealed that the presence of PpGS1b.2 is restricted to the genera Pinus and Picea and is not found in other conifers. Gene expression data suggest a putative role of PpGS1b.2 in plant development, similar to other GS1b genes from angiosperms, suggesting evolutionary convergence. The characterization of GS1b.1 and GS1b.2 at the structural, physicochemical, and kinetic levels has shown differences even though they have high sequence homology. GS1b.2 had a lower optimum pH (6 vs. 6.5) and was less thermally stable than GS1b.1. GS1b.2 exhibited positive cooperativity for glutamate and substrate inhibition for ammonium. However, GS1b.1 exhibited substrate inhibition behavior for glutamate and ATP. Alterations in the kinetic characteristics produced by site-directed mutagenesis carried out in this work strongly suggest an implication of amino acids at positions 264 and 267 in the active center of pine GS1b.1 and GS1b.2 being involved in affinity toward ammonium. Therefore, the amino acid differences between GS1b.1 and GS1b.2 would support the functioning of both enzymes to meet distinct plant needs.

Mitochondrial genome provides species-specific targets for the rapid detection of early invasive populations of Hylurgus ligniperda in China.

BACKGROUND: Hylurgus ligniperda, a major international forestry quarantine pest, was recently found to have invaded and posed a serious threat to the Pinus forests of the Jiaodong Peninsula in China. Continuous monitoring and vigilance of the early population is imperative, and rapid molecular detection technology is urgently needed. We focused on developing a single-gene-based species-specific PCR (SS-PCR) method. RESULTS: We sequenced and assembled the mitochondrial genome of H. ligniperda to identify suitable target genes. We identified three closely related species for detecting the specificity of SS-PCR through phylogenetic analysis based on 13 protein-coding genes (PCGs). Subsequently, we analyzed the evolution of 13 PCGs and selected four mitochondrial genes to represent slow-evolving gene (COI) and faster-evolving genes (e.g. ND2, ND4, and ND5), respectively. We developed four species-specific primers targeting COI, ND2, ND4, and ND5 to rapidly identify H. ligniperda. The results showed that the four species-specific primers exhibited excellent specificity and sensitivity in the PCR assays, with consistent performance across a broader range of species. This method demonstrates the ability to identify beetles promptly, even during their larval stage. The entire detection process can be completed within 2-3 h. CONCLUSIONS: This method is suitable for large-scale species detection in laboratory settings. Moreover, the selection of target genes in the SS-PCR method is not affected by the evolutionary rate. SS-PCR can be widely implemented at port and forestry workstations, significantly enhancing early management strategies and quarantine measures against H. ligniperda. This approach will help prevent the spread of the pest and effectively preserve the resources of Chinese pine forests.

Physiological responses to drought stress of three pine species and comparative transcriptome analysis of Pinus yunnanensis var. pygmaea.

Drought stress can significantly affect plant growth, development, and yield. Fewer comparative studies have been conducted between different species of pines, particularly involving Pinus yunnanensis var. pygmaea (P. pygmaea). In this study, the physiological indices, photosynthetic pigment and related antioxidant enzyme changes in needles from P. pygmaea, P. elliottii and P. massoniana under drought at 0, 7, 14, 21, 28 and 35 d, as well as 7 days after rehydration, were measured. The PacBio single-molecule real-time (SMRT) and Illumina RNA sequencing were used to uncover the gene expression differences in P. pygmaea under drought and rehydration conditions. The results showed that the total antioxidant capacity (TAOC) of P. pygmaea was significantly higher than P. massoniana and P. elliottii. TAOC showed a continuous increase trend across all species. Soluble sugar (SS), starch content and non-structural carbohydrate (NSC) of all three pines displayed a "W" pattern, declining initially, increasing, and then decreasing again. P. pygmaea exhibits stronger drought tolerance and greater recovery ability under prolonged drought conditions. Through the PacBio SMRT-seq, a total of 50,979 high-quality transcripts were generated, and 6,521 SSR and 5,561 long non-coding RNAs (LncRNAs) were identified. A total of 2310, 1849, 5271, 5947, 7710, and 6854 differentially expressed genes (DEGs) were identified compared to the control (Pp0D) in six pair-wise comparisons of treatment versus control. bHLH, NAC, ERF, MYB_related, C3H transcription factors (TFs) play an important role in drought tolerance of P. pygmaea. KEGG enrichment analysis and Gene set enrichment analysis (GSEA) analysis showed that P. pygmaea may respond to drought by enhancing metabolic processes such as ABA signaling pathway, alpha-linolenic acid. Weighted gene co-expression network analysis (WGCNA) revealed GST, CAT, LEC14B, SEC23 were associated with antioxidant enzyme activity and TAOC. This study provides a basis for further research on drought tolerance differences among coniferous species.

Integrated transcriptomic and metabolic analyses reveal the early response mechanism of Pinus tabulaeformis to pine wood nematodes.

Pine wilt disease (PWD) is a devastating disease of pine trees caused by the pine wood nematode (Bursapherenchus xylophilus, PWN). To study how Pinus tabulaeformis responds to PWD infection, we collected 3-year-old P. tabulaeformis seedlings at 2 days, 5 days, and 8 days after being infected with B. xylophilus. We identified genes and metabolites early responding to infection using transcriptome and metabolomic data obtained by high-throughput mRNA sequencing (RNA-seq) and liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based assays, respectively. The following results were obtained: (1) After inoculation with PWN, the average number of days taken for 3-year-old P. tabulaeformis seedlings to develop symptoms was 8 days. (2) Combined transcriptome and metabolome analysis revealed that phenylpropanoid biosynthesis and flavonoid biosynthesis are critically important pathways for P. tabulaeformis to respond to PWD. (3) The response of P. tabulaeformis to stress was mainly through positive regulation of gene expression, including some key genes related to plant hormones or transcription factors that have been widely studied. Genes related to pathways such as photosynthesis, plant-pathogen interactions, and DNA replication were downregulated. (4) Terpenoid biosynthesis genes involved during the development of pine wilt disease. This study demonstrated the defence and pathogenic mechanisms of P. tabulaeformis against PWD, providing a reference for the early diagnosis of PWD.

Response mechanism of carbon metabolism of Pinus massoniana to gradient high temperature and drought stress.

BACKGROUND: The carbon metabolism pathway is of paramount importance for the growth and development of plants, exerting a pivotal regulatory role in stress responses. The exacerbation of drought impacts on the plant carbon cycle due to global warming necessitates comprehensive investigation into the response mechanisms of Masson Pine (Pinus massoniana Lamb.), an exemplary pioneer drought-tolerant tree, thereby establishing a foundation for predicting future forest ecosystem responses to climate change. RESULTS: The seedlings of Masson Pine were utilized as experimental materials in this study, and the transcriptome, metabolome, and photosynthesis were assessed under varying temperatures and drought intensities. The findings demonstrated that the impact of high temperature and drought on the photosynthetic rate and transpiration rate of Masson Pine seedlings was more pronounced compared to individual stressors. The analysis of transcriptome data revealed that the carbon metabolic pathways of Masson Pine seedlings were significantly influenced by high temperature and drought co-stress, with a particular impact on genes involved in starch and sucrose metabolism. The metabolome analysis revealed that only trehalose and Galactose 1-phosphate were specifically associated with the starch and sucrose metabolic pathways. Furthermore, the trehalose metabolic heat map was constructed by integrating metabolome and transcriptome data, revealing a significant increase in trehalose levels across all three comparison groups. Additionally, the PmTPS1, PmTPS5, and PmTPPD genes were identified as key regulatory genes governing trehalose accumulation. CONCLUSIONS: The combined effects of high temperature and drought on photosynthetic rate, transpiration rate, transcriptome, and metabolome were more pronounced than those induced by either high temperature or drought alone. Starch and sucrose metabolism emerged as the pivotal carbon metabolic pathways in response to high temperature and drought stress in Masson pine. Trehalose along with PmTPS1, PmTPS5, and PmTPPD genes played crucial roles as metabolites and key regulators within the starch and sucrose metabolism.

Genome-wide analysis of phytochrome-interacting factor (PIF) families and their potential roles in light and gibberellin signaling in Chinese pine.

Phytochrome-interacting factors (PIFs) are a subgroup of transcription factors within the basic helix-loop-helix (bHLH) family, playing a crucial role in integrating various environmental signals to regulate plant growth and development. Despite the significance of PIFs in these processes, a comprehensive genome-wide analysis of PIFs in conifers has yet to be conducted. In this investigation, three PtPIF genes were identified in Chinese pine, categorized into three subgroups, with conserved motifs indicating the presence of the APA/APB motif and bHLH domain in the PtPIF1 and PtPIF3 proteins. Phylogenetic analysis revealed that the PtPIF1 and PtPIF3 proteins belong to the PIF7/8 and PIF3 groups, respectively, and were relatively conserved among gymnosperms. Additionally, a class of PIF lacking APA/APB motif was identified in conifers, suggesting its function may differ from that of traditional PIFs. The cis-elements of the PtPIF genes were systematically examined, and analysis of PtPIF gene expression across various tissues and under different light, temperature, and plant hormone conditions demonstrated similar expression profiles for PtPIF1 and PtPIF3. Investigations into protein-protein interactions and co-expression networks speculated the involvement of PtPIFs and PtPHYA/Bs in circadian rhythms and hormone signal transduction. Further analysis of transcriptome data and experimental validation indicated an interaction between PtPIF3 and PtPHYB1, potentially linked to diurnal rhythms. Notably, the study revealed that PtPIF3 may be involved in gibberellic acid (GA) signaling through its interaction with PtDELLAs, suggesting a potential role for PtPIF3 in mediating both light and GA responses. Overall, this research provides a foundation for future studies investigating the functions of PIFs in conifer growth and development.

Hybridization has localized effect on genetic variation in closely related pine species.

BACKGROUND: Hybridization is a known phenomenon in nature but its genetic impact on populations of parental species remains less understood. We investigated the evolutionary consequences of the interspecific gene flow in several contact zones of closely related pine species. Using a set of genetic markers from both nuclear and organellar genomes, we analyzed four hybrid zones (384 individuals) and a large panel of reference allopatric populations of parental taxa (2104 individuals from 96 stands). RESULTS: We observed reduced genetic diversity in maternally transmitted mitochondrial genomes of pure pine species and hybrids from contact zones compared to reference allopatric populations. The distribution of mtDNA haplotypes followed geographic rather than species boundaries. Additionally, no new haplotypes emerged in the contact zones, instead these zones contained the most common local variants. However, species diverged significantly at nuclear genomes and populations in contact zones exhibited similar or higher genetic diversity compared to the reference stands. There were no signs of admixture in any allopatric population, while clear admixture was evident in the contact zones, indicating that hybridization has a geographically localized effect on the genetic variation of the analyzed pine species. CONCLUSIONS: Our results suggest that hybrid zones act as sinks rather than melting pots of genetic diversity. Hybridization influences sympatric populations but is confined to contact zones. The spectrum of parental species ancestry in hybrids reflects the old evolutionary history of the sympatric populations. These findings also imply that introgression may play a crucial role in the adaptation of hybrids to specific environments.

Genetic parameter changes and age-age correlations in Pinus koraiensis growth over 40-year progeny testing.

BACKGROUND: Early selection in tree breeding could be achieved by addressing the longevity of tree improvement activities. Genetic parameter changes and age-age correlations are essential for determining the optimal timing of early selection. Practical tracking of genetic parameters of Pinus koraiensis, a major timber species with economic and ecological value, has become feasible as its progeny testing has entered the mid-term age in Korea. However, research on the age-age correlation of P. koraiensis as progeny trials approach rotation age is limited. This study aimed to investigate genetic parameter trends and age-age correlations in P. koraiensis progeny. P. koraiensis progeny were assessed at two sites using a linear mixed-effects model with two-dimensional spatial autoregressive structure. Height, diameter, and volume growth were measured in 11 assessments over 40 years. RESULTS: Genetic parameters, such as height and diameter, showed different patterns of change. The heritability ranged for the three growth traits in 0.083-0.710, 0.288-0.781, and 0.299-0.755 across the sites and age. Height heritability and its coefficient of variance decreased, whereas the diameter and volume estimates remained relatively constant. Correlations with Age 40 for phenotypic, genetic, and rank of breeding values ranged between 0.16 and 0.92, 0.594 and 0.988, and 0.412 and 0.965, respectively. These correlations generally increased as the age approached Age 40, with particularly high levels observed at Age 26 and Age 30. CONCLUSION: The observed genetic trends in P. koraiensis progeny testing offer valuable insights for early and precise selection. Notably, selecting superior genotypes at Ages 26-30 is supported by discernible genetic gains and robust correlations. Future research should integrate unbalanced data for selecting mother trees or families and conduct a comprehensive economic analysis of early selection to validate its practical benefits.

Genome-wide identification and expression analysis of ARF gene family in embryonic development of Korean pine (Pinus koraiensis).

BACKGROUND: The Auxin Responsive Factor (ARF) family plays a crucial role in mediating auxin signal transduction and is vital for plant growth and development. However, the function of ARF genes in Korean pine (Pinus koraiensis), a conifer species of significant economic value, remains unclear. RESULTS: This study utilized the whole genome of Korean pine to conduct bioinformatics analysis, resulting in the identification of 13 ARF genes. A phylogenetic analysis revealed that these 13 PkorARF genes can be classified into 4 subfamilies, indicating the presence of conserved structural characteristics within each subfamily. Protein interaction prediction indicated that Pkor01G00962.1 and Pkor07G00704.1 may have a significant role in regulating plant growth and development as core components of the PkorARFs family. Additionally, the analysis of RNA-seq and RT-qPCR expression patterns suggested that PkorARF genes play a crucial role in the development process of Korean pine. CONCLUSION: Pkor01G00962.1 and Pkor07G00704.1, which are core genes of the PkorARFs family, play a potentially crucial role in regulating the fertilization and developmental process of Korean pine. This study provides a valuable reference for investigating the molecular mechanism of embryonic development in Korean pine and establishes a foundation for cultivating high-quality Korean pine.

Pine wilt disease: what do we know from proteomics?

Pine wilt disease (PWD) is a devastating forest disease caused by the pinewood nematode (PWN), Bursaphelenchus xylophilus, a migratory endoparasite that infects several coniferous species. During the last 20 years, advances have been made for understanding the molecular bases of PWN-host trees interactions. Major advances emerged from transcriptomic and genomic studies, which revealed some unique features related to PWN pathogenicity and constituted fundamental data that allowed the development of postgenomic studies. Here we review the proteomic approaches that were applied to study PWD and integrated the current knowledge on the molecular basis of the PWN pathogenicity. Proteomics has been useful for understanding cellular activities and protein functions involved in PWN-host trees interactions, shedding light into the mechanisms associated with PWN pathogenicity and being promising tools to better clarify host trees PWN resistance/susceptibility.

Trunk distortion weakens the tree productivity revealed by half-sib progeny determination of Pinus yunnanensis.

Twisted trunks are not uncommon in trees, but their effects on tree growth are still unclear. Among coniferous tree species, the phenomenon of trunk distortion is more prominent in Pinus yunnanensis. To expand the germplasm of genetic resources, we selected families with excellent phenotypic traits to provide material for advanced generation breeding. The progeny test containing 93 superior families (3240 trees) was used as the research material. Phenotypic measurements and estimated genetic parameters (family heritability, realistic gain and genetic gain) were performed at 9, 15, and 18 years of age, respectively. The genetic evaluation yielded the following results (1) The intra-family variance component of plant height (PH) was greater than that of the inter-family, while the inter-family variance components of other traits (diameter at breast height (DBH), crown diameter (CD), height under branches (HUB), degree of stem-straightness (DS)) were greater than that of the intra-family, indicating that there was abundant variation among families and potential for selection. (2) At half rotation period (18 years old), there was a significant correlation among the traits. The proportion of trees with twisted trunks (level 1-3 straightness) reached 48%. The DS significantly affected growth traits, among which PH and DBH were the most affected. The volume loss rate caused by twisted trunk was 18.06-56.75%, implying that trunk distortion could not be completely eliminated after an artificial selection. (3) The influence of tree shape, crown width, and trunk on volume increased, and the early-late correlation between PH, DBH and volume was extremely significant. The range of phenotypic coefficient of variation, genetic variation coefficient and family heritability of growth traits (PH, DBH, and volume) were 44.29-127.13%, 22.88-60.87%, and 0.79-0.83, respectively. (4) A total of 21 superior families were selected by the method of membership function combined with independent selection. Compared with the mid-term selection (18 years old), the accuracy of early selection (9 years old) reached 77.5%. The selected families' genetic gain and realistic gain range were 5.79-19.82% and 7.12-24.27%, respectively. This study can provide some useful reference for the breeding of coniferous species.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers.

MAIN CONCLUSION: The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow.

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers with c . 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

Local adaptation of Pinus leiophylla under climate and land use change models in the Avocado Belt of Michoacán.

Climate change and land use change are two main drivers of global biodiversity decline, decreasing the genetic diversity that populations harbour and altering patterns of local adaptation. Landscape genomics allows measuring the effect of these anthropogenic disturbances on the adaptation of populations. However, both factors have rarely been considered simultaneously. Based on a set of 3660 SNPs from which 130 were identified as outliers by a genome-environment association analysis (LFMM), we modelled the spatial turnover of allele frequencies in 19 localities of Pinus leiophylla across the Avocado Belt in Michoacán state, Mexico. Then, we evaluated the effect of climate change and land use change scenarios, in addition to evaluating assisted gene flow strategies and connectivity metrics across the landscape to identify priority conservation areas for the species. We found that localities in the centre-east of the Avocado Belt would be more vulnerable to climate change, while localities in the western area are more threatened by land conversion to avocado orchards. Assisted gene flow actions could aid in mitigating both threats. Connectivity patterns among forest patches will also be modified by future habitat loss, with central and eastern parts of the Avocado Belt maintaining the highest connectivity. These results suggest that areas with the highest priority for conservation are in the eastern part of the Avocado Belt, including the Monarch Butterfly Biosphere Reserve. This work is useful as a framework that incorporates distinct layers of information to provide a more robust representation of the response of tree populations to anthropogenic disturbances.

Exploring interspecific hybridization dynamics in artificial forests of Pinus brutia and P. halepensis: Implications for sustainable afforestation.

Interspecific hybridization increases genetic diversity, which is essential for coping with changing environments. Hybrid zones, occurring naturally in overlapping habitats of closely related species, can be artificially established during afforestation. The resulting interspecific hybridization may promote sustainability in artificial forests, particularly in regions facing degradation due to climate change. Currently, there is limited evidence of hybridization during regeneration of artificial forests. Here, we studied the frequency of Pinus brutia Ten. × P. halepensis Mill. hybridization in five planted forests in Israel in three stages of forest regeneration: seeds before dispersal, emerged seedlings and recruited seedlings at the end of the dry season. We found hybrids on P. brutia, but not on P. halepensis trees due to asynchronous cone production phenology. Using 94 single-nucleotide polymorphism (SNP) markers, we found hybrids at all stages, most of which were hybrids of advanced generations. The hybrid proportions increased from 4.7 ± 2.1 to 8.2 ± 1.4 and 21.6 ± 6.4 per cent, from seeds to emerged seedlings and to recruited seedlings stages, respectively. The increased hybrid ratio implies an advantage of hybrids over P. brutia during forest regeneration. To test this hypothesis, we measured seedling growth rate and morphological traits under controlled conditions and found that the hybrid seedlings exhibited selected traits of the two parental species, which likely contributed to the fitness and survival of the hybrids during the dry season. This study highlights the potential contribution of hybrids to sustainable-planted forests and contributes to the understanding of genetic changes that occur during the regeneration of artificial forests.

Like mother like son: Transgenerational memory and cross-tolerance from drought to heat stress are identified in chloroplast proteome and seed provisioning in Pinus radiata.

How different stressors impact plant health and memory when they are imposed in different generations in wild ecosystems is still scarce. Here, we address how different environments shape heritable memory for the next generation in seeds and seedlings of Pinus radiata, a long-lived species with economic interest. The performance of the seedlings belonging to two wild clonal subpopulations (optimal fertirrigation vs. slightly stressful conditions) was tested under heat stress through physiological profiling and comparative time-series chloroplast proteomics. In addition, we explored the seeds conducting a physiological characterization and targeted transcriptomic profiling in both subpopulations. Our results showed differential responses between them, evidencing a cross-stress transgenerational memory. Seedlings belonging to the stressed subpopulation retained key proteins related to Photosystem II, chloroplast-to-nucleus signalling and osmoprotection which helped to overcome the applied heat stress. The seeds also showed a differential gene expression profile for targeted genes and microRNAs, as well as an increased content of starch and secondary metabolites, molecules which showed potential interest as biomarkers for early selection of primed plants. Thus, these finds not only delve into transgenerational cross-stress memory in trees, but also provide a new biotechnological tool for forest design.