Next Article in Journal
Research Progress in Skin Aging, Metabolism, and Related Products
Next Article in Special Issue
SMALL PLANT AND ORGAN 1 (SPO1) Encoding a Cellulose Synthase-like Protein D4 (OsCSLD4) Is an Important Regulator for Plant Architecture and Organ Size in Rice
Previous Article in Journal
Human Synovial Mesenchymal Stem Cells Expressed Immunoregulatory Factors IDO and TSG6 in a Context of Arthritis Mediated by Alphaviruses
Previous Article in Special Issue
Virome Analysis of Aconitum carmichaelii Reveals Infection by Eleven Viruses, including Two Potentially New Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 21
10.3390/ijms242115938
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress

by
Mengyang Zhang
,
Romaric Hippolyte Agassin
,
Zichen Huang
,
Dengbao Wang
,
Sheng Yao
and
Kongshu Ji
*
State Key Laboratory of Tree Genetics and Breeding, Key Open Laboratory of Forest Genetics and Gene Engineering of National Forestry and Grassland Administration, Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(21), 15938; https://doi.org/10.3390/ijms242115938
Submission received: 31 August 2023 / Revised: 21 October 2023 / Accepted: 24 October 2023 / Published: 3 November 2023
(This article belongs to the Collection Feature Papers in Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Pinus massoniana is an important coniferous tree species for barren mountain afforestation with enormous ecological and economic significance. It has strong adaptability to the environment. TEOSINTE BRANCHED 1/CYCLOIDEA/PCF (TCP) transcription factors (TFs) play crucial roles in plant stress response, hormone signal transduction, and development processes. At present, TCP TFs have been widely studied in multiple plant species, but research in P. massoniana has not been carried out. In this study, 13 PmTCP TFs were identified from the transcriptomes of P. massoniana. The phylogenetic results revealed that these PmTCP members were divided into two categories: Class I and Class II. Each PmTCP TF contained a conserved TCP domain, and the conserved motif types and numbers were similar in the same subgroup. According to the transcriptional profiling analysis under drought stress conditions, it was found that seven PmTCP genes responded to drought treatment to varying degrees. The qRT-PCR results showed that the majority of PmTCP genes were significantly expressed in the needles and may play a role in the developmental stage. Meanwhile, the PmTCPs could respond to several stresses and hormone treatments at different levels, which may be important for stress resistance. In addition, PmTCP7 and PmTCP12 were nuclear localization proteins, and PmTCP7 was a transcriptional suppressor. These results will help to explore the regulatory factors related to the growth and development of P. massoniana, enhance its stress resistance, and lay the foundation for further exploration of the physiological effects on PmTCPs.

1. Introduction

Transcriptional regulation is an important mechanism for regulating gene expression during plant growth, development, and metabolism [1]. Transcription factors (TFs) are regulatory proteins that can directly or indirectly interact with cis-acting elements in promoters to regulate the expression of target genes, thereby regulating plant development and responding to external stresses [2]. Previous studies have analyzed the conserved regions of TEOSINTE BRANCHED 1 (TB1) protein of Zea mays, CYCLOIDEA (CYC) protein of Antirrhinum majus, and PROLIFERATING CELL TRANSCRIPTION FACTORS 1/2 (PCF1/2) of Oryza sativa. It was found that they all contained the bHLH conserved domain; therefore, the TEOSINTE BRANCHED 1/CYCLOIDEA/PCF (TCP) TF family was defined [3]. Presently, TCP TFs have been studied in many species, for example, Arabidopsis thaliana [4], Gossypium raimondii [5], and Malus domestica [6].
The structural characteristics of TCP TFs are that they have a basic Helix–Loop–Helix (bHLH) structure, also known as the TCP domain, which mainly participates in protein interactions and DNA binding [7,8]. TCP TFs can be divided into two categories: Class I and Class II. The Class II protein members have four more conserved amino acid residues than Class I protein members in the basic region of the TCP domain [9]. Class I is also named the PCF subfamily. According to the differences in amino acid sequences, especially the differences in the basic regions, the Class II members can be further divided into CYC/TB1 and CIN subfamilies. In addition, a few members of the CIN subfamily contain an arginine-rich motif (R domain), while the majority of members of the CYC/TB1 subfamily contain the R domain [10,11]. The members of Class I can promote cell differentiation and plant growth [12]. CYC/TB1 subfamily (Class II) members influence the symmetry of plant floral patterns and inhibit the growth of axillary buds, thereby reducing the number of branches [13,14]. The members of the CIN subfamily (Class II) can regulate leaf development [4].
TCP TFs play important roles in various biological functions of plants, including seed germination, organ differentiation, hormone regulation, signal transduction, and stress response [15,16,17]. Relevant studies have shown that AtTCP4 can promote the biosynthesis of auxin by regulating the expression of genes related to auxin biosynthesis enzymes [18]. AtTCP4 can induce the expression of lipoxygenase-related genes, which promotes the synthesis of methyl jasmonate (MeJA) and the senescence of mature leaves [19]. PeTCP10 responds to abscisic acid (ABA) signals and enhances tolerance to drought stress [20]. OsPCF6 and OsTCP21 are target genes of microRNA319, which can significantly improve plant resistance to low temperatures [21]. The AtTCP14 gene can directly activate the cis-acting elements in vascular tissue that regulate axillary bud growth, thus mediating the expression of germination-related genes and promoting plant growth [22]. After Ziziphus jujuba is infected with phytoplasma, the expression levels of ZjTCP genes increase in leaves and phloem, indicating that ZjTCPs have important roles in resistance to phytoplasma [23].
Pinus massoniana is the pioneer tree of afforestation in barren mountains in south China. It is the main coniferous tree species and is generally used in construction, timber, and papermaking; it also plays crucial roles in environmental regulation and increasing carbon sequestration [24]. P. massoniana can secrete pine resin, which is a secondary metabolite mainly composed of terpenoids. Its components have significant functions on the defense systems [25,26,27]. However, P. massoniana is always affected by various unfavorable factors throughout its entire growth process. The distribution environment of P. massoniana is complex, and the uneven distribution of rainfall limits the growth of P. massoniana and limits its expansion in area [28]. Drought stress leads to stomatal closure of P. massoniana, which seriously reduces its photosynthetic activity [29]. To adjust to complex and changeable conditions, plants have evolved various regulatory mechanisms. Moreover, TCP TFs play crucial roles in numerous biological processes in plants [30,31]. For example, the AtTCP10 gene is involved in cell differentiation, leaf morphogenesis, and other biological processes [32]. Under bright light conditions, AtTCP15 can regulate anthocyanin accumulation and protect plant cells from damage [33]. Therefore, it is of great significance to study the functions of TCP TFs of P. massoniana. At present, there are still no available genome resources for P. massoniana. Therefore, we identified the members of the PmTCP TF family based on the transcriptomes, namely, CO2 stress transcriptome, tender shoots transcriptome, and drought stress transcriptome. Furthermore, identification and analysis of 13 PmTCP TFs using bioinformatics methods. In addition, PmTCP genes were used for expression pattern analysis under different conditions. These results not only laid a theoretical foundation for the study of PmTCP TFs but also provided a reference for the response mechanism to external stress of P. massoniana.

2. Results

2.1. Identification of TCP Family Proteins in P. massoniana

A total of 13 TCP proteins were identified from the P. massoniana transcriptomes after removing sequences without conserved structural domain and duplicate sequences. Then, the 13 PmTCP TFs were named PmTCP1 to PmTCP13, and their protein sequences were used for physicochemical properties analysis (Table S1). These protein sequences varied from 182 (PmTCP4) to 626 (PmTCP10) amino acids (aa) in length, with the molecular weight (MW) varying from 19,932.49 (PmTCP4) to 68,110.09 (PmTCP10) Da. The isoelectric point (pI) values ranged from 6.69 (PmTCP2) to 9.49 (PmTCP7), and the aliphatic index varied from 52.27 (PmTCP11) to 74.07 (PmTCP4). In addition, the instability index (II) of PmTCP TFs was all higher than 40, indicating that they were unstable proteins. The grand average of hydropathicity (GRAVY) values were all negative, indicating that PmTCPs were hydrophilic proteins.

2.2. Phylogenetic Analysis of the TCP Proteins

A phylogenetic evolutionary tree was constructed comprising 13 TCP protein sequences from P. massoniana, 24 TCP protein sequences from A. thaliana, 13 TCP protein sequences from Picea abies, and 4 TCP protein sequences from Pinus taeda. The clustering results (Figure 1) showed that all TCP TFs were classified into Class I PCF, Class II CIN, and Class II CYC/TB1. According to the evolutionary classification of TCP TFs in A. thaliana, the TCP TFs of the species were divided into two groups, namely, Class I and Class II. Class I—also known as the PCF subfamily—and Class II could be further divided into CYC/TB1 and CIN subfamilies (Figure 1). 7 PmTCP TF members were distributed in the PCF subfamily. The CIN subfamily contained six PmTCP TF members. In addition, Arabidopsis TCP members were observed in all subfamilies, while three gymnospermic TCP members were only found in PCF and CIN subfamilies, indicating that TCP TFs of P. massoniana were more closely related to P. abies and P. taeda in a phylogenetic relationship.

2.3. Conserved Structural Domain Analysis of PmTCP TFs

To confirm the conserved domain of the PmTCP TFs, multiple sequence alignments of conserved domains were performed using DNAMAN 6.0 software. The result is shown in Figure 2. All PmTCP proteins contained a complete TCP domain, including a Helix–Loop–Helix domain and a basic domain. In the basic domain, Class I members PmTCP2, PmTCP3, PmTCP4, PmTCP5, PmTCP6, PmTCP11, and PmTCP12 lacked four amino acid residues compared to Class II members. In the HLH region, conservative characteristic sites were also present with lines. The results indicated that PmTCP proteins were highly conserved.

2.4. Conserved Motif Identification of PmTCP TFs

A total of 10 conserved motifs were determined in the 13 PmTCP family TFs by the MEME 5.5.4 software, and they were named motifs 1–10. The sequences of motifs are listed in Table S2. The number of amino acids in the 10 motifs ranged from 11 to 50, and motif 1 was the characteristic domain of TCP TFs. The distribution pattern of the motifs is shown in Figure 3. All PmTCP proteins contain the conserved motif 1, indicating that the structures of PmTCP TFs were complete. In addition, members of the same subfamily generally contained similar motif types. For instance, all members of Class I contained motifs 1 and 2. In addition, there were some differences in the motif types of different subfamilies. For example, motifs 2, 3, and 10 only exist in Class I members; motifs 5, 7, and 4 were only found in Class II members. The results suggested that the distribution of motifs was consistent with the evolutionary classification, and the differences in motifs of different members may imply the differences in protein structure.

2.5. Subcellular Localization of PmTCP TFs

The Cell-PLoc 2.0 software predicted that all PmTCP TFs were nuclear localization proteins (Table S3). To further verify the subcellular localization characteristics of PmTCP TFs, PmTCP7, and PmTCP12 were randomly selected from the Class II (CIN) and Class I (PCF) subfamilies for instantaneous transformation experiments, respectively. They were highly expressed in the drought-stressed transcriptome. The fluorescence signal was observed in the instantaneously transformed tobacco leaves (Figure 4). The GFP signal was found to be distributed throughout the whole cell in the control; however, GFP fused with PmTCP7 and PmTCP12 only observed fluorescence signals in the nucleus. The results showed that PmTCP proteins were localized in the nucleus.

2.6. Transcriptional Analysis of PmTCPs Using Drought Transcriptome Data

To study the potential function of PmTCP genes in drought stress, the heatmap of PmTCP genes was constructed based on the drought stress transcriptome data (PRJNA595650). Under drought stress conditions, the expressions of several genes were too low to be detected, so only the expression levels of PmTCP6–PmTCP13 genes were given (Figure 5). The results demonstrated that the expression patterns of PmTCP genes were different. For example, the expression level of the PmTCP7 gene increased whether it was mild drought, moderate drought, or severe drought. The expression patterns of PmTCP8, PmTCP11, and PmTCP12 genes were similar; their expression levels slightly decreased under mild drought but increased under moderate and severe drought. The expression of PmTCP10 only significantly increased under moderate drought conditions. However, the expression levels of PmTCP6 and PmTCP13 genes were significantly down-regulated after moderate and severe drought treatment. PmTCP9 expression significantly decreased throughout the drought treatment. These results indicated that PmTCP genes could be responsive to drought and may play vital roles in drought conditions.

2.7. Expression Patterns of PmTCPs in Various Tissues

The expression levels of four selected PmTCP (PmTCP7, PmTCP8, PmTCP11, PmTCP12) genes were analyzed by the quantitative reverse transcription PCR (qRT-PCR) in eight tissues: terminal bud (TB), young stems (YS), old stems (OS), young needles (YN), old needles (ON), phloem (P), xylem (X), and roots (R). These genes were selected based on their high expression levels in the drought-stressed transcriptome. The results (Figure 6) indicated that PmTCP genes were expressed differently in these eight tissues. PmTCP7 was found to be expressed at higher levels in young and old needles than in other tissues. The expression levels in young needles and old needles were 46.2 and 29.7 times those in terminal buds, respectively. The expression of PmTCP8 in young needles, old needles, and roots was significantly higher than that in other tissues, which was 12.2, 13.5, and 12.9 times that in terminal buds, respectively. PmTCP11 was also highly expressed in young and old needles; its expression level was 3.1–3.7 times that of the terminal buds. The expression levels in young stems, old stems, and phloem were approximately 2.4–2.9 times that in terminal buds. PmTCP12 had the highest expression level in roots, approximately 4.1-fold of that in terminal bud. The expression levels in young stems, old stems, phloem, and xylem were about 2.4, 2.7, 2.6, and 3.0 times those in terminal buds, respectively. The expression level of PmTCP12 in young needles and old needles was relatively consistent with those in terminal buds.

2.8. Expression Levels of PmTCP Genes Respond to Different Treatments

The expression patterns of these four PmTCP genes under eight distinct treatments are presented in Figure 7. Under drought conditions, the expression of PmTCP7 and PmTCP8 significantly increased at 20 d. PmTCP11 expression increased at 12 d and 20 d and reached the highest expression level at 20 d. There was no significant difference in the expression of PmTCP12 throughout the process. Under PEG treatment conditions, the expression of PmTCP7 and PmTCP12 slightly decreased during the process. The expression level of PmTCP8 increased at 3 h and then decreased. The expression level of PmTCP11 also increased throughout the process. Under mechanical injury treatment, the expression level of PmTCP7 decreased at 3 h and 6 h, indicating that the expression of PmTCP7 was inhibited under mechanical injury treatment. The expression of PmTCP8 was upregulated at 3 h and 24 h. PmTCP11 reached the highest expression level at 6 h. PmTCP12 expression slightly decreased. After SA treatment, the expression level of PmTCP7 slightly decreased at 3 h. PmTCP8 expression increased at 3 h and then showed a decreased trend. The expression profile of PmTCP11 was higher than that of the control at 6 h and 24 h. The expression level of PmTCP12 decreased during the process. Under MeJA treatment, the expression of PmTCP7 decreased at 6 h, and the results showed that PmTCP7 expression was suppressed at a certain time point during MeJA treatment. PmTCP8 and PmTCP11 had the highest expression levels at 3 h and 12 h, respectively. The expression of PmTCP12 was downregulated. After ETH treatment, the expression pattern of PmTCP7 was lower than that of control at 6 h and reached its peak at 24 h. The expression levels of PmTCP8 and PmTCP12 were downregulated throughout the process. PmTCP11 expression pattern tended to increase slightly at 6 h. Under H2O2 treatment, the expression trends of PmTCP7 and PmTCP8 first increased, subsequently decreased, and then increased. The expression of PmTCP11 was upregulated at 12 h and 24 h. PmTCP12 expression pattern tends to stabilize, and the response to H2O2 was not significant. Under ABA treatment, the expression levels of PmTCP7 and PmTCP12 were downregulated and suppressed throughout the ABA process. The expression of PmTCP8 increased at 3 h. There was no significant difference in the expression of PmTCP11 throughout the whole process.

2.9. Transcriptional Activity Analysis of PmTCP7

Moreover, transcriptional activity of PmTCP7 was detected and analyzed (Figure 8). The PmTCP7 gene was successfully cloned. This gene had the highest expression in the drought-stressed transcriptome and responded to several hormonal treatments. In addition, its expression level is high in tissues, especially in needles. Yeast cells carrying the fusion expression vector of pGBKT7-PmTCP7 could grow on SD/-Trp yeast medium; however, it was unable to grow on yeast SD/-Trp/-Ade/-His selective media and did not display blue spots on SD/-Trp/-Ade/-His yeast media containing X-α-gal. This is depicted in various graphics. The findings indicated that PmTCP7 was a transcriptional inhibitor; this study laid the groundwork for further research on gene regulation.

3. Discussion

TCP TFs have been identified in many species, for example, 24 TCP TFs in A. thaliana [4], 22 TCP TFs in O. sativa [4], 33 TCP TFs in Populus euphratica [34], 9 TCP TF members in pineapple [35], and 18 VvTCPs were identified in Vitis vinifera [36]. The TCP family is widely involved in multiple biological processes in plants that have important effects on regulating plant growth, defense responses, and hormone signal transduction. Increasingly, studies have shown that TCP genes are promising candidates with regard to plant resistance and development [37,38,39]. In addition, P. massoniana has strong adaptability to the environment. Therefore, reports on the correlation between TCP TFs and P. massoniana stress and growth are meaningful. Therefore, this study provides a comprehensive analysis of the TCP TF members of P. massoniana, which combines PmTCPs with the aspects related to growth and stress, providing a basis for further understanding the biological functions of PmTCP TFs.
Thirteen PmTCP TF members were identified from the transcriptomes of P. massoniana and named PmTCP1–PmTCP13. According to the classification of A. thaliana TCP TFs, The 13 PmTCP members were divided into Class I and Class II, where Class I contained 7 PmTCP proteins and Class II contained 6 PmTCP TF members, with a ratio of about 1:1. Relevant studies have shown that there were 10 VvTCP TFs in Class I and 8 VvTCPs in Class II of V. vinifera [36]. In Prunus mume, Class I contained 10 TCP members, and Class II contained 9 TCP members [40], which was also close to 1:1. These results were consistent with the current results. Multiple sequence alignment analysis showed that these PmTCP TFs had bHLH functional conserved domains. The conserved domains of Class I proteins lacked four amino acid residues compared to Class II proteins, which may lead to differences in function. This phenomenon has also been observed in tobacco [41]. Motif 1, a characteristic motif of TCP protein, exists in all PmTCP amino acid sequences and is an important basis for identifying TCP members. In the same subgroup members, the number and type of conserved motifs were similar, and the distribution order was the same. However, different subgroup members had their own representative motifs, which strongly supported the phylogenetic relationship and classification of PmTCPs. The difference and consistency of motif distribution provided a foundation for the functional diversity and consistency of PmTCP TFs. This pattern was also observed in Cucumis sativus [42]. Overall, these results suggested that PmTCP TFs were structurally conservative.
The results of tobacco transient expression experiments and subcellular localization prediction implied that PmTCP TFs are nuclear localization proteins, which is consistent with the localization characteristics of TFs and the theory that TFs normally play roles in the nucleus.
To some extent, the expression patterns of genes are related to their functions. It has been reported that TCP TFs are involved in many aspects of plant biological processes. Drought resistance is an important response mechanism for plants to adapt to external stresses [28]. Therefore, the expression profile of PmTCP TFs was analyzed using the drought stress transcriptome data of P. massoniana. It was found that the expression levels of PmTCP7, PmTCP8, PmTCP11, and PmTCP12 genes significantly increased after moderate and severe drought treatments. The expression of the PmTCP10 gene was significantly upregulated after moderate drought stress. The expressions of PmTCP6, PmTCP9, and PmTCP13 genes were significantly downregulated after drought stress treatment. Additionally, there was no significant difference in the expression of most PmTCP genes compared to the control under mild drought, while the expression levels were obviously increased or decreased under moderate and severe drought conditions. It is speculated that PmTCPs could be responsive to drought stress, but their response mechanisms were complex and different. These findings laid a foundation for further research on the function of PmTCP genes in drought stress.
It has been observed that the function of the TCP genes can regulate the growth process of many tissues in plants, such as stems, leaves, flowers, and roots [32,43]. To investigate the expression characteristics of PmTCP genes in tissues, the expression of PmTCP members in different tissues of P. massoniana was analyzed. The results indicated that PmTCP genes could be expressed in various tissues. The expression levels of PmTCP7 and PmTCP8 genes in young needles and old needles were significantly higher than those in other tissues; in addition, PmTCP8 was also highly expressed in roots. According to the phylogenetic evolutionary relationships, PmTCP7 and PmTCP8 belong to Class II members. It has been reported that Class II members are specifically expressed in plant tissues and participate in the formation of leaf trichosomes and roots [44]. These results suggested that PmTCP7 and PmTCP8 may have regulation roles in needles and roots.
Relevant studies have shown that TCP TFs can respond to a variety of hormone and abiotic stresses [45,46]. The expression patterns of PmTCP genes were determined under different treatments. The expression of PmTCP7, PmTCP8, and PmTCP11 significantly increased under drought stress, and the expression of PmTCP8 and PmTCP11 was upregulated after PEG treatment. The results indicated that these genes had a positive response to drought and PEG stress, which may induce an increase in POD and SOD enzyme activities, as well as MDA and Pro contents, thereby enhancing osmotic regulation and reducing damage to P. massoniana [47]. The increased expression levels of PmTCP8 and PmTCP11 under mechanical injury suggested that they may be involved in biological processes, such as preventing the invasion of pathogens. SA, MeJA, ETH, H2O2, and ABA are important exogenous exciters that are induced by external stimuli and have an important role in response to stress [48]. Through gene expression levels, it was found that PmTCP7, PmTCP8, and PmTCP11 positively responded to H2O2. Meanwhile, PmTCP8 and PmTCP11 had a positive response to MeJA and SA. PmTCP8 and PmTCP7 could respond to ABA and ETH, respectively. The results showed that these genes were involved in signal transduction and stress perception, resulting in the adaptation changes of P. massoniana. The expression of TCP genes has been analyzed in many plants under these treatments; for example, AtTCP20 is involved in the biosynthesis of jasmonic acid [19]. Overexpression of the OsTCP19 gene can significantly improve drought resistance in mature plants and seedlings [49]. In general, TCP genes are widely involved in stress resistance and hormone signaling pathways in plants. However, further functional research is needed to determine how the TCP genes accurately respond to these signaling.
By observing the growth conditions of yeast cells containing pGBKT7–PmTCP7 fusion expression vector on the selected medium, it was impossible to detect the transcriptional self-activation activity of these cells. This study provides a basis for studying the regulatory network associated with TCP proteins.
In summary, this study comprehensively analyzed the TCP TFs family of P. massoniana. It not only contributes to screening suitable candidate genes for further investigation but also contributes to understanding the molecular regulatory mechanism and signal transduction of TCP genes in plant resistance and development. In the meantime, the functional characterization of PmTCPs provides insights into the function of TCP TFs in P. massoniana.

4. Materials and Methods

4.1. Identification and Sequence Analysis of TCP Family Members of P. massoniana

A hidden Markov model (HMM) of the TCP conserved domain (PF03634) was acquired through the Pfam database (http://pfam.xfam.org/ (accessed on 17 April 2023)) [50], and TCP protein members were searched using the HMMER 3.0 software from the P. massoniana transcriptomes with its defined threshold of E < 10−5. The transcriptomes of P. massoniana originated from the previously determined drought stress transcriptome (PRJNA595650), CO2 stress transcriptome (PRJNA561037), and tender shoots transcriptome (PRJNA655997). PmTCP proteins were identified from these transcriptomes. The SMART online website (http://smart.embl-heidelberg.de/, accessed on 18 April 2023) was used to predict the conserved domain, and proteins without structural domain were deleted. Then, the protein sequences were aligned, and sequences with similarity exceeding 97% were deleted. The protein sequences are listed in Table S4. The physicochemical properties, such as instability index (II), grand average of hydropathicity (GRAVY), isoelectric point (pI), and molecular weight (MW) of the screened TCP proteins were analyzed with ExPASy (http://web.expasy.org/protparam/ (accessed on 19 April 2023)).

4.2. Phylogenetic Evolutionary of TCP TFs

The TCP protein sequences of A. thaliana, P. taeda, and P. abies were downloaded from the Plant Transcription Factor Database (http://planttfdb.cbi.pku.edu.cn/ (accessed on 23 April 2023)). Their protein sequences and those of P. massoniana were aligned with the ClustalW 2.1 program. A phylogenetic evolutionary tree was constructed using the Neighbor-Joining (NJ) method with MEGA 7.0 software; the bootstrap value was 1000, and other parameters were default values [51]. Subsequently, the Evolview online tool (http://evolgenius.info//evolview-v2/ (accessed on 27 April 2023)) was used to embellish the phylogenetic tree.

4.3. Multiple Sequence Alignments and Motif Analysis

DNAMAN 6.0 software was used to align the protein sequences of PmTCP members and coTBtools their conserved domains. The Multiple Em for Motif Elicitation (MEME) (http://meme.nbcr.net/meme/intro.html (accessed on 30 April 2023)) online website was used to obtain the conserved motifs organization of the 13 PmTCP proteins, and 10 motifs were searched using default parameters [52]. Then, the TBtools-II v2.019 platform was used to analyze their conserved structures and produce diagrams [53].

4.4. Subcellular Localization Analysis

The Cell-PLoc 2.0 (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/ (accessed on 1 July 2023)) online software was used to predict subcellular localization of PmTCP proteins. PmTCP7 and PmTCP12 were randomly selected for transient transformation experiments. The primers for constructing vectors are shown in Table S5. The open reading frame (ORF) regions of PmTCP7 and PmTCP12 without stop codon were linked with pBI121-GFP vector, the 35S::PmTCP7-GFP and 35S::PmTCP12-GFP expression vectors were constructed by recombinase. Then, the vectors were transferred to Agrobacterium GV3101. These strains, along with P19 (RNA Silencing Inhibitor) Agrobacterium strains, were cultured in LB media for 48 h at 28 °C. Then, they were suspended in a solution containing 150 μM acetosyringone, 10 mM MgCl2, and 10 mM 2-(N-morpholino)ethanesulfonic acid (MES). The suspension cells were mixed with p19 at a ratio of 1:1, and the Nicotiana benthamiana leaves were injected with the mixed solution. Subsequently, the infiltrated N. benthamiana plants were maintained for 3 d under the conditions of 16 h light and 8 h dark photoperiod. The LSM710 confocal microscope (Zeiss, Jena, Germany) was used to capture the GFP signal.

4.5. Transcriptional Profile Analysis of the PmTCP Genes

To study the expression of PmTCP genes, we utilized the original drought-stressed transcriptome (PRJNA595650) of P. massoniana to obtain the RNA sequencing (RNA-seq) data. The soil moisture content for the growth of P. massoniana was set to four gradients: normal (CK), mild (T1), medium (T2), and serious (T3). The corresponding water holding capacity was CK (80 ± 5)%, T1 (65 ± 5)%, T2 (50 ± 5)% and T3 (35 ± 5)%, respectively. They were placed at 75% humidity for 60 d and subsequently sequenced. Fragments per kilobase of exon model per million mapped fragments (FPKM) values (Table S6) were computed to assess the expression level of PmTCPs. Then, a heat map was generated based on log2(FPKM+1) values using TBtools-II v2.019 software.

4.6. Plant Materials and Treatments

Two-year-old P. massoniana seedlings were planted in pots containing a soil mixture. The volume ratio of vermiculite, perlite, and peat of the soil mixture was 3:1:1. Then, seedlings were placed under the conditions with a temperature of 24 °C and photoperiod of 16 h light and 8 h dark for growth. When P. massoniana seedlings grew for about a month and remained in stable conditions, the seedlings with consistent growth status were selected for subsequent treatments.
To investigate the expression levels of PmTCP genes, eight tissues were sampled from the seedlings, including young needles, old needles, young stems, old stems, phloem, xylem, terminal buds, and roots. Furthermore, the seedlings were treated with these treatments: 1 mM salicylic acid (SA), 10 mM methyl jasmonate (MeJA), 50 µM ethephon (ETH), 100 µM abscisic acid (ABA), 10 mM hydrogen peroxide (H2O2), 15% polyethylene glycol (PEG6000), mechanical injury, and drought. The mechanical injury was carried out by cutting off the upper half of the seedling; 15%PEG treatment involved immersing seedlings in this solution, causing osmotic stress. Drought treatment was performed with watering at 0 d and then naturally evaporated for 20 d. The remaining treatments were all sprayed on the surface of the seedlings. Afterward, samples were collected at 0 h, 3 h, 6 h, 12 h, and 24 h, except for drought. The drought-treated samples were collected at 0 d, 3 d, 7 d, 12 d and 20 d. The expression levels of the untreated samples were used as controls. Three biological replicates were used for each treatment.

4.7. RNA Extraction and qRT-PCR Analysis

Total RNA was extracted using a FastPure Plant Total RNA Isolation Kit (RC401-01, Vazyme Biotech, Nanjing, China). Then, a NanoDrop2000 (Thermo Fisher Scientific, Waltham, MA, USA) was used to detect the purity and concentration of RNA. Afterward, 1000 ng total RNA was reverse-transcribed to synthesize 20 μL of cDNA using the first-strand cDNA synthesis kit (AT311, TransGen Biotech, Beijing, China). Primers for qRT-PCR were designed by Primer 5.0, and α-tubulin (TUA) was used as an internal reference gene [54]. These primers are displayed in Table S6. qRT-PCR was carried out under the StepOne Plus program (Applied Biosystems, Foster City, CA, USA) using SYBR Green Real-time PCR Master Mix (QPK-201, Toyobo Bio-Technology, Shanghai, China). Each PCR mixture was a 10 μL system that contained 5 μL of SYBR Green Real-time PCR Master Mix, 1 μL of 20-fold diluted cDNA, 3.2 μL of sterile water, and 0.4 μL of each primer. The amplification procedure was as follows: predenaturation at 95 °C for 2 min, denaturation at 95 °C for 10 s, annealing and extension at 60 °C for 30 s; this process was repeated 40 times. The remaining procedures were the program default parameters. Each reaction was detected with three independent replicates. The relative expression levels of PmTCP genes were computed according to the 2−ΔΔCT method [55]. One-way ANOVA and multiple comparison tests were analyzed by GraphPad Prism 8.0 software.

4.8. Transcription Self-Activation Detection

The open reading frame (ORF) of the PmTCP7 gene was fused with the pGBKT7 vector cut by EcoRI and BamHI double enzymes. The primers for vector construction are provided in Table S5. Then, the correctly sequenced recombinant vector was transformed into AH109 yeast strains (YC1010, Weidi Biotech, Shanghai, China), which were subsequently cultured on the SD/-Trp medium at 28 °C for 3 d. Single yeast colonies were selected for PCR validation, the successfully transformed yeast single colonies were collected in SD/-Trp liquid medium, and 5 μL of diluent was inoculated onto the surface of yeast medium SD/-Trp and SD/-Trp/-Ade/-His, respectively. Meanwhile, 5 μL of diluent was inoculated on the SD/-Trp/-Ade/-His solid medium containing X-α-gal. Subsequently, the growth status of yeast cells was observed and recorded.

5. Conclusions

In this study, we identified 13 PmTCP TFs from transcriptomes of P. massoniana. These PmTCPs were divided into Class I and Class II categories, and different subgroup members had different motif compositions. Members of the same subgroup had similar conserved motif composition, indicating that they may perform similar functions. According to the transcriptome of drought stress, seven PmTCPs were found to be responsive to drought stress. Meanwhile, the expression patterns of the selected four PmTCP genes under different conditions suggested that they may have important functions on several stress responses and play multiple roles in different stages of plant growth and development. These findings provide a theoretical basis for the study of TCP TFs, are conducive to further understanding the functions of PmTCPs, supply potential strategies for breeding in P. massoniana, and contribute to the study of the mechanism resistance ability of other plants.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms242115938/s1.

Author Contributions

Conceptualization, M.Z.; software, M.Z. and D.W.; methods, M.Z. and S.Y.; analysis, M.Z.; investigation, M.Z. and Z.H.; resources, M.Z.; writing—original draft preparation, M.Z.; writing—reviewing and editing, M.Z., R.H.A. and K.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key R&D Program of China (2022YFD2200202), and the project was funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, M.M.; Wang, M.M.; Yang, J.; Wen, J.; Guo, P.C.; Wu, Y.W.; Ke, Y.Z.; Li, P.F.; Li, J.N.; Du, H. Evolutionary and comparative expression analyses of TCP transcription factor gene family in land plants. Int. J. Mol. Sci. 2019, 20, 3591. [Google Scholar] [CrossRef]
  2. Perez, M.; Guerringue, Y.; Ranty, B.; Pouzet, C.; Aldon, D. Specific TCP transcription factors interact with and stabilize PRR2 within different nuclear sub-domains. Plant Sci. 2019, 287, 110197. [Google Scholar] [CrossRef]
  3. Cubas, P.; Lauter, N.; Doebley, J.; Coen, E. The TCP domain: A motif found in proteins regulating plant growth and development. Plant J. 1999, 18, 215–222. [Google Scholar] [CrossRef]
  4. Yao, X.; Ma, H.; Wang, J.; Zhang, D.B. Genome-wide comparative analysis and expression pattern of TCP gene families in Arabidopsis thaliana and Oryza sativa. J. Integr. Plant Biol. 2007, 49, 885–897. [Google Scholar] [CrossRef]
  5. Ma, J.; Wang, Q.L.; Sun, R.R.; Xie, F.L.; Jones, D.C.; Zhang, B.H. Genome-wide identification and expression analysis of TCP transcription factors in Gossypium raimondii. Sci. Rep. 2014, 4, 6645. [Google Scholar] [CrossRef]
  6. Xu, R.R.; Sun, P.; Jia, F.J.; Lu, L.T.; Li, Y.Y.; Zhang, S.Z.; Huang, J.G. Genome-wide analysis of TCP transcription factor gene family in Malus domestica. J. Genet. 2014, 93, 733–746. [Google Scholar] [CrossRef]
  7. Manassero, N.G.U.; Viola, I.L.; Welchen, E.; Gonzalez, D.H. TCP transcription factors: Architectures of plant form. Biomol. Concepts 2013, 4, 111–127. [Google Scholar] [CrossRef]
  8. Kieffer, M.; Master, V.; Waites, R.; Davies, B. TCP14 and TCP15 affect internode length and leaf shape in Arabidopsis. Plant J. 2011, 68, 147–158. [Google Scholar] [CrossRef]
  9. Martin-Trillo, M.; Cubas, P. TCP genes: A family snapshotten years later. Trends Plant Sci. 2010, 15, 31–39. [Google Scholar] [CrossRef]
  10. Howarth, D.; Donoghue, M. Phylogenetic analysis of the “ECE” (CYC/TB1) clade reveals duplications predating the core eudicots. Proc. Natl. Acad. Sci. USA 2006, 103, 9101–9106. [Google Scholar] [CrossRef]
  11. Navaud, O.; Dabos, P.; Carnus, E.; Tremousaygue, D.; Hervé, C. TCP transcription factors predate the emergence of land plants. J. Mol. Evol. 2007, 65, 23–33. [Google Scholar] [CrossRef]
  12. Kosugi, S.; Ohashi, Y. PCF1 and PCF2 specifically bind to cis elements in the rice proliferating cell nuclear antigen gene. Plant Cell 1997, 9, 1607–1619. [Google Scholar]
  13. Broholm, S.K.; Tahtiharju, S.; Laitinen, R.A.E.; Albert, V.A.; Teeri, T.H.; Elomaa, P. A TCP domain transcription factor controls flower type specification along the radial axis of the Gerbera (Asteraceae) inflorescence. Proc. Natl. Acad. Sci. USA 2008, 105, 9117–9122. [Google Scholar] [CrossRef]
  14. Finlayson, S.A. Arabidopsis TEOSINTE BRANCHED1-LIKE 1 regulates axillary bud outgrowth and is homologous to monocot TEOSINTE BRANCHED1. Plant Cell Physiol. 2007, 48, 667–677. [Google Scholar] [CrossRef]
  15. Parapunova, V.; Busscher, M.; Busscher-Lange, J.; Lammers, M.; Karlova, R.; Bovy, A.D.; Angenent, G.C.; de Maagd, R.A. Identification, cloning and characterization of the tomato TCP transcription factor family. BMC Plant Biol. 2014, 14, 157. [Google Scholar] [CrossRef]
  16. Danisman, S. TCP Transcription factors at the interface between environmental challenges and the plant’s growth responses. Front. Plant Sci. 2016, 7, 1930. [Google Scholar] [CrossRef]
  17. Bao, S.; Zhang, Z.X.; Lian, Q.; Sun, Q.H.; Zhang, R.F. Evolution and expression of genes encoding TCP transcription factors in Solanum tuberosum reveal the involvement of StTCP23 in plant defence. BMC Genet. 2019, 20, 91. [Google Scholar] [CrossRef]
  18. Challa, K.R.; Aggarwal, P.; Nath, U. Activation of YUCCA5 by the transcription factor TCP4 integrates developmental and environmental signals to promote hypocotyl elongation in Arabidopsis. Plant Cell 2016, 28, 2117–2130. [Google Scholar]
  19. Danisman, S.; Van der Wal, F.; Dhondt, S.; Waites, R.; de Folter, S.; Bimbo, A.; van Dijk, A.J.; Muino, J.M.; Cutri, L.; Dornelas, M.C.; et al. Arabidopsis Class I and Class II TCP transcription factors regulate jasmonic acid metabolism and leaf development antagonistically. Plant Physiol. 2012, 159, 1511–1523. [Google Scholar] [CrossRef]
  20. Liu, H.L.; Gao, Y.M.; Wu, M.; Shi, Y.N.; Wang, H.; Wu, L.; Xiang, Y. TCP10, a TCP transcription factor in moso bamboo (Phyllostachys edulis), confers drought tolerance to transgenic plants. Environ. Exp. Bot. 2020, 172, 104002. [Google Scholar] [CrossRef]
  21. Wang, S.T.; Sun, X.L.; Hoshino, Y.; Yu, Y.; Jia, B.; Sun, Z.W.; Sun, M.Z.; Duan, X.B.; Zhu, Y.M. MicroRNA319 positively regulates cold tolerance by targeting OsPCF6 and OsTCP21 in rice (Oryza sativa L.). PLoS ONE 2014, 9, e91357. [Google Scholar] [CrossRef]
  22. Tatematsu, K.; Nakabayashi, K.; Kamiya, Y.; Nambara, E. Transcription factor AtTCP14 regulates embryonic growth potential during seed germination in Arabidopsis thaliana. Plant J. 2008, 53, 42–52. [Google Scholar] [CrossRef]
  23. Chen, P.; Li, J.D.; Ye, X.; Tan, B.; Zheng, X.B.; Cheng, J.; Wang, W.; Wang, H.Y.; Gu, L.Y.; Feng, J.C. Genomewide identification of Ziziphus jujuba TCP transcription factors and their expression in response to infection with jujube witches’ broom phytoplasma. Acta Physiol. Plant. 2019, 41, 86. [Google Scholar] [CrossRef]
  24. Maleki, S.S.; Mohammadi, K.; Ji, K.S. Study on factors influencing transformation efficiency in Pinus massoniana using Agrobacterium tumefaciens. Plant Cell Tissue Organ Cult. 2018, 133, 437–445. [Google Scholar] [CrossRef]
  25. Bai, Q.S.; He, B.X.; Cai, Y.L.; Lian, H.M.; Zhang, Q. Transcriptomic and metabolomic analyses reveal several critical metabolic pathways and candidate genes involved in resin biosynthesis in Pinus massoniana. Mol. Genet. Genomics 2020, 295, 327–341. [Google Scholar] [CrossRef]
  26. Karanikas, C.; Walker, V.; Scaltsoyiannes, A.; Comte, G.; Bertrand, C. High vs. low yielding oleoresin Pinus halepensis Mill. trees GC terpenoids profiling as diagnostic tool. Ann. For. Sci. 2010, 67, 412. [Google Scholar] [CrossRef]
  27. Tholl, D. Biosynthesis and biological functions of terpenoids in plants. Adv. Biochem. Eng. Biotechnol. 2015, 148, 63–106. [Google Scholar]
  28. Fang, Y.J.; Xiong, L.Z. General mechanisms of drought response and their application in drought resistance improvement in plants. Cell. Mol. Life Sci. 2015, 72, 673–689. [Google Scholar] [CrossRef]
  29. Martin-St Paul, N.; Delzon, S.; Cochard, H. Plant resistance to drought depends on timely stomatal closure. Ecol. Lett. 2017, 20, 1437–1447. [Google Scholar] [CrossRef]
  30. Zhou, H.Y.; Hwarari, D.; Ma, H.Y.; Xu, H.B.; Yang, L.M.; Luo, Y.M. Genomic survey of TCP transcription factors in plants: Phylogenomics, evolution and their biology. Front. Genet. 2022, 13, 1060546. [Google Scholar] [CrossRef]
  31. Yin, Z.J.; Li, Y.; Zhu, W.D.; Fu, X.Q.; Han, X.L.; Wang, J.J.; Lin, H.; Ye, W.W. Identification, characterization, and expression patterns of TCP genes and microRNA319 in Cotton. Int. J. Mol. Sci. 2018, 19, 3655. [Google Scholar] [CrossRef] [PubMed]
  32. Koyama, T.; Mitsuda, N.; Seki, M.; Shinozaki, K.; Ohme-Takagi, M. TCP transcription factors regulate the activities of ASYMMETRIC LEAVES1 and miR164, as Well as the auxin response, during differentiation of leaves in Arabidopsis. Plant Cell 2010, 22, 3574–3588. [Google Scholar] [CrossRef]
  33. Viola, I.L.; Camoirano, A.; Gonzalez, D.H. Redox-dependent modulation of anthocyanin biosynthesis by the TCP transcription factor TCP15 during exposure to high light intensity conditions in Arabidopsis. Plant Physiol. 2016, 170, 74–85. [Google Scholar] [CrossRef] [PubMed]
  34. Ma, X.D.; Ma, J.C.; Fan, D.; Li, C.F.; Jiang, Y.Z.; Luo, K.M. Genome-wide identification of TCP family transcription factors from Populus euphratica and their involvement in leaf shape regulation. Sci. Rep. 2016, 6, 32795. [Google Scholar] [CrossRef] [PubMed]
  35. Li, Z.Q.; Ouyang, Y.W.; Pan, X.L.; Zhang, X.H.; Zhao, L.; Wang, C.; Xu, R.; Zhang, H.N.; Wei, Y.Z. TCP transcription factors in pineapple: Genome-wide characterization and expression profile analysis during flower and fruit development. Horticulturae 2023, 9, 799. [Google Scholar] [CrossRef]
  36. Leng, X.P.; Wei, H.R.; Xu, X.Z.; Ghuge, S.A.; Jia, D.J.; Liu, G.S.; Wang, Y.Z.; Yuan, Y.B. Genome-wide identification and transcript analysis of TCP transcription factors in grapevine. BMC Genom. 2019, 20, 786. [Google Scholar] [CrossRef] [PubMed]
  37. Takeda, T.; Suwa, Y.; Kitano, H.; Ueguchi-Tanaka, M.; Ashikari, M.; Matsuoka, M.; Ueguchi, C. The OsTB1 gene negatively regulates lateral branching in rice. Plant Cell Physiol. 2003, 33, 513–520. [Google Scholar] [CrossRef]
  38. Guan, P.Z.; Ripoll, J.J.; Wang, R.H.; Vuong, L.; Bailey-Steinitz, L.J.; Ye, D.N.; Crawford, N.M. Interacting TCP and NLP transcription factors control plant responses to nitrate availability. Proc. Natl. Acad. Sci. USA 2017, 114, 2419–2424. [Google Scholar] [CrossRef]
  39. Yang, M.F.; He, G.D.; Hou, Q.D.; Fan, Y.; Duan, L.L.; Li, K.Y.; Wei, X.L.; Qiu, Z.L.; Chen, E.J.; He, T.B. Systematic analysis and expression profiles of TCP gene family in Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) revealed the potential function of FtTCP15 and FtTCP18 in response to abiotic stress. BMC Genom. 2022, 23, 415. [Google Scholar] [CrossRef]
  40. Zhou, Y.Z.; Xu, Z.D.; Zhao, K.; Yang, W.R.; Cheng, T.R.; Wang, J.; Zhang, Q.X. Genome-wide identification, characterization and expression analysis of the TCP Gene family in Prunus mume. Front. Plant Sci. 2016, 7, 1301. [Google Scholar] [CrossRef]
  41. Chen, L.; Chen, Y.Q.; Ding, A.M.; Chen, H.; Xia, F.; Wang, W.F.; Sun, Y.H. Genome-wide analysis of TCP family in tobacco. Genet. Mol. Res. 2016, 15, 1507728. [Google Scholar] [CrossRef] [PubMed]
  42. Wen, H.F.; Chen, Y.; Du, H.; Zhang, L.Y.; Zhang, K.Y.; He, H.L.; Pan, J.S.; Cai, R.; Wang, G. Genome-wide identification and characterization of the TCP gene family in cucumber (Cucumis sativus L.) and their transcriptional responses to different treatments. Genes 2020, 11, 1379. [Google Scholar] [CrossRef] [PubMed]
  43. Chai, W.B.; Jiang, P.F.; Huang, G.Y.; Jiang, H.Y.; Li, X.Y. Identification and expression profiling analysis of TCP family genes involved in growth and development in maize. Physiol. Mol. Biol. Plants 2017, 23, 779–791. [Google Scholar] [CrossRef] [PubMed]
  44. Qi, X.Y.; Qu, Y.X.; Gao, R.; Jiang, J.F.; Fang, W.M.; Guan, Z.Y.; Zhang, F.; Zhao, S.; Chen, S.M.; Chen, F.D. The heterologous expression of a Chrysanthemum nankingense TCP transcription factor blocks cell division in yeast and Arabidopsis thaliana. Int. J. Mol. Sci. 2019, 20, 4848. [Google Scholar] [CrossRef] [PubMed]
  45. Wei, W.Q.; He, J.Z.; Luo, Y.W.; Yang, Z.; Xia, X.Y.; Han, Y.X.; He, N.J. Root-growth-related MaTCP transcription factors responsive to drought stress in mulberry. Forests 2023, 14, 143. [Google Scholar] [CrossRef]
  46. Liu, S.; Yin, X.C.; Feng, T.H.; Kang, Z.P.; Zhang, X.M.; Dong, J.N.; Liang, Z.S. Genome-wide identification and expression analysis of the TCP genes in Senna tora reveal the regulatory mechanism of their response to MeJA. Ind. Crop Prod. 2022, 187, 115511. [Google Scholar] [CrossRef]
  47. Nguyen-Queyrens, A.; Bouchet-Lannat, F. Osmotic adjustment in three-year-old seedlings of five provenances of maritime pine (Pinus pinaster) in response to drought. Tree Physiol. 2003, 23, 397–404. [Google Scholar] [CrossRef]
  48. Šņepste, I.; Krivmane, B.; Šķipars, V.; Zaluma, A.; Ruņģis, D.E. Induction of defense responses in Pinus sylvestris seedlings by methyl jasmonate and response to Heterobasidion annosum and Lophodermium seditiosum inoculation. Forests 2021, 12, 628. [Google Scholar] [CrossRef]
  49. Mukhopadhyay, P.; Tyagi, A.K. OsTCP19 influences developmental and abiotic stress signaling by modulating ABI4-mediated pathways. Sci. Rep. 2015, 5, 9998. [Google Scholar] [CrossRef]
  50. Finn, R.D.; Mistry, J.; Schuster-Böckler, B.; Griffiths-Jones, S.; Hollich, V.; Lassmann, T.; Moxon, S.; Marshall, M.; Khanna, A.; Durbin, R.; et al. Pfam: Clans, web tools and services. Nucleic Acids Res. 2006, 34, D247–D251. [Google Scholar] [CrossRef]
  51. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  52. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME Suite: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  54. Zhu, P.H.; Ma, Y.Y.; Zhu, L.Z.; Chen, Y.; Li, R.; Ji, K.S. Selection of suitable reference genes in Pinus massoniana Lamb. under different abiotic stresses for qPCR normalization. Forests 2019, 10, 632–649. [Google Scholar] [CrossRef]
  55. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29, 9. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 15938 g001
Figure 1. Phylogenetic tree of the TCP protein family members of P. massoniana, A. thaliana, P. abies, and P. taeda. The different color backgrounds and surrounding letters denote different groups. The magenta circles, purple circles, yellow circles, and cyan circles represent TCP TFs of P. massoniana, A. thaliana, P. abies, and P. taeda, respectively.
Figure 1. Phylogenetic tree of the TCP protein family members of P. massoniana, A. thaliana, P. abies, and P. taeda. The different color backgrounds and surrounding letters denote different groups. The magenta circles, purple circles, yellow circles, and cyan circles represent TCP TFs of P. massoniana, A. thaliana, P. abies, and P. taeda, respectively.
Ijms 24 15938 g001
Ijms 24 15938 g002
Figure 2. Multiple sequence alignments of conserved domains of 13 TCP family proteins in P. massoniana. The brown, light red, light blue, and dark blue backgrounds imply protein identities of 33%, 50%, 75%, and 100%, respectively. The HLH region and basic region are marked with black lines.
Figure 2. Multiple sequence alignments of conserved domains of 13 TCP family proteins in P. massoniana. The brown, light red, light blue, and dark blue backgrounds imply protein identities of 33%, 50%, 75%, and 100%, respectively. The HLH region and basic region are marked with black lines.
Ijms 24 15938 g002
Ijms 24 15938 g003
Figure 3. Phylogenetic tree, Conserved motif distribution, and structural domain analysis of the TCP proteins in P. massoniana. The motif structures were acquired by MEME analysis. Ten conserved motifs are presented of PmTCP TFs, and the different colors indicate different kinds of motifs.
Figure 3. Phylogenetic tree, Conserved motif distribution, and structural domain analysis of the TCP proteins in P. massoniana. The motif structures were acquired by MEME analysis. Ten conserved motifs are presented of PmTCP TFs, and the different colors indicate different kinds of motifs.
Ijms 24 15938 g003
Ijms 24 15938 g004
Figure 4. Subcellular localization experiments of PmTCP7 and PmTCP12 proteins. Transient expression of GFP (control), PmTCP7-GFP, and PmTCP12-GFP in tobacco leaves. The scale bar in the images of GFP is 50 μM, the scale bar in the images of PmTCP7-GFP is 5 μM, and the scale bar in the images of PmTCP12-GFP is 20 μM.
Figure 4. Subcellular localization experiments of PmTCP7 and PmTCP12 proteins. Transient expression of GFP (control), PmTCP7-GFP, and PmTCP12-GFP in tobacco leaves. The scale bar in the images of GFP is 50 μM, the scale bar in the images of PmTCP7-GFP is 5 μM, and the scale bar in the images of PmTCP12-GFP is 20 μM.
Ijms 24 15938 g004
Ijms 24 15938 g005
Figure 5. Transcriptional analysis of PmTCP genes in P. massoniana under four drought stresses: CK (80 ± 5)%, T1 (65 ± 5)%, T2 (50 ± 5)%, and T3 (35 ± 5)%. A heatmap was generated using the log2(FPKM+1) value, and the color scale denotes the relative expression level.
Figure 5. Transcriptional analysis of PmTCP genes in P. massoniana under four drought stresses: CK (80 ± 5)%, T1 (65 ± 5)%, T2 (50 ± 5)%, and T3 (35 ± 5)%. A heatmap was generated using the log2(FPKM+1) value, and the color scale denotes the relative expression level.
Ijms 24 15938 g005
Ijms 24 15938 g006
Figure 6. The expression levels of PmTCPs in eight different tissues in P. massoniana. TB: terminal bud; YS: young stems; OS: old stems; YN: young needles; ON: old needles; P: phloem; X: xylem; R: roots. The relative expression level is indicated as the mean ± standard deviation (SD). Asterisks show significant differences in the expression level between the treated groups and the control group (0 h) (** p < 0.01, *** p < 0.001, **** p < 0.0001).
Figure 6. The expression levels of PmTCPs in eight different tissues in P. massoniana. TB: terminal bud; YS: young stems; OS: old stems; YN: young needles; ON: old needles; P: phloem; X: xylem; R: roots. The relative expression level is indicated as the mean ± standard deviation (SD). Asterisks show significant differences in the expression level between the treated groups and the control group (0 h) (** p < 0.01, *** p < 0.001, **** p < 0.0001).
Ijms 24 15938 g006
Ijms 24 15938 g007aIjms 24 15938 g007b
Figure 7. The expression patterns of PmTCP genes under different treatments. (a) drought; (b) PEG; (c) mechanical injury; (d) H2O2; (e) ABA; (f) SA; (g) MeJA; (h) ETH. The relative expression is indicated as the mean ± standard deviation (SD). Asterisks show significant differences in the expression level between the treated groups and the control group (0 h) (* p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001).
Figure 7. The expression patterns of PmTCP genes under different treatments. (a) drought; (b) PEG; (c) mechanical injury; (d) H2O2; (e) ABA; (f) SA; (g) MeJA; (h) ETH. The relative expression is indicated as the mean ± standard deviation (SD). Asterisks show significant differences in the expression level between the treated groups and the control group (0 h) (* p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001).
Ijms 24 15938 g007aIjms 24 15938 g007b
Ijms 24 15938 g008
Figure 8. Transcriptional activity assay of PmTCP7. (a) PCR amplification of PmTCP7 gene; (b) Transcriptional activity analysis, pGBKT7, was used as a negative control, and pGBKT7-GAL4 was used as a positive control.
Figure 8. Transcriptional activity assay of PmTCP7. (a) PCR amplification of PmTCP7 gene; (b) Transcriptional activity analysis, pGBKT7, was used as a negative control, and pGBKT7-GAL4 was used as a positive control.
Ijms 24 15938 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, M.; Agassin, R.H.; Huang, Z.; Wang, D.; Yao, S.; Ji, K. Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress. Int. J. Mol. Sci. 2023, 24, 15938. https://doi.org/10.3390/ijms242115938

AMA Style

Zhang M, Agassin RH, Huang Z, Wang D, Yao S, Ji K. Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress. International Journal of Molecular Sciences. 2023; 24(21):15938. https://doi.org/10.3390/ijms242115938

Chicago/Turabian Style

Zhang, Mengyang, Romaric Hippolyte Agassin, Zichen Huang, Dengbao Wang, Sheng Yao, and Kongshu Ji. 2023. "Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress" International Journal of Molecular Sciences 24, no. 21: 15938. https://doi.org/10.3390/ijms242115938

APA Style

Zhang, M., Agassin, R. H., Huang, Z., Wang, D., Yao, S., & Ji, K. (2023). Transcriptome-Wide Identification of TCP Transcription Factor Family Members in Pinus massoniana and Their Expression in Regulation of Development and in Response to Stress. International Journal of Molecular Sciences, 24(21), 15938. https://doi.org/10.3390/ijms242115938

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop