Next Article in Journal
Exosomal miRNA Changes Associated with Restoration to Sinus Rhythm in Atrial Fibrillation Patients
Next Article in Special Issue
Intra-Phenotypic and -Genotypic Variations of Beauveria bassiana (Bals.) Vuill. Strains Infecting Aedes aegypti L. Adults
Previous Article in Journal
Targeting Vascular Impairment, Neuroinflammation, and Oxidative Stress Dynamics with Whole-Body Cryotherapy in Multiple Sclerosis Treatment
Previous Article in Special Issue
Genome-Wide Identification and Expression Pattern Analysis of BAHD Acyltransferase Family in Taxus mairei
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 7
10.3390/ijms25073859
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

Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis

by
Hang Shen
†,
Ying Hou
†,
Xiaorong Wang
,
Yaru Li
,
Jiasheng Wu
* and
Heqiang Lou
*
State Key Laboratory of Subtropical Silviculture, Zhejiang A&F University, Hangzhou 311300, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(7), 3859; https://doi.org/10.3390/ijms25073859
Submission received: 3 February 2024 / Revised: 10 March 2024 / Accepted: 21 March 2024 / Published: 29 March 2024
(This article belongs to the Special Issue Advances in Genetics and Phylogenomics of Tree)
Browse Figures
Review Reports Versions Notes

Abstract

:
The multidrug and toxin efflux (MATE) family participates in numerous biological processes and plays important roles in abiotic stress responses. However, information about the MATE family genes in Torreya grandis remains unclear. In this study, our genome-wide investigation identified ninety MATE genes in Torreya grandis, which were divided into five evolutionary clades. TgMATE family members are located on eleven chromosomes, and a total of thirty TgMATEs exist in tandem duplication. The promoter analysis showed that most TgMATEs contain the cis-regulatory elements associated with stress and hormonal responses. In addition, we discovered that most TgMATE genes responded to abiotic stresses (aluminum, drought, high temperatures, and low temperatures). Weighted correlation network analysis showed that 147 candidate transcription factor genes regulated the expression of 14 TgMATE genes, and it was verified through a double-luciferase assay. Overall, our findings offer valuable information for the characterization of the TgMATE gene mechanism in responding to abiotic stress and exhibit promising prospects for the stress tolerance breeding of Torreya grandis.

1. Introduction

Torreya grandis is an evergreen conifer species of the Taxaceae family. It is an endangered plant species and only grows in a few mountainous areas of some southern provinces of China [1]. As the native economic nut tree species, the fruit of Torreya grandis has been documented to have many medicinal properties and has been used as food for thousands of years in China [2]. Torreya grandis seed oil is abundant in unsaturated fatty acids and bioactive compounds, rendering it of excellent nutritional and medical value [3]. Due to the high-value components, Torreya grandis seed oil has been widely employed in medicine, food, and cosmetics in Japan and China. [4,5]. Nut crops are an integral component of food production systems, which are vital to nutrition and food security [6]. However, the frequent occurrence of abiotic stress environments in recent years has shown adverse effects on crop yields and thus poses a challenge to food security [7,8]. In recent years, acid rain and long-term improper fertilization have resulted in increased soil acidification and aluminum (Al) toxicity stress in Torreya grandis plantations, which adversely affects the yield and quality of Torreya grandis [9,10]. In addition, severe drought occurrences are common in Torreya grandis planting areas from the summer to autumn, causing Torreya grandis seedling growth and development to be hampered by water stress, which inhibits their growth and development [2]. Therefore, it is important to improve the resistance of Torreya grandis to various abiotic stresses. Nowadays, genetic engineering is widely used for crop improvement and provides a means to address the challenges of food security associated with environmental stress [11]. Numerous studies have shown that molecular breeding methods, such as genetic modification, can enhance plant stress tolerance traits [12,13]. Consequently, the breeding of Torreya grandis with greater abiotic stress tolerance is important for ensuring food security.
Over a long period of evolution, plants have evolved plenty of strategies to adapt to abiotic stress. When subjected to external stress, the stress-specific signal transduction pathways in plants will be triggered, thus activating relevant regulatory mechanisms to adapt themselves to the stress [14]. One of the regulatory mechanisms is the differential expression of stress-related genes [15]. The molecular processes of plant responses to abiotic challenges have been widely explored, and many genes critical for improving tolerance to various abiotic pressures have been found. The MATE gene family, which encodes multidrug and toxic compound extrusion transporters, is believed to be relevant to numerous abiotic stresses [16].
The multidrug and toxic compound extrusion (MATE) family is a newly classified family of multidrug efflux transporters that are ubiquitous in both prokaryotes and eukaryotes [16,17]. In plants, MATE transporters are engaged in a variety of biological processes [16]. They were discovered to have an essential role in the transport of metabolites, including flavonoids and alkaloids [17]. For example, Arabidopsis TT12 encoding the MATE protein was found to be able to transport flavonoids, which are involved in the vacuolar accumulation of proanthocyanidin precursors in the seed [18]. In addition, MATE transporters have been proven to be involved in the endogenous and exogenous mechanisms of detoxification [19] and the transport of plant hormones, such as efflux abscisic acid (ABA) [20] and salicylic acid (SA) [21]. At present, MATE family members have been discovered in many angiosperms, such as Oryza sativa [22], Solanum tuberosum [23], and Arabidopsis thaliana [24]. Studies in these species have revealed that plant MATE genes play a key role in the response to various abiotic stresses. In Arabidopsis thaliana, MATE family member AtDTX50 can promote cellular ABA efflux and consequently impact drought tolerance [20]. The barley MATE gene, HvAACT1, could accelerate the citrate efflux and improve the Al tolerance of the transgene plants [25]. However, the MATE family has not been discovered in gymnosperms to date. The structural features of MATE genes in Torreya grandis and their effects in the response to abiotic stress need to be further studied. Thankfully, the whole genome sequencing of Torreya grandis has been completed recently [26], which provides the basis for analyzing the MATE family in gymnosperms.
In this investigation, we aimed to conduct a comprehensive genome-wide analysis of the Torreya grandis MATE gene family. We systematically analyzed their phylogenetic relationships, gene structure, chromosome localization, potential cis-acting elements, and expression profiles in Al, drought, high-temperature, and low-temperature responses. In addition, we investigated potential transcription factors that regulate TgMATE members’ expression under abiotic stress through weighted correlation network analysis (WGCNA) and validated their reliability through a dual-luciferase assay. In summary, the response of TgMATE genes to abiotic stress was initially investigated in this research, which is expected to be helpful for improving the stress tolerance of Torreya grandis through genetic engineering.

2. Results

2.1. Identification and Characterization of MATE Genes in Torreya grandis

In our research, there are 90 MATE family members that were identified in Torreya grandis through HMM search and conserved domain analysis. We named those genes TgMATE1-TgMATE90 according to their chromosomal locations. To comprehend the physicochemical characteristics of TgMATE proteins, a bioinformatics evaluation was carried out. The outcomes indicated that TgMATE proteins have an amino acid count ranging from 75 to 755 aa. The molecular weight of these proteins varies from 8.253 to 83.36 kDa, while the isoelectric point (pI) ranges from 4.74 to 10.42. According to the subcellular localization prediction results, 69 TgMATE proteins are located on the plasma membrane, while less than a quarter of the TgMATEs are distributed on other organelles (Table S2).

2.2. Chromosomal Localization of TgMATEs

TBtools was used to process the genome annotation file and visualize the chromosome distribution of TgMATE genes [27]. The results showed that the TgMATEs gene was distributed on all 11 chromosomes of Torreya grandis. Chromosome 5 contains the largest number of TgMATE genes, which has 19 TgMATEs, while the least number of TgMATE genes distributed is chromosome 9, which has just 1 TgMATE gene (Figure 1). The processes of gene replication are crucial to the evolution of organisms [28]. Tandem duplication, one of the gene replication processes, is defined by the coexistence of family members in the same or nearby intergenic areas, and it might encourage the growth of gene families [29]. The tandem duplicates in the TgMATE gene family were investigated through the MCScanX tool [30], and those genes were marked with red arcs (Figure 1). Finally, there are a total of 30 tandem repeat TgMATE genes that were found, which constituted 20 tandem duplication events. The outcomes indicate that there were numerous tandem duplication events that occurred, which fueled the growth of the TgMATE genes in Torreya grandis.

2.3. Phylogenetic Analysis of TgMATE Proteins

A phylogenetic tree was constructed based on the protein sequences of 18 different species to evaluate the evolutionary relationship of the TgMATE gene family. These MATE proteins were categorized into five clades (I–V), which are illustrated in Figure 2. The TgMATE genes were distributed in all clades, and clade Ⅴ has the largest number of TgMATEs (49 TgMATEs). In clades I, Ⅱ, and Ⅲ, there are 24, 6, and 4 TgMATE genes, respectively. Clade Ⅳ (seven TgMATEs) contains the most abundant species, with a total of fourteen plant species of MATE proteins, indicating that they may have important functions and are relatively conservative in evolution. Interestingly, we found that many homologous MATE genes in clade Ⅳ have been reported to be associated with tolerance to Al stress, including TaMATE1B [31], SbMATE [32], VuMATE1 [33], VuMATE2 [34], OsFRDL4 [35], HvAACT1 [36], PtrMATE1, and PtrMATE2 [37]. These results imply that the TgMATE genes in clade Ⅳ may also have resistance to Al stress.

2.4. Investigation of Protein Motifs and Gene Structure of TgMATEs

We created the exon–intron structure diagram of TgMATE genes according to the Torreya grandis genome annotation file to investigate the structural variety of the TgMATE family members (Figure 3B). The exon number of the TgMATE gene ranges from 2 to 14. The bulk of the TgMATE genes in clades II, Ⅲ, and IV have quite lengthy intron sequences, whereas the majority of TgMATE genes in clades Ⅰ and V have practically short intron sequences. Noticeably, the TgMATE genes with close genetic distances have comparable exon–intron structures, whereas the structural patterns of TgMATE genes on different phylogenetic branches are distinct. Additionally, the evolutionary characteristics of TgMATE proteins were further explored through the conserved domain analysis in MEME, and 10 conserved motifs have been discovered (Figure 3C). Figure 3A shows the distribution of conserved motifs on TgMATE proteins. TgMATE proteins with nearly genetic distances have similar conserved motif categories and distributions. The protein motif composition of clades III and IV differed dramatically from that of clades I, II, and V. Members of clades III and IV, for example, have just 1–2 motifs in common, while TgMATE proteins in clades I, II, and V have 2–10 motifs in common, and most of them have 10 conserved motifs. Our findings are compatible with phylogenetic analysis, which supports the categorization of MATE family members.

2.5. Cis-Acting Elements Analysis of TgMATE Gene Promoters

After the prediction in PlantCARE, a total of 67 cis-acting elements with known functions were identified in the TgMATE gene promoters. According to their main functional characteristics, these elements were classified into six categories, encompassing the environmental stress-related, hormone-responsive, light-responsive, development-related, site-binding-related, and promoter-related elements (Table S3 and Figure 4). As the results demonstrate, the light-responsive elements were the most abundant cis-acting elements in TgMATE promoters, accounting for 31 different types, followed by the hormone-responsive elements, accounting for 11 cis-elements. There are seven categories of environmental stress-related elements that were detected, which were DRE, LTR, ARE, GC-motif, TC-rich repeats, MBS, and WUN-motif in the TgMATE promoters. The development-related elements include eight important types, such as cell cycle regulation, circadian control, etc. Notably, all the TgMATE promoters contain light-responsive cis-elements. In addition, the stress-response cis-elements exist in almost all TgMATE promoters, except for the promoters of TgMATE44, TgMATE69, TgMATE70, and TgMATE87 (Table S3). Moreover, 83 TgMATE promoters have hormone-responsive elements, and 59 TgMATE promoters have development-related elements (Table S3).

2.6. Analysis of the Expression Profiles of TgMATEs under Various Abiotic Stresses

To obtain some insights into how TgMATE genes respond to abiotic stress, we analyzed the transcript profile of TgMATEs based on the transcriptome data under Al, drought, high-temperature, and low-temperature stress. A total of 63 TgMATEs were found to have differential expression (Fold Change > 1.5 or <0.8) in the four types of abiotic stress treatment (Table S4). As shown in Figure 5A, the number of up-regulated and down-regulated TgMATE genes was on the same order of magnitude in Al and low-temperature stress (around 20 TgMATEs), while the quantity of down-regulated TgMATE genes is significantly greater than that of up-regulated TgMATE genes under drought and high-temperature stress. Next, we contrasted the overlapping relationships of the differentially expressed genes (DEGs) under different stress situations through the Venn diagram (Figure 5B). It is interesting to note that 18 TgMATE genes responded to Al, drought, high-temperature, and low-temperature stress, and 25 TgMATE genes responded to at least two stresses (Figure 5B), suggesting TgMATE family genes are extensively implicated in abiotic stress response processes in Torreya grandis.
Further, we constructed the gene expression heat maps to illustrate the expression profile of TgMATE genes under different stresses (Figure 5C–F). During Al stress, the patterns of these TgMATE genes’ expression can be divided into four groups (patterns Al_1–4) (Figure 5C). Expression pattern Al_1 exists in 12 TgMATE genes, where their expressions reached the highest levels at 8 h but decreased at 24 h after AlCl3 treatment (Figure 5C). Pattern Al_2 comprises 13 TgMATEs whose expression levels increased during the Al treatment and peaked at 24 h (Figure 5C). Pattern Al_3 contains nine TgMATEs, and their expression was lowest at 8 h (Figure 5C). Pattern Al_4 (14 TgMATEs) displayed the opposite expression pattern to pattern Al_2, with relatively high expression levels at 0 h of Al treatment but decreasing after 8 h and 24 h of Al treatment (Figure 5C). During drought stress, pattern Dr_1 (13 TgMATEs) exhibited comparatively high expression at 60 d compared to 0 d and 40 d (Figure 5D). Pattern Dr_2 (18 TgMATEs) was generally down-regulated at 40 d and 60 d (Figure 5D). Pattern Dr_3 (6 TgMATEs) was increased at 40 d but decreased at 60 d (Figure 5D). For high-temperature conditions, pattern HT_1 (17 TgMATEs) exhibited an up-regulation trend, while a majority of genes in pattern HT_2 showed significant down-regulation at 1 d, 2 d, and 4 d of high-temperature treatment (Figure 5E). Moreover, under low-temperature stress, pattern LT_1 consisting of 20 TgMATEs, was up-regulated at nearly all the time points (Figure 5F). Most genes in pattern LT_2 were down-regulated at 2 d, 4 d, and 6 d of low-temperature treatment, while in contrast, a few of them were up-regulated at 6 d of low-temperature treatment (Figure 5F). Pattern LT_3 (6 TgMATEs) positively responded to low-temperature stress in the early phases but was down-regulated during the later phases (Figure 5F). These findings highlight the possibility of these TgMATEs involved in Al, drought, high-temperature, and low-temperature stress responses and imply the different contributions of TgMATE family members to abiotic stress tolerance in Torreya grandis.

2.7. Expression Pattern of TgMATEs in Different Phylogenetic Clades under Abiotic Stresses

We also examined the expression patterns of each phylogenic clade of TgMATE family genes in the four types of abiotic stress. Clade I has a total of 24 TgMATEs, among which 13, 6, 9, and 9 genes respond in Al, drought, high-temperature, and low-temperature, with proportions of 54.17%, 25%, 37.5%, and 37.5%, respectively (Figure 6). There are two TgMATEs (33.33%) in clade II responding to Al, drought, and low-temperature stress, and five TgMATEs (83.33%) responding to high-temperature treatment (Figure 6). Clade III contains only four TgMATEs, with one gene (25.0%) responding in Al and low-temperature stress, three genes (75.0%) in drought stress, and all of them responding in high-temperature stress (Figure 6). Clade IV has more than half of TgMATEs that respond to each type of stress, with five (71.43%), four (57.14%), seven (100.0%), and seven (100.0%) genes in Al, drought, high-temperature, and low-temperature stress, respectively (Figure 6). Clade V contains the largest number of TgMATE members (49 in total), and 17, 20, 26, and 24 TgMATEs in this clade responded under the four treatments, accounting for 34.69%, 40.82%, 53.06%, and 48.98%, respectively (Figure 6). According to our results, different evolutionary branches of the TgMATE family exhibit varying proportions of response genes to abiotic stress, indicating that these branches may contribute differently to the abiotic stress response. Interestingly, clade IV has the highest proportion of TgMATEs responding in Al, high-temperature, and low-temperature treatments of the five polygenetic branches, implying that TgMATE members in this clade are widely involved in Torreya grandis resistance to abiotic stress. Additionally, clade Ⅲ showed the highest percentage of response genes to drought stress, suggesting that this clade may be more resistant to drought stress. Our findings demonstrate the broad involvement of the TgMATE family in the response to abiotic stress and imply that genes from various evolutionary branches may contribute differently to facing abiotic stress, which offers possible candidate genes for Torreya grandis stress-tolerant breeding.

2.8. QRT-PCR Validation

To validate the results of transcriptome analysis, we evaluated the expressions of the TgMATE family members in clade Ⅳ under Al treatment by qRT-PCR. Unsurprisingly, the expression patterns of these TgMATEs obtained from qRT-PCR almost matched the outcomes of the transcriptome study. During the Al stress treatment, the expressions of TgMATE1, TgMATE2, and TgMATE4 dramatically rose, reaching their peak levels at 24 h (Figure 7). TgMATE87 experienced a significant upregulation at 8 h while exhibiting a downregulation at 24 h (Figure 7). These results demonstrated the reliability of the expression pattern of TgMATEs under various abiotic stresses obtained through transcriptome data.

2.9. Identification of the Potential Transcription Factors by WGCNA

To investigate the reactions between TgMATEs and the transcript factors in Torreya grandis, a WGCNA was completed using the transcriptome data of Torreya grandis. A cluster tree was constructed according to the correlation between gene expression levels, and 11 different modules were generated by the dynamic cutting method (Figure 8A,B). Correlation weights calculated by WGCNA were used to measure the interaction degree between transcription factors and TgMATE genes. We chose a weight threshold of 0.01 as a criterion for identifying transcription factors that might interact with TgMATEs and used Cytoscape to visualize their co-expression relationship. According to the WGCNA analysis, a total of 147 transcription factors were found to have an association with 14 TgMATE family members (Figure 8C and Table S5). The transcription factors were arranged into 28 families, and the majority of these transcription factors belong to WRKY, bHLH, HD-ZIP, bZIP, NAC, and ERF families (Table S5), implying that these transcription factor families may play an important role in regulating TgMATE gene expression.

2.10. Verification of WGCNA Results through the Dual-Luciferase Reporter Assay

To verify the prediction results of WGCNA, TgMATE2 and the transcription factor predicted to be associated with it in the turquoise module were selected for the dual-luciferase reporter assay. Through PCR reaction, we successfully amplified the TgMATE2 promoter fragment and four transcription factors (evm.TU.PTG000304L.6, evm.TU.PTG009519L.1, evm.TU.PTG000725L.65, and evm.TU.PTG016833L.4) from Torreya sinensis DNA and cDNA. Their sequences were cloned into the reporters and the effectors, respectively, followed by the transient co-expressed in tobacco leaves through Agrobacterium tumefaciens. The empty 62-SK plasmid was used as the negative control. LUC and Rluc values were measured, and the interaction intensity of the TgMATE2 gene promoter region and the transcription factor proteins was reflected by the ratio of LUC to Rluc. As shown in Figure 9, the transcription factors in the co-expression network had different activation capacities for TgMATE2 promoters. Notably, transcript factor evm.TU.PTG016833L.4 from the HD-ZIP family, had a significantly higher LUC/Rluc ratio than the negative control (Figure 9B). In addition, the fluorescence of evm.TU.PTG016833L.4 detected through Tanon 5200 is brighter than that of the control (Figure 9C), implying this HDZIP transcription factor could activate the expression of TgMATE2. These findings support the validity of the approach for identifying transcription factors that control TgMATE members via WGCNA.

3. Discussion

3.1. Torreya Grandis Genome Possessed a Large TgMATE Gene Family

Our research discovered 90 TgMATE family members by searching the Torreya grandis genome with the hidden Markov model of the MATE protein. According to earlier studies, there are 56 MATE family genes in Arabidopsis [24], 46 in Oryza sativa [22], 42 in Cucumis melo [38], 48 in Solanum tuberosum [39], and 42 in Capsicum annuum [23]. Compared to these reported species, Torreya grandis has a substantially higher number of MATE family genes, suggesting their diversity and important role in plants. Gene replication is one of the main driving forces for the evolution and expansion of gene families [40]. Previous studies have shown that tandem and fragment replication mainly promote the replication of the MATE gene family [22,23,41]. The collinearity analysis of the Torreya grandis genome revealed there are 30 (33.4%) tandem repeat genes in the TgMATEs, which was close to the number of tandem repeat MATE genes in pepper and potato (38.1% and 35%, respectively) [23]. Therefore, we consider that tandem duplication may also promote the expansion of the MATE gene family in Torreya grandis. The results of subcellular localization prediction showed that the majority of TgMATEs, 69 in number, were located on the plasma membrane. This is compatible with the function of toxic compound efflux.
The length of MATE family proteins varies in different species. Arabidopsis MATE proteins contain 414 to 539 residues [24], and OsMATE proteins range in length from 370 to 598 aa [22]. MATE proteins in Gossypium arboreum range from 153 to 722 aa [42]. Our research discovered that the amino acid length of the TgMATE protein varied from 75 to 755, implying that the TgMATE family has higher diversity and complexity. According to phylogenetic analysis, MATE proteins from eighteen species could be divided into five clades, which is consistent with the MATE gene family previously reported in Capsicum annuum [23]. Further analysis reveals that the conserved motif distribution and gene structure of TgMATEs in different clusters were significantly distinct, while similar patterns were found among TgMATE genes in the same phylogenetic branch (Figure 3), which supports the classification results of MATE protein in the phylogenetic analysis. In summary, the differences among TgMATE family members suggest their diverse gene functions. Exon and intron play an important role in the differentiation of gene structure and function [42]. The divergence of gene structure could be primarily caused by the exonization of intronic sequences or pseudoexonization of exonic sequences [43]. We noticed TgMATE genes have very long intron lengths, especially in clades II, III, and IV. In tobacco, the intron length of MATE family genes is within 30 kb, and some members are intron-less [44]. While in Torreya grandis, all TgMATE genes have introns, and some are even more than 100 kb in length. This suggests that the expansion of TgMATE genes may be controlled by changes in the number and length of introns.

3.2. Members of the TgMATE Family Play an Essential Part to Response the ABIOTIC Stress in Torreya grandis

Nowadays, abiotic stressors, including Al, drought, high-temperature, and low-temperature stresses, are recognized as significant barriers to plant growth and agricultural productivity. Understanding the molecular response mechanism of plants under stress is helpful for the development of plant resistance breeding [45]. It is generally believed that one of the means for plants to cope with Al stress is to form complexes with extracellular aluminum ions through the secretion of citrate and other organic acid anions in the roots to alleviate Al toxicity stress [46]. Under drought stress, plants will accumulate a large amount of ABA and regulate the expression of downstream-related genes to resist drought [47]. In addition, plants will produce many reactive oxygen species under abiotic stress, resulting in tissue oxidative damage [48]. Numerous studies have shown that MATE transporters can respond to a variety of abiotic stresses through these mechanisms [15,31,32,49]. Nevertheless, little is known about how the TgMATE genes in Torreya grandis respond to abiotic stress. To date, transcriptome analysis of gene expression profiles has become one of the most important means of biological research and is generally considered to provide important clues for the exploration of gene functions [40]. Therefore, we investigated the expression pattern of TgMATE genes in response to Al, drought, high-temperature, and low-temperature stress based on the transcriptomic data. We found that more than half of TgMATE genes (70%) were differentially expressed under abiotic stress. The number of up-regulated genes in Al and low-temperature stress is almost equal to that of down-regulated genes, whereas the number of up-regulated genes is less than that of down-regulated genes in drought and high-temperature stress, suggesting that TgMATE genes tend to be down-regulated in response to drought and high-temperature stress in contrast to Al and low-temperature stress. Under Al stress, we found that 19 TgMATEs were up-regulated and 22 TgMATEs were down-regulated. There is some evidence showing that MATE transporters are responsible for the secretion of organic anion citrate, which affects the Al tolerance of plants. For example, the introgression of bread wheat TaMATE1B into durum wheat enables better growth in the soil with high Al3+ [31]. The transgenic sugarcanes with the SbMATE gene overexpressed have an improved tolerance to Al [32]. These data suggested that TgMATE family genes might be involved in Al stress response. The Arabidopsis homologous gene DTX50 was found to encode a transporter with ABA effector function and thus play a significant role in response to drought stress [48]. In our research, 15 TgMATEs were up-regulated and 22 TgMATEs were down-regulated. Under drought treatment, suggesting the possibility of these TgMATEs responding to drought stress. high-temperature and low-temperature stress are common abiotic stresses that accumulate reactive oxygen species and cause damage to plants. Previous studies have shown that the Cotton DTX/MATE gene accelerates the expression of antioxidant genes, which reduces the content of reactive oxygen species and enhances abiotic stress tolerance in transgenic Arabidopsis [15]. In Torreya grandis, there are 11 and 24 genes that were up-regulated, and 44 and 22 genes were down-regulated under high-temperature and low-temperature stress, respectively, indicating their important roles in responding to extreme temperature. Notably, 25 TgMATEs responded to at least two stresses, and 18 TgMATEs responded to all four types of stress, indicating that TgMATE family members actively participate in the abiotic stress response mechanisms in Torreya grandis. Moreover, hierarchical clustering results showed that TgMATEs had different expression patterns under four different stresses, implying that the TgMATE genes may contribute differently to cope with these abiotic stresses. In summary, these findings suggest that the TgMATE genes are extensively associated with the abiotic stress response in Torreya grandis and could provide the molecular foundation for the breeding of Torreya grandis cultivars with improved abiotic stress tolerance.

3.3. Expression Profiles of Different TgMATE Subfamilies under Abiotic Stress

During evolution, different phylogenic clades within a gene family typically acquire distinct functions to enhance their adaptability to changing environmental conditions [50]. In Arabidopsis and rice, the subfamily MATE III has distinct physiochemical properties that are different from the other subfamilies, as Type-II divergence only occurred in the MATE III subfamily, while there were no Type-II functional divergence sites in any other subfamilies [51]. Nevertheless, little research has been conducted on how MATE family members from different phylogenic clades react to diverse abiotic stresses. Here, we analyzed the expression patterns of TgMATE family members during abiotic stress to explore their functional differentiation in abiotic stress responses. TgMATE family members from five evolutionary clades have different numbers of responses across the four stress types. For example, clade IV has the highest percentage of TgMATE genes responding in Al, high-temperature, and low-temperature stress, which suggests that TgMATEs in this clade are extensively implicated in plant resistance to abiotic stress. Notably, clade Ⅳ also combines many MATE homolog genes that have been confirmed in other plants to be tolerant to Al resistance, e.g., VuMATE1 [33], etc., which highlights the important roles of TgMATEs in clade IV in the abiotic stress. Moreover, we found clade Ⅲ had the highest proportion of response genes under drought stress, indicating this clade may tend to resist drought stress. Our findings suggest that different TgMATE family evolution branches may contribute differentially to the abiotic stress tolerance in Torreya grandis.

3.4. Mining of Candidate Transcription Factors That Interact with TgMATEs

Transcription factors control gene expression at the transcriptional level by attaching to cis-acting elements on target promoters [52]. In plants, transcription factors are essential in numerous biological activities such as the regulation of development, the induction of defense, and the stress response [53]. Recent studies have reported that transcription factors regulate the expression of the MATE gene in other plants. In rice, Yokosho et al. [54] discovered that a c2h2 type zinc finger Al tolerance transcription factor ART1 controls the expression of the rice MATE gene (OsFRDL4). However, the transcription factors that regulate the family members of TgMATEs in Torreya grandis have not yet been discovered. In this study, transcription factors co-expressed with TgMATE genes were predicted through WGCNA analysis. A total of 147 transcription factors were discovered to potentially regulate the expression of TgMATE genes and their relationship network was constructed. Many of these transcription factors belong to the WRKY, bHLH, HDZIP, bZIP, NAC, and ERF families, which were believed to be widely involved in plant responses to stress. For instance, wheat TaWRKY10 could regulate the expression of a series of stress-related genes to improve the stress resistance of tobacco [55]. In Vigna umbellate, an NAC-type transcription VuNAR1 could regulate the Al resistance through interaction with the cell wall-associated receptor gene promoter to regulate the metabolism of cell wall pectin [56]. The two HDZIP family members, AtHB7 and AtHB12, oppositely control Al resistance by affecting the root cell wall Al accumulation in Arabidopsis thaliana [57]. Therefore, we speculated that the transcription factors screened by WGCNA analysis could respond to abiotic stress and regulate the expression of TgMATE genes, thus affecting the stress tolerance of Torreya grandis. Furthermore, the dual-luciferase reporter assay showed that a transcription factor (evm.TU.PTG016833L.4) from the HDZIP family in the co-expression network can positively regulate the expression of TgMATE2, suggesting the credibility in mining the transcription factors of TgMATEs through WGCNA. The potential transcription factors predicted by WGCNA can provide a molecular basis for further exploring the mechanism by which TgMATE genes respond to abiotic stress. Above all, the study of transcription factors related to TgMATE genes has important biological significance; it could not only enrich our knowledge of the mechanism of Torreya grandis in response to environmental stress but also contribute to providing new theoretical guidance for the genetic improvement of Torreya grandis in stress conditions.

4. Methods

4.1. Identification and Characterization of the MATE Family Members in Torreya grandis

The entire genome sequence of Torreya grandis and the annotation file were downloaded from Figshare (https://doi.org/10.6084/m9.figshare.21089869.v1, accessed on 29 April 2023) [58]. To find the potential TgMATE genes, we obtained the Hidden Markov model file of MATE protein (PF01554) in PFAM (https://pfam-legacy.xfam.org, accessed on 29 April 2023) and used the HMMER3.0 software to employ an HMM search (E-value < 1 × 10−5) on the Torreya grandis proteomes. The sequences of those putative TgMATE proteins were validated further based on their MATE domain, which was discovered through the NCBI web tool CDD-Search (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi, accessed on 29 April 2023). The isoelectric point (PI) and molecular weight (MW) were determined via the ExPASY online tools (https://web.expasy.org/protparam/, accessed on 8 May 2023) [59]. We used the TargetP 2.0 web tool (https://services.healthtech.dtu.dk/services/TargetP-2.0, accessed on 8 May 2023) [60] to forecast the protein subcellular location of each TgMATE.

4.2. Phylogenetic Analysis

Homologous MATE protein sequences in other plants were obtained from NCBI. The proteins of 18 species were used for the phylogenetic analysis, including 11 dicotyledonous plants and 6 monocotyledons in addition to Torreya grandis. We aligned all protein sequences through ClustalW [61], and the phylogenetic tree was constructed based on the Maximin Likelihood method with the bootstrap value 1000 in IQTREE software v1.6.12 [62]. The phylogenetic relationship of these MATE proteins was visualized through ITOL (https://www.itol.org/, accessed on 11 July 2023) [63].

4.3. Chromosomal Location, Gene Structure, and Conserved Motifs Analysis

We used TBtools software v2.042 to extract the exons, introns, and UTR information of TgMATEs from the Torreya grandis genome annotation file [27]. The 10 most conserved motifs of TgMATE proteins were detected through MEME SUITE (https://meme-suite.org/meme/, accessed on 13 May 2023) [64]. The chromosome distribution records of TgMATEs were derived via the genome annotation file. Finally, these outcomes were visualized via the TBtools software v2.042 [27].

4.4. Analysis of the Promoter Cis-Acting Elements of TgMATEs

Based on the genomic annotation file, the promoter sequence located in the upper 2000 bp of the coding sequence of the TgMATEs was extracted. The cis-acting elements of their promoters were predicted through PlantCARE websites (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 9 May 2023) [65], and the results were visualized through TBtools software v2.042 [27].

4.5. Transcriptomic Analyses of TgMATE Genes

Our lab has previously constructed a transcriptome database of 1-year-old Torreya grandis at various abiotic stresses. Here, we examined the expression of TgMATE family members under various abiotic stressors using these transcriptomes of Torreya grandis under Al stress (50 µM AlCl3 treatment for 0, 8, and 24 h), high-temperature stress (40 °C high-temperature treatment for 0, 1, 2, and 3 days), low-temperature stress (2 °C low-temperature treatment for 0, 1, 2, 4, and 6 days), and drought stress (drought without watering for 0, 40, and 60 days). The Fragments per Kilobase Million (FPKM) data of TgMATE family genes from the transcriptomes were used to generate a gene expression matrix of TgMATEs. The up-regulated gene was defined as having a Fold Change value greater than 1.5, and the Fold Change of the down-regulated gene was less than 0.8. Finally, TBtools software v2.042 was utilized to visualize the gene expression heat map of TgMATEs [27].

4.6. Plant Materials and Al Stress Treatment

One-year-old Torreya grandis seedlings (bought from Zhuji City, China) were rinsed with sterile water and subsequently transferred to a container filled with 1/5 of Hoagland’s nutrients [66]. During the adaptation period, the seedlings were cultured for 15 days, and Hoagland’s solution was changed every 3 days. After that, they were moved to a 1/30 Hoagland nutritional solution [66] with 50 μM AlCl3 and 1 mM CaCl2 (pH 5.0). Root tips of 0–1 cm length of Torreya grandis seedlings were taken after treatment at 0, 8, and 24 h. For each treatment, three biological replicates were carried out.

4.7. RNA Extraction and qRT-PCR Validation

An RNAprep Pure Plant Plus Kit (Tiangen Biotech, Beijing, China) was used to extract total plant RNA from Torreya grandis seedling samples, which was then reverse transcribed into cDNA via a FastQuant RT Kit (with gDNase) (Tiangen Biotech, Beijing, China). Based on the findings of phylogenetic analysis, 7 TgMATE genes in clade Ⅳ that may have the ability to resist Al toxicity were chosen for qRT-PCR validation. The fluorescent quantitative primers used for the qRT-PCR are detailed in Table S1, and the reaction steps were performed according to the method of Luo et al. [67]. There are three technical duplications that were conducted for each sample, and the relative expressions were calculated according to a 2−ΔΔCt algorithm [68].

4.8. Weighted Gene Co-Expression Network Analysis

Transcription factors in Torreya grandis were forecasted through the PlantTFD (http://planttfdb.gao-lab.org/, accessed on 26 June 2023) [69]. To establish the gene co-expression networks, the WGCNA shiny program in TBtools was used [27]. The gene set is first filtered to ensure that every gene participating in the WGCNA has an expression level greater than 1 in at least 90% of the samples. A soft threshold power value of 14 was chosen to build a network. The minimum module size was chosen as 50, and the threshold for similar module merging was set to 0.35. Finally, the Cytoscape software v3.7.2 was used for the visualization of the interested module [70].

4.9. Dual-Luciferase Reporter Assay

A transient dual-luciferase experiment was performed to examine the transactivation activity of the co-expression transcription factors in the promoters of TgMATE2. Whole coding sequences of transcription factors were cloned and placed in the pGreenII 62-SK plasmid as the effector, while the promoter region of TgMATE2 was placed in the pGreenII 0800-LUC plasmid as the reporter. Primers used in the PCR process are detailed in Table S1. The effector and reporter were then converted into the Agrobacterium tumefaciens receptive bacteria strain GV3101 (with psoup) (Weidi, Shanghai, China) by the chemical transformation methods. After overnight incubation (28°, 200 rpm) in Luria–Bertani medium containing kanamycin and rifampicin, the target strains were collected and then re-suspended with infection buffers (10 mM morph ethanesulfonic acid, 10 mM MgCl2, 50 mM acetyl syringone, and pH 5.6) and the optical density was adjusted to 0.8. Consequently, the effector and the reporter were injected into the 5-week-old tobacco leaves at a ratio of 10 to 1. After 2 days of culture in the artificial climate chamber, intensities of firefly luciferase (LUC) and renila luciferase (Rluc) were detected according to the methods as demonstrated before [53]. A fully automated Tanon 5200 chemiluminescent image analyzer was used to observe fluorescence.

5. Conclusions

In this study, a comprehensive analysis of the MATE gene family in Torreya grandis was conducted. A total of ninety MATE genes were identified in the Torreya grandis genome, which can be phylogenetically classified into five clades. The conserved motifs and gene structure analysis revealed significant conservation of TgMATEs within each clade, which confirmed our classification scheme. The tandem duplication plays an important role in the expansion of the TgMATE gene family. Promoter analysis revealed that the majority of TgMATEs contain hormone-response and environmental stress-related elements. The expression profiles of MATE family members in Torreya grandis under Al, drought, high-temperature, and low-temperature stress confirmed that TgMATEs played critical roles in response to abiotic stress. Furthermore, potential transcription factors interacting with TgMATEs were identified by WGCNA analysis, and the dual-luciferase assay validated their credibility. Taken together, these findings provide a comprehensive and systematic characterization of the MATE gene family in Torreya grandis and provide a molecular basis for the stress-tolerant breeding of Torreya grandis.

Supplementary Materials

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

Author Contributions

Conceptualization, J.W. and H.L.; Methodology, H.S., Y.H., X.W., and Y.L.; Data curation and Formal analysis, H.S.; Resources, J.W. and H.L.; Writing—original draft preparation, H.S.; Writing—review and editing, H.S., Y.H., and H.L.; Visualization, H.S.; Supervision, H.L.; Funding acquisition, J.W. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Breeding of New Varieties of Torreya grandis Program (2021C02066-11), and the Scientific Research Startup Fund Project of Zhejiang A&F University (2018FR028).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MATE, multidrug and toxic compound extrusion transporters; PCR, Polymerase Chain Reaction; qRT-PCR quantitative real-time PCR; Al, aluminum; Dr, drought; HT, high temperature; LT, low temperature; WGCNA, weighted correlation network analysis; LUC, firefly luciferase; Rluc, renila luciferase.

References

  1. Zeng, J.; Chen, J.; Kou, Y.; Wang, Y. Application of EST-SSR markers developed from the transcriptome of Torreya grandis (Taxaceae), a threatened nut-yielding conifer tree. PeerJ 2018, 6, e5606. [Google Scholar] [CrossRef] [PubMed]
  2. Shen, C.H.; Hu, Y.Y.; Du, X.H.; Li, T.T.; Tang, H.; Wu, J.S. Salicylic acid induces physiological and biochemical changes in Torreya grandis cv. Merrillii seedlings under drought stress. Trees 2014, 28, 961–970. [Google Scholar] [CrossRef]
  3. He, Z.; Zhu, H.; Li, W.; Zeng, M.; Wu, S.; Chen, S.; Qin, F.; Chen, J. Chemical components of cold pressed kernel oils from different Torreya grandis cultivars. Food Chem. 2016, 209, 196–202. [Google Scholar] [CrossRef] [PubMed]
  4. Feng, T.; Hu, Z.; Song, S.; Yao, L.; Sun, M.; Zhu, X.; Lu, J. The antioxidant and tyrosinase inhibition properties of essential oil from the peel of Chinese Torreya grandis Fort. RSC Adv. 2019, 9, 42360–42366. [Google Scholar] [CrossRef] [PubMed]
  5. Shi, L.K.; Mao, J.H.; Zheng, L.; Zhao, C.W.; Jin, Q.Z.; Wang, X.G. Chemical characterization and free radical scavenging capacity of oils obtained from Torreya grandis Fort. ex. Lindl. and Torreya grandis Fort. var. Merrillii: A comparative study using chemometrics. Ind. Crops Prod. 2018, 115, 250–260. [Google Scholar] [CrossRef]
  6. Savadi, S.; Mangalassery, S.; Sandesh, M.S. Advances in genomics and genome editing for breeding next generation of fruit and nut crops. Genomics 2021, 113, 3718–3734. [Google Scholar] [CrossRef] [PubMed]
  7. Rahman, M.; Zulfiqar, S.; Raza, M.A.; Ahmad, N.; Zhang, B. Engineering Abiotic Stress Tolerance in Crop Plants through CRISPR Genome Editing. Cells 2022, 11, 3590. [Google Scholar] [CrossRef] [PubMed]
  8. Kumar, S. Abiotic Stresses and Their Effects on Plant Growth, Yield and Nutritional Quality of Agricultural Produce. Int. J. Food Sci. Agric. 2020, 4, 367–378. [Google Scholar] [CrossRef]
  9. Xia, L.; Shao, C.; Zhang, N.; Wu, A.; Xie, J.; Qiu, Y.; He, X.; Pei, J.; Wang, X.; Wang, Y. Improved Tolerance of Mycorrhizal Torreya grandis Seedlings to Sulfuric Acid Rain Related to Phosphorus and Zinc Contents in Shoots. J. Fungi 2021, 7, 296. [Google Scholar] [CrossRef]
  10. Zhou, Z.; Yrjälä, K.; Chen, J.; Yu, C.; Shi, W.; Qin, H.; Yu, W.; Dai, W.; Hu, Y.; Wu, J. Organic amendments combined with biochar for improving soil and plant quality in a Torreya grandis plantation. J. Soils Sediments 2022, 22, 1080–1094. [Google Scholar] [CrossRef]
  11. Ouyang, B.; Gu, X.; Holford, P. Plant genetic engineering and biotechnology: A sustainable solution for future food security and industry. Plant Growth Regul. 2017, 83, 171–173. [Google Scholar] [CrossRef]
  12. James, D.; Borphukan, B.; Fartyal, D.; Ram, B.; Singh, J.; Manna, M.; Sheri, V.; Panditi, V.; Yadav, R.; Achary, V.M.M.; et al. Concurrent Overexpression of OsGS1;1 and OsGS2 Genes in Transgenic Rice (Oryza sativa L.): Impact on Tolerance to Abiotic Stresses. Front. Plant Sci. 2018, 9, 786. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, Y.; Xue, Y.; Li, J. Towards molecular breeding and improvement of rice in China. Trends Plant Sci. 2005, 10, 610–614. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, Y.; Xu, J.; Li, R.; Ge, Y.; Li, Y.; Li, R. Plants’ Response to Abiotic Stress: Mechanisms and Strategies. Int. J. Mol. Sci. 2023, 24, 10915. [Google Scholar] [CrossRef] [PubMed]
  15. Lu, P.; Magwanga, R.O.; Kirungu, J.N.; Hu, Y.; Dong, Q.; Cai, X.; Zhou, Z.; Wang, X.; Zhang, Z.; Hou, Y.; et al. Overexpression of Cotton a DTX/MATE Gene Enhances Drought, Salt, and Cold Stress Tolerance in Transgenic Arabidopsis. Front. Plant Sci. 2019, 10, 299. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, J.; Li, Y.; Wang, W.; Gai, J.; Li, Y. Genome-wide analysis of MATE transporters and expression patterns of a subgroup of MATE genes in response to aluminum toxicity in soybean. BMC Genom. 2016, 17, 223. [Google Scholar] [CrossRef] [PubMed]
  17. Takanashi, K.; Shitan, N.; Yazaki, K. The multidrug and toxic compound extrusion (MATE) family in plants. Plant Biotechnol. 2014, 31, 417–430. [Google Scholar] [CrossRef]
  18. Marinova, K.; Pourcel, L.; Weder, B.; Schwarz, M.; Barron, D.; Routaboul, J.M.; Debeaujon, I.; Klein, M. The Arabidopsis MATE Transporter TT12 Acts as a Vacuolar Flavonoid/H+-Antiporter Active in Proanthocyanidin-Accumulating Cells of the Seed Coat. Plant Cell 2007, 19, 2023–2038. [Google Scholar] [CrossRef] [PubMed]
  19. Xu, L.; Shen, Z.L.; Chen, W.; Si, G.Y.; Meng, Y.; Guo, N.; Sun, X.; Cai, Y.P.; Lin, Y.; Gao, J.S. Phylogenetic analysis of upland cotton MATE gene family reveals a conserved subfamily involved in transport of proanthocyanidins. Mol. Biol. Rep. 2019, 46, 161–175. [Google Scholar] [CrossRef]
  20. Zhang, H.; Zhu, H.; Pan, Y.; Yu, Y.; Luan, S.; Li, L. A DTX/MATE-Type Transporter Facilitates Abscisic Acid Efflux and Modulates ABA Sensitivity and Drought Tolerance in Arabidopsis. Mol. Plant 2014, 7, 1522–1532. [Google Scholar] [CrossRef]
  21. Serrano, M.; Wang, B.; Aryal, B.; Garcion, C.; Abou-Mansour, E.; Heck, S.; Geisler, M.; Mauch, F.; Nawrath, C.; Métraux, J.P. Export of Salicylic Acid from the Chloroplast Requires the Multidrug and Toxin Extrusion-Like Transporter EDS5. Plant Physiol. 2013, 162, 1815–1821. [Google Scholar] [CrossRef] [PubMed]
  22. Du, Z.; Su, Q.; Wu, Z.; Huang, Z.; Bao, J.; Li, J.; Tu, H.; Zeng, C.; Fu, J.; He, H. Genome-wide characterization of MATE gene family and expression profiles in response to abiotic stresses in rice (Oryza sativa). BMC Ecol. Evol. 2021, 21, 141. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, Q.; Wang, L.; Liu, D.; Ma, S.; Dai, Y.; Zhang, X.; Wang, Y.; Hu, T.; Xiao, M.; Zhou, Y.; et al. Identification and Expression of the Multidrug and Toxic Compound Extrusion (MATE) Gene Family in Capsicum annuum and Solanum tuberosum. Plants 2020, 9, 1448. [Google Scholar] [CrossRef] [PubMed]
  24. Li, L.; He, Z.; Pandey, G.K.; Tsuchiya, T.; Luan, S. Functional Cloning and Characterization of a Plant Efflux Carrier for Multidrug and Heavy Metal Detoxification. J. Biol. Chem. 2002, 277, 5360–5368. [Google Scholar] [CrossRef] [PubMed]
  25. Furukawa, J.; Yamaji, N.; Wang, H.; Mitani, N.; Murata, Y.; Sato, K.; Katsuhara, M.; Takeda, K.; Ma, J.F. An Aluminum-Activated Citrate Transporter in Barley. Plant Cell Physiol. 2007, 48, 1081–1091. [Google Scholar] [CrossRef] [PubMed]
  26. Lou, H.Q.; Song, L.L.; Li, X.L.; Zi, H.L.; Chen, W.J.; Gao, Y.D.; Zheng, S.; Fei, Z.J.; Sun, X.P.; Wu, J.S. The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis. Nat. Commun. 2023, 14, 1315. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  28. Kong, H.; Landherr, L.L.; Frohlich, M.W.; Leebens-Mack, J.; Ma, H.; DePamphilis, C.W. Patterns of gene duplication in the plant SKP1 gene family in angiosperms: Evidence for multiple mechanisms of rapid gene birth. Plant J. 2007, 50, 873–885. [Google Scholar] [CrossRef]
  29. Ramamoorthy, R.; Jiang, S.Y.; Kumar, N.; Venkatesh, P.N.; Ramachandran, S. A Comprehensive Transcriptional Profiling of the WRKY Gene Family in Rice Under Various Abiotic and Phytohormone Treatments. Plant Cell Physiol. 2008, 49, 865–879. [Google Scholar] [CrossRef]
  30. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef]
  31. Pooniya, V.; Palta, J.A.; Chen, Y.; Delhaize, E.; Siddique, K.H.M. Impact of the TaMATE1B gene on above and below-ground growth of durum wheat grown on an acid and Al3+ toxic soil. Plant Soil 2020, 447, 73–84. [Google Scholar] [CrossRef]
  32. Ribeiro, A.P.; Vinecky, F.; Duarte, K.E.; Santiago, T.R.; das Chagas Noqueli Casari, R.A.; Hell, A.F.; da Cunha, B.A.D.B.; Martins, P.K.; da Cruz Centeno, D.; de Oliveira Molinari, P.A.; et al. Enhanced aluminum tolerance in sugarcane: Evaluation of SbMATE overexpression and genome-wide identification of ALMTs in Saccharum spp. BMC Plant Biol. 2021, 21, 300. [Google Scholar] [CrossRef] [PubMed]
  33. Liu, M.Y.; Xu, J.M.; Lou, H.Q.; Fan, W.; Yang, J.L.; Zheng, S.J. Characterization of VuMATE1 Expression in Response to Iron Nutrition and Aluminum Stress Reveals Adaptation of Rice Bean (Vigna umbellata) to Acid Soils through Cis Regulation. Front. Plant Sci. 2016, 7, 511. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, M.Y.; Lou, H.Q.; Chen, W.W.; Piñeros, M.A.; Xu, J.M.; Fan, W.; Kochian, L.V.; Zheng, S.J.; Yang, J.L. Two citrate transporters coordinately regulate citrate secretion from rice bean root tip under aluminum stress. Plant Cell Environ. 2018, 41, 809–822. [Google Scholar] [CrossRef] [PubMed]
  35. Yokosho, K.; Yamaji, N.; Fujii-Kashino, M.; Ma, J.F. Retrotransposon-Mediated Aluminum Tolerance through Enhanced Expression of the Citrate Transporter OsFRDL4. Plant Physiol. 2016, 172, 2327–2336. [Google Scholar] [CrossRef] [PubMed]
  36. Zhou, G.; Delhaize, E.; Zhou, M.; Ryan, P.R. The barley MATE gene, HvAACT1, increases citrate efflux and Al3+ tolerance when expressed in wheat and barley. Ann. Bot. 2013, 112, 603–612. [Google Scholar] [CrossRef] [PubMed]
  37. Li, N.; Meng, H.; Xing, H.; Liang, L.; Zhao, X.; Luo, K. Genome-wide analysis of MATE transporters and molecular characterization of aluminum resistance in Populus. J. Exp. Bot. 2017, 68, 5669–5683. [Google Scholar] [CrossRef] [PubMed]
  38. Shah, I.H.; Manzoor, M.A.; Sabir, I.A.; Ashraf, M.; Haq, F.; Arif, S.; Abdullah, M.; Niu, Q.; Zhang, Y. Genome-wide identification and comparative analysis of MATE gene family in Cucurbitaceae species and their regulatory role in melon (Cucumis melo) under salt stress. Hortic. Environ. Biotechnol. 2022, 63, 595–612. [Google Scholar] [CrossRef]
  39. Li, Y.; He, H.; He, L.F. Genome-wide analysis of the MATE gene family in potato. Mol. Biol. Rep. 2019, 46, 403–414. [Google Scholar] [CrossRef]
  40. Zhu, H.; Wu, J.; Jiang, Y.; Jin, J.; Zhou, W.E.I.; Wang, Y.U.; Han, G.; Zhao, Y.; Cheng, B. Genomewide analysis of MATE-type gene family in maize reveals microsynteny and their expression patterns under aluminum treatment. J. Genet. 2016, 95, 691–704. [Google Scholar] [CrossRef]
  41. dos Santos, A.L.; Chaves-Silva, S.; Yang, L.; Maia, L.G.S.; Chalfun-Júnior, A.; Sinharoy, S.; Zhao, J.; Benedito, V.A. Global analysis of the MATE gene family of metabolite transporters in tomato. BMC Plant Biol. 2017, 17, 185. [Google Scholar] [CrossRef] [PubMed]
  42. Lu, P.; Magwanga, R.O.; Guo, X.; Kirungu, J.N.; Lu, H.; Cai, X.; Zhou, Z.; Wei, Y.; Wang, X.; Zhang, Z.; et al. Genome-Wide Analysis of Multidrug and Toxic Compound Extrusion (MATE) Family in Gossypium raimondii and Gossypium arboreum and Its Expression Analysis Under Salt, Cadmium, and Drought Stress. G3 Genes Genomes Genet. 2018, 8, 2483–2500. [Google Scholar] [CrossRef] [PubMed]
  43. Fan, K.; Wang, M.; Miao, Y.; Ni, M.; Bibi, N.; Yuan, S.; Li, F.; Wang, X. Molecular Evolution and Expansion Analysis of the NAC Transcription Factor in Zea mays. PLoS ONE 2014, 9, e111837. [Google Scholar] [CrossRef] [PubMed]
  44. Gani, U.; Sharma, P.; Tiwari, H.; Nautiyal, A.K.; Kundan, M.; Wajid, M.A.; Kesari, R.; Nargotra, A.; Misra, P. Comprehensive genome-wide identification, characterization, and expression profiling of MATE gene family in Nicotiana tabacum. Gene 2021, 783, 145554. [Google Scholar] [CrossRef] [PubMed]
  45. Jang, M.G.; Kim, Y.J.; Jang, G.H.; Sukweenadhi, J.; Kwon, W.S.; Yang, D.C. Ectopic overexpression of the aluminum-induced protein gene from Panax ginseng enhances heavy metal tolerance in transgenic Arabidopsis. Plant Cell Tissue Organ Cult. 2014, 119, 95–106. [Google Scholar] [CrossRef]
  46. Kochian, L.V.; Piñeros, M.A.; Liu, J.; Magalhaes, J.V. Plant Adaptation to Acid Soils: The Molecular Basis for Crop Aluminum Resistance. Annu. Rev. Plant Biol. 2015, 66, 571–598. [Google Scholar] [CrossRef] [PubMed]
  47. Yoshida, T.; Mogami, J.; Yamaguchi-Shinozaki, K. ABA-dependent and ABA-independent signaling in response to osmotic stress in plants. Curr. Opin. Plant Biol. 2014, 21, 133–139. [Google Scholar] [CrossRef] [PubMed]
  48. Zandalinas, S.I.; Mittler, R.; Balfagón, D.; Arbona, V.; Gómez-Cadenas, A. Plant adaptations to the combination of drought and high temperatures. Physiol. Plant. 2018, 162, 2–12. [Google Scholar] [CrossRef] [PubMed]
  49. Remy, E.; Cabrito, T.R.; Baster, P.; Batista, R.A.; Teixeira, M.C.; Friml, J.; Sá-Correia, I.; Duque, P. A Major Facilitator Superfamily Transporter Plays a Dual Role in Polar Auxin Transport and Drought Stress Tolerance in Arabidopsis. Plant Cell 2013, 25, 901–926. [Google Scholar] [CrossRef]
  50. Qin, J.; Jiang, Y.; Lu, Y.; Zhao, P.; Wu, B.; Li, H.; Wang, Y.; Xu, S.; Sun, Q.; Liu, Z. Genome-wide identification and transcriptome profiling reveal great expansion of SWEET gene family and their wide-spread responses to abiotic stress in wheat (Triticum aestivum L.). J. Integr. Agric. 2020, 19, 1704–1720. [Google Scholar] [CrossRef]
  51. Wang, L.; Bei, X.; Gao, J.; Li, Y.; Yan, Y.; Hu, Y. The similar and different evolutionary trends of MATE family occurred between rice and Arabidopsis thaliana. BMC Plant Biol. 2016, 16, 207. [Google Scholar] [CrossRef] [PubMed]
  52. Zeng, D.; da Silva, J.A.; Zhang, M.; Yu, Z.; Si, C.; Zhao, C.; Dai, G.; He, C.; Duan, J. Genome-Wide Identification and Analysis of the APETALA2 (AP2) Transcription Factor in Dendrobium officinale. Int. J. Mol. Sci. 2021, 22, 5221. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, W.; Zhang, J.; Zheng, S.; Wang, Z.; Xu, C.; Zhang, Q.; Wu, J.; Lou, H. Metabolite profiling and transcriptome analyses reveal novel regulatory mechanisms of melatonin biosynthesis in hickory. Hortic. Res. 2021, 8, 196. [Google Scholar] [CrossRef] [PubMed]
  54. Yokosho, K.; Yamaji, N.; Ma, J.F. An Al-inducible MATE gene is involved in external detoxification of Al in rice. Plant J. 2011, 68, 1061–1069. [Google Scholar] [CrossRef] [PubMed]
  55. Wang, C.; Deng, P.; Chen, L.; Wang, X.; Ma, H.; Hu, W.; Yao, N.; Feng, Y.; Chai, R.; Yang, G.; et al. A Wheat WRKY Transcription Factor TaWRKY10 Confers Tolerance to Multiple Abiotic Stresses in Transgenic Tobacco. PLoS ONE 2013, 8, e65120. [Google Scholar] [CrossRef] [PubMed]
  56. Lou, H.Q.; Fan, W.; Jin, J.F.; Xu, J.M.; Chen, W.W.; Yang, J.L.; Zheng, S.J. A NAC-type transcription factor confers aluminium resistance by regulating cell wall-associated receptor kinase 1 and cell wall pectin. Plant Cell Environ. 2020, 43, 463–478. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, Y.; Xu, J.; Guo, S.; Yuan, X.; Zhao, S.; Tian, H.; Dai, S.; Kong, X.; Ding, Z. AtHB7/12 Regulate Root Growth in Response to Aluminum Stress. Int. J. Mol. Sci. 2020, 21, 4080. [Google Scholar] [CrossRef] [PubMed]
  58. Sun, X. Genome assembly and annotation of Torreya grandis. Available online: https://doi.org/10.6084/m9.figshare.21089869.v1 (accessed on 15 September 2022).
  59. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef]
  60. Armenteros, J.J.A.; Salvatore, M.; Emanuelsson, O.; Winther, O.; von Heijne, G.; Elofsson, A.; Nielsen, H. Detecting sequence signals in targeting peptides using deep learning. Life Sci. Alliance 2019, 2, e201900429. [Google Scholar] [CrossRef]
  61. Thompson, J.D.; Gibson, T.J.; Higgins, D.G. Multiple Sequence Alignment Using ClustalW and ClustalX. Curr. Protoc. Bioinf. 2003, 2.3.1–2.3.22. [Google Scholar]
  62. Chernomor, O.; von Haeseler, A.; Minh, B.Q. Terrace Aware Data Structure for Phylogenomic Inference from Supermatrices. Syst. Biol. 2016, 65, 997–1008. [Google Scholar] [CrossRef]
  63. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  64. Bailey, T.L.; Johnson, J.; Grant, C.E.; Noble, W.S. The MEME Suite. Nucleic Acids Res. 2015, 43, W39–W49. [Google Scholar] [CrossRef] [PubMed]
  65. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef] [PubMed]
  66. Waheed, H.; Javaid, M.M.; Shahid, A.; Ali, H.H.; Nargis, J.; Mehmood, A. Impact of foliar-applied Hoagland’s nutrient solution on growth and yield of mash bean (Vigna mungo L.) under different growth stages. J. Plant Nutr. 2019, 42, 1133–1141. [Google Scholar] [CrossRef]
  67. Luo, Z.; Wang, L.; Wang, Y.; Zhang, W.; Guo, Y.; Shen, Y.; Jiang, L.; Wu, Q.; Zhang, C.; Cai, Y.; et al. Identifying and characterizing SCRaMbLEd synthetic yeast using ReSCuES. Nat. Commun. 2018, 9, 1930. [Google Scholar] [CrossRef] [PubMed]
  68. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  69. Guo, A.Y.; Chen, X.; Gao, G.; Zhang, H.; Zhu, Q.H.; Liu, X.C.; Zhong, Y.F.; Gu, X.; He, K.; Luo, J. PlantTFDB: A comprehensive plant transcription factor database. Nucleic Acids Res. 2008, 36, D966–D969. [Google Scholar] [CrossRef]
  70. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
Ijms 25 03859 g001
Figure 1. The distribution of the TgMATE gene family on the chromosomes and the occurrences of tandem duplication events (presented in red brackets).
Figure 1. The distribution of the TgMATE gene family on the chromosomes and the occurrences of tandem duplication events (presented in red brackets).
Ijms 25 03859 g001
Ijms 25 03859 g002
Figure 2. The phylogenetic tree of the expansion MATE genes in Oryza sativa, Nicotiana tabacum, Arabidopsis thaliana, Ricinus communis, Brassica oleracea, Brassica rapa subsp. Oleifera, Eucalyptus camaldulensis, Glycine max, Hordeum vulgare subsp. Vulgare, Medicago truncatula, Populus trichocarpa, Secale cereale, Sorghum bicolor, Torreya grandis, Triticum aestivum, Vigna umbellate, Zea mays, and Vitis vinifera. The Roman numerals I through V represent the five evolutionary groups of the MATE family.
Figure 2. The phylogenetic tree of the expansion MATE genes in Oryza sativa, Nicotiana tabacum, Arabidopsis thaliana, Ricinus communis, Brassica oleracea, Brassica rapa subsp. Oleifera, Eucalyptus camaldulensis, Glycine max, Hordeum vulgare subsp. Vulgare, Medicago truncatula, Populus trichocarpa, Secale cereale, Sorghum bicolor, Torreya grandis, Triticum aestivum, Vigna umbellate, Zea mays, and Vitis vinifera. The Roman numerals I through V represent the five evolutionary groups of the MATE family.
Ijms 25 03859 g002
Ijms 25 03859 g003
Figure 3. Protein motifs and gene structure of TgMATEs. (A) The conserved motif information of the TgMATE proteins. The Maximin Likelihood approach was used to build the phylogenetic tree of the TgMATE proteins. (B) Gene structure of the TgMATE genes. Yellow boxes are used to symbolize the untranslated regions (UTRs), and the green boxes represent the CDS areas. The introns are represented by the scale bar. (C) Conserved motifs of TgMATE proteins. The different font colors represent different amino acids.
Figure 3. Protein motifs and gene structure of TgMATEs. (A) The conserved motif information of the TgMATE proteins. The Maximin Likelihood approach was used to build the phylogenetic tree of the TgMATE proteins. (B) Gene structure of the TgMATE genes. Yellow boxes are used to symbolize the untranslated regions (UTRs), and the green boxes represent the CDS areas. The introns are represented by the scale bar. (C) Conserved motifs of TgMATE proteins. The different font colors represent different amino acids.
Ijms 25 03859 g003
Ijms 25 03859 g004
Figure 4. The different categories of cis-acting elements in TgMATE gene promoters. The Roman numerals I through V represent the five evolutionary groups of the TgMATE family.
Figure 4. The different categories of cis-acting elements in TgMATE gene promoters. The Roman numerals I through V represent the five evolutionary groups of the TgMATE family.
Ijms 25 03859 g004
Ijms 25 03859 g005
Figure 5. The expression patterns of the TgMATE family genes in response to aluminum (Al) stress, drought (Dr) stress, high-temperature (HT) stress, and low-temperature (LT) stress in Torreya grandis. (A) The numbers of up-regulated and down-regulated TgMATE genes under Al, Dr, HT, and LT stress. (B) Venn diagrams show the overlap of differentially expressed TgMATEs in response to Al (green), Dr (blue), HT (red), and LT (yellow). (C–F) Heat maps show the expression patterns of the stress-responsive TgMATEs under Al (C), Dr (D), HT (E), and LT stress (F). The FPKM values are log-normalized and clustered in rows.
Figure 5. The expression patterns of the TgMATE family genes in response to aluminum (Al) stress, drought (Dr) stress, high-temperature (HT) stress, and low-temperature (LT) stress in Torreya grandis. (A) The numbers of up-regulated and down-regulated TgMATE genes under Al, Dr, HT, and LT stress. (B) Venn diagrams show the overlap of differentially expressed TgMATEs in response to Al (green), Dr (blue), HT (red), and LT (yellow). (C–F) Heat maps show the expression patterns of the stress-responsive TgMATEs under Al (C), Dr (D), HT (E), and LT stress (F). The FPKM values are log-normalized and clustered in rows.
Ijms 25 03859 g005
Ijms 25 03859 g006
Figure 6. The proportion of differentially expressed genes (DEGs) and unresponsive genes in different polygenetic clades of TgMATEs under Al stress (Al), drought stress (Dr), high-temperature stress (HT), and low-temperature stress (LT).
Figure 6. The proportion of differentially expressed genes (DEGs) and unresponsive genes in different polygenetic clades of TgMATEs under Al stress (Al), drought stress (Dr), high-temperature stress (HT), and low-temperature stress (LT).
Ijms 25 03859 g006
Ijms 25 03859 g007
Figure 7. qRT-PCR analysis of the 7 TgMATE genes from clade IV under Al stress. The gene expression is normalized relative to the Torreya grandis reference gene, TgActin. The data present the average of three biological replicates, and ANOVA was employed to test the hypothesis. Different letters represent statistically significant p-values that are less than 0.05.
Figure 7. qRT-PCR analysis of the 7 TgMATE genes from clade IV under Al stress. The gene expression is normalized relative to the Torreya grandis reference gene, TgActin. The data present the average of three biological replicates, and ANOVA was employed to test the hypothesis. Different letters represent statistically significant p-values that are less than 0.05.
Ijms 25 03859 g007
Ijms 25 03859 g008
Figure 8. Relationships between TgMATEs and possible transcription factors predicted by WGCNA methods. (A) Determination of soft threshold power. (B) Eigengene adjacency heatmap of WGCNA. (C) Co-expression network represents the relationships between TgMATEs and possible transcription factors predicted by WGCNA methods. The blue ovals represent the TgMATE genes, and the yellow ovals represent transcription factors. Different background colors represent different modules.
Figure 8. Relationships between TgMATEs and possible transcription factors predicted by WGCNA methods. (A) Determination of soft threshold power. (B) Eigengene adjacency heatmap of WGCNA. (C) Co-expression network represents the relationships between TgMATEs and possible transcription factors predicted by WGCNA methods. The blue ovals represent the TgMATE genes, and the yellow ovals represent transcription factors. Different background colors represent different modules.
Ijms 25 03859 g008
Ijms 25 03859 g009
Figure 9. Dual-luciferase experiment between TgMATE2 and its candidate transcription factors. (A) Structure schematic of reporter and effector. (B) Transcriptional activation ability of 4 transcription factors for TgMATE2 in tobacco leaves. Each data presents the average of three biological replicates. A quadruple asterisk denotes the significant difference (p < 0.0001) assessed by ANOVA. (C) Fluorescence is released as a result of the reaction between the substrate and firefly luciferase. The scale bar corresponds to 1 cm.
Figure 9. Dual-luciferase experiment between TgMATE2 and its candidate transcription factors. (A) Structure schematic of reporter and effector. (B) Transcriptional activation ability of 4 transcription factors for TgMATE2 in tobacco leaves. Each data presents the average of three biological replicates. A quadruple asterisk denotes the significant difference (p < 0.0001) assessed by ANOVA. (C) Fluorescence is released as a result of the reaction between the substrate and firefly luciferase. The scale bar corresponds to 1 cm.
Ijms 25 03859 g009
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

Shen, H.; Hou, Y.; Wang, X.; Li, Y.; Wu, J.; Lou, H. Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis. Int. J. Mol. Sci. 2024, 25, 3859. https://doi.org/10.3390/ijms25073859

AMA Style

Shen H, Hou Y, Wang X, Li Y, Wu J, Lou H. Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis. International Journal of Molecular Sciences. 2024; 25(7):3859. https://doi.org/10.3390/ijms25073859

Chicago/Turabian Style

Shen, Hang, Ying Hou, Xiaorong Wang, Yaru Li, Jiasheng Wu, and Heqiang Lou. 2024. "Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis" International Journal of Molecular Sciences 25, no. 7: 3859. https://doi.org/10.3390/ijms25073859

APA Style

Shen, H., Hou, Y., Wang, X., Li, Y., Wu, J., & Lou, H. (2024). Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis. International Journal of Molecular Sciences, 25(7), 3859. https://doi.org/10.3390/ijms25073859

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