High-yield hybrid breeding of Camellia oleifolia based on ISSR molecular markers.

BACKGROUND: C. Oleifera is among the world's largest four woody plants known for their edible oil production, yet the contribution rate of improved varieties is less than 20%. The species traditional breeding is lengthy cycle (20-30 years), occupation of land resources, high labor cost, and low accuracy and efficiency, which can be enhanced by molecular marker-assisted selection. However, the lack of high-quality molecular markers hinders the species genetic analysis and molecular breeding. RESULTS: Through quantitative traits characterization, genetic diversity assessment, and association studies, we generated a selection population with wide genetic diversity, and identified five excellent high-yield parental combinations associated with four reliable high-yield ISSR markers. Early selection criteria were determined based on kernel fresh weight and cultivated 1-year seedling height, aided by the identification of these 4 ISSR markers. Specific assignment of selected individuals as paternal and maternal parents was made to capitalize on their unique attributes. CONCLUSIONS: Our results indicated that molecular markers-assisted breeding can effectively shorten, enhance selection accuracy and efficiency and facilitate the development of a new breeding system for C. oleifera.

Regulatory networks of senescence-associated gene-transcription factors promote degradation in Moso bamboo shoots.

Bamboo cultivation, particularly Moso bamboo (Phyllostachys edulis), holds significant economic importance in various regions worldwide. Bamboo shoot degradation (BSD) severely affects productivity and economic viability. However, despite its agricultural consequences, the molecular mechanisms underlying BSD remain unclear. Consequently, we explored the dynamic changes of BSD through anatomy, physiology and the transcriptome. Our findings reveal ruptured protoxylem cells, reduced cell wall thickness and the accumulation of sucrose and reactive oxygen species (ROS) during BSD. Transcriptomic analysis underscored the importance of genes related to plant hormone signal transduction, sugar metabolism and ROS homoeostasis in this process. Furthermore, BSD appears to be driven by the coexpression regulatory network of senescence-associated gene transcription factors (SAG-TFs), specifically PeSAG39, PeWRKY22 and PeWRKY75, primarily located in the protoxylem of vascular bundles. Yeast one-hybrid and dual-luciferase assays demonstrated that PeWRKY22 and PeWRKY75 activate PeSAG39 expression by binding to its promoter. This study advanced our understanding of the molecular regulatory mechanisms governing BSD, offering a valuable reference for enhancing Moso bamboo forest productivity.

Phylogenetic conservatism and coordination in traits of Chinese woody endemic flora.

Range-limited endemic species, often labeled as endangered due to their low adaptability to climate change, exhibit unclear evolutionary mechanisms influencing their distribution. This study explores the relationship between leaf length, maximum height, and seed diameter and their linkage to phylogeny and climate in the macroecology of 1,370 woody endemics. Using Bayesian analytical method that allows partitioning phylogenetic and environmental variances and covariance, we revealed moderate to high phylogenetic signals in these traits, indicating evolutionary constraints potentially impacting climate change adaptability. The study uncovered a phylogenetically conserved coordination between height and leaf length which showed to be independent of macroecological patterns of temperature and precipitation. These findings emphasize the role of phylogenetic ancestry in shaping the distribution of woody endemics, highlighting the need for prioritized in-situ conservation and providing insights for ex situ conservation strategies.

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.

Spatial Distribution and Ecological Determinants of Coexisting Hybrid Oak Species: A Study in Yushan's Mixed Forest.

Ecological niche partitioning is crucial in reducing interspecific competition, fostering species coexistence, and preserving biodiversity. Our research, conducted in a hybrid mixed oak forest in Yushan, Jiangsu, China, focuses on Quercus acutissima, Q. variabilis, Q. fabri, and Q. serrata var. brevipetiolata. Using Point Pattern Analysis, we investigated the spatial relationships and ecological trait autocorrelation, including total carbon (TC), nitrogen (TN), phosphorus (TP), potassium (TK), and breast height diameter (DBH). Our findings show aggregated distribution patterns within the oak populations. The Inhomogeneous Poisson Point model highlights the impact of environmental heterogeneity on Q. variabilis, leading to distinct distribution patterns, while other species showed wider dispersion. This study reveals aggregated interspecific interactions, with a notable dispersal pattern between Q. acutissima and Q. variabilis. We observed significant variability in nutrient elements, indicating distinct nutrient dynamics and uptake processes. The variations in total carbon (TC), nitrogen (TN), phosphorus (TP), and potassium (TK) suggest distinct nutrient dynamics, with TK showing the highest variability. Despite variations in TC, TK, and TP, the species did not form distinct classes, suggesting overlapping nutritional strategies and environmental adaptations. Furthermore, spatial autocorrelation analysis indicates strong positive correlations for DBH, TC, and TP, whereas TK and TN correlations are non-significant. The results suggest habitat filtering as a key driver in intraspecific relationships, with a finer spatial scale of ecological niche division through TC and TP, which is crucial for maintaining coexistence among these oak species.

Integration of Physiological, Transcriptomic, and Metabolomic Analyses Reveal Molecular Mechanisms of Salt Stress in Maclura tricuspidata.

Salt stress is a universal abiotic stress that severely affects plant growth and development. Understanding the mechanisms of Maclura tricuspidate's adaptation to salt stress is crucial for developing salt-tolerant plant varieties. This article discusses the integration of physiology, transcriptome, and metabolome to investigate the mechanism of salt adaptation in M. tricuspidata under salt stress conditions. Overall, the antioxidant enzyme system (SOD and POD) of M. tricuspidata exhibited higher activities compared with the control, while the content of soluble sugar and concentrations of chlorophyll a and b were maintained during salt stress. KEGG analysis revealed that deferentially expressed genes were primarily involved in plant hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, alkaloids, and MAPK signaling pathways. Differential metabolites were enriched in amino acid metabolism, the biosynthesis of plant hormones, butanoate, and 2-oxocarboxylic acid metabolism. Interestingly, glycine, serine, and threonine metabolism were found to be important both in the metabolome and transcriptome-metabolome correlation analyses, suggesting their essential role in enhancing the salt tolerance of M. tricuspidata. Collectively, our study not only revealed the molecular mechanism of salt tolerance in M. tricuspidata, but also provided a new perspective for future salt-tolerant breeding and improvement in salt land for this species.

Temporal dynamics of genetic architecture governing leaf development in Populus.

Leaf development is a multifaceted and dynamic process orchestrated by a myriad of genes to shape the proper size and morphology. The dynamic genetic network underlying leaf development remains largely unknown. Utilizing a synergistic genetic approach encompassing dynamic genome-wide association study (GWAS), time-ordered gene co-expression network (TO-GCN) analyses and gene manipulation, we explored the temporal genetic architecture and regulatory network governing leaf development in Populus. We identified 42 time-specific and 18 consecutive genes that displayed different patterns of expression at various time points. We then constructed eight TO-GCNs that covered the cell proliferation, transition, and cell expansion stages of leaf development. Integrating GWAS and TO-GCN, we postulated the functions of 27 causative genes for GWAS and identified PtoGRF9 as a key player in leaf development. Genetic manipulation via overexpression and suppression of PtoGRF9 revealed its primary influence on leaf development by modulating cell proliferation. Furthermore, we elucidated that PtoGRF9 governs leaf development by activating PtoHB21 during the cell proliferation stage and attenuating PtoLD during the transition stage. Our study provides insights into the dynamic genetic underpinnings of leaf development and understanding the regulatory mechanism of PtoGRF9 in this dynamic process.

N6-methyladenosine mRNA methylation positively regulated the response of poplar to salt stress.

As the most abundant form of methylation modification in messenger RNA (mRNA), the distribution of N6-methyladenosine (m6A) has been preliminarily revealed in herbaceous plants under salt stress, but its function and mechanism in woody plants were still unknown. Here, we showed that global m6A levels increased during poplar response to salt stress. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) revealed that m6A significantly enriched in the coding sequence region and 3'-untranslated regions in poplar, by recognising the conserved motifs, AGACU, GGACA and UGUAG. A large number of differential m6A transcripts have been identified, and some have been proved involving in salt response and plant growth and development. Further combined analysis of MeRIP-seq and RNA-seq revealed that the m6A hypermethylated and enrich in the CDS region preferred to positively regulate expression abundance. Writer inhibitor, 3-deazaneplanocin A treatment increased the sensitivity of poplar to salt stress by reducing mRNA stability to regulate the expression of salt-responsive transcripts PagMYB48, PagGT2, PagNAC2, PagGPX8 and PagARF2. Furthermore, we verified that the methyltransferase PagFIP37 plays a positively role in the response of poplar to salt stress, overexpressed lines have stronger salt tolerance, while RNAi lines were more sensitive to salt, which relied on regulating mRNA stability in an m6A manner of salt-responsive transcripts PagMYB48, PagGT2, PagNAC2, PagGPX8 and PagARF2. Collectively, these results revealed the regulatory role of m6A methylation in poplar response to salt stress, and revealed the importance and mechanism of m6A methylation in the response of woody plants to salt stress for the first time.

Transcriptomic and Phenotypic Analyses Reveal the Molecular Mechanism of Dwarfing in Tetraploid Robinia pseudoacacia L.

Polyploid breeding techniques aid in the cultivation of new forestry cultivars, thus expanding the suite of strategies for the improvement of arboreal traits and innovation within the field of forestry. Compared to diploid Robinia pseudoacacia L. (black locust) 'D26-5①' (2×), its dwarfed homologous tetraploid 'D26-5②' (4×) variety has better application prospects in garden vegetation guardrails and urban landscape. However, the molecular mechanism of the generation and growth of this dwarf variety is still unclear. Here, plant growth and development as well as histological differences between the diploid and its autotetraploid were investigated. Levels of endogenous hormones at three different developmental stages (20, 40, and 70 days) of 2× and homologous 4× tissue culture plantlets were assessed, and it was found that the brassinosteroid (BR) contents of the former were significantly higher than the latter. Transcriptome sequencing data analysis of 2× and homologous 4× showed that differentially expressed genes (DEGs) were significantly enriched in plant hormone synthesis and signal transduction, sugar and starch metabolism, and the plant circadian rhythm pathway, which are closely related to plant growth and development. Therefore, these biological pathways may be important regulatory pathways leading to dwarfism and slow growth in tetraploids. Additionally, utilizing weighted gene coexpression network analysis (WGCNA), we identified three crucial differentially expressed genes (DEGs)-PRR5, CYP450, and SPA1-that potentially underlie the observed ploidy variation. This study provides a new reference for the molecular mechanism of dwarfism in dwarfed autotetraploid black locusts. Collectively, our results of metabolite analysis and comparative transcriptomics confirm that plant hormone signaling and the circadian rhythm pathway result in dwarfism in black locusts.

Efficient genomics-based 'end-to-end' selective tree breeding framework.

Since their initiation in the 1950s, worldwide selective tree breeding programs followed the recurrent selection scheme of repeated cycles of selection, breeding (mating), and testing phases and essentially remained unchanged to accelerate this process or address environmental contingencies and concerns. Here, we introduce an "end-to-end" selective tree breeding framework that: (1) leverages strategically preselected GWAS-based sequence data capturing trait architecture information, (2) generates unprecedented resolution of genealogical relationships among tested individuals, and (3) leads to the elimination of the breeding phase through the utilization of readily available wind-pollinated (OP) families. Individuals' breeding values generated from multi-trait multi-site analysis were also used in an optimum contribution selection protocol to effectively manage genetic gain/co-ancestry trade-offs and traits' correlated response to selection. The proof-of-concept study involved a 40-year-old spruce OP testing population growing on three sites in British Columbia, Canada, clearly demonstrating our method's superiority in capturing most of the available genetic gains in a substantially reduced timeline relative to the traditional approach. The proposed framework is expected to increase the efficiency of existing selective breeding programs, accelerate the start of new programs for ecologically and environmentally important tree species, and address climate-change caused biotic and abiotic stress concerns more effectively.

Genome-wide identification of the Pinus tabuliformis CONSTANS-like gene family and their potential roles in reproductive cone development.

The CONSTANS-like (COL) genes, as a core transcription factor in the photoperiod regulation pathway, play a key role in plant reproduction development. However, their molecular characterization has rarely been studied in Pinus tabuliformis. Here, 10 PtCOL genes were identified in the P. tabuliformis genome and multiple sequence alignments have indicated that the PtCOL proteins contained highly conserved B-BOX1 and CCT domains. Sequence similarity analysis showed that PtCOL1 and PtCOL3 had the higher similarity with Norway spruce COLs (PaCOL2 and PaCOL1) and Arabidopsis COLs (AtCOL3, 4 and 5), respectively. Phylogeny and gene structure analyses revealed that PtCOLs were divided into three subgroups, each with identical or similar distributions of exons, introns, and motifs. Moreover, 10 PtCOLs were distributed on 6 chromosomes and PtCOL9 has syntenic gene pairs in both Ginkgo biloba and Sequoiadendron giganteum. Interestingly, in transcriptome profiles, most PtCOLs exhibited a diurnal oscillation pattern under both long (LD) and short (SD) day conditions. Additionally, PtCOLs were highly expressed in needles and female cones, and showed different spatial expression patterns. Among the ten PtCOLs, PtCOL1/3 heterologous overexpression Arabidopsis displayed a delayed-flowering phenotype under SD, indicating that they are likely to play a crucial role in the reproductive development. Additionally, PtCOL1 and PtCOL3 were not only capable of interacting with each other, but they were each capable of interacting with themselves. Furthermore, PtCOL1 and PtCOL3 were also involved in the MADS-box protein-protein interaction (PPI) network in P. tabuliformis cone development. Direct interactions of PtDAL11 with PtCOL1/3 impeded PtCOL1/3 translocation into the nucleus. In summary, this study provided comprehensive understanding for the functions of the PtCOL gene family and revealed their biological roles in the photoperiod-dependent P. tabuliformis cone development.

Allelic variations of WAK106-E2Fa-DPb1-UGT74E2 module regulate fibre properties in Populus tomentosa.

Wood formation, intricately linked to the carbohydrate metabolism pathway, underpins the capacity of trees to produce renewable resources and offer vital ecosystem services. Despite their importance, the genetic regulatory mechanisms governing wood fibre properties in woody plants remain enigmatic. In this study, we identified a pivotal module comprising 158 high-priority core genes implicated in wood formation, drawing upon tissue-specific gene expression profiles from 22 Populus samples. Initially, we conducted a module-based association study in a natural population of 435 Populus tomentosa, pinpointing PtoDPb1 as the key gene contributing to wood formation through the carbohydrate metabolic pathway. Overexpressing PtoDPb1 led to a 52.91% surge in cellulose content, a reduction of 14.34% in fibre length, and an increment of 38.21% in fibre width in transgenic poplar. Moreover, by integrating co-expression patterns, RNA-sequencing analysis, and expression quantitative trait nucleotide (eQTN) mapping, we identified a PtoDPb1-mediated genetic module of PtoWAK106-PtoDPb1-PtoE2Fa-PtoUGT74E2 responsible for fibre properties in Populus. Additionally, we discovered the two PtoDPb1 haplotypes that influenced protein interaction efficiency between PtoE2Fa-PtoDPb1 and PtoDPb1-PtoWAK106, respectively. The transcriptional activation activity of the PtoE2Fa-PtoDPb1 haplotype-1 complex on the promoter of PtoUGT74E2 surpassed that of the PtoE2Fa-PtoDPb1 haplotype-2 complex. Taken together, our findings provide novel insights into the regulatory mechanisms of fibre properties in Populus, orchestrated by PtoDPb1, and offer a practical module for expediting genetic breeding in woody plants via molecular design.

Inferring the Regulatory Network of miRNAs on Terpene Trilactone Biosynthesis Affected by Environmental Conditions.

As a medicinal tree species, ginkgo (Ginkgo biloba L.) and terpene trilactones (TTLs) extracted from its leaves are the main pharmacologic activity constituents and important economic indicators of its value. The accumulation of TTLs is known to be affected by environmental stress, while the regulatory mechanism of environmental response mediated by microRNAs (miRNAs) at the post-transcriptional levels remains unclear. Here, we focused on grafted ginkgo grown in northwestern, southwestern, and eastern-central China and integrally analyzed RNA-seq and small RNA-seq high-throughput sequencing data as well as metabolomics data from leaf samples of ginkgo clones grown in natural environments. The content of bilobalide was highest among detected TTLs, and there was more than a twofold variation in the accumulation of bilobalide between growth conditions. Meanwhile, transcriptome analysis found significant differences in the expression of 19 TTL-related genes among ginkgo leaves from different environments. Small RNA sequencing and analysis showed that 62 of the 521 miRNAs identified were differentially expressed among different samples, especially the expression of miRN50, miR169h/i, and miR169e was susceptible to environmental changes. Further, we found that transcription factors (ERF, MYB, C3H, HD-ZIP, HSF, and NAC) and miRNAs (miR319e/f, miRN2, miRN54, miR157, miR185, and miRN188) could activate or inhibit the expression of TTL-related genes to participate in the regulation of terpene trilactones biosynthesis in ginkgo leaves by weighted gene co-regulatory network analysis. Our findings provide new insights into the understanding of the regulatory mechanism of TTL biosynthesis but also lay the foundation for ginkgo leaves' medicinal value improvement under global change.

Characterization of pollen tube development in distant hybridization of Chinese cork oak (Quercus variabilis L.).

MAIN CONCLUSION: This work mainly found that the stigma and style of Q. variabilis did not completely lose the specific recognition towards heterologous pollen, a fact which is different from previous studies. Quercus is the foundation species in the Northern Hemisphere, with extreme prevalence for interspecific hybridization. It is not yet entirely understood whether or how the pollen tube-female tissue interaction contributes to the "extensive hybridization" in oaks. Pollen storage conditions correlate with distant hybridization. We conducted hybridization experiments with Q. variabilis as female and Q. variabilis and Q. mongolica as male parents. And the differences in pollen tube (PT) development between intra- and distant interspecific hybridization were studied by fluorescence microscopy and scanning electron microscopy (SEM). Our results showed that -20 °C allowed pollen of both species to maintain some viability. Both Q. variabilis and Q. mongolica pollen germinated profusely on the stigmas. SEM results indicated that in the intraspecific hybridization, Q. variabilis pollen started to germinate at 6 h after pollination (hap), PTs elongated significantly at 12 hap, and entered the stigma at 24 hap. By contrast, Q. mongolica pollen germinated at 15 hap, and the PTs entered the stigma at 27 hap. By fluorescence microscopical studies it was observed that some PTs of Q. variabilis gathered at the style-joining at 96 hap, unlike the Q. mongolica which reached the style junction at 144 hap. The above results indicate that the abundant germination of heterologous pollen (HP) on the stigma and the "Feeble specificity recognition" of the stigma and transmitting tract to HP may create opportunities for the "extensive hybridization" of oaks. This work provides a sexual developmental reference for clarifying the causes of Quercus "extensive hybridization".

Role of soil nutrient elements transport on Camellia oleifera yield under different soil types.

BACKGROUND: Most of Camellia oleifera forests have low fruit yield and poor oil quality that are largely associated with soil fertility. Soil physical and chemical properties interact with each other affecting soil fertility and C. oleifera growing under different soil conditions produced different yield and oil composition. Three main soil types were studied, and redundancy, correlation, and double-screening stepwise regression analysis were used for exploring the relationships between C. oleifera nutrients uptake and soil physical and chemical properties, shedding light on the transport law of nutrient elements from root, leaves, and kernel, and affecting the regulation of fruit yield and oil composition. RESULTS: In the present study, available soil elements content of C. oleifera forest were mainly regulated by water content, pH value, and total N, P and Fe contents. Seven elements (N, P, K, Mg, Cu, Mn and C) were key for kernel's growth and development, with N, P, K, Cu and Mn contents determining 74.0% the yield traits. The transport characteristics of these nutrients from root, leaves to the kernel had synergistic and antagonistic effects. Increasing oil production and unsaturated fatty acid content can be accomplished in two ways: one through increasing N, P, Mg, and Zn contents of leaves by applying corresponding N, P, Mg, Zn foliar fertilizers, while the other through maintaining proper soil moisture content by applying Zn fertilizer in the surface layer and Mg and Ca fertilizer in deep gully. CONCLUSION: Soil type controlled nutrient absorption by soil pH, water content and total N, P and Fe content. There were synergistic and antagonistic effects on the inter-organ transport of nutrient elements, ultimately affecting N, P, K, Cu and Mn contents in kernel, which determined the yield and oil composition of C. oleifera.

Genome-wide identification of late embryogenesis abundant protein family and their key regulatory network in Pinus tabuliformis cold acclimation.

Cold acclimation is a crucial biological process that enables conifers to overwinter safely. The late embryogenesis abundant (LEA) protein family plays a pivotal role in enhancing freezing tolerance during this process. Despite its importance, the identification, molecular functions and regulatory networks of the LEA protein family have not been extensively studied in conifers or gymnosperms. Pinus tabuliformis, a conifer with high ecological and economic values and with high-quality genome sequence, is an ideal candidate for such studies. Here, a total of 104 LEA genes were identified from P. tabuliformis, and we renamed them according to their subfamily group: PtLEA1-PtLEA92 (group LEA1-LEA6), PtSMP1-PtSMP6 (group seed maturation protein) and PtDHN1-PtDHN6 (group Dehydrin). While the sequence structure of P. tabuliformis  LEA genes are conserved, their physicochemical properties exhibit unique characteristics within different subfamily groupings. Notably, the abundance of low-temperature responsive elements in PtLEA genes was observed. Using annual rhythm and temperature gradient transcriptome data, PtLEA22 was identified as a key gene that responds to low-temperature induction while conforming to the annual cycle of cold acclimation. Overexpression of PtLEA22 enhanced Arabidopsis freezing tolerance. Furthermore, several transcription factors potentially co-expressed with PtLEA22 were validated using yeast one-hybrid and dual-luciferase assays, revealing that PtDREB1 could directly bind PtLEA22 promoter to positively regulate its expression. These findings reveal the genome-wide characterization of P. tabuliformis  LEA genes and their importance in the cold acclimation, while providing a theoretical basis for studying the molecular mechanisms of cold acclimation in conifers.

Natural variation in the prolyl 4-hydroxylase gene PtoP4H9 contributes to perennial stem growth in Populus.

Perennial trees must maintain stem growth throughout their entire lifespan to progressively increase in size as they age. The overarching question of the molecular mechanisms that govern stem perennial growth in trees remains largely unanswered. Here we deciphered the genetic architecture that underlies perennial growth trajectories using genome-wide association studies (GWAS) for measures of growth traits across years in a natural population of Populus tomentosa. By analyzing the stem growth trajectory, we identified PtoP4H9, encoding prolyl 4-hydroxylase 9, which is responsible for the natural variation in the growth rate of diameter at breast height (DBH) across years. Quantifying the dynamic genetic contribution of PtoP4H9 loci to stem growth showed that PtoP4H9 played a pivotal role in stem growth regulation. Spatiotemporal expression analysis showed that PtoP4H9 was highly expressed in cambium tissues of poplars of various ages. Overexpression and knockdown of PtoP4H9 revealed that it altered cell expansion to regulate cell wall modification and mechanical characteristics, thereby promoting stem growth in Populus. We showed that natural variation in PtoP4H9 occurred in a BASIC PENTACYSTEINE transcription factor PtoBPC1-binding promoter element controlling PtoP4H9 expression. The geographic distribution of PtoP4H9 allelic variation was consistent with the modes of selection among populations. Altogether, our study provides important genetic insights into dynamic stem growth in Populus, and we confirmed PtoP4H9 as a potential useful marker for breeding or genetic engineering of poplars.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa.

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes.

BACKGROUND: Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species. RESULTS: Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification. CONCLUSION: We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Allelic variation in transcription factor PtoWRKY68 contributes to drought tolerance in Populus.

Drought stress limits woody species productivity and influences tree distribution. However, dissecting the molecular mechanisms that underpin drought responses in forest trees can be challenging due to trait complexity. Here, using a panel of 300 Chinese white poplar (Populus tomentosa) accessions collected from different geographical climatic regions in China, we performed a genome-wide association study (GWAS) on seven drought-related traits and identified PtoWRKY68 as a candidate gene involved in the response to drought stress. A 12-bp insertion and/or deletion and three nonsynonymous variants in the PtoWRKY68 coding sequence categorized natural populations of P. tomentosa into two haplotype groups, PtoWRKY68hap1 and PtoWRKY68hap2. The allelic variation in these two PtoWRKY68 haplotypes conferred differential transcriptional regulatory activities and binding to the promoters of downstream abscisic acid (ABA) efflux and signaling genes. Overexpression of PtoWRKY68hap1 and PtoWRKY68hap2 in Arabidopsis (Arabidopsis thaliana) ameliorated the drought tolerance of two transgenic lines and increased ABA content by 42.7% and 14.3% compared to wild-type plants, respectively. Notably, PtoWRKY68hap1 (associated with drought tolerance) is ubiquitous in accessions in water-deficient environments, whereas the drought-sensitive allele PtoWRKY68hap2 is widely distributed in well-watered regions, consistent with the trends in local precipitation, suggesting that these alleles correspond to geographical adaptation in Populus. Moreover, quantitative trait loci analysis and an electrophoretic mobility shift assay showed that SHORT VEGETATIVE PHASE (PtoSVP.3) positively regulates the expression of PtoWRKY68 under drought stress. We propose a drought tolerance regulatory module in which PtoWRKY68 modulates ABA signaling and accumulation, providing insight into the genetic basis of drought tolerance in trees. Our findings will facilitate molecular breeding to improve the drought tolerance of forest trees.