Characterization of a High Hierarchical Regulator, PtrGATA12, Functioning in Differentially Regulating Secondary Wall Component Biosynthesis in Populus trichocarpa

Overexpression of the PtrNF-YA6 gene inhibits secondary cell wall thickening in poplar

The NF-Y gene family in plants plays a crucial role in numerous biological processes, encompassing hormone response, stress response, as well as growth and development. In this study, we first used bioinformatics techniques to identify members of the NF-YA family that may function in wood formation. We then used molecular biology techniques to investigate the role and molecular mechanism of PtrNF-YA6 in secondary cell wall (SCW) formation in Populus trichocarpa. We found that PtrNF-YA6 protein was localized in the nucleus and had no transcriptional activating activity. Overexpression of PtrNF-YA6 had an inhibitory effect on plant growth and development and significantly suppressed hemicellulose synthesis and SCW thickening in transgenic plants. Yeast one-hybrid and ChIP-PCR assays revealed that PtrNF-YA6 directly regulated the expression of hemicellulose synthesis genes (PtrGT47A-1, PtrGT8C, PtrGT8F, PtrGT43B, PtrGT47C, PtrGT8A and PtrGT8B). In conclusion, PtrNF-YA6 can inhibit plant hemicellulose synthesis and SCW thickening by regulating the expression of downstream SCW formation-related target genes.

GATA family transcription factors in alga Chlamydomonas reinhardtii

GATA family transcription factors (GATA-TFs) are metalloproteins that regulate many metabolic pathways. These conserved proteins recognize the consensus sequence (A/T)GATA(A/G) in the promoter regions of many genes and regulate their transcription in response to environmental signals. Currently, the study of GATA-TFs is of increasing interest. GATA genes and their proteins are most actively studied in vascular plants and fungi. Based on the results of numerous studies, it has been shown that GATA factors regulate the metabolic pathways of nitrogen and carbon, and also play a major role in the processes induced by light and circadian rhythms. In algae, GATA-TFs remain poorly studied, and information about them is scattered. In this work, all known data on GATA-TFs in the unicellular green alga Chlamydomonas reinhardtii has been collected and systematized. The genome of this alga contains 12 GATA coding genes. Using the phylogenetic analysis, we identified three classes of GATA factors in C. reinhardtii according to the structure of the zinc finger domain and showed their difference from the classification of GATA factors developed on vascular plants.

Spatiotemporal miRNA and transcriptomic network dynamically regulate the developmental and senescence processes of poplar leaves

Poplar is an important afforestation and urban greening species. Poplar leaf development occurs in stages, from young to mature and then from mature to senescent; these are accompanied by various phenotypic and physiological changes. However, the associated transcriptional regulatory network is relatively unexplored. We first used principal component analysis to classify poplar leaves at different leaf positions into two stages: developmental maturity (the stage of maximum photosynthetic capacity); and the stage when photosynthetic capacity started to decline and gradually changed to senescence. The two stages were then further subdivided into five intervals by gene expression clustering analysis: young leaves, the period of cell genesis and functional differentiation (L1); young leaves, the period of development and initial formation of photosynthetic capacity (L3-L7); the period of maximum photosynthetic capacity of functional leaves (L9-L13); the period of decreasing photosynthetic capacity of functional leaves (L15-L27); and the period of senescent leaves (L29). Using a weighted co-expression gene network analysis of regulatory genes, high-resolution spatiotemporal transcriptional regulatory networks were constructed to reveal the core regulators that regulate leaf development. Spatiotemporal transcriptome data of poplar leaves revealed dynamic changes in genes and miRNAs during leaf development and identified several core regulators of leaf development, such as GRF5 and MYB5. This in-depth analysis of transcriptional regulation during leaf development provides a theoretical basis for exploring the biological basis of the transcriptional regulation of leaf development and the molecular design of breeding for delaying leaf senescence.

Comprehensive Analysis and Characterization of the GATA Gene Family, with Emphasis on the GATA6 Transcription Factor in Poplar

GATA transcription factors are ubiquitously present in eukaryotic organisms and play a crucial role in multiple biological processes, such as plant growth, stress response, and hormone signaling. However, the study of GATA factors in poplar is currently limited to a small number of proteins, despite their evident functional importance. In this investigation, we utilized the most recent genome annotation and stringent criteria to identify 38 GATA transcription factor genes in poplar. Subsequently, we conducted a comprehensive analysis of this gene family, encompassing phylogenetic classification, protein characterization, analysis of promoter cis-acting elements, and determination of chromosomal location. Our examination of gene duplication events indicated that both tandem and segmental duplications have contributed to the expansion of the GATA gene family in poplar, with segmental duplication potentially being a major driving force. By performing collinearity analysis of genes across six different species, we identified 74 pairs of co-linear genes, which provide valuable insights for predicting gene functions from a comparative genomics perspective. Furthermore, through the analysis of gene expression patterns, we identified five GATA genes that exhibited differential expression in leaf-stem-root tissues and eight genes that were responsive to salt stress. Of particular interest was GATA6, which displayed strong induction by salt stress and overlapped between the two gene sets. We discovered that GATA6 encodes a nuclear-localized protein with transcription activation activity, which is continuously induced by salt stress in leaf and root tissues. Moreover, we constructed a co-expression network centered around GATA6, suggesting the potential involvement of these genes in the growth, development, and response to abiotic stress processes in poplar through cell transport systems and protein modification mechanisms, such as vesicle-mediated transport, intracellular transport, ubiquitination, and deubiquitination. This research provides a foundation for further exploration of the functions and mechanisms of GATA transcription factors in poplar.

Construction of a genetic map and QTL mapping of seed size traits in soybean

Soybean seed size and seed shape traits are closely related to plant yield and appearance quality. In this study, 186 individual plants of the F2 generation derived from crosses between Changjiang Chun 2 and JiYu 166 were selected as the mapping population to construct a molecular genetic linkage map, and the phenotypic data of hundred-grain weight, seed length, seed width, and seed length-to-width ratio of soybean under three generations of F2 single plants and F2:3 and F2:4 lines were combined to detect the QTL (quantitative trait loci) for the corresponding traits by ICIM mapping. A soybean genetic map containing 455 markers with an average distance of 6.15 cM and a total length of 2799.2 cM was obtained. Forty-nine QTLs related to the hundred-grain weight, seed length, seed width, and seed length-to-width ratio of soybean were obtained under three environmental conditions. A total of 10 QTLs were detected in more than two environments with a phenotypic variation of over 10%. Twelve QTL clusters were identified on chromosomes 1, 2, 5, 6, 8, 13, 18, and 19, with the majority of the overlapping intervals for hundred-grain weight and seed width. These results will lay the theoretical and technical foundation for molecularly assisted breeding in soybean seed weight and seed shape. Eighteen candidate genes that may be involved in the regulation of soybean seed size were screened by gene functional annotation and GO enrichment analysis.

The putative vacuolar processing enzyme gene TaVPE3cB is a candidate gene for wheat stem pith-thickness

KEY MESSAGE

The vacuolar processing enzyme gene TaVPE3cB is identified as a candidate gene for a QTL of wheat pith-thickness on chromosome 3B by BSR-seq and differential expression analyses. The high pith-thickness (PT) of the wheat stem could greatly enhance stem mechanical strength, especially the basal internodes which support the heavier upper part, such as upper stems, leaves and spikes. A QTL for PT in wheat was previously discovered on 3BL in a double haploid population of 'Westonia' × 'Kauz'. Here, a bulked segregant RNA-seq analysis was applied to identify candidate genes and develop associated SNP markers for PT. In this study, we aimed at screening differentially expressed genes (DEGs) and SNPs in the 3BL QTL interval. Sixteen DEGs were obtained based on BSR-seq and differential expression analyses. Twenty-four high-probability SNPs in eight genes were identified by comparing the allelic polymorphism in mRNA sequences between the high PT and low PT samples. Among them, six genes were confirmed to be associated with PT by qRT-PCR and sequencing. A putative vacuolar processing enzyme gene TaVPE3cB was screened out as a potential PT candidate gene in Australian wheat 'Westonia'. A robust SNP marker associated with TaVPE3cB was developed, which can assist in the introgression of TaVPE3cB.b in wheat breeding programs. In addition, we also discussed the function of other DEGs which may be related to pith development and programmed cell death (PCD). A five-level hierarchical regulation mechanism of stem pith PCD in wheat was proposed.

PtrWOX13A Promotes Wood Formation and Bioactive Gibberellins Biosynthesis in Populus trichocarpa

WUSCHEL-related homeobox (WOX) genes are plant-specific transcription factors (TFs) involved in multiple processes of plant development. However, there have hitherto no studies on the WOX TFs involved in secondary cell wall (SCW) formation been reported. In this study, we identified a Populus trichocarpa WOX gene, PtrWOX13A, which was predominantly expressed in SCW, and then characterized its functions through generating PtrWOX13A overexpression poplar transgenic lines; these lines exhibited not only significantly enhanced growth potential, but also remarkably increased SCW thicknesses, fiber lengths, and lignin and hemicellulose contents. However, no obvious change in cellulose content was observed. We revealed that PtrWOX13A directly activated its target genes through binding to two cis-elements, ATTGATTG and TTAATSS, in their promoter regions. The fact that PtrWOX13A responded to the exogenous GAs implies that it is responsive to GA homeostasis caused by GA inactivation and activation genes (e.g., PtrGA20ox4, PtrGA2ox1, and PtrGA3ox1), which were regulated by PtrWOX13A directly or indirectly. Since the master switch gene of SCW formation, PtrWND6A, and lignin biosynthesis regulator, MYB28, significantly increased in PtrWOX13A transgenic lines, we proposed that PtrWOX13A, as a higher hierarchy TF, participated in SCW formation through controlling the genes that are components of the known hierarchical transcription regulation network of poplar SCW formation, and simultaneously triggering a gibberellin-mediated signaling cascade. The discovery of PtrWOX13A predominantly expressed in SCW and its regulatory functions in the poplar wood formation has important implications for improving the wood quality of trees via genetic engineering.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analyses Provide Comprehensive Insight Into Molecular Mechanisms of Cellulose and Lignin Biosynthesis in Cunninghamia lanceolata

Cunninghamia lanceolata is an essential timber species that provide 20%-30% raw materials for China's timber industry. Although a few transcriptomes have been published in C. lanceolata, full-length mRNA transcripts and regulatory mechanisms behind the cellulose and lignin biosynthesis have not been thoroughly investigated. Here, PacBio Iso-seq and RNA-seq analyses were adapted to identify the full-length and differentially expressed transcripts along a developmental gradient from apex to base of C. lanceolata shoots. A total of 48,846 high-quality full-length transcripts were obtained, of which 88.0% are completed transcriptome based on benchmarking universal single-copy orthologs (BUSCO) assessment. Along stem developmental gradient, 18,714 differentially expressed genes (DEGs) were detected. Further, 28 and 125 DEGs were identified as enzyme-coding genes of cellulose and lignin biosynthesis, respectively. Moreover, 57 transcription factors (TFs), including MYB and NAC, were identified to be involved in the regulatory network of cellulose and lignin biosynthesis through weighted gene co-expression network analysis (WGCNA). These TFs are composed of a comparable regulatory network of secondary cell wall formation in angiosperms, revealing a similar mechanism may exist in gymnosperms. Further, through qRT-PCR, we also investigated eight specific TFs involved in compression wood formation. Our findings provide a comprehensive and valuable source for molecular genetics breeding of C. lanceolata and will be beneficial for molecular-assisted selection.

Genome-Wide Identification of the Eucalyptus urophylla GATA Gene Family and Its Diverse Roles in Chlorophyll Biosynthesis

GATA transcription factors have been demonstrated to play key regulatory roles in plant growth, development, and hormonal response. However, the knowledge concerning the evolution of GATA genes in Eucalyptus urophylla and their trans-regulatory interaction is indistinct. Phylogenetic analysis and study of conserved motifs, exon structures, and expression patterns resolved the evolutionary relationships of these GATA proteins. Phylogenetic analysis showed that EgrGATAs are broadly distributed in four subfamilies. Cis-element analysis of promoters revealed that EgrGATA genes respond to light and are influenced by multiple hormones and abiotic stresses. Transcriptome analysis revealed distinct temporal and spatial expression patterns of EgrGATA genes in various tissues of E. urophylla S.T.Blake, which was confirmed by real-time quantitative PCR (RT-qPCR). Further research revealed that EurGNC and EurCGA1 were localized in the nucleus, and EurGNC directly binds to the cis-element of the EurGUN5 promoter, implying its potential roles in the regulation of chlorophyll synthesis. This comprehensive study provides new insights into the evolution of GATAs and could help to improve the photosynthetic assimilation and vegetative growth of E. urophylla at the genetic level.

Comprehensive Transcriptome Analysis of Stem-Differentiating Xylem Upon Compression Stress in Cunninghamia Lanceolata

Compression wood (CW) in gymnosperm brings great difficulties to wood industry using wood as raw materials since CW presents special wood structure and have different physical and chemical properties from those of normal wood (NW). Chinese fir (Cunninghamia lanceolata) is widely distributed in China. However, global transcriptome profiling of coding and long non-coding RNA in response to compression stress has not been reported in the gymnosperm species. In this study, we revealed that CW in Chinese fir exhibited distinct morphology and cytology properties compared with those of NW, including high lignin content, thick and round tracheid cells. Furthermore, we combined both PacBio long-read SMRT sequencing (Iso-Seq) and Illumina short-read RNA-Seq to reveal the transcriptome in stem-differentiating xylem (SDX) under different time points (2, 26, and 74 h) upon compression stress in NW, CW, and OW (opposite wood), respectively. Iso-Seq was successfully assembled into 41,253 de-novo full-length transcriptome reference (average length 2,245 bp). Moreover, there were striking differences in expression upon compression stress, which were involved 13 and 7 key enzyme genes in the lignin and cellulose synthesis, respectively. Especially, we revealed 11 secondary growth-related transcription factors show differential expression under compression stress, which was further validated by qRT-PCR. Finally, the correlation between 6,533 differentially expressed coding genes and 372 differentially expressed long non-coding RNAs (lncRNAs) indicates that these lncRNAs may affect cell wall biogenesis and xyloglucan metabolism. In conclusion, our results provided comprehensive cytology properties and full-length transcriptome profiling of wood species upon compression stress. Especially we explored candidate genes, including both coding and long non-coding genes, and provided a theoretical basis for further research on the formation mechanism of CW in gymnosperm Chinese fir.