1. Introduction
In plants, the calcium (Ca
2+) ion is an important element for maintaining the intracellular ion homeostasis of plant cells, as well as for regulating structural processes, metabolism, and signaling [
1,
2]. The major function of the Ca
2+ ion is to enhance the plant’s tolerance to stabilize cell walls and membranes against various plant stresses [
2,
3,
4]. The Ca
2+ ion regulates the biochemical reactions in the plant cell by controlling the activity of ion-channel-related proteins, such as cyclic nucleotide-gated channels (CNGCs), ionotropic glutamate receptors, two-pore channel 1 (TPC1), annexins, and mechanosensitive channels [
2]. The Ca
2+ ion is a versatile secondary messenger in biotic and abiotic stress responses [
3,
4]. Increased Ca
2+ concentrations through various abiotic stresses are recognized, amplified, and transmitted through Ca
2+-binding proteins [
4]. These Ca
2+-binding proteins play a role as Ca
2+ sensors, including calmodulin (CaM) and calmodulin-like proteins (CMLs), calcineurin B-like proteins (CBLs), and calcium-dependent protein kinases (CDPKs) [
5]. These signal cascades in plants physiologically regulate the division and elongation of cells, stomatal movement, and developmental processes in plant tolerance or adaptation responses to abiotic stresses [
5,
6].
The CaMs with one or more EF-hand Ca
2+-binding motifs are associated with a multitude of biochemical reactions, such as the control of transcriptional expression and enzyme activity [
6,
7]. Although the Ca
2+ ion lacks cellular levels to activate CaM under natural conditions, when exposed to abiotic stresses, the Ca
2+ levels are rapidly increased, generating the Ca
2+ signals, and there is a boost in the activity of CaM [
8]. CaM is the best characterized of the Ca
2+ receptors, and the Ca
2+–CaM complex regulates the ion channels, transporters, transcription factors (TFs), and various enzymes, including protein kinase, phosphatases, and metabolic enzymes [
9,
10]. The CaM-driven TFs such as calmodulin-binding transcription activators (CAMTAs), DREB, MYB, WRKY, NAC, bZIP, bHLH, and MADS-box play a pivotal role in multiple cellular processes, including plant stress responses and developmental processes [
5,
10].
A total of 465
CAMTA genes have been identified from 112 plant species, including dicots, monocots, basal angiosperms, gymnosperms, lycophyta, monilophytes, bryophytes, and chlorophytes [
10]. Of the total, 302
CAMTAs were obtained from 49 dicots among plant species; in addition, 15 CAMTAs are the best isolated from soybean (
Glycine max L.), respectively [
10]. CAMTA proteins in soybean have been known to comprise conserved functional domains, including the DNA-binding (CG-1) domain, transcription-associated immunoglobulin (TIG) domain, calmodulin-binding (CaMB) domain, ankyrin (ANK) repeats, and IQ motifs [
11]. The CG-1 domain recognizes and binds to specific DNA sequences of the CGCG box for activating gene transcription [
12,
13]. The ANK repeats and the TIG domain with binding to the non-specific DNA are involved in protein–protein interactions with either homo- or heterodimers [
12,
14]. A diverse number of IQ motifs in the C-terminal regions of CAMTAs contain amino acid sequences, such as IQXXXRGXXX, as Ca
2+-independent calmodulin-binding (CaMB) sites [
5,
12,
15]. The CaMB domain with its α-helix structure is an important Ca
2+-dependent CaMB domain for determining characteristics of CAMTAs in plants [
12,
16].
The expression levels of
CAMTAs show significant differences in the extensive developmental stages and plant stress responses [
5,
10]. In
Arabidopsis, the expression of
AtCAMTA6 (
AT3G16940) shows a high level in mature leaves and siliques compared with that in other tissues [
10]. In addition, the expression levels of
AtCAMTA2 (
AT5G64220) in seeds and young leaves are significantly lower than in other tissues [
10]. In tomato (
Solanum lycopersicum),
PtCAMTA4 (
JN566050) shows high levels of transcriptional expression in turning and orange fruit stages, while
PtCAMTA2 (
JN566047) is highly expressed in all developmental fruit stages except the mature green and breaker stages [
16]. In tea (
Camellia sinensis) plant, the expression of
CsCAMTA1 (
TEA025813.1) is at higher levels in old leaves than that in both other tissues and other
CsCAMTA genes [
17]. In black cottonwood (
Populus trichocarpa),
PtCAMTA1 (
POPTR_0001s13700) is highly expressed in the mature leaves, while
PtCAMTA2 (
POPTR_0005s07660) and
PtCAMTA3 (
POPTR_0007s05410) are more weakly expressed in the roots than in different tissues [
18]. In rice (
Oryza sativa), the expression levels of all seven
OsCAMTA genes are displayed more in leaves after flowering compared with those before flowering [
10]. The
MuCAMTA1 (
Musa acuminate;
LOC1039776) and
phavuCAMTA1 (
Phaseolus vulgaris;
PHAVU_006G206400g) genes are highly expressed in plant responses under drought stress conditions [
19,
20]. In the root of maize (
Zea mays L.), the expression of
ZmCAMTA1 (
GRMZM2G171600) is significantly decreased to cold stress, but increases to salt stress and jasmonic acid (JA) treatment [
21]. In
Medicago truncatula, the expression of
MtCAMTA7 (
Medtr8g090205) is induced by ABA treatment compared with other
MtCAMTA genes [
22]. Therefore, the expression of
CAMTA genes presents substantial specificity in the developmental stages and tissues, as well as in stress responses.
CAMTAs in plants are known to be positive or negative regulators in plant development processes, as well as in biotic and abiotic stress responses [
23,
24,
25].
AtCAMTA1 (
AT5G09410) positively regulates the ABA-dependent response for enhancing plant tolerance to drought stress [
23].
atcamta1 mutants showed hypersensitivity to drought stress by regulating stress-responsive genes, such as
RD26,
ERD7,
RAB18,
LTPs, and
COR78 [
23].
AtCAMTA3 (
AT2G22300) functions as a negative regulator in plant defense responses against bacteria [
26] and fungi [
27].
atcamta3 mutants showed enhanced resistance to
Pseudomonas syringae pv.
tomato DC3000 (
PstDC3000) and
Botrytis cinerea by the high expression of
pathogenesis-related (
PR) genes [
26,
27]. In contrast,
AtCAMTA3 acts as a positive regulator in plant tolerance to freezing stress [
24]. A mutation of
atcamta3 weakened the freezing tolerance caused by reducing the cold-induced accumulation of
CBF genes transcript levels during exposure to low temperature [
24].
AtCAMTA1 and
AtCAMTA5 (
AT4G16150) positively regulate the pollen developmental process by increasing the expression of
Arabidopsis V-PPase 1 (
AVP1) [
28]. Transcriptome profiling using a
camta6 mutant under salt stress showed that 638 up-regulated and 1,242 down-regulated genes were involved in
AtCAMTA6-dependent manners [
29]. In rice,
OsCBT/
OsCAMTA5 (
LOC_Os07g30774) plays a negative role in the plant defense response against fungal pathogens by the high expression of
PR genes, such as
PR1,
PR4,
PR10, and
PBZ [
30]. Rice
oscbt-1 mutants showed enhanced plant disease resistance against both rice blast fungus (
Magnaporthe grisea) and bacteria (
Xanthomonas oryzae) [
30]. Soybean
GmCAMTA12 (
Glyma.17G031900) regulated the expression of its regulatory network-related genes with the
CGCG/
CGTTG motif under drought stress in
Arabidopsis and soybean hairy roots [
31]. In peach (
Prunus persica L. Batsch),
PpCAMTA1 (
Prupe.1G108700) in
atcamta2/3 double mutants suppressed salicylic acid (SA) biosynthesis and the expression of SA-related genes in the plant defense response against
PstDC3000 [
32]. Thus, CAMTAs in plants play crucial roles as transcription factors in both biotic and abiotic stress responses.
Various environmental changes have many implications for growth, developmental processes, and yield in crop plants. Drought stress is an important restricting factor for crop growth and the developmental tissue processes, holding both expansion and elongation down [
33]. Since drought stress triggers low turgor pressure in the cell, the osmotic maintenance of cell turgor in crops is of great importance for plant growth and crop survival [
33]. In many crops, soybean is highly sensitive to drought stress responses through various morphological changes [
34]. When exposed to drought stress, soybean shows momentous changes in most of the tissue morphology, such as a reduction in new branches and trifoliate leaves per plant, growth inhibition of initiation leaves, and reductions in leaf area via repression of lamina expansion [
34]. Drought stress has a negative effect on the developmental processes of stems and leaves by decreasing their elongation and expansion [
33].
Soybean (
Glycine max L.) is one of the most economical crops with a nutritional balance of protein, fat, and carbohydrate. It is mainly used as a portion in important foods, feeds, and bio-feedstocks [
35]. The Williams 82 cultivar is a typical soybean breed that was used to be established the reference genome sequence [
36]. The 15
GmCAMTAs were identified in soybean by the SoyBase database (
http://soybase.org/), and the analysis of gene expression showed that the expression patterns of
GmCAMTAs were altered during abiotic stress (salt, drought, cold, and oxidative) and hormone (ABA, SA, and JA) responses (
Supplementary Table S1) [
11]. Recently, it was reported that GmCAMTA12 plays an important role in the tolerance responses of
Arabidopsis and soybean to drought stress [
31]. However, GmCAMTAs still do not have enough study for the functional characterization in plant development processes and stress tolerance responses. In this study,
GmCAMTA2 and
GmCAMTA8 characterized the substantial specificity and reliability in the developmental stages and abiotic stress responses. The results of our study can be effective for further functional analysis of GmCAMTAs in soybean.
3. Discussion
The characterization of various
CAMTA genes, as core transcription factors with CaM binding sites, has been reported in recent years by genome-wide analysis and expression profiles [
10,
11,
16,
17,
18,
19,
20,
21,
22,
25]. Although 15
GmCAMTA genes in soybean have been identified, the functions of
GmCAMTAs remained almost unknown, except for
GmCAMTA12 [
31]. In this study, we identified 15
GmCAMTA genes and established that
GmCAMTA transcripts were controlled by circadian rhythms (
Figure 1). Previous studies of functions of
CAMTA genes were little known about circadian regulation in crops. Plant cellular responses, such as gene expression, development, and metabolic processes, were circadian regulated by circadian oscillators [
40,
41]. Circadian responses in plants were attributable to cellular oscillations in the intracellular Ca
2+ concentration ([Ca
2+]) in the surrounding environment, including photoperiod and light intensity [
42,
43]. The calmodulin-like 24 (CML24), one of the Ca
2+-binding proteins, regulated the circadian period in Ca
2+-dependent signaling [
41]. The most important properties of CAMTAs in plants have a high correlation with CaM proteins, as a Ca
2+ sensor to regulate the gene expression in Ca
2+-dependent plant responses [
5]. These facts were important to explain that the circadian regulation of
GmCAMTAs transcripts was mediated by the cellular oscillation of [Ca
2+].
Many
CAMTAs in crops have been directly or indirectly associated with various developmental and cellular processes, whether through autonomously regulating a major transcription factor or targeting other transcription factors [
5]. The fifteen
GmCAMTAs in soybean [
11], seven
PtCAMTAs in
Populus trichocarpa [
18], nine
ZmCAMTAs in maize [
21], seven
MtCAMTAs in
Medicago truncatula [
22], and fifteen
VvCAMTAs in
Vitis vinifera [
44] were endowed with substantial specificity in the developmental stages by highly tissue-specific expression. Root-specific-induced
PtCAMTA2 and
PtCAMTA3 genes were negatively regulated under cold stress in a short time [
18]. In this study,
GmCAMTA2,
GmCAMTA8, and
GmCAMTA12 were especially highly expressed in stem tissues compared to other
GmCAMTAs (
Figure 3). The growth of plant tissues, such as flowers, stems, and leaves, were greatly affected by abiotic stress [
43]. The specificity in tissues is an important characteristic for investigating the functions of GmCMATAs in diverse molecular processes against environmental stresses. Since plants are constantly exposed to environmental stimuli during their lifetime, they have developed various adaptive mechanisms between the developmental process and stress responses [
44,
45].
Drought stress is an important restricting factor for crop growth and developmental tissue processes, holding both expansion and elongation down [
46]. Since drought stress triggers low turgor pressure in the cell, the osmotic maintenance of cell turgor in crops is of great importance for plant growth and crop survival [
46]. In many crops, soybean is highly sensitive to drought stress responses through various morphological changes [
47]. When exposed to drought stress, soybean showed momentous changes in most of the tissue morphology, such as a reduction in new branches and trifoliate leaves per plant, growth inhibition of initiation leaves, and reduction in leaf area through repression of lamina expansion [
47]. Drought stress has a negative effect on the developmental processes of stems and leaves by decreasing their elongation and expansion [
46]. In particular, drought stress often leads to the expression of common morphogenesis- and stress-responsive genes through molecular mechanisms, including reactive oxygen species (ROS)- and phytohormone-dependent signaling [
45].
ATP-dependent metalloprotease 4 (
AtFtsH4) played a tissue-specific role in maintaining the stem cell activity of shoot apical meristem (SAM) stages by ROS accumulation under high-temperature conditions [
48].
MscS-like 2 (
MSL2) and
MSL3 were mediated to cell proliferation by regulation of
WUSCHEL (
WUS) gene expression in osmotic stress responses [
49]. The mutation of
NADPH genes showed unstable shoot development through decreasing the stem cell population [
50]. In case of
CAMTA genes, the relationships between development and stress responses have been studied more in fruit species than in crops. The tomato
SISR2 and
SISR3L transcripts were affected by ethylene signaling during fruit development and ripening [
16]. The expression of peach
PpCAMTA1 was induced in cold stress, but
PpCAMTA3 was repressed with UV-B irradiation during fruit development [
32]. Our results indicated that
GmCAMTA2,
GmCAMTA8, and
GmCAMTA12 were significantly expressed in both stem tissues and drought stress responses (
Figure 2 and 3). When
GmCAMTA2 and
GmCAMTA8 were overexpressed in
Arabidopsis,
GmCAMTA2-OX and
GmCAMTA8-OX displayed the sensitive phenotype under drought stress conditions by losing faster water than the WT plant (
Figure 4). These results suggest that
GmCAMTA2 and
GmCAMTA8 act as negative regulators in drought stress responses. Considering the tolerant phenotype produced by
GmCAMTA12 overexpression [
31] and sensitive phenotypes by
GmCAMTA2 and
GmCAMTA8 overexpression (
Figure 4) under drought stress conditions,
GmCAMTA2,
GmCAMTA8, and
GmCAMTA12 may be associated with drought tolerance by different regulating mechanisms. Our results suggest that GmCAMTA2 and GmCAMTA8 play an important role in both the development process and drought stress responses by circadian regulation.