Identification of Homeobox Genes Associated with Lignification and Their Expression Patterns in Bamboo Shoots

Impacts of in-vitro zebularine treatment on genome-wide DNA methylation and transcriptomic profiles in Salix purpurea L

Renewable energy resources such as biomass are crucial for a sustainable global society. Trees are a major source of lignocellulosic biomass, which can vary in response to different environmental factors owing to epigenetic regulation, such as DNA C-methylation. To investigate the effects of DNA methylation on plant development and wood formation, and its impacts on gene expression, with a focus on secondary cell wall (SCW)-associated genes, Salix purpurea plantlets were cloned from buds derived from a single hybrid tree for both treatment and control conditions. For the treatment condition, buds were exposed to 50 μM zebularine in vitro and a combined strategy of whole-genome bisulfite sequencing (WGBS) and RNA-seq was employed to examine the methylome and transcriptome profiles of different tissues collected at various time points under both conditions. Transcriptomic and methylome data revealed that most of the promoter and gene body demethylation had no marked effects on the expression profiles of genes. Nevertheless, gene expression tended to decrease with the increased methylation levels of genes with highly methylated promoters. Results indicated that demethylation is less evident in centromeric regions and sex chromosomes. Promoters of secondary cell wall-associated genes, such as 4-coumarate-CoA ligase-like and Rac-like GTP-binding protein RHO, were differentially methylated in the secondary xylem samples collected from two-month potted treated plants compared to control samples. Our results provide novel insights into DNA methylation and gene expression landscapes and a basis for investigating the epigenetic regulation of wood formation in S. purpurea as a model plant for bioenergy species.

Expression Profiling and MicroRNA Regulatory Networks of Homeobox Family Genes in Sugarcane Saccharum spontaneum L

Homeobox (HB) genes play important roles in plant growth and development processes, particularly in the formation of lateral organs. Thus, they could influence leaf morphogenesis and biomass formation in plants. However, little is known about HBs in sugarcane, a crucial sugar crop, due to its complex genetic background. Here, 302 allelic sequences for 104 HBs were identified and divided into 13 subfamilies in sugarcane Saccharum spontaneum. Comparative genomics revealed that whole-genome duplication (WGD)/segmental duplication significantly promoted the expansion of the HB family in S. spontaneum, with SsHB26, SsHB63, SsHB64, SsHB65, SsHB67, SsHB95, and SsHB96 being retained from the evolutionary event before the divergence of dicots and monocots. Based on the analysis of transcriptome and degradome data, we speculated that SsHB15 and SsHB97 might play important roles in regulating sugarcane leaf morphogenesis, with miR166 and SsAGO10 being involved in the regulation of SsHB15 expression. Moreover, subcellular localization and transcriptional activity detection assays demonstrated that these two genes, SsHB15 and SsHB97, were functional transcription factors. This study demonstrated the evolutionary relationship and potential functions of SsHB genes and will enable the further investigation of the functional characterization and the regulatory mechanisms of SsHBs.

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.

Transcriptomic analysis of tuberous root in two sweet potato varieties reveals the important genes and regulatory pathways in tuberous root development

Background

Tuberous root formation and development is a complex process in sweet potato, which is regulated by multiple genes and environmental factors. However, the regulatory mechanism of tuberous root development is unclear.

Results

In this study, the transcriptome of fibrous roots (R0) and tuberous roots in three developmental stages (Rl, R2, R3) were analyzed in two sweet potato varieties, GJS-8 and XGH. A total of 22,914 and 24,446 differentially expressed genes (DEGs) were identified in GJS-8 and XGH respectively, 15,920 differential genes were shared by GJS-8 and XGH. KEGG pathway enrichment analysis showed that the DEGs shared by GJS-8 and XGH were mainly involved in “plant hormone signal transduction” “starch and sucrose metabolism” and “MAPK signal transduction”. Trihelix transcription factor (Tai6.25300) was found to be closely related to tuberous root enlargement by the comprehensive analysis of these DEGs and weighted gene co-expression network analysis (WGCNA).

Conclusion

A hypothetical model of genetic regulatory network for tuberous root development of sweet potato is proposed, which emphasizes that some specific signal transduction pathways like “plant hormone signal transduction” “Ca²⁺signal” “MAPK signal transduction” and metabolic processes including “starch and sucrose metabolism” and “cell cycle and cell wall metabolism” are related to tuberous root development in sweet potato. These results provide new insights into the molecular mechanism of tuberous root development in sweet potato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08670-x.

Genome-wide identification, expression analyses of Wuschel-related homeobox (WOX) genes in Brachypodium distachyon and functional characterization of BdWOX12

As one kind of plant-specific transcription factors (TFs), WOX (Wuschel-related homeobox) plays an essential role in plant growth and development. In this study, 21 WOX TFs were identified in Brachypodium distachyon. They were divided into ancient, intermediate, and WUS clades based on phylogenetic analysis. These 21 BdWOX genes are mapped on 5 chromosomes unevenly. In the promoters, the most abundant cis-elements are ABRE, TGACG-motif, and G-box. qRT-PCR results showed that most BdWOX genes are expressed in vegetative and reproductive organs. Meanwhile, the expression of 14, 12, and 15 BdWOX genes are up-regulated by exogenous 6-BA, NAA, and GA, respectively. These results indicated that BdWOX genes participate in hormone signaling and regulate plant growth and development. Overexpression of BdWOX12 in Arabidopsis improved the root system, further indicating the functions of BdWOX genes in growth and development. This study provided a basis for the functional elucidation of BdWOX genes.

Genome-Wide Identification, Expansion, and Evolution Analysis of Homeobox Gene Family Reveals TALE Genes Important for Secondary Cell Wall Biosynthesis in Moso Bamboo (Phyllostachys edulis)

The TALE gene family is a subfamily of the homeobox gene family and has been implicated in regulating plant secondary growth. However, reports about the evolutionary history and function of the TALE gene family in bamboo are limited. Here, the homeobox gene families of moso bamboo Olyra latifolia and Bonia amplexicaulis were identified and compared. Many duplication events and obvious expansions were found in the TALE family of woody bamboo. PhTALEs were found to have high syntenies with TALE genes in rice. Through gene co-expression analysis and quantitative real-time PCR analysis, the candidate PhTALEs were thought to be involved in regulating secondary cell wall development of moso bamboo during the fast-growing stage. Among these candidate PhTALEs, orthologs of OsKNAT7, OSH15, and SH5 in moso bamboo may regulate xylan synthesis by regulating the expression of IRX-like genes. These results suggested that PhTALEs may participate in the secondary cell wall deposition in internodes during the fast-growing stage of moso bamboo. The expansion of the TALE gene family may be implicated in the increased lignification of woody bamboo when divergent from herbaceous bamboos.

Genome-Wide Identification, Expression, and Interaction Analysis of BEL-Like Homeodomain Gene Family in Peach

BEL1-like homeodomain (BLH) family genes as homeodomain transcription factors are found ubiquitously in plants to play important regulatory roles in reproductive development, morphological development, and stress response. Although BLH proteins have been reported in some species, there is little information about BLH genes in peach. In this study, we identified 11 peach PpBLH genes based on the conserved domain. Phylogenetic analysis suggested that the PpBLH proteins could be divided into five groups, which might be involved in different aspects of morphogenesis. Genomics structure analysis revealed that there were four exons in the PpBLH gene, and the length of the third exon was 61 bp. Chromosomal location analysis showed that the PpBLH genes were not distributed uniformly on six chromosomes. Promoter analysis showed that the promoter sequences of six PpBLH genes contained multiple cis-acting elements for hormones and stress. Six PpBLH genes were cloned by RT-PCR, and PpBLH1, PpBLH4, and PpBLH7 showed different expression patterns in the tested fruits under common temperature and high temperature. Y2H results indicated that PpBLH7 andPpBLH10 interacted with the PpOFP6 protein, and PpBLH1 interacted with the PpOFP1, PpOFP2, PpOFP4, and PpOFP13 proteins. These results provide new insight for further study of PpBLH genes, and construction of regulatory networks of PpBLH proteins in the growth, development, and stress response of peach.

Transcriptomic, Proteomic, and Metabolic Profiles of Catalpa bungei Tension Wood Reveal New Insight Into Lignin Biosynthesis Involving Transcription Factor Regulation

Lignin is a complex polymer in plant cell walls whose proportion is second only to that of cellulose and plays an important role in the mechanical properties of wood and stress resistance of plants. Here, we induced tension wood (TW) formation in Catalpa bungei by artificial bending and analyzed the lignin metabolism of the TW. LC-MS analysis showed that a significantly higher content of coniferyl aldehyde was observed in the TW cell wall than in the opposite wood (OW) and normal wood (NW) cell walls. TW had significantly lower contents of coniferyl alcohol than OW and NW. Raman spectroscopy results indicated that TW had lower total lignin than OW and NW. The transcription and translation levels of most of the differentially expressed genes (DEGs) involved in lignin monomer biosynthesis indicated upregulation in TW/OW and TW/NW. We found no significant difference in the transcription levels of three collision gases (CADs) between TW and OW or between NW, but their translation levels were significantly downregulated in TW, suggesting post-transcriptional control for CAD. We predicted and analyzed transcription factors that could target DEGs involved in lignin monomer biosynthesis in TW. Based on the analysis of the relationships of targeting and coexpression, we found that NAC (evm.model.group1.695) could potentially target 4CLs and CCoAOMT, that HD-Zip (evm.model.group7.1157) had potential targeting relationships with CCoAOMT, F5H, and CCR, and that their expression levels were significantly positive. It is speculated that the upregulation of NAC and HD-ZIP transcription factors activates the expression of downstream target genes, which leads to a significant increase in coniferyl aldehyde in TW. However, the decrease in total lignin in TW may be caused by the significant downregulation of CAD translation and the significant decrease in precursors (coniferyl alcohol). Whether the expression of CAD genes is regulated by post-transcriptional control and affects TW lignin metabolism needs further study.

Genome-Wide Analysis of the Homeobox Gene Family and Identification of Drought-Responsive Members in Populus trichocarpa

Homeobox (HB) genes play critical roles in the regulation of plant morphogenesis, growth and development. Here, we identified a total of 156 PtrHB genes from the Populus trichocarpa genome. According to the topologies and taxonomy of the phylogenetic tree constructed by Arabidopsis thaliana HB members, all PtrHB proteins were divided into six subgroups, namely HD-ZIP, ZF-HD, HB-PHD, TALE, WOX and HB-OTHERS. Multiple alignments of conserved homeodomains (HDs) revealed the conserved loci of each subgroup, while gene structure analysis showed similar exon–intron gene structures, and motif analysis indicated the similarity of motif number and pattern in the same subgroup. Promoter analysis indicated that the promoters of PtrHB genes contain a series of cis-acting regulatory elements involved in responding to various abiotic stresses, indicating that PtrHBs had potential functions in these processes. Collinearity analysis revealed that there are 96 pairs of 127 PtrHB genes mainly distributing on Chromosomes 1, 2, and 5. We analyzed the spatio-temporal expression patterns of PtrHB genes, and the virus-induced gene silencing (VIGS) of PtrHB3 gene resulted in the compromised tolerance of poplar seedlings to mannitol treatment. The bioinformatics on PtrHB family and preliminary exploration of drought-responsive genes can provide support for further study of the family in woody plants, especially in drought-related biological processes. It also provides a direction for developing new varieties of poplar with drought resistance. Overall, our results provided significant information for further functional analysis of PtrHB genes in poplar and demonstrated that PtrHB3 is a dominant gene regulating tolerance to water stress treatment in poplar seedlings.

Genome-wide characterization and expression and co-expression analysis suggested diverse functions of WOX genes in bread wheat

WUSCHEL-related homeobox (WOX) genes belong to the homeobox superfamily, are plant-specific and play vital functions in the growth and development. Herein, we identified a total of 43 TaWOX genes in the allohexaploid (AABBDD) genome of Triticum aestivum L. These genes were distributed on the various chromosomes of each subgenome (A, B and D). The phylogenetic analysis showed the clustering of TaWOXs into three clades: ancient, intermediate and modern or WUS. The gene and protein structures including exon/intron organization, intron phases, and domain and motif distribution were found to be conserved in each phylogenetic clade. The subcellular localization was predicted as nuclear. The Ka/Ks analyses suggested the purifying selection of paralogous genes. The differential expression profiling of various TaWOXs in numerous tissue developmental stages and different layers of grains suggested their role in growth and development. Moreover, a few genes exhibited modulated expression during abiotic and biotic stress conditions, which revealed their roles in stress response. The occurrence of various cis-acting regulatory elements further confirmed their role in plant development and stress tolerance. The co-expression analyses suggested the interactions of these genes with other genes, involved in various processes including plant development, signalling and stress responses. The present study reported several characteristic features of TaWOXs genes that can be useful for further characterization in future studies.

Transcriptome Reveals the Specificity of Phyllostachys edulis ‘Pachyloen’ Shoots at Different Developmental Stages

Phyllostachys edulis ‘Pachyloen’ can have a stalk wall thickness of up to 2.5 cm at a height of 1.3 m, which is 1.8 times that of normal Moso bamboo (Phyllostachys edulis); this serves as an excellent cultivar, comprising both wood and bamboo shoots. We collected bamboo shoot samples of Phyllostachys edulis ‘Pachyloen’ and Moso bamboo on a monthly basis from September to April and used transcriptome sequencing to explore the differences in their development. The results showed that there were 666–1839 Phyllostachys edulis ‘Pachyloen’-specific genes at different developmental stages enriched in 20 biological processes, 15 cellular components, 12 molecular functions, and 137 metabolic pathways, 52 of which were significant. Among these, 27 metabolic pathways such as tyrosine metabolism and their uniquely expressed genes were found to play important roles in the thickening of Phyllostachys edulis ‘Pachyloen’. This study provides insights into the mechanisms underlying the thickening of the culm wall of Phyllostachys edulis ‘Pachyloen’.