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.

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.

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.

Identification of key gene networks controlling organic acid and sugar metabolism during star fruit (Averrhoa carambola) development.

The sugar and organic acid content significantly impacts the flavor quality of star fruit, and it undergoes dynamic changes during development. However, the metabolic network and molecular mechanisms governing the formation of sugar and organic acid in star fruit remain unclear. In this study, 23 of 743 components were detected by metabonomic analysis. The highest metabolites contents were organic acids and derivatives. The highest sugar content in the fruit was fructose and glucose, followed by sucrose, which proved that A. carambola is a hexose accumulation type fruit. Genome identification preliminarily screened 141 genes related to glucose metabolism and 67 genes related to acid metabolism. A total of 7,881 unigenes were found in transcriptome data, 6,124 differentially expressed genes were screened, with more up-regulated than down-regulated genes. Transcriptome and metabolome association analysis screened seven core candidate genes related to glucose metabolism and 17 core genes highly related to organic acid pathway, and eight differentially expressed sugar and acid genes were selected for qRT-PCR verification. In addition, 29 bHLHs and eight bZIPs transcription factors were predicted in the glucose metabolism pathway, and 23 MYBs, nine C2H2s transcription factors and one GRAS transcription factor was predicted in the acid metabolism pathway, and transcription factors have both positive and negative regulatory effects on sugar and acid structure genes. This study increased our understanding of A. carambola fruit flavor and provided basic information for further exploring the ornamental and edible values of star fruit.

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.

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.

Enhanced genome-wide association reveals the role of YABBY11-NGATHA-LIKE1 in leaf serration development of Populus.

Leaf margins are complex plant morphological features that contribute to leaf shape diversity, which affects plant structure, yield, and adaptation. Although several leaf margin regulators have been identified to date, the genetic basis of their natural variation has not been fully elucidated. In this study, we profiled two distinct leaf morphology types (serrated and smooth) using the persistent homology mathematical framework (PHMF) in two poplar species (Populus tomentosa and Populus simonii, respectively). A combined genome-wide association study (GWAS) and expression quantitative trait nucleotide (eQTN) mapping were applied to create a leaf morphology control module using data from P. tomentosa and P. simonii populations. Natural variation in leaf margins was associated with YABBY11 (YAB11) transcript abundance in poplar. In P. tomentosa, PtoYAB11 carries a premature stop codon (PtoYAB11PSC), resulting in the loss of its positive regulation of NGATHA-LIKE1 (PtoNGAL-1) and RIBULOSE BISPHOSPHATE CARBOXYLASE LARGE SUBUNIT (PtoRBCL). Overexpression of PtoYAB11PSC promoted serrated leaf margins, enlarged leaves, enhanced photosynthesis, and increased biomass. Overexpression of PsiYAB11 in P. tomentosa promoted smooth leaf margins, higher stomatal density, and greater light damage repair ability. In poplar, YAB11-NGAL1 is sensitive to environmental conditions, acts as a positive regulator of leaf margin serration, and may also link environmental signaling to leaf morphological plasticity.

Growth-regulating factor 15-mediated gene regulatory network enhances salt tolerance in poplar.

Soil salinity is an important determinant of crop productivity and triggers salt stress response pathways in plants. The salt stress response is controlled by transcriptional regulatory networks that maintain regulatory homeostasis through combinations of transcription factor (TF)-DNA and TF-TF interactions. We investigated the transcriptome of poplar 84 K (Populus alba × Populus glandulosa) under salt stress using samples collected at 4- or 6-h intervals within 2 days of salt stress treatment. We detected 24,973 differentially expressed genes, including 2,231 TFs that might be responsive to salt stress. To explore these interactions and targets of TFs in perennial woody plants, we combined gene regulatory networks, DNA affinity purification sequencing, yeast two-hybrid-sequencing, and multi-gene association approaches. Growth-regulating factor 15 (PagGRF15) and its target, high-affinity K+ transporter 6 (PagHAK6), were identified as an important regulatory module in the salt stress response. Overexpression of PagGRF15 and PagHAK6 in transgenic lines improved salt tolerance by enhancing Na+ transport and modulating H2O2 accumulation in poplar. Yeast two-hybrid assays identified more than 420 PagGRF15-interacting proteins, including ETHYLENE RESPONSE FACTOR TFs and a zinc finger protein (C2H2) that are produced in response to a variety of phytohormones and environmental signals and are likely involved in abiotic stress. Therefore, our findings demonstrate that PagGRF15 is a multifunctional TF involved in growth, development, and salt stress tolerance, highlighting the capability of a multifaceted approach in identifying regulatory nodes in plants.

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.

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.

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.

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.

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.

Inbreeding in Chinese fir: Insights into the adaptive growth traits of selfed progeny from mRNA, miRNA, and copy number variation.

PREMISE: The impact of inbreeding on biological processes is well documented in individuals with severe inbreeding depression. However, the biological processes influencing the adaptive growth of normal selfed individuals are unknown. Here, we aimed to investigate how inbreeding affects gene expression for adaptive growth of normal selfed seedlings from a self-fertilizing parent in Chinese fir (Cunninghamia lanceolata). METHODS: Using RNA-seq data from above- and underground tissues of abnormal and normal selfed seedlings, we analyzed GO biological processes network. We also sequenced small RNAs in the aboveground tissues and measured the copy number variations (CNV) of the hub genes. RESULTS: Phenotypic fitness analysis revealed that the normal seedlings were better adapted than their abnormal counterparts. Upregulated differentially expressed genes (DEGs) were associated with development processes, and downregulated DEGs were mainly enriched in fundamental metabolism and stress response. Results of mRNA-miRNA parallel sequencing revealed that upregulated target genes were predominantly associated with development, highlighting their crucial role in phosphorylation in signal transduction networks. We also discovered a moderate correlation (0.1328 < R2 < 0.6257) between CNV and gene expression levels for three hub genes (TMKL1, GT2, and RHY1A). CONCLUSIONS: We uncovered the key biological processes underpinning the growth of normal selfed seedlings and established the relationship between CNV and the expression levels of hub genes in selfed seedlings. Understanding the candidate genes involved in the growth of selfed seedlings will help us comprehend the genetic mechanisms behind inbreeding depression in the evolutionary biology of plants.

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.

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.

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.

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".

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.

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata.

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.