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.

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.

Evidence that Ophiostomatoid Fungal Symbionts of Mountain Pine Beetle Do Not Play a Role in Overcoming Lodgepole Pine Defenses During Mass Attack.

Mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins) is a devastating forest insect pest that has killed millions of hectares of pines in western North America over the past two decades. Like other bark beetles, MPB vectors ophiostomatoid fungal species, some of which are pathogenic to host pine species. The phytopathogenicity of these fungal symbionts has sparked considerable debate regarding their role in facilitating MPB attack success. We tested the hypothesis that MPB ophiostomatoid fungal associates like Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield contribute to overwhelming host defenses during MPB mass attack. We compared responses of mature lodgepole pine (Pinus contorta Dougl. ex Loud. var. latifolia Engelm.) trees growing in natural stands that were mass attacked by MPB with those inoculated with G. clavigera by examining host defense hormones, secondary metabolites, and gene expression profiles. The jasmonate and ethylene signatures of necrotrophic pathogen-triggered response were identified in G. clavigera-inoculated trees, but only the jasmonate signature of a herbivore-triggered response was measured in MPB-attacked trees. Several G. clavigera-induced changes in pine phenolic metabolite profiles and phenolic biosynthesis gene expression patterns were absent in MPB-attacked pines. These findings indicate that ophiostomatoid fungi like G. clavigera are not a major factor in overwhelming host defenses during MPB mass attack. Instead, fungal pathogenicity likely is more important in aiding MPB colonization and development within the host tree. Phenolics appear to play a larger role in the host response to G. clavigera than to MPB, although phenolics may also influence MPB feeding and behavior. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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.

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.

Strong Purifying Selection in Haploid Tissue-Specific Genes of Scots Pine Supports the Masking Theory.

The masking theory states that genes expressed in a haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in a diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intraspecific competition. Evidence from other plant groups is scarce, and results show no support for the masking theory. Here, we have used a gymnosperm Scots pine megagametophyte, a maternally derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue and that this signal of strong selection is not an artifact driven by high expression levels.

A delayed response in phytohormone signaling and production contributes to pine susceptibility to Fusarium circinatum.

BACKGROUND: Fusarium circinatum is the causal agent of pine pitch canker disease, which affects Pinus species worldwide, causing significant economic and ecological losses. In Spain, two Pinus species are most affected by the pathogen; Pinus radiata is highly susceptible, while Pinus pinaster has shown moderate resistance. In F. circinatum-Pinus interactions, phytohormones are known to play a crucial role in plant defense. By comparing species with different degrees of susceptibility, we aimed to elucidate the fundamental mechanisms underlying resistance to the pathogen. For this purpose, we used an integrative approach by combining gene expression and metabolomic phytohormone analyses at 5 and 10 days post inoculation. RESULTS: Gene expression and metabolite phytohormone contents suggested that the moderate resistance of P. pinaster to F. circinatum is determined by the induction of phytohormone signaling and hormone rearrangement beginning at 5 dpi, when symptoms are still not visible. Jasmonic acid was the hormone that showed the greatest increase by 5 dpi, together with the active gibberellic acid 4 and the cytokinin dehydrozeatin; there was also an increase in abscisic acid and salicylic acid by 10 dpi. In contrast, P. radiata hormonal changes were delayed until 10 dpi, when symptoms were already visible; however, this increase was not as high as that in P. pinaster. Indeed, in P. radiata, no differences in jasmonic acid or salicylic acid production were found. Gene expression analysis supported the hormonal data, since the activation of genes related to phytohormone synthesis was observed earlier in P. pinaster than in the susceptible P. radiata. CONCLUSIONS: We determine that the moderate resistance of P. pinaster to F. circinatum is in part a result of early and strong activation of plant phytohormone-based defense responses before symptoms become visible. We suggest that jasmonic acid signaling and production are strongly associated with F. circinatum resistance. In contrast, P. radiata susceptibility was attributed to a delayed response to the fungus at the moment when symptoms were visible. Our results contribute to a better understanding of the phytohormone-based defense mechanism involved in the Pinus-F. circinatum interactions and provide insight into the development of new strategies for disease mitigation.

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.

Adaptive potential of maritime pine under contrasting environments.

BACKGROUND: Predicting the adaptability of forest tree populations under future climates requires a better knowledge of both the adaptive significance and evolvability of measurable key traits. Phenotypic plasticity, standing genetic variation and degree of phenotypic integration shape the actual and future population genetic structure, but empirical estimations in forest tree species are still extremely scarce. We analysed 11 maritime pine populations covering the distribution range of the species (119 families and 8 trees/family, ca. 1300 trees) in a common garden experiment planted at two sites with contrasting productivity. We used plant height as a surrogate of fitness and measured five traits (mean and plasticity of carbon isotope discrimination, specific leaf area, needle biomass, Phenology growth index) related to four different strategies (acquisitive economics, photosynthetic organ size, growth allocation and avoidance of water stress). RESULTS: Estimated values of additive genetic variation would allow adaptation of the populations to future environmental conditions. Overall phenotypic integration and selection gradients were higher at the high productivity site, while phenotypic integration within populations was higher at the low productivity site. Response to selection was related mainly to photosynthetic organ size and drought-avoidance mechanisms rather than to water use efficiency. Phenotypic plasticity of water use efficiency could be maladaptive, resulting from selection for height growth. CONCLUSIONS: Contrary to the expectations in a drought tolerant species, our study suggests that variation in traits related to photosynthetic organ size and acquisitive investment of resources drive phenotypic selection across and within maritime pine populations. Both genetic variation and evolvability of key adaptive traits were considerably high, including plasticity of water use efficiency. These characteristics would enable a relatively fast micro-evolution of populations in response to the ongoing climate changes. Moreover, differentiation among populations in the studied traits would increase under the expected more productive future Atlantic conditions.

Transcriptome profiles reveal response mechanisms and key role of PsNAC1 in Pinus sylvestris var. mongolica to drought stress.

BACKGROUND: Drought stress severely impedes plant growth, and only a limited number of species exhibit long-term resistance to such conditions. Pinus sylvestris var. mongolica, a dominant tree species in arid and semi-arid regions of China, exhibits strong drought resistance and plays a crucial role in the local ecosystem. However, the molecular mechanisms underlying this resistance remain poorly understood. RESULTS: Here, we conducted transcriptome sequence and physiological indicators analysis of needle samples during drought treatment and rehydration stages. De-novo assembly yielded approximately 114,152 unigenes with an N50 length of 1,363 bp. We identified 6,506 differentially expressed genes (DEGs), with the majority being concentrated in the heavy drought stage (4,529 DEGs). Functional annotation revealed enrichment of drought-related GO terms such as response to water (GO:0009415: enriched 108 genes) and response to water deprivation (GO:0009414: enriched 106 genes), as well as KEGG categories including MAPK signaling pathway (K04733: enriched 35 genes) and monoterpenoid biosynthesis (K21374: enriched 27 genes). Multiple transcription factor families and functional protein families were differentially expressed during drought treatment. Co-expression network analysis identified a potential drought regulatory network between cytochrome P450 genes (Unigene4122_c1_g1) and a core regulatory transcription factor Unigene9098_c3_g1 (PsNAC1) with highly significant expression differences. We validated PsNAC1 overexpression in Arabidopsis and demonstrated enhanced drought resistance. CONCLUSIONS: These findings provide insight into the molecular basis of drought resistance in P. sylvestris var. mongolica and lay the foundation for further exploration of its regulatory network.

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.

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.

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.

Scots pine - panmixia and the elusive signal of genetic adaptation.

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Physiological trait coordination and variability across and within three Pinus species.

Studies have explored how traits separate plants ecologically and the trade-offs that underpin this separation. However, uncertainty remains as to the taxonomic scale at which traits can predictably separate species. We studied how physiological traits separated three Pinus (Pinus banksiana, Pinus resinosa, and Pinus strobus) species across three sites. We collected traits from four common leaf and branch measurements (light-response curves, CO2-response curves, pressure-volume curves, and hydraulic vulnerability curves) across each species and site. While common, these measurements are not typically measured together due to logistical constraints. Few traits varied across species and sites as expected given the ecological preferences of the species and environmental site characteristics. Some trait trade-offs present at broad taxonomic scales were observed across the three species, but most were absent within species. Certain trade-offs contrasted expectations observed at broader scales but followed expectations given the species' ecological preferences. We emphasize the need to both clarify why certain traits are being studied, as variation in unexpected but ecologically meaningful ways often occurs and certain traits might not vary substantially within a given lineage (e.g. hydraulic vulnerability in Pinus), highlighting the role a trait selection in trait ecology.

Local carbon reserves are insufficient for phloem terpene induction during drought in Pinus edulis in response to bark beetle-associated fungi.

Stomatal closure during drought inhibits carbon uptake and may reduce a tree's defensive capacity. Limited carbon availability during drought may increase a tree's mortality risk, particularly if drought constrains trees' capacity to rapidly produce defenses during biotic attack. We parameterized a new model of conifer defense using physiological data on carbon reserves and chemical defenses before and after a simulated bark beetle attack in mature Pinus edulis under experimental drought. Attack was simulated using inoculations with a consistent bluestain fungus (Ophiostoma sp.) of Ips confusus, the main bark beetle colonizing this tree, to induce a defensive response. Trees with more carbon reserves produced more defenses but measured phloem carbon reserves only accounted for c. 23% of the induced defensive response. Our model predicted universal mortality if local reserves alone supported defense production, suggesting substantial remobilization and transport of stored resin or carbon reserves to the inoculation site. Our results show that de novo terpene synthesis represents only a fraction of the total measured phloem terpenes in P. edulis following fungal inoculation. Without direct attribution of phloem terpene concentrations to available carbon, many studies may be overestimating the scale and importance of de novo terpene synthesis in a tree's induced defense response.

A pioneer nematode effector suppresses plant reactive oxygen species burst by interacting with the class III peroxidase.

Bursaphelenchus xylophilus is the pathogen of pine wilt disease, which can devastate the pine forest ecosystem. Usually, plant cells generate reactive oxygen species (ROS) as a defensive substance or signalling molecules to resist the infection of nematodes. However, little is known about how B. xylophilus effectors mediate the plant ROS metabolism. Here, we identified a pioneer B. xylophilus Prx3-interacting effector 1 (BxPIE1) expressed in the dorsal gland cells and the intestine. Silencing of the BxPIE1 gene resulted in reduced nematode reproduction and a delay in disease progression during parasitic stages, with the upregulation of pathogenesis-related (PR) genes PtPR-3 (class Ⅳ chitinase) and PtPR-9 (peroxidase). The protein-protein interaction assays further demonstrated that BxPIE1 interacts with a Pinus thunbergii class III peroxidase (PtPrx3), which produces H2O2 under biotic stress. The expression of BxPIE1 and PtPrx3 was upregulated during the infection stage. Furthermore, BxPIE1 effectively inhibited H2O2 generating from class III peroxidase and ascorbate can recover the virulence of siBxPIE1-treated B. xylophilus by scavenging H2O2. Taken together, BxPIE1 is an important virulence factor, revealing a novel mechanism utilized by nematodes to suppress plant immunity.

Deciphering the Role of SVP-Like Genes and Their Key Regulation Networks During Reproductive Cone Development in Pinus tabuliformis.

Reproductive development plays an essential role in the perpetuation of genetic material and environmental adaptation. In angiosperms, the Short Vegetative Phase (SVP) serves as a flowering repressor, influencing the development of floral organs. In this study, heterologous transformation of Arabidopsis thaliana with SVP-like genes (PtSVL1 and PtSVL2) derived from Pinus tabuliformis significantly impacted stamen formation and pollen fertility, without altering flowering time. Gene co-expression networks revealed that SVP-like and SOC1-like genes function as key coregulatory transcription factors during the initial stages of cone development in P. tabuliformis. Interestingly, the regulatory module of SOC1 regulated by SVP in angiosperms is absent in conifers and conifer SVP-like exercises its function in a form that is physically bound to SOC1-like. Furthermore, combining the yeast one-hybrid scanning with co-expression network analysis, revealed that SPLs and TPSs were the principal downstream target genes of PtSVL1. Notably, the PtSPL16 promoter is positively regulated by PtSVL1, and overexpression of PtSPL16 results in delayed flowering in Arabidopsis, suggesting that the PtSVL1-PtSPL16 module plays a crucial role in regulating reproductive development in conifers. Collectively, these findings enhance our understanding of the roles of SVP-like genes in conifers and the key regulatory networks centred on PtSVL1 during reproductive cone development.