One of the largest families of transcription factors in plants, the MYB transcription factors family (Myeloblastosis, MYB TF), plays a vital role in regulating plant biochemical and physiological processes. The role of MYB TF in coping with stresses, such as drought, salt and cold, has been reported. Unfortunately, a comprehensive identification of R2R3-MYB TF in Masson pine (Pinus massoniana) has not been achieved. In this study, a total of 49 sequences were identified as R2R3-MYB TF. The structure, function and phylogenetic relationships of the conserved structural domains of Masson pine R2R3-MYB TF and Populus trichocarpa Torr. & A.Gray ex Hook. TFs were compared using bioinformatics tools. The results showed that Masson pine R2R3-MYB TF was divided into 24 groups, mainly located in the nucleus, and mostly lacking signal peptides and transmembrane structural domains with multiple phosphorylation sites. The drought stress-responsive R2R3-MYB gene, PmMYB4, was selected from the drought stress transcriptome based on analysis of the expression pattern and tissue specificity of PmMYB4 gene under abiotic stress using qPCR. The results showed that PmMYB4 can be involved in drought stress treatment through ABA signaling, as well as in multiple stress responses such as salt stress, and there were significant differences in the expression of PmMYB4 in the eight tissues. These results provide a reference scheme for the functional identification of R2R3-MYB transcription factors, which may be involved in plant responses to multiple stresses such as drought, and enrich our understanding of the functions of R2R3-MYB transcription factors in plants.

Transcription factors (TFs) are a class of proteins that play an important regulatory role in controlling the expression of plant target genes by interacting with downstream regulatory genes. The lateral organ boundary (LOB) structural domain (LBD) genes are a family of genes encoding plant-specific transcription factors that play important roles in regulating plant growth and development, nutrient metabolism, and environmental stresses. However, the LBD gene family has not been systematically identified in Pinus massoniana, one of the most important conifers in southern China. Therefore, in this study, we combined cell biology and bioinformatics approaches to identify the LBD gene family of P. massoniana by systematic gene structure and functional evolutionary analysis. We obtained 47 LBD gene family members, and all PmLBD members can be divided into two subfamilies, (Class I and Class II). By treating the plants with abiotic stress and growth hormone, etc., under qPCR-based analysis, we found that the expression of PmLBD genes was regulated by growth hormone and abiotic stress treatments, and thus this gene family in growth and development may be actively involved in plant growth and development and responses to adversity stress, etc. By subcellular localization analysis, PmLBD is a nuclear protein, and two of the genes, PmLBD44 and PmLBD45, were selected for functional characterization; secondly, yeast self-activation analysis showed that PmLBD44, PmLBD45, PmLBD46 and PmLBD47 had no self-activating activity. This study lays the foundation for an in-depth study of the role of the LBD gene family in other physiological activities of P. massoniana.

Cellulose synthase (CESA) is a key enzyme in the synthesis of cellulose, which plays an important role in cell wall construction and plant growth and development. In this study, seven CesA genes of P. massoniana were identified by searching the transcriptome data. Bioinformatics analysis showed that the putative CESA proteins were composed of 984–1101 amino acids, each containing the typical motifs of CESA proteins. Phylogenetic analysis showed that Transcript4609, Tran-script2643 and Transcript1263 were clustered into three groups with proteins related to regulating secondary wall synthesis, while Transcript691, Transcript1283, Transcript418 and Transcript556 were categorized into three clades with those associated with the formation of the primary cell walls. RT-qPCR analysis showed that the CesA genes were differentially expressed in different tissues, and most of the genes were induced by different abiotic stress and hormones. Transcript4609, Tran-script2643 and Transcript1263 were mainly expressed in the xylem and could respond to drought and salt stress induced by ABA, MeJA, ETH and SA hormones, indicating that these three CesA genes may play an important role in the response to abiotic stress in P. massoniana. This study revealed the possible biochemical and physiological functions of the CesA gene in P. massoniana, which can provide a basis for further exploration of the function of the CesA gene in cell wall formation and the response to external stress.