Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Post-Transcriptional Induction of the Antiviral Host Factor GILT/IFI30 by Interferon Gamma
Previous Article in Journal
Novel Clinical, Immunological, and Metabolic Features Associated with Persistent Post-Acute COVID-19 Syndrome
 
 
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 17
10.3390/ijms25179662
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

Identification and Analysis of Aluminum-Activated Malate Transporter Gene Family Reveals Functional Diversification in Orchidaceae and the Expression Patterns of Dendrobium catenatum Aluminum-Activated Malate Transporters

by
Fu-Cheng Peng
,
Meng Yuan
,
Lin Zhou
,
Bao-Qiang Zheng
and
Yan Wang
*
Key Laboratory of Tree Breeding and Cultivation of National Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(17), 9662; https://doi.org/10.3390/ijms25179662
Submission received: 21 August 2024 / Revised: 30 August 2024 / Accepted: 3 September 2024 / Published: 6 September 2024
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Review Reports Versions Notes

Abstract

:
Aluminum-activated malate transporter (ALMT) genes play an important role in aluminum ion (Al3+) tolerance, fruit acidity, and stomatal movement. Although decades of research have been carried out in many plants, there is little knowledge about the roles of ALMT in Orchidaceae. In this study, 34 ALMT genes were identified in the genomes of four orchid species. Specifically, ten ALMT genes were found in Dendrobium chrysotoxum and D. catenatum, and seven were found in Apostasia shenzhenica and Phalaenopsis equestris. These ALMT genes were further categorized into four clades (clades 1–4) based on phylogenetic relationships. Sequence alignment and conserved motif analysis revealed that most orchid ALMT proteins contain conserved regions (TM1, GABA binding motif, and WEP motif). We also discovered a unique motif (19) belonging to clade 1, which can serve as a specifically identified characteristic. Comparison with the gene structure of AtALMT genes (Arabidopsis thaliana) showed that the gene structure of ALMT was conserved across species, but the introns were longer in orchids. The promoters of orchid ALMT genes contain many light-responsive and hormone-responsive elements, suggesting that their expression may be regulated by light and phytohormones. Chromosomal localization and collinear analysis of D. chrysotoxum indicated that tandem duplication (TD) is the main reason for the difference in the number of ALMT genes in these orchids. D. catenatum was chosen for the RT-qPCR experiment, and the results showed that the DcaALMT gene expression pattern varied in different tissues. The expression of DcaALMT1-9 was significantly changed after ABA treatment. Combining the circadian CO2 uptake rate, titratable total acid, and RT-qPCR data analysis, most DcaALMT genes were highly expressed at night and around dawn. The result revealed that DcaALMT genes might be involved in photosynthate accumulation. The above study provides more comprehensive information for the ALMT gene family in Orchidaceae and a basis for subsequent functional analysis.

1. Introduction

Aluminum-activated malate transporters (ALMT) are a gene family encoding anion channel protein [1,2,3]. ALMT genes share the common feature of containing a conserved structural domain, PF11744, which is comprised of several transmembrane structural domains (TMDs) [4]. ALMT genes are widely present in terrestrial plants and encode channel proteins that can transport organic anions (e.g., malate) and inorganic anions (e.g., Cl−, SO4−, NO3−) in cells [5,6]. A large number of studies have demonstrated that the functions of ALMT gene family members are involved in plant Al3+ tolerance [7,8,9,10], stomatal movement [11,12,13,14,15], mineral nutrition [16,17], and fruit acidity [18,19]. In the C4 and CAM pathways, ALMT is mainly involved in the carboxylation phase of malate, which transports malate into the vacuole [20,21]. Furthermore, transmembrane transport of malate across the vacuole membrane has osmoregulatory effects, especially in guard cells, which can drive stomatal movement [22].
The ALMT gene family has been identified in several species, and the functions of its single members have been extensively studied, including model plants such as A. thaliana and rice [10,13,23] and non-model plants such as apple and soybean [24,25]. Fourteen AtALMT genes were identified in A. thaliana, which can be divided into four clades. Nine OsALMT genes were identified in rice, and clade 5 appeared, containing two members: OsALMT1 and OsALMT2. Subsequent studies demonstrated that its function is related to cellular osmoregulation and the maintenance of electroneutrality [23]. A total of 34 GmALMT genes were identified in the soybean genome. In a P-deficient environment, GmALMT5 improved the utilization of insoluble phosphorus sources and then improved phosphorus efficiency in soybeans [25]. Twenty-five ALMT genes were identified in the domesticated apple (Malus × domestica Borkh.), and whole-genome duplication (WGD) played an important role in the expansion of the MdALMT gene family [24]. In addition, eight candidate ALMT genes, including three ALMT9 homologous genes (AcoALMT9-1-3) and five ALMT1 homologous genes (AcoALMT1-1-5), were identified in the CAM plant pineapple (Ananas comosus). AcoALMT1 had higher transcript levels under drought treatment [26]. The functions of single ALMT genes have also been studied. The expression levels of ALMT genes related to Al3+ tolerance (such as TaALMT1 and AtALMT1) were affected by conditions of Al3+ and low pH [27]. The overexpression of StALMT6 and StALMT10(Solanum tuberosum) in A. thaliana enhanced tolerance to Al3+ toxicity [28]. SlALMT4 (Solanum lycopersicum) and SlALMT5 are mainly expressed during fruit development. Their overexpression significantly increases the malic and citric acid content of seeds in the fruit [29]. AtALMT12 was expressed mainly in leaf guard cells. The almt12 mutant plants have impaired stomatal closure induced by darkness, CO2, and ABA, resulting in accelerated plant wilting [14]. AtALMT6 is also preferentially expressed in leaf guard cells and is involved in stomatal opening [15].
The CAM pathway opens their stomata at night to absorb CO2 and close them during the day. This behavior reduces the transpiration loss of water. CAM plants, which have higher water use efficiency (WUE), are better adapted to arid environments compared to the C3 and C4 plants [30,31]. ALMTs, as the main bearer of malate transport and stomatal movement, have attracted the interest of many researchers. The diel (diel means diurnal or day/night) expression of the putative ALMT6 (Kaladp0062s0038) gene in the typical CAM plant Kalanchoë fedtschenkoi has a distinct circadian rhythm [32]. It is mainly expressed at night and may be related to vacuole malate transport. Two genes homologous to AtALMT9 (Aco003023.1 and Aco010725.1) are members of the ALMT gene family in pineapple [26]. They exhibit the highest expression in leaf photosynthetic tissues and may be prime candidates for the malate influx pathway. Aco003023.1 is predicted to be a target of miR172d-5p, which is involved in the regulation of photosynthesis (predicted by psRNATarget and TAPIR). The functions of these genes have not been validated in CAM species; therefore, molecular identification and functional characterization of ALMT genes in CAM plants are crucial goals.
CAM plants occupy a large proportion of Bromeliaceae, Crassulaceae, Orchidaceae, and Asclepiadaceae [33]. Orchidaceae is the second-largest family of angiosperms, with at least 112 genera showing strong CAM behavior [34]. Orchids are not only prized for their economic and ornamental value but also play a significant role in evolutionary and ecological studies [35]. Whole-genome sequencing of more orchids provides a platform for gene family analyses [36,37,38,39]. Currently, most of the studies on the orchid ALMT genes have focused on the transcriptome level, lacking analyses such as structural and evolutionary analyses at the genome level [26,32]. There is a lack of studies that delve into the expression patterns of orchid ALMT genes. In this study, we mainly identified the members of the ALMT gene family in four orchids, namely A. shenzhenica, P. equestris, D. catenatum, and D. chrysotoxum. Next, we analyzed the characteristics, structure, phylogenetic relationship, and cis-acting elements of these members. D. catenatum was selected as a representative for expression pattern analysis under the tissue, circadian, and ABA treatments. Our results provide a reference for further clarification of the function of the orchid ALMT gene family.

2. Results

2.1. Identification and Characterization of the ALMT Gene Family in Orchidaceae

Thirty-four ALMT genes were identified from four orchid genomes (ten in D. catenatum and D. chrysotoxum, seven in A. shenzhenica and P. equestris) (Table 1). The protein sequences of these ALMT ranged from 253 to 887 aa. The molecular weights ranged from 27,270.69 to 96,227.92 Da, and the theoretical isoelectric points ranged from 5.7 to 9.39. Seventy percent of the ALMT proteins possessed high isoelectric points (pI > 7.0). The majority of them had instability indices below 40 (II), which indicates good stability. The average aliphatic index (AI) of the 34 ALMT proteins was 99.87, indicating high thermal stability. In addition, the average hydrophilicity index (GRAVY) of all orchid ALMT proteins ranged between 0.5 and −0.5, indicating that they are amphoteric amino acids. Transmembrane structural domains (TMDs) were also predicted, and all ALMT gene family members were detected at the N-terminus with numbers ranging from two to seven. Subcellular localization predictions showed that the majority (79.4%) were localized to the plasma membrane. The detailed orchid ALMT protein sequences are displayed in Supplementary Table S1.
The gene ontology analysis showed significant enrichment of orchid ALMT genes in various biological processes (Figure S1). The molecular function was focused on “transmembrane transporter activity”; the most abundant cellular component was the “plant-type vacuole membrane”, and the most abundant bioprocess was the “C4-dicarboxylate transport”. These results suggest that orchid ALMT genes are closely related to C4 acid transport and may be involved in dicarboxylate transport in the photosynthetic pathway.

2.2. Chromosomal Localization of DchALMT Genes

To reveal the distribution of ALMT genes on orchid chromosomes, we mapped the location of ALMT using genome annotation information (Figure 1). Since only the D. chrysotoxum genome was assembled to the chromosome level, the positions of the A. shenzhenica, P. equestris, and D. catenatum genes on the scaffold were used only for reference (Figure S2). We found that the DchALMT gene was unevenly distributed on seven chromosomes. Among them, chromosomes 08 and 16 have two ALMT genes at the same locus. We further performed collinearity analysis of DchALMT genes using the TBtools MCScanX (version 2.119) tool. The results showed that there were tandem duplications (TD) in two pairs of genes: DchALMT2/7 and DchALMT8/10. Their Ka/Ks values were 0.249092716 and 0.49015544 (Table S2), suggesting that the ALMT genes were subjected to purifying selection during the evolutionary process. Interestingly, even though the D. catenatum genome is unassembled to the chromosome, we still found genes located at the same locus: DcaALMT9/4 and DcaALMT2/8/5 at the scaffold position. The same situation was not found in A. shenzhenica and P. equestris.

2.3. Multiple Sequence Alignment and Phylogenetic Analysis of ALMT Genes

Previous studies have shown that ALMT-conserved domains contain transmembrane domains (TMD), Gamma-aminobutyric acid (GABA) binding motifs, and WEP motifs [40,41,42]. Multiple sequence alignment of orchid ALMT protein revealed that almost all ALMTs contain similar conserved sequences (Figure 2). The TMD1 sequence consists of the highly conserved TVVVVVFE sequence, which is preceded by other conserved residues. The transmembrane domain, along with the cytoplasmic helical domain (CHD), constitutes the structural foundation of the ion channel [40]. In previous studies, amino acid residues have been associated with channel activity [43,44]. As a signaling substance, GABA is involved in the regulation of ALMT, while ALMT can act as a channel for GABA binding and transport [45]. Mutations in residues F213 and F215 can weaken the interactions between the two [46]. We found that some of the F residues are replaced with L in the orchid ALMT GABA motif, and the effect of this site on recognition and transport needs to be further investigated. Glutamate E is important for ALMT channel activity and was found in TaALMT1 (E284), AtALMT1 (E256), and AtALMT12 (E276) [43,44]. We found that all ALMTs contain a conserved WEP motif except PeALMT7. Meanwhile, PeALMT7 lacks multiple conserved residues such as PW, Y, and R. We speculate that PeALMT7 may have lost channel activity.
To investigate the phylogenetic relationships of the ALMT family, phylogenetic trees were constructed using full-length ALMT protein sequences from Orchidaceae (34 genes), arabidopsis (13 genes), rice (9 genes), and other species (14 genes) (Figure 3). A total of 70 ALMT proteins were classified into four clades, consistent with previous classification in AtALMT [3]. Clade 1 had the highest number, and clade 4 had the lowest number of members. There were 10, 7, 8, and 9 orchid ALMT members in clades 1–4, respectively. The functions of most genes in clades 1–3 have been characterized, except clade 4.

2.4. Gene Structure and Motif Analysis of ALMT Genes

The conserved domains of ALMT proteins were searched using NCBI-CD search tools. All orchid ALMT genes contain PF11744 structural domains, including some truncated structural domains, such as DcaALMT5, AsALMT4, PeAlMT7, DcaALMT10, and DchALMT4 (Figure 4C). Previous studies have found that truncation of domains leads to altered function [2]. It was also found that DchALMT7 has a repetitive structure with two almost the same structural domains, like the ALMT identified in grape [4].
To understand the structure of orchid ALMT genes, 34 ALMT proteins were analyzed using the MEME online tool, setting the motif search limit to 20. Most of the ALMT proteins had the same conserved motifs in each branch, and the number of ALMT motifs ranged from 6 to 29 (Figure 4B, Table S3). The results showed that the distribution of most ALMT motifs is relatively conserved, with motifs 6/15/1/10/9/3/5/17/4/2/8/7 being typical. The TM region corresponds to motif 1, the GABA binding motif corresponds to motif 5, and the WEP motif corresponds to motif 4. In addition, motif19 was found to be present only in clade 1, which may be a unique feature of the orchid ALMT clade 1 compared with Arabidopsis (Figure S3). The motif arrangement in clades 2 and 3 was more different, which may predict the functional differentiation in the clade.
Meanwhile, we analyzed the intron-exon structure of the orchid ALMT genes (Figure 4D). The results showed that all ALMTs contained 4–12 exons and 3–11 introns. We found that compared with Arabidopsis ALMT, orchid ALMT is a typical long intron type (Figure S3). This has also been found in the P. equestris and D. chrysotoxum genomes. This has been regarded as a general feature of the orchid genome [37,39]. The longer intron sequences have a higher probability of homologous recombination [47]. This could mean that orchids can produce rich phenotypes in the face of natural selection during evolution.

2.5. Cis-Elements in the Promoter Regions of ALMT Genes

To explore the regulatory roles of orchid ALMTs, we searched the 2000 bp region upstream of the ALMT gene to identify potential cis-acting elements (Figure 5 and Supplementary Tables S3 and S4). We identified a total of 667 cis-acting elements, including 21 types and 14 response functions, in addition to the common TAAT-box and CAAT box (associated with transcription initiation and promotion). The cis-acting elements included responses to phytohormones such as gibberellins, abscisic acid, salicylic acid, auxin, ethylene, and methyl jasmonate; abiotic stresses such as drought, low temperature, defense, and stress responses; and plant growth and development processes such as light response and anaerobic sensing. We found binding sites for WRKY, ERF, and MYB transcription factors, indicating that orchid ALMT genes are regulated by multiple transcription factors. Each ALMT gene contains multiple elements, of which the most abundant is the light-responsive element (Box 4), followed by ethylene-responsive element (ERE), MYB-binding site, and abscisic acid-responsive element (ABRE). The information on cis-elements is displayed in Tables S4 and S5.

2.6. Tissue-Specific Expression of DcaALMT

To understand the function of orchid ALMT genes, DcaALMT genes were selected to determine the expression pattern in four different tissues (roots, stems, leaves, and flowers) by RT-qPCR. The results showed that all DcaALMT were expressed in the four tissues except DcaALMT10 (Figure 6). We found that DcaALMT1/3/5/6 were highly expressed in leaves, and DcaALMT2/7/8 were expressed highly in flowers. DcaALMT4/9/10 were highly expressed in roots. Previous studies have shown that Zea mays ZmALMT1 is specifically expressed in roots and plays a role in mineral nutrient acquisition and transport [16]. In this study, DcaALMT4/9/10 were highly expressed in roots, suggesting that they may play a role in root nutrient uptake.

2.7. Diel Expression of DcaALMT

To further investigate the expression pattern of DcaALMTs, the diel expression patterns of ALMT genes in D. catenatum leaves were determined and analyzed in combination with the rate of CO2 uptake and the content of titratable total acid in the leaves. The results showed that the diurnal changes in CO2 uptake rate and titratable total acid in D. catenatum were consistent with the characteristics of the CAM pathway (Figure 7A,B). We found that the highest expression of DcaALMT1/2/3/4/8/9 was at night and dawn, which may be related to CO2 uptake and total acid accumulation (Figure 7C). The peak expression of DcaALMT1/2/4/8 was at night, which was already reduced to a certain level before dawn, and then the expression was lower throughout the photoperiod. The expression of DcaALMT3/9 was elevated immediately at dawn and reached a peak after 1–2 h (6:00 or 7:00). These results suggest that the expression level of DcaALMTs is affected by circadian rhythms and is likely to be involved in the organic acid accumulation of photosynthesis pathway.

2.8. Expression of DcaALMT under ABA Treatment

To investigate whether DcaALMT is involved in the response to ABA, we used RT-qPCR to study the expression of DcaALMT in D. catenatum leaves under ABA treatment. The results showed that all DcaALMT were sensitive to ABA treatment (Figure 8). Among them, the expression of DcaALMT3/4/5/8 was up-regulated after ABA treatment, and they reached a peak at 1.5 h or 3 h after treatment and then down-regulated. DcaALMT5 was the most significantly up-regulated, which was 3.5-fold higher than 0 h. In contrast, DcaALMT1/2/6/7/9 showed a significant down-regulation after 6 h of treatment, with the most significant down-regulation (0.04) in the expression level of DcaALMT9.

3. Discussion

ALMT was first identified in wheat roots and named after its ability to enhance plant Al3+ tolerance [10]. Subsequent studies of ALMT over the last decade have revealed that its function is not limited to Al3+ resistance but is also associated with stomatal movement, mineral nutrition, fruit acidity, seed development, and GABA signaling [7,11,14,19,23,45,46]. Large-scale plant genome sequencing has laid the foundation for analyzing this family. The ALMT gene family has been characterized in the Brassicaceae [42], Poaceae [48,49], Rosaceae [50], and Fabaceae [51,52]. In Orchidaceae, the second largest family among angiosperms, we know little about the ALMT family. In this study, 34 ALMT genes were identified from four orchid species. The results revealed that each orchid genome contained seven to ten ALMT genes, which were fewer than those of A. thaliana (14), Brassica napus (39), and apple (25). This may be due to the WGD events that occurred in the evolutionary history of Orchidaceae, most of which were accompanied by gene loss after the WGD event [53]. Chromosome distribution showed that 10 DchALMTs were distributed on seven chromosomes. DchALMT10/8 and DchALMT2/7 were two pairs of tandem duplicated genes. Tandem duplications (TDs) are important for the adaptation of plants to rapidly changing environments [54]. KaKs calculations indicated that these two pairs of genes were subjected to purifying selection during evolution. This is the main reason for the difference in the number of ALMT genes between orchids and A. thaliana. The same situation occurs in the soybean MYB family. MYB genes in the soybean show stronger purifying selection, with the appearance of pseudogenes and functionally redundant genes [55]. Purification selection can help plants remove deleterious mutations and retain important genes.
Phylogenetic analyses showed that ALMT genes from several species, including Orchidaceae, can be divided into four clades, which is consistent with the ALMT family in A. thaliana and Hevea brasiliensis [56], confirming the reliability of grouping. Phylogenetic, gene structure, and conserved motif analyses of orchid ALMT genes revealed that most ALMT genes contain conserved N-segment structural domains, which contain 5–7 TMDs (Table 1). The transmembrane region of the N-terminal structural domains is mainly responsible for the transporter activity of the protein [57]. We also found members containing truncated ALMT-conserved domains, such as AsALMT4, DcaALMT5, AsALMT2, and PeALMT7. AtALMT11 represents the same structure that may signify a loss of function within these genes. It is also possible that this is due to incomplete genome annotation, which frequently occurs in the orchid genome [36,39]. In further studies of gene structure, we found that the presence of motif12 only in clade 1 provides more referability to the branching of the orchid ALMT phylogeny. The gene structures and motifs of clades 1 and 2 were more conserved within groups, while clades 2 and 3 showed large differences within groups. We speculate that functional divergence may have occurred in clades 2 and 3. In clades 2 and 3, a large number of members (60%) detected novel motifs, mostly concentrated at the C-terminus, with a small concentration at the N-terminus. The C-terminus tends to be associated with independent activity rather than Al3+ resistance [58]. Combined with multispecies phylogenetic analyses, AtALMT4/6/9 in clade 2 and AtALMT12 in clade 3 were associated with stomatal movement [12,14]. The function of these two clades may be related to the unique stomatal behavior of the CAM pathway. Compared to Arabidopsis, the number of introns and exons of the orchid ALMT genes is more similar to that of the AtALMTs, but the intron length is significantly longer than that of AtALMTs. This has become an important feature of orchid species [37]. Longer introns are conducive to the occurrence of homologous recombination, which forms the different combinations of protein domains and enhances phenotypic polymorphism [47].
Gene expression is regulated by upstream transcription factors, making the prediction of cis-acting elements of gene promoters essential. A variety of cis-acting elements were identified in the promoter region of orchid ALMT genes, among which light-responsive elements were the most abundant, followed by ethylene, MYB-binding sites, and ABA-responsive elements. These results suggest that light and hormones may affect the expression of the orchid ALMT gene. The effects of ethylene on ALMT were mainly in Al3+ tolerance and fruit acidity regulation [59,60]. In wheat, ethylene can negatively regulate TaALMT1, reducing malate efflux from the root and making wheat more sensitive to Al3+ toxicity stress [59]. In apples, MdESE3 can bind to the ethylene-response element (ERE) located in the MdMa11 promoter and thus activate its expression, promoting malate accumulation [60]. Under drought stress, both AtALMT4 and AtALMT12 are regulated by ABA, which in turn closes stomata to reduce water dissipation [61]. MYB is an important transcription factor regulating malic acid accumulation. The R2R3-type MdMYB73 protein directly binds to the promoter of MdALMT9 to enhance its expression, resulting in the accumulation of malic acid in the vacuole [62]. The above results suggest that Al3+ tolerance, stomatal movement, and malate accumulation in the vacuole may be the main direction of the functional study of orchid ALMT genes, while the link between MYB and orchid ALMT needs to be further explored.
Previous studies have shown that ALMT has different expression patterns in different plant tissues. For example, MsALMT1 (Medicago sativa) is mainly expressed in roots, which is involved in the response to heavy metal ion stress. PtaALMTs (Pinus tabuliformis) are specifically expressed in flowers and roots, which may be related to reproductive organs and nutrient uptake. AtALMT6 and AtALMT9 are highly expressed in the mesophyll tissue of leaves, and their functions are related to stomatal opening and closing. In our study, the 10 DcaALMT genes had different tissue expression patterns. Four of them (DcaALMT1/3/5/6) were highly expressed in leaves, three (DcaALMT2/7/8) in flowers, and three (DcaALMT4/9/10) in roots. In addition, DcaALMT10 was specifically expressed in roots and flowers. These studies suggest that DcaALMT genes may perform different functions in different tissues. The expression of the BnALMT genes was up-regulated in roots under phosphorus-deficient conditions. Therefore, we need to further analyze the expression pattern of DcaALMT genes to determine the tissue expression pattern under different stress treatments.
An important feature of the CAM pathway is its circadian rhythms [63,64], and the expression of DcaALMT genes may also be regulated by the biological clock. We further analyzed the circadian expression pattern of DcaALMT genes. The results showed that the expression of DcaALMT1/2/3/4/8/9 was affected by circadian rhythms and corresponded to diurnal fluctuations in the rate of CO2 uptake and titratable total acid. This is consistent with previous reports in the CAM plant K. fedtschenkoi. KfALMT6 and KfALMT1 are highly expressed at night and dawn [32]. The role of AtALMT12, which is located in the same clade as DcaALMT3, in stomatal movement has been verified several times [5,14]. DcaALMT3 was expressed abundantly at dawn, which may lead to states of stomatal closure at a later stage. These states could correspond to the CO2 uptake rate curves and total acid curves.
ABA is an important hormone for stomatal closure in response to drought stress in plants [65,66]. Plants regulate the expression of AtALMT12 through ABA signaling via OST1 and ABI1, leading to stomatal closure [67,68]. Later studies show that xylem sulfate can also directly activate AtALMT12 expression to accomplish stomatal behavior under drought [69]. Recent studies have shown that AtALMT4 is also regulated by ABA and mediates stomatal closure by phosphorylation [64]. In this study, the DcaALMT genes showed different expression patterns under ABA treatment. DcaALMT3/4/5/8 reached a peak expression within 1.5–3 h after treatment, whereas DcaALMT1/2/6/7/9 showed decreased expression after treatment. These results suggest that they may play different functions in the ABA pathway. Stomatal movement is a complex process involving multiple signals and transporters. The Ca2+ signal in the guard cell can also regulate stomatal opening and closing [70]. Cytoplasmic GABA signaling can negatively regulate AtALMT9 to promote stomatal opening [71]. We need to study the expression and function of ALMT genes under various signaling molecules to understand how they regulate stomatal movement in CAM plants.
Integration of the CAM pathway into C3 plant crops has been desired as a potential strategy to improve water use efficiency in plants [21]. Previous studies have found that transforming genes from the carboxylation module of the CAM pathway in Mesembryanthemum crystallinum into Arabidopsis significantly increased plant rosette diameter, leaf area, and leaf fresh weight, and some of the genes also led to an increase in stomatal conductance and titratable total acid accumulation, whereas genes from the decarboxylation module led to a decrease in stomatal conductance and titratable total acid content [20]. Overexpression of the Populus euphratica xyloglucan endotransglucosylase/hydrolase gene (PeXTH) in Nicotiana tabacum decreased intracellular air space in palisade tissues while increasing leaf water content and cellular accumulation, resulting in improved salt tolerance [72]. Our findings that the orchid ALMT genes may play a role in stomatal movement and malate transport will help us to better study the stomatal movement module of the CAM pathway. The study of orchid ALMT genes in the CAM pathway can help us in the future to accelerate the progress of C3-to-CAM biosystems. This will help to improve the performance of non-CAM plants under stressful environments.

4. Materials and Methods

4.1. Identification and Bioinformatic Analysis of ALMT Genes in Orchid

Genome data of A. shenzhenica [36], P. equestris [39], D. catenatum [38], and D. chrysotoxum [37] were downloaded from the NCBI (https://www.ncbi.nlm.nih.gov/) (accessed on 11 August 2023). The AtALMT protein sequence was obtained from the TAIR (https://www.arabidopsis.org/) (accessed on 11 August 2023). The 14 AtALMTs protein sequences were used as query sequences, performing a blast on the genomic protein sequence. Next, an online blast was conducted using the NCBI Blastp tool (select the Swissport database). Duplicate sequences were removed after merging the above two results. NCBI CDD search tool and HMMTOP software (version 2.0) were used to identify the conserved domain and transmembrane domains. Combining the results, candidate members that do not contain ALMT domains and TMDs were removed. The basic information of ALMT protein was predicted using the ExPASY website ComputepI/Mwtool and ProtParam tools (https://web.expasy.org/compute_pi/) (accessed on 11 August 2023). Subcellular prediction is accomplished through Plant mPLo (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/) (accessed on 11 August 2023). Gene ontology analysis was performed on the eggNOG mapper (http://eggnog-mapper.embl.de/) (accessed on 11 August 2023), and the annotation results were visualized using TBtools [73].

4.2. Multiple Sequence Alignment and Evolutionary Tree Construction

Multiple comparisons of orchid ALMT protein sequences were performed using the Muscle program in MEGA7 software (version 7.0.14) [74]. The results were visualized and embellished using Jalview software (version 2.10.5) [75], and conserved amino acid sites and structural domains were annotated. The phylogenetic tree was constructed using MEGA7, the method was Maximum Likelihood, and the parameters were set: bootstrap value was set to 1000 times, and the reference model was partial deletion. The evolutionary tree was beautified using Evolview (http://www.evolgenius.info/evolview/) (accessed on 17 August 2023) [76].

4.3. Gene Structure and Motif of ALMT Genes

The gene structure (including exons, introns, and UTRs) information was obtained from genome annotation information and visualized using Bio-sequence Structure Illustrator in TBtools [73]. Motif analysis was performed through the online website MEME (http://meme-suite.org/tools/ (accessed on 28 August 2023)) with the following parameters: the number of motifs was set to 20, the site distribution was set to “Zero or one occurrence per sequence”, and the other parameters were kept as default.

4.4. Prediction of ALMT Genes Cis-Acting Elements

The 2000 bp upstream of orchid ALMTs was extracted by TBtools, and the online site PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (accessed on 28 August 2023) was used to identify cis-acting elements in the promoter region [77]. The number of elements and heatmap were visualized using TBtools [73].

4.5. Chromosome Localization and Collinearity Analysis

D. chrysotoxum with chromosome-level genome was selected as the subject of analysis. The distribution of ALMT genes on chromosomes was obtained from genome annotation files; tandem duplicate sequences were searched for using the MCScanX program [78]; visualization was performed using TBtools. The KaKs values of tandem duplicated genes were calculated using the TBtools Simple Kaks calculator [79].

4.6. Plant Materials and Growth Conditions

The seed source of D. catenatum was introduced from the South China Botanical Garden, Chinese Academy of Sciences, and cultivated in the greenhouse of the Chinese Academy of Forestry.
D. catenatum plants were uniformly concentrated in an artificial climate chamber for cultivation and management under photoperiodic dimensions of 12 h darkness and 12 h light, temperature of 26 ± 2 °C, and relative humidity of 60 ± 5%.
The leaves were treated with 100 µmol/L ABA and harvested from the same position after 0, 1.5, 3, and 6 h of the treatment for three biological replicates. The sampling time was from 18:00 to 18:00 the next day. The sampling interval was one hour for 24 h of consecutive sampling and three biological replicates for each sampling site. These samples were placed in liquid nitrogen and stored in a −80 °C refrigerator.

4.7. Physiological Index

The Li-6400XT was used to measure the diel CO2 exchange rate of D. catenatum leaves for 24 h. The instrument was equipped with a 5 L air buffer bottle with a constant gas flow rate (200 μms). Three plants were randomly measured at each sampling point. The readings of one sample point per plant were taken three times with an interval of 30 s.
We conducted the diel titratable total acid of D. catenatum leaves for 24 h. Three biological replicates were set up for sampling at each time point, with sample weights of approximately 0.3–0.6 g. The titratable acid content of the leaves was determined by acid-base titration.

4.8. Total RNA Extraction and RT-qPCR

RT-qPCR was used to analyze the expression of the DcaALMT gene with three biological replicates. The extraction of total RNA was used by the Polysaccharide Polyphenol Plant Total RNA Extraction Kit (TIANGEN; DP441). Easy Script One-Step gDNA Removal and cDNA Synthesis Super Mix (TRANS; AE311) were used for cDNA synthesis. The RT-qPCR was carried out using TB Green Premix Ex Taq (Takara; RR420A) with Roche LightCycler®480 Real-Time PCR system. Relative quantification of genes was calculated using 2-ΔΔCT. Primers for RT-qPCR were designed by Primer6 software (version 6.24) (Supplementary Table S6) and validated using primer-blast in NCBI.
ANNOVA multiple comparisons were used to analyze the qPCR expression levels between different samples.

5. Conclusions

In this study, 34 ALMT family members were identified in four orchid genomes and analyzed for their protein basic information and functional annotation. Phylogenetic analyses showed that orchid ALMT genes can be divided into four clades. Protein structural domains, conserved motifs, and gene structures analysis showed that clade 1 and clade 4 were conserved across species, while clade 2 and clade 3 were more variable across species, suggesting that functional differentiation may have occurred in clade 2 and clade 3. In D. chrysotoxum, tandem duplications were an important driver of ALMT evolution and a major cause of differences in gene number between species. Cis-acting element analysis proved that light and hormones are the main factors regulating ALMT genes. In addition, the expression patterns of different tissues of DcaALMT genes supported the functional differentiation of clade 2 and clade 3. The diel and ABA treatment expression patterns of DcaALMT genes predicted a role for DcaALMT in the CAM and ABA pathway. In this study, we comprehensively analyzed the basic information and expression patterns of ALMT genes in orchids. These results provide insight into the function of ALMT genes in the CAM pathway. The next task focuses on the functional validation of ALMT genes to improve richer genetic materials for CAM genetic engineering.

Supplementary Materials

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

Author Contributions

Y.W., L.Z. and B.-Q.Z. participated in the conceptualization and design of this study. Material preparation, data collection, and analysis were performed by F.-C.P. and M.Y. The first draft of the manuscript was written by F.-C.P. All authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Basic Research Fund of the Chinese Academy of Forestry, CAF (Grant No. CAFYBB2019ZB001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data are provided within this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liu, J.; Zhou, M. The ALMT Gene Family Performs Multiple Functions in Plants. Agronomy 2018, 8, 20. [Google Scholar] [CrossRef]
  2. Dabravolski, S.A.; Isayenkov, S.V. Recent Updates on ALMT Transporters’ Physiology, Regulation, and Molecular Evolution in Plants. Plants 2023, 12, 3167. [Google Scholar] [CrossRef]
  3. Medeiros, D.B.; Fernie, A.R.; Araújo, W.L. Discriminating the Function(s) of Guard Cell ALMT Channels. Trends Plant Sci. 2018, 23, 649–651. [Google Scholar] [CrossRef]
  4. Sharma, T.; Dreyer, I.; Kochian, L.; Piñeros, M.A. The ALMT Family of Organic Acid Transporters in Plants and Their Involvement in Detoxification and Nutrient Security. Front. Plant Sci. 2016, 7, 1488. [Google Scholar] [CrossRef] [PubMed]
  5. Dreyer, I.; Gomez-Porras, J.L.; Riaño-Pachón, D.M.; Hedrich, R.; Geiger, D. Molecular Evolution of Slow and Quick Anion Channels (SLACs and QUACs/ALMTs). Front. Plant Sci. 2012, 3, 263. [Google Scholar] [CrossRef]
  6. Eisenach, C.; De Angeli, A. Ion Transport at the Vacuole during Stomatal Movements. Plant Physiol. 2017, 174, 520–530. [Google Scholar] [CrossRef] [PubMed]
  7. Hoekenga, O.A.; Maron, L.G.; Piñeros, M.A.; Cançado, G.M.A.; Shaff, J.; Kobayashi, Y.; Ryan, P.R.; Dong, B.; Delhaize, E.; Sasaki, T.; et al. AtALMT1, Which Encodes a Malate Transporter, Is Identified as One of Several Genes Critical for Aluminum Tolerance in Arabidopsis. Proc. Natl. Acad. Sci. USA 2006, 103, 9738–9743. [Google Scholar] [CrossRef]
  8. Wang, J.; Yu, X.; Ding, Z.J.; Zhang, X.; Luo, Y.; Xu, X.; Xie, Y.; Li, X.; Yuan, T.; Zheng, S.J.; et al. Structural Basis of ALMT1-Mediated Aluminum Resistance in Arabidopsis. Cell Res. 2022, 32, 89–98. [Google Scholar] [CrossRef]
  9. 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]
  10. Sasaki, T.; Yamamoto, Y.; Ezaki, B.; Katsuhara, M.; Ahn, S.J.; Ryan, P.R.; Delhaize, E.; Matsumoto, H. A Wheat Gene Encoding an Aluminum-Activated Malate Transporter. Plant J. 2004, 37, 645–653. [Google Scholar] [CrossRef] [PubMed]
  11. De Angeli, A.; Zhang, J.; Meyer, S.; Martinoia, E. AtALMT9 Is a Malate-Activated Vacuolar Chloride Channel Required for Stomatal Opening in Arabidopsis. Nat. Commun. 2013, 4, 1804. [Google Scholar] [CrossRef]
  12. Gruber, B.D.; Ryan, P.R.; Richardson, A.E.; Tyerman, S.D.; Ramesh, S.; Hebb, D.M.; Howitt, S.M.; Delhaize, E. HvALMT1 from Barley Is Involved in the Transport of Organic Anions. J. Exp. Bot. 2010, 61, 1455–1467. [Google Scholar] [CrossRef]
  13. Kovermann, P.; Meyer, S.; Hörtensteiner, S.; Picco, C.; Scholz-Starke, J.; Ravera, S.; Lee, Y.; Martinoia, E. The Arabidopsis Vacuolar Malate Channel Is a Member of the ALMT Family. Plant J. 2007, 52, 1169–1180. [Google Scholar] [CrossRef]
  14. Meyer, S.; Mumm, P.; Imes, D.; Endler, A.; Weder, B.; Al-Rasheid, K.A.S.; Geiger, D.; Marten, I.; Martinoia, E.; Hedrich, R. AtALMT12 Represents an R-type Anion Channel Required for Stomatal Movement in Arabidopsis Guard Cells. Plant J. 2010, 63, 1054–1062. [Google Scholar] [CrossRef]
  15. Meyer, S.; Scholz-Starke, J.; Angeli, A.D.; Kovermann, P.; Burla, B.; Gambale, F.; Martinoia, E. Malate Transport by the Vacuolar AtALMT6 Channel in Guard Cells Is Subject to Multiple Regulation. Plant J. 2011, 67, 247–257. [Google Scholar] [CrossRef]
  16. Piñeros, M.A.; Cançado, G.M.A.; Maron, L.G.; Lyi, S.M.; Menossi, M.; Kochian, L.V. Not All ALMT1-type Transporters Mediate Aluminum-activated Organic Acid Responses: The Case of ZmALMT1—An Anion-selective Transporter. Plant J. 2007, 53, 352–367. [Google Scholar] [CrossRef]
  17. Ligaba, A.; Maron, L.; Shaff, J.; Kochian, L.; Piñeros, M. Maize ZmALMT2 Is a Root Anion Transporter That Mediates Constitutive Root Malate Efflux. Plant Cell Environ. 2012, 35, 1185–1200. [Google Scholar] [CrossRef]
  18. De Angeli, A.; Baetz, U.; Francisco, R.; Zhang, J.; Chaves, M.M.; Regalado, A. The Vacuolar Channel VvALMT9 Mediates Malate and Tartrate Accumulation in Berries of Vitis Vinifera. Planta 2013, 238, 283–291. [Google Scholar] [CrossRef]
  19. Bai, Y.; Dougherty, L.; Li, M.; Fazio, G.; Cheng, L.; Xu, K. A Natural Mutation-Led Truncation in One of the Two Aluminum-Activated Malate Transporter-like Genes at the Ma Locus Is Associated with Low Fruit Acidity in Apple. Mol. Genet. Genom. 2012, 287, 663–678. [Google Scholar] [CrossRef] [PubMed]
  20. Lim, S.D.; Lee, S.; Choi, W.-G.; Yim, W.C.; Cushman, J.C. Laying the Foundation for Crassulacean Acid Metabolism (CAM) Biodesign: Expression of the C4 Metabolism Cycle Genes of CAM in Arabidopsis. Front. Plant Sci. 2019, 10, 101. [Google Scholar] [CrossRef] [PubMed]
  21. Yuan, G.; Hassan, M.M.; Liu, D.; Lim, S.D.; Yim, W.C.; Cushman, J.C.; Markel, K.; Shih, P.M.; Lu, H.; Weston, D.J.; et al. Biosystems Design to Accelerate C 3 -to-CAM Progression. BioDesign Res. 2020, 2020, 3686791. [Google Scholar] [CrossRef]
  22. Meyer, S.; De Angeli, A.; Fernie, A.R.; Martinoia, E. Intra- and Extra-Cellular Excretion of Carboxylates. Trends Plant Sci. 2010, 15, 40–47. [Google Scholar] [CrossRef]
  23. Liu, J. Biology and Function of the OsALMT1 Gene in Rice (Oryza sativa L.). Ph.D. Thesis, University of Tasmania, Hobart, Australia, 2016. [Google Scholar]
  24. Ma, B.; Yuan, Y.; Gao, M.; Qi, T.; Li, M.; Ma, F. Genome-Wide Identification, Molecular Evolution, and Expression Divergence of Aluminum-Activated Malate Transporters in Apples. Int. J. Mol. Sci. 2018, 19, 2807. [Google Scholar] [CrossRef]
  25. Peng, W.; Wu, W.; Peng, J.; Li, J.; Lin, Y.; Wang, Y.; Tian, J.; Sun, L.; Liang, C.; Liao, H. Characterization of the Soybean GmALMT Family Genes and the Function of GmALMT5 in Response to Phosphate Starvation. J. Integr. Plant Biol. 2018, 60, 216–231. [Google Scholar] [CrossRef]
  26. Chen, L.-Y.; Xin, Y.; Wai, C.M.; Liu, J.; Ming, R. The Role of Cis-Elements in the Evolution of Crassulacean Acid Metabolism Photosynthesis. Hortic. Res. 2020, 7, 5. [Google Scholar] [CrossRef]
  27. Kobayashi, Y.; Hoekenga, O.A.; Itoh, H.; Nakashima, M.; Saito, S.; Shaff, J.E.; Maron, L.G.; Piñeros, M.A.; Kochian, L.V.; Koyama, H. Characterization of AtALMT1 Expression in Aluminum-Inducible Malate Release and Its Role for Rhizotoxic Stress Tolerance in Arabidopsis. Plant Physiol. 2007, 145, 843–852. [Google Scholar] [CrossRef]
  28. Zhang, F.; Jiang, S.; Li, Q.; Song, Z.; Yang, Y.; Yu, S.; Nie, Z.; Chu, M.; An, Y. Identification of the ALMT Gene Family in the Potato (Solanum tuberosum L.) and Analysis of the Function of StALMT6/10 in Response to Aluminum Toxicity. Front. Plant Sci. 2023, 14, 1274260. [Google Scholar] [CrossRef]
  29. Sasaki, T.; Tsuchiya, Y.; Ariyoshi, M.; Nakano, R.; Ushijima, K.; Kubo, Y.; Mori, I.C.; Higashiizumi, E.; Galis, I.; Yamamoto, Y. Two Members of the Aluminum-Activated Malate Transporter Family, SlALMT4 and SlALMT5, Are Expressed during Fruit Development, and the Overexpression of SlALMT5 Alters Organic Acid Contents in Seeds in Tomato (Solanum lycopersicum). Plant Cell Physiol. 2016, 57, 2367–2379. [Google Scholar] [CrossRef]
  30. Sunagawa, H.; Cushman, J.; Agarie, S. Crassulacean Acid Metabolism May Alleviate Production of Reactive Oxygen Species in a Facultative CAM Plant, the Common Ice Plant Mesembryanthemum crystallinum L. Plant Prod. Sci. 2010, 13, 256–260. [Google Scholar] [CrossRef]
  31. Borland, A.M.; Griffiths, H.; Hartwell, J.; Smith, J.A.C. Exploiting the Potential of Plants with Crassulacean Acid Metabolism for Bioenergy Production on Marginal Lands. J. Exp. Bot. 2009, 60, 2879–2896. [Google Scholar] [CrossRef]
  32. Yang, X.; Hu, R.; Yin, H.; Jenkins, J.; Shu, S.; Tang, H.; Liu, D.; Weighill, D.A.; Cheol Yim, W.; Ha, J.; et al. The Kalanchoë Genome Provides Insights into Convergent Evolution and Building Blocks of Crassulacean Acid Metabolism. Nat. Commun. 2017, 8, 1899. [Google Scholar] [CrossRef]
  33. Silvera, K.; Neubig, K.M.; Whitten, W.M.; Williams, N.H.; Winter, K.; Cushman, J.C. Evolution along the Crassulacean Acid Metabolism Continuum. Funct. Plant Biol. 2010, 37, 995–1010. [Google Scholar] [CrossRef]
  34. Gamisch, A.; Winter, K.; Fischer, G.A.; Comes, H.P. Evolution of Crassulacean Acid Metabolism (CAM) as an Escape from Ecological Niche Conservatism in Malagasy Bulbophyllum (Orchidaceae). New Phytol. 2021, 231, 1236–1248. [Google Scholar] [CrossRef]
  35. Bolaños-Villegas, P.; Chen, F.-C. Advances and Perspectives for Polyploidy Breeding in Orchids. Plants 2022, 11, 1421. [Google Scholar] [CrossRef]
  36. Zhang, G.-Q.; Liu, K.-W.; Li, Z.; Lohaus, R.; Hsiao, Y.-Y.; Niu, S.-C.; Wang, J.-Y.; Lin, Y.-C.; Xu, Q.; Chen, L.-J.; et al. The Apostasia Genome and the Evolution of Orchids. Nature 2017, 549, 379–383. [Google Scholar] [CrossRef]
  37. Zhang, Y.; Zhang, G.-Q.; Zhang, D.; Liu, X.-D.; Xu, X.-Y.; Sun, W.-H.; Yu, X.; Zhu, X.; Wang, Z.-W.; Zhao, X.; et al. Chromosome-Scale Assembly of the Dendrobium Chrysotoxum Genome Enhances the Understanding of Orchid Evolution. Hortic. Res. 2021, 8, 183. [Google Scholar] [CrossRef] [PubMed]
  38. Zhang, G.-Q.; Xu, Q.; Bian, C.; Tsai, W.-C.; Yeh, C.-M.; Liu, K.-W.; Yoshida, K.; Zhang, L.-S.; Chang, S.-B.; Chen, F.; et al. The Dendrobium Catenatum Lindl. Genome Sequence Provides Insights into Polysaccharide Synthase, Floral Development and Adaptive Evolution. Sci. Rep. 2016, 6, 19029. [Google Scholar] [CrossRef] [PubMed]
  39. Cai, J.; Liu, X.; Vanneste, K.; Proost, S.; Tsai, W.-C.; Liu, K.-W.; Chen, L.-J.; He, Y.; Xu, Q.; Bian, C.; et al. The Genome Sequence of the Orchid Phalaenopsis Equestris. Nat. Genet. 2015, 47, 65–72. [Google Scholar] [CrossRef] [PubMed]
  40. Qin, L.; Tang, L.; Xu, J.; Zhang, X.; Zhu, Y.; Zhang, C.; Wang, M.; Liu, X.; Li, F.; Sun, F.; et al. Cryo-EM Structure and Electrophysiological Characterization of ALMT from Glycine Max Reveal a Previously Uncharacterized Class of Anion Channels. Sci. Adv. 2022, 8, eabm3238. [Google Scholar] [CrossRef]
  41. Xu, B.; Long, Y.; Feng, X.; Zhu, X.; Sai, N.; Chirkova, L.; Betts, A.; Herrmann, J.; Edwards, E.J.; Okamoto, M.; et al. GABA Signalling Modulates Stomatal Opening to Enhance Plant Water Use Efficiency and Drought Resilience. Nat. Commun. 2021, 12, 1952. [Google Scholar] [CrossRef]
  42. Din, I.; Ullah, I.; Wang, W.; Zhang, H.; Shi, L. Genome-Wide Analysis, Evolutionary History and Response of ALMT Family to Phosphate Starvation in Brassica napus. Int. J. Mol. Sci. 2021, 22, 4625. [Google Scholar] [CrossRef] [PubMed]
  43. Mumm, P.; Imes, D.; Martinoia, E.; Al-Rasheid, K.A.S.; Geiger, D.; Marten, I.; Hedrich, R. C-Terminus-Mediated Voltage Gating of Arabidopsis Guard Cell Anion Channel QUAC1. Mol. Plant 2013, 6, 1550–1563. [Google Scholar] [CrossRef]
  44. Ligaba, A.; Dreyer, I.; Margaryan, A.; Schneider, D.J.; Kochian, L.; Piñeros, M. Functional, Structural and Phylogenetic Analysis of Domains Underlying the A l Sensitivity of the Aluminum-activated Malate/Anion Transporter, T a ALMT 1. Plant J. 2013, 76, 766–780. [Google Scholar] [CrossRef] [PubMed]
  45. Domingos, P.; Dias, P.N.; Tavares, B.; Portes, M.T.; Wudick, M.M.; Konrad, K.R.; Gilliham, M.; Bicho, A.; Feijó, J.A. Molecular and Electrophysiological Characterization of Anion Transport in Arabidopsis thaliana Pollen Reveals Regulatory Roles for PH, Ca2+ and GABA. New Phytol. 2019, 223, 1353–1371. [Google Scholar] [CrossRef]
  46. Ramesh, S.A.; Kamran, M.; Sullivan, W.; Chirkova, L.; Okamoto, M.; Degryse, F.; McLaughlin, M.; Gilliham, M.; Tyerman, S.D. Aluminum-Activated Malate Transporters Can Facilitate GABA Transport. Plant Cell 2018, 30, 1147–1164. [Google Scholar] [CrossRef] [PubMed]
  47. Jo, B.-S.; Choi, S.S. Introns: The Functional Benefits of Introns in Genomes. Genom. Inform. 2015, 13, 112–118. [Google Scholar] [CrossRef]
  48. Fontecha, G.; Silva-Navas, J.; Benito, C.; Mestres, M.A.; Espino, F.J.; Hernández-Riquer, M.V.; Gallego, F.J. Candidate Gene Identification of an Aluminum-Activated Organic Acid Transporter Gene at the Alt4 Locus for Aluminum Tolerance in Rye (Secale cereale L.). Theor. Appl. Genet. 2006, 114, 249–260. [Google Scholar] [CrossRef]
  49. Ryan, P.R.; Raman, H.; Gupta, S.; Sasaki, T.; Yamamoto, Y.; Delhaize, E. The Multiple Origins of Aluminium Resistance in Hexaploid Wheat Include Aegilops Tauschii and More Recent Cis Mutations to TaALMT1: Evolution of Aluminium Resistance in Wheat. Plant J. 2010, 64, 446–455. [Google Scholar] [CrossRef]
  50. Linlin, X.; Xin, Q.; Mingyue, Z.; Shaoling, Z. Genome-Wide Analysis of Aluminum-Activated Malate Transporter Family Genes in Six Rosaceae Species, and Expression Analysis and Functional Characterization on Malate Accumulation in Chinese White Pear. Plant Sci. 2018, 274, 451–465. [Google Scholar] [CrossRef]
  51. Liang, C.; Piñeros, M.A.; Tian, J.; Yao, Z.; Sun, L.; Liu, J.; Shaff, J.; Coluccio, A.; Kochian, L.V.; Liao, H. Low pH, Aluminum, and Phosphorus Coordinately Regulate Malate Exudation through GmALMT1 to Improve Soybean Adaptation to Acid Soils. Plant Physiol. 2013, 161, 1347–1361. [Google Scholar] [CrossRef]
  52. Chen, Q.; Zhang, X.D.; Wang, S.S.; Wang, Q.F.; Wang, G.Q.; Nian, H.J.; Li, K.Z.; Yu, Y.X.; Chen, L.M. Transcriptional and Physiological Changes of Alfalfa in Response to Aluminium Stress. J. Agric. Sci. 2011, 149, 737–751. [Google Scholar] [CrossRef]
  53. Zhang, T.; Huang, W.; Zhang, L.; Li, D.-Z.; Qi, J.; Ma, H. Phylogenomic Profiles of Whole-Genome Duplications in Poaceae and Landscape of Differential Duplicate Retention and Losses among Major Poaceae Lineages. Nat. Commun. 2024, 15, 3305. [Google Scholar] [CrossRef]
  54. Zhang, W.; Wu, J.; He, J.; Liu, C.; Yi, W.; Xie, J.; Wu, Y.; Xie, T.; Ma, J.; Zhong, Z.; et al. AcMYB266, a Key Regulator of the Red Coloration in Pineapple Peel: A Case of Subfunctionalization in Tandem Duplicated Genes. Hortic. Res. 2024, 11, uhae116. [Google Scholar] [CrossRef]
  55. Xun, H.; Zhang, Z.; Zhou, Y.; Qian, X.; Dong, Y.; Feng, X.; Pang, J.; Wang, S.; Liu, B. Identification and Functional Characterization of R3 MYB Transcription Factor Genes in Soybean. J. Plant Biol. 2018, 61, 85–96. [Google Scholar] [CrossRef]
  56. Ma, X.; An, F.; Wang, L.; Guo, D.; Xie, G.; Liu, Z. Genome-Wide Identification of Aluminum-Activated Malate Transporter (ALMT) Gene Family in Rubber Trees (Hevea brasiliensis) Highlights Their Involvement in Aluminum Detoxification. Forests 2020, 11, 142. [Google Scholar] [CrossRef]
  57. Zhang, J.; Baetz, U.; Krügel, U.; Martinoia, E.; De Angeli, A. Identification of a Probable Pore-Forming Domain in the Multimeric Vacuolar Anion Channel AtALMT9. Plant Physiol. 2013, 163, 830–843. [Google Scholar] [CrossRef] [PubMed]
  58. Li, C.; Dougherty, L.; Coluccio, A.E.; Meng, D.; El-Sharkawy, I.; Borejsza-Wysocka, E.; Liang, D.; Piñeros, M.A.; Xu, K.; Cheng, L. Apple ALMT9 Requires a Conserved C-Terminal Domain for Malate Transport Underlying Fruit Acidity. Plant Physiol. 2020, 182, 992–1006. [Google Scholar] [CrossRef] [PubMed]
  59. Tian, Q.; Zhang, X.; Ramesh, S.; Gilliham, M.; Tyerman, S.D.; Zhang, W.-H. Ethylene Negatively Regulates Aluminium-Induced Malate Efflux from Wheat Roots and Tobacco Cells Transformed with TaALMT1. J. Exp. Bot. 2014, 65, 2415–2426. [Google Scholar] [CrossRef]
  60. Zheng, L.; Ma, W.; Liu, P.; Song, S.; Wang, L.; Yang, W.; Ren, H.; Wei, X.; Zhu, L.; Peng, J.; et al. Transcriptional Factor MdESE3 Controls Fruit Acidity by Activating Genes Regulating Malic Acid Content in Apple. Plant Physiol. 2024, 196, kiae282. [Google Scholar] [CrossRef]
  61. Eisenach, C.; Baetz, U.; Huck, N.V.; Zhang, J.; De Angeli, A.; Beckers, G.J.M.; Martinoia, E. ABA-Induced Stomatal Closure Involves ALMT4, a Phosphorylation-Dependent Vacuolar Anion Channel of Arabidopsis. Plant Cell 2017, 29, 2552–2569. [Google Scholar] [CrossRef]
  62. Hu, D.; Li, Y.; Zhang, Q.; Li, M.; Sun, C.; Yu, J.; Hao, Y. The R2R3-MYB Transcription Factor Md MYB 73 Is Involved in Malate Accumulation and Vacuolar Acidification in Apple. Plant J. 2017, 91, 443–454. [Google Scholar] [CrossRef] [PubMed]
  63. Davies, B.N.; Griffiths, H. Competing Carboxylases: Circadian and Metabolic Regulation of Rubisco in C3 and CAM Mesembryanthemum crystallinum L. Plant Cell Environ. 2012, 35, 1211–1220. [Google Scholar] [CrossRef] [PubMed]
  64. Armarego-Marriott, T. Stop the Clock: Optimized Carbon Fixation and Circadian Rhythm in a CAM Plant. Plant Cell 2017, 29, 2314–2315. [Google Scholar] [CrossRef]
  65. Merilo, E.; Yarmolinsky, D.; Jalakas, P.; Parik, H.; Tulva, I.; Rasulov, B.; Kilk, K.; Kollist, H. Stomatal VPD Response: There Is More to the Story Than ABA. Plant Physiol. 2018, 176, 851–864. [Google Scholar] [CrossRef] [PubMed]
  66. Cai, S.; Chen, G.; Wang, Y.; Huang, Y.; Marchant, D.B.; Wang, Y.; Yang, Q.; Dai, F.; Hills, A.; Franks, P.J.; et al. Evolutionary Conservation of ABA Signaling for Stomatal Closure. Plant Physiol. 2017, 174, 732–747. [Google Scholar] [CrossRef]
  67. Zhang, J.; Chen, X.; Song, Y.; Gong, Z. Integrative Regulatory Mechanisms of Stomatal Movements under Changing Climate. J. Integr. Plant Biol. 2024, 66, 368–393. [Google Scholar] [CrossRef]
  68. Geiger, D.; Scherzer, S.; Mumm, P.; Stange, A.; Marten, I.; Bauer, H.; Ache, P.; Matschi, S.; Liese, A.; Al-Rasheid, K.A.S.; et al. Activity of Guard Cell Anion Channel SLAC1 Is Controlled by Drought-Stress Signaling Kinase-Phosphatase Pair. Proc. Natl. Acad. Sci. USA 2009, 106, 21425–21430. [Google Scholar] [CrossRef]
  69. Malcheska, F.; Ahmad, A.; Batool, S.; Müller, H.M.; Ludwig-Müller, J.; Kreuzwieser, J.; Randewig, D.; Hänsch, R.; Mendel, R.R.; Hell, R.; et al. Drought-Enhanced Xylem Sap Sulfate Closes Stomata by Affecting ALMT12 and Guard Cell ABA Synthesis. Plant Physiol. 2017, 174, 798–814. [Google Scholar] [CrossRef]
  70. Ward’, J.M. Calcium-Activated K+ Channels and Calcium-Lnduced Calcium Release by Slow Vacuolar Lon Channels in Guard Cell Vacuoles Lmplicated in the Control of Stomatal Closure. Plant Cell 1994, 6, 669–683. [Google Scholar] [CrossRef]
  71. Xu, B.; Sai, N.; Gilliham, M. The Emerging Role of GABA as a Transport Regulator and Physiological Signal. Plant Physiol. 2021, 187, 2005–2016. [Google Scholar] [CrossRef]
  72. Han, Y.; Wang, W.; Sun, J.; Ding, M.; Zhao, R.; Deng, S.; Wang, F.; Hu, Y.; Wang, Y.; Lu, Y.; et al. Populus Eu-phratica XTH Overexpression Enhances Salinity Tolerance by the Development of Leaf Succulence in Transgenic Tobacco Plants. J. Exp. Bot. 2013, 64, 4225–4238. [Google Scholar] [CrossRef]
  73. 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]
  74. 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]
  75. Troshin, P.V.; Procter, J.B.; Barton, G.J. Java Bioinformatics Analysis Web Services for Multiple Sequence Alignment—JABAWS: MSA. Bioinformatics 2011, 27, 2001–2002. [Google Scholar] [CrossRef] [PubMed]
  76. He, Z.; Zhang, H.; Gao, S.; Lercher, M.J.; Chen, W.-H.; Hu, S. Evolview v2: An Online Visualization and Management Tool for Customized and Annotated Phylogenetic Trees. Nucleic Acids Res. 2016, 44, W236–W241. [Google Scholar] [CrossRef]
  77. Lescot, M. 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]
  78. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; 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]
  79. Wang, D.; Zhang, Y.; Zhang, Z.; Zhu, J.; Yu, J. KaKs_Calculator 2.0: A Toolkit Incorporating Gamma-Series Methods and Sliding Window Strategies. Genom. Proteom. Bioinform. 2010, 8, 77–80. [Google Scholar] [CrossRef]
Ijms 25 09662 g001
Figure 1. Chromosomal localization of DchALMT genes. The red line represents tandem duplicated genes.
Figure 1. Chromosomal localization of DchALMT genes. The red line represents tandem duplicated genes.
Ijms 25 09662 g001
Ijms 25 09662 g002
Figure 2. Multiple sequence alignments and typical motifs of orchid ALMT protein.
Figure 2. Multiple sequence alignments and typical motifs of orchid ALMT protein.
Ijms 25 09662 g002
Ijms 25 09662 g003
Figure 3. Phylogenetic trees of ALMT proteins based on 34 orchid ALMT proteins, 13 AtALMT proteins, 9 OsALMT proteins, and 14 other ALMT proteins were constructed by Maximum Likelihood methods. The ALMT gene family was classified into four clades. A. shenzhenica, D. catenatum, D. chrysotoxum, P. equestris, A. thaliana, Oryza sativa, Aegilops tauschii, Citrus sinensis, Glycine max, Holcus lanatus, Malus domestica, Medicago sativa, Triticum aestivum, Zea mays, Brassica napus, Brassica oleracea, Hordeum vulgare, and Secale cereale are labeled as Ash, Dca, Dch, Pe, At, Os, Aet, Cs, Gm, Hl, Ma, Ms, Ta, Zm, Bn, Bo, Hv, and Sc, respectively.
Figure 3. Phylogenetic trees of ALMT proteins based on 34 orchid ALMT proteins, 13 AtALMT proteins, 9 OsALMT proteins, and 14 other ALMT proteins were constructed by Maximum Likelihood methods. The ALMT gene family was classified into four clades. A. shenzhenica, D. catenatum, D. chrysotoxum, P. equestris, A. thaliana, Oryza sativa, Aegilops tauschii, Citrus sinensis, Glycine max, Holcus lanatus, Malus domestica, Medicago sativa, Triticum aestivum, Zea mays, Brassica napus, Brassica oleracea, Hordeum vulgare, and Secale cereale are labeled as Ash, Dca, Dch, Pe, At, Os, Aet, Cs, Gm, Hl, Ma, Ms, Ta, Zm, Bn, Bo, Hv, and Sc, respectively.
Ijms 25 09662 g003
Ijms 25 09662 g004
Figure 4. Gene structure, conserved motifs, and domains of ALMTs. (A) The ML tree contains 34 orchid ALMTs. (B) Squares of different colors represent conserved motifs of ALMTs. (C) Squares of different colors represent the conserved domains of ALMTs. (D) Squares of different colors represent the gene structures.
Figure 4. Gene structure, conserved motifs, and domains of ALMTs. (A) The ML tree contains 34 orchid ALMTs. (B) Squares of different colors represent conserved motifs of ALMTs. (C) Squares of different colors represent the conserved domains of ALMTs. (D) Squares of different colors represent the gene structures.
Ijms 25 09662 g004
Ijms 25 09662 g005
Figure 5. Cis-acting elements in the 2k bp of upstream and downstream regions of orchid ALMT genes. (A) Elements with similar regulatory functions are displayed in the same color. (B) Numbers of each type of element.
Figure 5. Cis-acting elements in the 2k bp of upstream and downstream regions of orchid ALMT genes. (A) Elements with similar regulatory functions are displayed in the same color. (B) Numbers of each type of element.
Ijms 25 09662 g005
Ijms 25 09662 g006
Figure 6. Expression analysis of DcaALMT genes at four tissues. The error bars indicate three RT-qPCR biological replicates. Statistical analysis using ANNOVA multiple comparisons for qPCR expression levels between each tissue. The asterisk indicates the P value in the significance test (* p < 0.05, *** p < 0.001, **** p < 0.0001).
Figure 6. Expression analysis of DcaALMT genes at four tissues. The error bars indicate three RT-qPCR biological replicates. Statistical analysis using ANNOVA multiple comparisons for qPCR expression levels between each tissue. The asterisk indicates the P value in the significance test (* p < 0.05, *** p < 0.001, **** p < 0.0001).
Ijms 25 09662 g006
Ijms 25 09662 g007
Figure 7. (A) The CO2 uptake rate of D. catenatum leaves over a day and night. (B) The fluctuation of titratable total acid over day and night in D. catenatum leaves. (C) The expression profiles of DcaALMT genes over day and night. The black line represents night.
Figure 7. (A) The CO2 uptake rate of D. catenatum leaves over a day and night. (B) The fluctuation of titratable total acid over day and night in D. catenatum leaves. (C) The expression profiles of DcaALMT genes over day and night. The black line represents night.
Ijms 25 09662 g007
Ijms 25 09662 g008
Figure 8. Expression analysis of DcaALMT genes at 0 h, 1.5 h, 3 h, and 6 h after ABA treatment. The error bars indicate three RT-qPCR biological replicates. Statistical analysis using ANNOVA multiple comparisons for qPCR expression levels between each time point. The asterisk indicates the P value in the significance test (* p < 0.05, *** p < 0.001, **** p < 0.0001).
Figure 8. Expression analysis of DcaALMT genes at 0 h, 1.5 h, 3 h, and 6 h after ABA treatment. The error bars indicate three RT-qPCR biological replicates. Statistical analysis using ANNOVA multiple comparisons for qPCR expression levels between each time point. The asterisk indicates the P value in the significance test (* p < 0.05, *** p < 0.001, **** p < 0.0001).
Ijms 25 09662 g008
Table 1. Information on ALMT genes in four orchid genomes.
Table 1. Information on ALMT genes in four orchid genomes.
SpeciesNameGene IDAA(aa)MW (Da)pIIIAIGRAVYTMDLocalization
D. catenatumDcaALMT1Dca00704746350,348.469.0238.1599.330.1796PM
DcaALMT2Dca00664748252,841.688.4237101.20.2356PM
DcaALMT3Dca02471652759,656.098.8533.396.760.046PM
DcaALMT4Dca00040548654,431.845.7740.58102.940.0196PM
DcaALMT5Dca00664928749,094.168.631.297.580.2365PM
DcaALMT6DcaN0360252759,414.266.336.7990.340.0865Vacuole
DcaALMT7Dca00073954961,017.746.2338.94103.330.166ER
DcaALMT8Dca00664846350,712.198.933.7899.050.195PM
DcaALMT9Dca00040443548,375.728.9932.78106.920.2547PM
DcaALMT10Dca01327244148,376.828.4534.65114.720.3295PM
D. chrysotoxumDchALMT1Maker6680752759,567.048.8535.8896.56−0.046PM
DchALMT2Maker5353344548,443.58.4437.598.270.2266PM
DchALMT3Maker6709446350,479.826.7644.694.160.1996PM
DchALMT4Maker9516739542,849.99.3942.92100.230.1456PM
DchALMT5Maker8023636840,776.738.3425103.320.3756PM
DchALMT6Maker6616752859,246.116.2637.5689.03−0.0937ER
DchALMT7Maker5359688796,227.928.4733.4197.630.2312Extracellular
DchALMT8Maker5815549354,475.685.739.9398.560.0155PM
DchALMT9Maker3959645449,362.178.8440.5696.580.1146PM
DchALMT10Maker5817243348,090.469.232.621090.2425PM
P. equestrisPeALMT1XM_020717192.145549,782.377.6344.2293.370.0726PM
PeALMT2XM_020720250.145649,218.339.1335.64100.810.2386PM
PeALMT3XM_020726555.147952,949.236.1245.1498.350.046PM
PeALMT4XM_020728361.146650,732.388.2336.3896.90.2317PM
PeALMT5XM_020727176.152359,234.728.933.8395.63−0.0466PM
PeALMT6XM_020730006.159065,819.36.4236.63100.150.0696ER
PeALMT7XM_020743275.128230,544.028.8226.07123.810.5336PM
A. shenzhenicaAsALMT1Ash00975851358,212.418.6433.5497.29−0.0836PM
AsALMT2Ash01892553259,665.426.3748.4993.07−0.1474ER
AsALMT3Ash01141844748,104.868.0941.7101.790.2336PM
AsALMT4Ash00772325327,270.698.2437.2898.730.325PM
AsALMT5Ash02023343148,267.479.1237.49105.230.1697PM
AsALMT6Ash01125245749,292.448.3733.41102.980.2566PM
AsALMT7Ash00526352058,377.166.5240.5491.93−0.0566ER
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

Peng, F.-C.; Yuan, M.; Zhou, L.; Zheng, B.-Q.; Wang, Y. Identification and Analysis of Aluminum-Activated Malate Transporter Gene Family Reveals Functional Diversification in Orchidaceae and the Expression Patterns of Dendrobium catenatum Aluminum-Activated Malate Transporters. Int. J. Mol. Sci. 2024, 25, 9662. https://doi.org/10.3390/ijms25179662

AMA Style

Peng F-C, Yuan M, Zhou L, Zheng B-Q, Wang Y. Identification and Analysis of Aluminum-Activated Malate Transporter Gene Family Reveals Functional Diversification in Orchidaceae and the Expression Patterns of Dendrobium catenatum Aluminum-Activated Malate Transporters. International Journal of Molecular Sciences. 2024; 25(17):9662. https://doi.org/10.3390/ijms25179662

Chicago/Turabian Style

Peng, Fu-Cheng, Meng Yuan, Lin Zhou, Bao-Qiang Zheng, and Yan Wang. 2024. "Identification and Analysis of Aluminum-Activated Malate Transporter Gene Family Reveals Functional Diversification in Orchidaceae and the Expression Patterns of Dendrobium catenatum Aluminum-Activated Malate Transporters" International Journal of Molecular Sciences 25, no. 17: 9662. https://doi.org/10.3390/ijms25179662

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