Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation

Macadamia nuts are one of the most important economic food items in the world. Pericarp thickness and flavonoid composition are the key quality traits of Macadamia nuts, but the underlying mechanism of pericarp formation is still unknown. In this study, three varieties with significantly different pericarp thicknesses, namely, A38, Guire No.1, and HAES 900, at the same stage of maturity, were used for transcriptome analysis, and the results showed that there were significant differences in their gene expression profile. A total of 3837 new genes were discovered, of which 1532 were functionally annotated. The GO, COG, and KEGG analysis showed that the main categories in which there were significant differences were flavonoid biosynthesis, phenylpropanoid biosynthesis, and the cutin, suberine, and wax biosynthesis pathways. Furthermore, 63 MiMYB transcription factors were identified, and 56 R2R3-MYB transcription factors were clustered into different subgroups compared with those in Arabidopsis R2R3-MYB. Among them, the S4, S6, and S7 subgroups were involved in flavonoid biosynthesis and pericarp formation. A total of 14 MiMYBs' gene expression were verified by RT-qPCR analysis. These results provide fundamental knowledge of the pericarp formation regulatory mechanism in macadamia nuts.

Identification and characterization of the WOX Gene Family revealed two WUS Clade Members associated with embryo development in Cunninghamia lanceolata

The WUSCHEL-related homeobox (WOX) gene family is vital for plant development and stress response. In this study, we conducted a comprehensive analysis of WOX genes in Cunninghamia lanceolata (C. lanceolata) and subsequently explored the potential roles of two ClWOX genes within the WUS clade. In total, six ClWOX genes were identified through a full-length transcriptome analysis. These genes, exhibiting conserved structural and functional motifs, were assigned to the ancient clade and Modern/WUS clade, respectively, through a phylogenetic analysis. Our expression analysis indicated that these ClWOX genes were highly expressed in the middle and late developmental stages of zygotic embryos in C. lanceolata. Moreover, only ClWOX5 and ClWOX6 within the Modern/WUS clade exhibited transcriptional activity, and their expressions were also induced in response to auxin and wounding. Overexpression of ClWOX5 and ClWOX6 in Arabidopsis caused a partially sterile phenotype, resulting in a very low seed setting rate. Transcriptomic analysis revealed that expressions of many embryo-defective (EMB) genes, phytohormone-related genes, and transcription factors (TFs) were dramatically altered in ClWOX5 and ClWOX6 transgenic plants, which suggested that ClWOX5 and ClWOX6 may play specific important roles in embryo development via complex gene networks. In addition, overexpression of ClWOX5 and ClWOX6 in leaf segments promoted shoot regeneration in tobacco, indicating that ClWOX5 and ClWOX6 can promote plant regeneration and could be used to improve genetic transformation. In conclusion, these results help to elucidate the function of the WOX gene and provide a valuable basis for future studies of the developmental regulation and applications of WOX genes in C. lanceolata.

An R2R3-MYB network modulates stem strength by regulating lignin biosynthesis and secondary cell wall thickening in herbaceous peony

Stem strength is an important agronomic trait affecting plant lodging, and plays an essential role in the quality and yield of plants. Thickened secondary cell walls in stems provide mechanical strength that allows plants to stand upright, but the regulatory mechanism of secondary cell wall thickening and stem strength in cut flowers remains unclear. In this study, first, a total of 11 non-redundant Paeonia lactiflora R2R3-MYBs related to stem strength were identified and isolated from cut-flower herbaceous peony, among which PlMYB43, PlMYB83 and PlMYB103 were the most upregulated differentially expressed genes. Then, the expression characteristics revealed that these three R2R3-MYBs were specifically expressed in stems and acted as transcriptional activators. Next, biological function verification showed that these P. lactiflora R2R3-MYBs positively regulated stem strength, secondary cell wall thickness and lignin deposition. Furthermore, yeast-one-hybrid and dual luciferase reporter assays demonstrated that they could bind to the promoter of caffeic acid O-methyltransferase gene (PlCOMT2) and/or laccase gene (PlLAC4), two key genes involved in lignin biosynthesis. In addition, the function of PlLAC4 in increasing lignin deposition was confirmed by virus-induced gene silencing and overexpression. Moreover, PlMYB83 could also act as a transcriptional activator of PlMYB43. The findings of the study propose a regulatory network of R2R3-MYBs modulating lignin biosynthesis and secondary cell wall thickening for improving stem lodging resistance, and provide a resource for molecular genetic engineering breeding of cut flowers.

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.