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.

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.

Binding of PtoRAP2.12 to demethylated and accessible chromatin regions in the PtoGntK promoter stimulates growth of poplar.

DNA methylation is an essential epigenetic modification for gene regulation in plant growth and development. However, the precise mechanisms of DNA methylation remain poorly understood, especially in woody plants. We employed whole-genome bisulfite sequencing (WGBS), assays for transposase-accessible chromatin using sequencing (ATAC-seq), and RNA-Seq to investigate epigenetic regulatory relationships in Populus tomentosa treated with DNA methylation inhibitor 5-azacitidine. Expression-quantitative trait methylation analysis (eQTM), epigenome-wide association study (EWAS), and joint linkage-linkage disequilibrium mapping were used to explore the epigenetic regulatory genes, and using CRISPR/Cas9 to identify the role of candidate genes. Plant developmental abnormalities occurred when DNA methylation levels were substantially reduced. DNA methylation regulated 112 expressed genes via chromatin accessibility, of which 61 genes were significantly influenced by DNA methylation variation at the population level. One DNA methylation-regulated gene, PtoGntK, was located in a major quantitative trait locus (QTL) for poplar growth. Overexpression and CRISPR/Cas9 of PtoGntK revealed it affected poplar height and stem diameter. The PtoRAP2.12 was found to bind to the demethylated accessible region in the PtoGntK promoter, thereby promoting growth in poplar. This study identified key genes with epigenetic regulation for plant growth and provides insights into epigenetic regulation mechanisms in woody plants.

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.

Rare variations within the serine/arginine-rich splicing factor PtoRSZ21 modulate stomatal size to determine drought tolerance in Populus.

Rare variants contribute significantly to the 'missing heritability' of quantitative traits. The genome-wide characteristics of rare variants and their roles in environmental adaptation of woody plants remain unexplored. Utilizing genome-wide rare variant association study (RVAS), expression quantitative trait loci (eQTL) mapping, genetic transformation, and molecular experiments, we explored the impact of rare variants on stomatal morphology and drought adaptation in Populus. Through comparative analysis of five world-wide Populus species, we observed the influence of mutational bias and adaptive selection on the distribution of rare variants. RVAS identified 75 candidate genes correlated with stomatal size (SS)/stomatal density (SD), and a rare haplotype in the promoter of serine/arginine-rich splicing factor PtoRSZ21 emerged as the foremost association signal governing SS. As a positive regulator of drought tolerance, PtoRSZ21 can recruit the core splicing factor PtoU1-70K to regulate alternative splicing (AS) of PtoATG2b (autophagy-related 2). The rare haplotype PtoRSZ21hap2 weakens binding affinity to PtoMYB61, consequently affecting PtoRSZ21 expression and SS, ultimately resulting in differential distribution of Populus accessions in arid and humid climates. This study enhances the understanding of regulatory mechanisms that underlie AS induced by rare variants and might provide targets for drought-tolerant varieties breeding in Populus.

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.

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.

Local diversity of drought resistance and resilience in Populus tomentosa correlates with the variation of DNA methylation.

Little information is known about DNA methylation variation in shaping environment-specific drought resistance and resilience for tree adaptation. In this study, we leveraged RNA sequencing and whole-genome bisulfite sequencing data to dissect the distinction of epigenetic regulation under drought stress and rewater condition of Populus tomentosa accessions from three geographical regions. We demonstrated low resistance and high resilience for accessions from South. Non-CG methylation levels in promoter regions of Southern accessions were lower than accessions from higher latitudes and negatively regulated gene expression. CHH context methylation was more sensitive to drought stress, and the geographical-specific differentially methylated regions were scarcely changed by environmental fluctuation. We identified 60 conserved hub genes within the co-expression networks that correlate with photosynthetic and stomatal morphological traits. Epigenome-wide association studies and genome-wide association studies of these 60 hub genes revealed the interdependency between genetic and epigenetic variation in GATA9 and LECRK-VIII.2, which was associated with stomatal morphology and chlorophyll content. The natural epigenetic variation in GATA9 was also faithfully transmitted to progenies in two family-based F1 populations. This study indicates a functional relationship of DNA methylation diversity with drought resistance and resilience which offers new insights into plants' local adaptation to a stressful environment.

Leaf physiology variations are modulated by natural variations that underlie stomatal morphology in Populus.

Stomata are essential for photosynthesis and abiotic stress tolerance. Here, we used multiomics approaches to dissect the genetic architecture and adaptive mechanisms that underlie stomatal morphology in Populus tomentosa juvenile natural population (303 accessions). We detected 46 candidate genes and 15 epistatic gene-pairs, associated with 5 stomatal morphologies and 18 leaf development and photosynthesis traits, through genome-wide association studies. Expression quantitative trait locus mapping revealed that stomata-associated gene loci were significantly associated with the expression of leaf-related genes; selective sweep analysis uncovered significant differentiation in the allele frequencies of genes that underlie stomatal variations. An allelic regulatory network operating under drought stress and adequate precipitation conditions, with three key regulators (DUF538, TRA2 and AbFH2) and eight interacting genes, was identified that might regulate leaf physiology via modulation of stomatal shape and density. Validation of candidate gene variations in drought-tolerant and F1 hybrid populations of P. tomentosa showed that the DUF538, TRA2 and AbFH2 loci cause functional stabilisation of spatiotemporal regulatory, whose favourable alleles can be faithfully transmitted to offspring. This study provides insights concerning leaf physiology and stress tolerance via the regulation of stomatal determination in perennial plants.

Genome-wide association studies reveal the coordinated regulatory networks underlying photosynthesis and wood formation in Populus.

Photosynthesis and wood formation underlie the ability of trees to provide renewable resources and perform ecological functions; however, the genetic basis and regulatory pathways coordinating these two linked processes remain unclear. Here, we used a systems genetics strategy, integrating genome-wide association studies, transcriptomic analyses, and transgenic experiments, to investigate the genetic architecture of photosynthesis and wood properties among 435 unrelated individuals of Populus tomentosa, and unravel the coordinated regulatory networks resulting in two trait categories. We detected 222 significant single-nucleotide polymorphisms, annotated to 177 candidate genes, for 10 traits of photosynthesis and wood properties. Epistasis uncovered 74 epistatic interactions for phenotypes. Strikingly, we deciphered the coordinated regulation patterns of pleiotropic genes underlying phenotypic variations for two trait categories. Furthermore, expression quantitative trait nucleotide mapping and coexpression analysis were integrated to unravel the potential transcriptional regulatory networks of candidate genes coordinating photosynthesis and wood properties. Finally, heterologous expression of two pleiotropic genes, PtoMYB62 and PtoMYB80, in Arabidopsis thaliana demonstrated that they control regulatory networks balancing photosynthesis and stem secondary cell wall components, respectively. Our study provides insights into the regulatory mechanisms coordinating photosynthesis and wood formation in poplar, and should facilitate genetic breeding in trees via molecular design.

Genetic Architecture Underlying the Metabolites of Chlorogenic Acid Biosynthesis in Populus tomentosa.

Chlorogenic acid (CGA) plays a crucial role in defense response, immune regulation, and the response to abiotic stress in plants. However, the genetic regulatory network of CGA biosynthesis pathways in perennial plants remains unclear. Here, we investigated the genetic architecture for CGA biosynthesis using a metabolite-based genome-wide association study (mGWAS) and expression quantitative trait nucleotide (eQTN) mapping in a population of 300 accessions of Populus tomentosa. In total, we investigated 204 SNPs which were significantly associated with 11 metabolic traits, corresponding to 206 genes, and were mainly involved in metabolism and cell growth processes of P. tomentosa. We identified 874 eQTNs representing 1066 genes, in which the expression and interaction of causal genes affected phenotypic variation. Of these, 102 genes showed significant signatures of selection in three geographical populations, which provided insights into the adaptation of CGA biosynthesis to the local environment. Finally, we constructed a genetic network of six causal genes that coordinately regulate CGA biosynthesis, revealing the multiple regulatory patterns affecting CGA accumulation in P. tomentosa. Our study provides a multiomics strategy for understanding the genetic basis underlying the natural variation in the CGA biosynthetic metabolites of Populus, which will enhance the genetic development of abiotic-resistance varieties in forest trees.

Genetic Architecture and Genome-Wide Adaptive Signatures Underlying Stem Lenticel Traits in Populus tomentosa.

The stem lenticel is a highly specialized tissue of woody plants that has evolved to balance stem water retention and gas exchange as an adaptation to local environments. In this study, we applied genome-wide association studies and selective sweeping analysis to characterize the genetic architecture and genome-wide adaptive signatures underlying stem lenticel traits among 303 unrelated accessions of P. tomentosa, which has significant phenotypic and genetic variations according to climate region across its natural distribution. In total, we detected 108 significant single-nucleotide polymorphisms, annotated to 88 candidate genes for lenticel, of which 9 causative genes showed significantly different selection signatures among climate regions. Furthermore, PtoNAC083 and PtoMYB46 showed significant association signals and abiotic stress response, so we overexpressed these two genes in Arabidopsis thaliana and found that the number of stem cells in all three overexpression lines was significantly reduced by PtoNAC083 overexpression but slightly increased by PtoMYB46 overexpression, suggesting that both genes are involved in cell division and expansion during lenticel formation. The findings of this study demonstrate the successful application of an integrated strategy for dissecting the genetic basis and landscape genetics of complex adaptive traits, which will facilitate the molecular design of tree ideotypes that may adapt to future climate and environmental changes.

Genetic interactions among Pto-miR319 family members and their targets influence growth and wood properties in Populus tomentosa.

MicroRNAs (miRNAs) play crucial roles in all aspects of plant growth and development, but the genetic interactions of miRNAs and their target genes in woody plants are largely unknown. Here, we integrated association genetics and expression profiling to decipher the allelic variations and interactions of the Pto-MIR319 family of miRNAs and 12 putative Pto-miR319 target genes related to wood formation in 435 unrelated individuals of Populus tomentosa Carrière (Chinese white poplar). Expression pattern analysis showed that among all pairings between expressions of pre-miRNA of Pto-MIR319 members and targets, 70.0% showed negative correlation of expression levels (r = - 0.944 to 0.674, P < 0.01) in eight tissues and organs of poplar, suggesting that Pto-miR319 may participate in the regulatory network of wood formation. Single SNP-based association studies identified 137 significant associations (P < 0.01, Q < 0.1), representing 126 unique SNPs from Pto-MIR319 members and their targets, with 10 tree growth traits, revealing that these genetic factors have common roles related to wood formation. Epistasis analysis uncovered 105 significant SNP-SNP associations (P < 0.01) influencing the 10 traits, demonstrating the close genetic interactions between Pto-MIR319 family members and the 12 Pto-miR319 target genes. Notably, one common SNP, in the precursor region of Pto-MIR319e, affected the stability of Pto-MIR319e's secondary structure by altering the stem-loop structure and minimum free energy, contributing to variations in the expression of Pto-MIR319e and Pto-miR319e target genes. This study enriches the understanding of the functions of miR319 family miRNAs in poplar and exemplifies a feasible approach to exploring the genetic effects underlying miRNA-mRNA interactions related to complex traits in trees.

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus.

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein-encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue-specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis-related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker-assisted breeding efforts to genetically improve woody plants.

Linkage-linkage disequilibrium dissection of the epigenetic quantitative trait loci (epiQTLs) underlying growth and wood properties in Populus.

Increasing evidence indicates that DNA methylation is heritable and serves as an essential marker contributing to phenotypic variation. Linkage-linkage disequilibrium mapping was used to decipher the epigenetic architecture underlying nine growth and wood property traits in a linkage population (550 F1 progeny) and a natural population (435 unrelated individuals) of Populus using methylation-sensitive amplification polymorphism (MSAP)-based analysis. The interactions between genetic and epigenetic variants in the causative genes was further unveiled using expression quantitative trait methylation (eQTM) and nucleotide (eQTN) mapping strategies. A total of 163 epigenetic quantitative trait loci (epiQTLs; LOD ≥ 3.0), explaining 1.7-44.5% of phenotypic variations, were mapped to a high-resolution epigenetic map with 19 linkage groups, which was supported by the significant MSAP associations (P < 0.001) in the two populations. There were 23 causal genes involved in growth regulation and wood formation, whose markers were located in epiQTLs and associated with the same traits in both populations. Further eQTN and eQTM mapping showed that causal genetic and epigenetic variants within the 23 candidate genes may interact more in trans in gene expression and phenotype. The present study provides strategies for investigating epigenetic architecture and the interaction between genetic and epigenetic variants modulating complex traits in forest trees.

Indole-3-acetic acid has long-term effects on long non-coding RNA gene methylation and growth in Populus tomentosa.

DNA methylation and long non-coding RNAs (lncRNAs) regulate plant growth and development, but their relationship and effect on responses to the auxin phytohormone indole-3-acetic acid (IAA) remain largely unknown, particularly in woody plants such as poplar (Populus tomentosa). Following treatment of 1-year-old clonal plants with 100 µM IAA, key poplar lncRNA genes showed changes in methylation, but whole-genome methylation levels showed no significant change. Moreover, 100 µM IAA inhibited growth of the 1-year-old poplar clones, possibly through the suppression of photosynthesis. This inhibition had a long-term effect, persisting at 1 month after removal of the exogenous IAA. Transcriptome analysis identified two candidate lncRNA genes that show changes in expression following IAA treatment, TCONS_00003480 and TCONS_00004832. TCONS_00003480 contains the same microRNA target sites of ptc-miR6464 as the 4-coumarate: CoA ligase 2 transcript, which encodes a lignin biosynthesis enzyme. And TCONS_00004832 shares the same target sites of ptc-miR6437a with the Photosystem II reaction center protein D and Cytochrome C Oxidase 17 transcripts, which are related to photosynthesis. The two lncRNAs as the mimics to corresponding target genes of miRNAs to prevent them from degrading. Examination of lncRNA gene expression and methylation revealed a negative relationship (r = - 0.29, P < 0.05); moreover, hypermethylation of the two candidate lncRNA genes remained 1 month after IAA treatment, suggesting that changes in methylation might be involved in the long-term effects of plant hormones. Therefore, our study reveals a long-term effect of IAA on the growth of P. tomentosa, possibly via methylation-mediated epigenetic changes in lncRNA gene expression and the interaction with corresponding miRNAs, leading to regulation of genes related to photosynthesis and growth.

Genetic architecture underlying the lignin biosynthesis pathway involves noncoding RNAs and transcription factors for growth and wood properties in Populus.

Lignin provides structural support in perennial woody plants and is a complex phenolic polymer derived from phenylpropanoid pathway. Lignin biosynthesis is regulated by coordinated networks involving transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). However, the genetic networks underlying the lignin biosynthesis pathway for tree growth and wood properties remain unknown. Here, we used association genetics (additive, dominant and epistasis) and expression quantitative trait nucleotide (eQTN) mapping to decipher the genetic networks for tree growth and wood properties in 435 unrelated individuals of Populus tomentosa. We detected 124 significant associations (P ≤ 6.89E-05) for 10 growth and wood property traits using 30 265 single nucleotide polymorphisms from 203 lignin biosynthetic genes, 81 TF genes, 36 miRNA genes and 71 lncRNA loci, implying their common roles in wood formation. Epistasis analysis uncovered 745 significant pairwise interactions, which helped to construct proposed genetic networks of lignin biosynthesis pathway and found that these regulators might affect phenotypes by linking two lignin biosynthetic genes. eQTNs were used to interpret how causal genes contributed to phenotypes. Lastly, we investigated the possible functions of the genes encoding 4-coumarate: CoA ligase and cinnamate-4-hydroxylase in wood traits using epistasis, eQTN mapping and enzymatic activity assays. Our study provides new insights into the lignin biosynthesis pathway in poplar and enables the novel genetic factors as biomarkers for facilitating genetic improvement of trees.

Time-specific and pleiotropic quantitative trait loci coordinately modulate stem growth in Populus.

In perennial woody plants, the coordinated increase of stem height and diameter during juvenile growth improves competitiveness (i.e. access to light); however, the factors underlying variation in stem growth remain unknown in trees. Here, we used linkage-linkage disequilibrium (linkage-LD) mapping to decipher the genetic architecture underlying three growth traits during juvenile stem growth. We used two Populus populations: a linkage mapping population comprising a full-sib family of 1,200 progeny and an association mapping panel comprising 435 unrelated individuals from nearly the entire natural range of Populus tomentosa. We mapped 311 quantitative trait loci (QTL) for three growth traits at 12 timepoints to 42 regions in 17 linkage groups. Of these, 28 regions encompassing 233 QTL were annotated as 27 segmental homology regions (SHRs). Using SNPs identified by whole-genome re-sequencing of the 435-member association mapping panel, we identified significant SNPs (P ≤ 9.4 × 10-7 ) within 27 SHRs that affect stem growth at nine timepoints with diverse additive and dominance patterns, and these SNPs exhibited complex allelic epistasis over the juvenile growth period. Nineteen genes linked to potential causative alleles that have time-specific or pleiotropic effects, and mostly overlapped with significant signatures of selection within SHRs between climatic regions represented by the association mapping panel. Five genes with potential time-specific effects showed species-specific temporal expression profiles during the juvenile stages of stem growth in five representative Populus species. Our observations revealed the importance of considering temporal genetic basis of complex traits, which will facilitate the molecular design of tree ideotypes.

Genetic variants in microRNA biogenesis genes as novel indicators for secondary growth in Populus.

MicroRNAs (miRNAs) function as key regulators of complex traits, but how genetic alterations in miRNA biogenesis genes (miRBGs) affect quantitative variation has not been elucidated. We conducted transcript analyses and association genetics to investigate how miRBGs, miRNA genes (MIRNAs) and their respective targets contribute to secondary growth in a natural population of 435 Populus tomentosa individuals. This analysis identified 29 843 common single-nucleotide polymorphisms (SNPs; frequency > 0.10) within 682 genes (80 miRBGs, 152 MIRNAs, and 457 miRNA targets). Single-SNP association analysis found SNPs in 234 candidate genes exhibited significant additive/dominant effects on phenotypes. Among these, specific candidates that associated with the same traits produced 791 miRBG-MIRNA-target combinations, suggesting possible genetic miRBG-MIRNA and MIRNA-target interactions, providing an important clue for the regulatory mechanisms of miRBGs. Multi-SNP association found 4672 epistatic pairs involving 578 genes that showed significant associations with traits and identified 106 miRBG-MIRNA-target combinations. Two multi-hierarchical networks were constructed based on correlations of miRBG-miRNA and miRNA-target expression to further probe the mechanisms of trait diversity underlying changes in miRBGs. Our study opens avenues for the investigation of miRNA function in perennial plants and underscored miRBGs as potentially modulating quantitative variation in traits.

Allelic variation in PtoPsbW associated with photosynthesis, growth, and wood properties in Populus tomentosa.

Photosynthesis is one of the most important reactions on earth. PsbW, a nuclear-encoded subunit of photosystem II (PSII), stabilizes PSII structure and plays an important role in photosynthesis. Here, we used candidate gene-based linkage disequilibrium (LD) mapping to detect significant associations between allelic variations of PtoPsbW and traits related to photosynthesis, growth, and wood properties in Populus tomentosa. PtoPsbW showed the highest expression in leaves and it increased during the development of these leaves, suggesting that PtoPsbW may play an important role in plant growth and development. Analysis of nucleotide diversity and LD revealed that PtoPsbW has low single-nucleotide polymorphism (SNP) diversity (π tot = 0.0048 and θ w = 0.0050) and relatively low average value of LD (0.1500), indicating that PtoPsbW is conserved due to its indispensable function. Using single-SNP associations in an association population of 435 individuals, we identified five significant associations at the threshold of P ≤ 0.05, explaining 3.28-15.98 % of the phenotypic variation. Haplotype-based association analyses indicated that 13 haplotypes (P ≤ 0.05) from six blocks were associated with photosynthesis, growth, and wood properties. Our work shows that identifying allelic variation and LD can help to decipher the genetic basis of photosynthesis and could potentially be applied for molecular marker-assisted selection in Populus.