Transcriptome sequencing of active buds from Populus deltoides CL. and Populus × zhaiguanheibaiyang reveals phytohormones involved in branching

Genome-wide identification and analysis of GH1-containing H1 histones among poplar species

Histone H1s are basic nuclear proteins, which played key role in the binding of DNA and nucleosome, eventually the stability of eukaryotic chromatin. In most species, H1s possess an evolutionarily conserved nucleosome-DNA binding globular domain (GH1), which is conserved between species, especially in mammals. However, there is limited information on the phylogeny, structure and function of H1s in poplar. In the present research, 21 GH1-containing proteins found in Populus trichocarpa were classified into three subgroups (H1s, Myb (SANK) GH1 and AT-hook GH1) based on their domains. The Populus H1 proteins contained lysine-rich N-, C-terminal tails and a conserved GH1 domain, particularly the characteristic amino acids in the helix and strand structures of the five H1 subtypes. The phylogenetic and structure diversity analysis of GH1 proteins across different Populus species and model plants revealed three conserved subgroups with characteristic amino acids. The variation in the number of members across the five subtypes was consistent with the evolutionary relationships among Populus species. The conserved characteristic amino acids among same Populus subtype can be served as markers for subtype identification. Furthermore, the abundance analysis of H1s in Populus indicated their unique functions in young tissues and stages, which may be related to DNA methylation. The consistent expression pattern of H1 across Populus species was in accordance with collinearity pairs. Present analyses provided valuable information on the diversity and evolution of H1s in Populus, advocating further research of H1s in plants.

Deletion of the PtrDJ1C gene leads to increased branching in poplar

The PtrDJ1C gene is essential for poplar growth and early chloroplast development. Disruption of PtrDJ1C expression results in an albino leaf phenotype and increased branching. However, the underlying mechanism for the increased branching remains unknown. In this study, we employed integrated approaches to investigate the function of PtrDJ1C in the branch-increasing phenotype. Our results revealed that levels of indole-3-acetic acid (IAA), gibberellin (GA), and abscisic acid (ABA) were significantly reduced in ptrdj1c mutants, while cytokinin (CK) levels were slightly increased. Transcriptomic and proteomic analyses identified several key genes and proteins involved in hormone regulation and branching development that were differentially expressed. Specifically, the expression levels of TAA, ZEP, and GA20ox-genes involved in IAA, GA, and ABA biosynthesis-were significantly reduced in ptrdj1c, while IPT and LOG, which regulate CK synthesis, were upregulated. Moreover, immunoblot analysis further validated reduced levels of key biosynthetic enzymes for IAA, GA, and ABA, alongside increased levels of IPT and LOG enzymes. Interestingly, our findings suggest that hormone signaling pathways act in concert with the transcription factor WUSCHEL (WUS) to synergistically promote branching development. These results provide novel insight into the regulatory role of PtrDJ1C in hormone balance and its downstream effects on poplar branching.

A protocol for identifying universal reference genes within a genus based on RNA-Seq data: a case study of poplar stem gene expression

Real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) plays a crucial role in relative gene expression analysis, and accurate normalization relies on suitable reference genes (RGs). In this study, a pipeline for identifying candidate RGs from publicly available stem-related RNA-Seq data of different Populus species under various developmental and abiotic stress conditions is presented. DESeq2's median of ratios yielded the smallest coefficient of variance (CV) values in a total of 292 RNA-Seq samples and was therefore chosen as the method for sample normalization. A total of 541 stably expressed genes were retrieved based on the CV values with a cutoff of 0.3. Universal gene-specific primer pairs were designed based on the consensus sequences of the orthologous genes of each Populus RG candidate. The expression levels of 12 candidate RGs and six reported RGs in stems under different abiotic stress conditions or in different Populus species were assessed by RT-qPCR. The expression stability of selected genes was further evaluated using ΔCt, geNorm, NormFinder, and BestKeeper. All candidate RGs were stably expressed in different experiments and conditions in Populus. A test dataset containing 117 RNA-Seq samples was then used to confirm the expression stability, six candidate RGs and three reported RGs met the requirement of CV ≤ 0.3. In summary, this study was to propose a systematic and optimized protocol for the identification of constitutively and stably expressed genes based on RNA-Seq data, and Potri.001G349400 (CNOT2) was identified as the best candidate RG suitable for gene expression studies in poplar stems.

Integrated transcriptomic and gibberellin analyses reveal genes related to branch development in Eucalyptus urophylla

Tree branches affect the planting density and basal scab, which act as important attributes in the yield and quality of trees. Eucalyptus urophylla is an important pioneer tree with characteristics of strong adaptability, fast growth, short rotation period, and low disease and pest pressures. In this study, we collected ZQUC14 and LDUD26 clones and compared their transcriptomes and metabolomes from mature xylem, phloem, and developing tissues to identify factors that may influence branch development. In total, 32,809 differentially expressed genes (DEGs) and 18 gibberellin (GA) hormones were detected in the five sampled tissues. Searches of the kyoto Encyclopedia of Genes and Genomes pathways identified mainly genes related to diterpenoid biosynthesis, plant MAPK signaling pathways, plant hormone signal transduction, glycerolipid metabolism, peroxisome, phenylpropanoid biosynthesis, ABC transporters, and brassinosteroid biosynthesis. Furthermore, gene expression trend analysis and weighted gene co-expression network analysis revealed 13 genes likely involved in diterpenoid biosynthesis, including five members of the 2OG-Fe(II) oxygenase superfamily, four cytochrome P450 genes, and four novel genes. In GA signal transduction pathways, 24 DEGs were found to positively regulate branch formation. These results provide a comprehensive analysis of branch development based on the transcriptome and metabolome, and help clarify the molecular mechanisms of E. urophylla.

Roles of auxin in the inhibition of shoot branching in 'Dugan' fir

Shoot branching substantially impacts vegetative and reproductive growth as well as wood characteristics in perennial woody species by shaping the shoot system architecture. Although plant hormones have been shown to play a fundamental role in shoot branching in annual species, their corresponding actions in perennial woody plants are largely unknown, in part due to the lack of branching mutants. Here, we demonstrated the role of plant hormones in bud dormancy transition toward activation and outgrowth in woody plants by comparing the physiological and molecular changes in the apical shoot stems of 'Yangkou' 020 fir and 'Dugan' fir, two Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.) clones with normal and completely abolished branching phenotypes, respectively. Our studies showed that the defect in bud outgrowth was the cause of failed shoot branching in 'Dugan' fir whereas apically derived signals acted as triggers of this ectopic bud activity. Further studies indicated that auxin played a key role in inhibiting bud outgrowth in 'Dugan' fir. During bud dormancy release, the differential auxin resistant 1/Like AUX1 (AUX1/LAX) and PIN-formed (PIN) activity resulted in an ectopic auxin/indole-3-acetic acid (IAA) accumulation in the apical shoot stem of 'Dugan' fir, which could inhibit the cell cycle in the axillary meristem by decreasing cytokinin (CK) biosynthesis but increasing abscisic acid (ABA) production and response through the signaling pathway. In contrast, during bud activation and outgrowth, the striking increase in auxin biosynthesis and PIN activity in the shoot tip of 'Dugan' fir may trigger the correlative inhibition of axillary buds by modulating the polar auxin transport stream (PATS) and connective auxin transport (CAT) in shoots, and by influencing the biosynthesis of secondary messengers, including CK, gibberellin (GA) and ABA, thereby inducing the paradormancy of axillary buds in 'Dugan' fir by apical dominance under favorable conditions. The findings of this study provide important insights into the roles of plant hormones in bud outgrowth control in perennial woody plants.

Transcriptome and metabolite analysis related to branch development in two genotypes of Eucalyptus urophylla

Branching in long-lived plants can cause scarring at the base and affect wood density, which greatly inhibits wood yield and quality. Eucalyptus urophylla is one of the most important commercial forest tree species in South China, with diverse branch number and branch angles under different genetic backgrounds. However, the main elements and regulatory mechanisms associated with different branching traits in E. urophylla remain unclear. To identify the factors that may influence branching, the transcriptome and metabolome were performed on the shoot apex (SA), lateral shoot apex (LSA), and stem segment at the 5th axillary bud from the shoot apex (S1) in lines ZQUC14 (A) and LDUD26 (B), with A exhibiting a smaller Ba than B. A total of 307.3 million high-quality clean reads and nine hormones were identified from six libraries. Several differentially expressed regulatory factors were identified between the two genotypes of E. urophylla. The Kyoto Encyclopedia of Genes and Genomes pathways were enriched in plant hormone signal transduction, plant hormone biosynthesis and their transport pathways. Furthermore, gene expression pattern analysis identified genes that were significantly downregulated or upregulated in S1 relative to the SA and LSA segments, and the plant hormone signal transduction pathway was constructed to explain branching development. This study clarified the main plant hormones and genes underlying branch numbers and angles of E. urophylla, confirmed that ABA and SA could promote a larger branch angle and smaller branch number, while IAA has an opposite function. Numbers of key candidate genes involved in plant hormone signal transduction were found in the positive regulation of branch formation. These novel findings should aid molecular breeding of branching in Eucalyptus.

Molecular mechanism of lateral bud differentiation of Pinus massoniana based on high-throughput sequencing

Knot-free timber cultivation is an important goal of forest breeding, and lateral shoots affect yield and stem shape of tree. The purpose of this study was to analyze the molecular mechanism of lateral bud development by removing the apical dominance of Pinus massoniana young seedlings through transcriptome sequencing and identify key genes involved in lateral bud development. We analyzed hormone contents and transcriptome data for removal of apical dominant of lateral buds as well as apical and lateral buds of normal development ones. Data were analyzed using an comprehensive approach of pathway- and gene-set enrichment analysis, Mapman visualization tool, and gene expression analysis. Our results showed that the contents of auxin (IAA), Zea and strigolactone (SL) in lateral buds significantly increased after removal of apical dominance, while abscisic acid (ABA) decreased. Gibberellin (GA) metabolism, cytokinin (CK), jasmonic acid, zeatin pathway-related genes positively regulated lateral bud development, ABA metabolism-related genes basically negatively regulated lateral bud differentiation, auxin, ethylene, SLs were positive and negative regulation, while only A small number of genes of SA and BRASSINOSTEROID, such as TGA and TCH4, were involved in lateral bud development. In addition, it was speculated that transcription factors such as WRKY, TCP, MYB, HSP, AuxIAA, and AP2 played important roles in the development of lateral buds. In summary, our results provided a better understanding of lateral bud differentiation and lateral shoot formation of P. massoniana from transcriptome level. It provided a basis for molecular characteristics of side branch formation of other timber forests, and contributed to knot-free breeding of forest trees.