Next Article in Journal
Rosavin Alleviates LPS-Induced Acute Lung Injure by Modulating the TLR-4/NF-κB/MAPK Singnaling Pathways
Next Article in Special Issue
OsLRR-RLP2 Gene Regulates Immunity to Magnaporthe oryzae in Japonica Rice
Previous Article in Journal
Construction of a g-C3N4/Bi(OH)3 Heterojunction for the Enhancement of Visible Light Photocatalytic Antibacterial Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

OsACA9, an Autoinhibited Ca2+-ATPase, Synergically Regulates Disease Resistance and Leaf Senescence in Rice

College of Life Sciences, Zhejiang Normal University, Jinhua 321004, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(3), 1874; https://doi.org/10.3390/ijms25031874
Submission received: 2 January 2024 / Revised: 27 January 2024 / Accepted: 30 January 2024 / Published: 3 February 2024

Abstract

:
Calcium (Ca2+) is a versatile intracellular second messenger that regulates several signaling pathways involved in growth, development, stress tolerance, and immune response in plants. Autoinhibited Ca2+-ATPases (ACAs) play an important role in the regulation of cellular Ca2+ homeostasis. Here, we systematically analyzed the putative OsACA family members in rice, and according to the phylogenetic tree of OsACAs, OsACA9 was clustered into a separated branch in which its homologous gene in Arabidopsis thaliana was reported to be involved in defense response. When the OsACA9 gene was knocked out by CRISPR/Cas9, significant accumulation of reactive oxygen species (ROS) was detected in the mutant lines. Meanwhile, the OsACA9 knock out lines showed enhanced disease resistance to both rice bacterial blight (BB) and bacterial leaf streak (BLS). In addition, compared to the wild-type (WT), the mutant lines displayed an early leaf senescence phenotype, and the agronomy traits of their plant height, panicle length, and grain yield were significantly decreased. Transcriptome analysis by RNA-Seq showed that the differentially expressed genes (DEGs) between WT and the Osaca9 mutant were mainly enriched in basal immune pathways and antibacterial metabolite synthesis pathways. Among them, multiple genes related to rice disease resistance, receptor-like cytoplasmic kinases (RLCKs) and cell wall-associated kinases (WAKs) genes were upregulated. Our results suggest that the Ca2+-ATPase OsACA9 may trigger oxidative burst in response to various pathogens and synergically regulate disease resistance and leaf senescence in rice.

1. Introduction

In plants, Ca2+ signaling is involved in all aspects of life, including growth regulation, development, reproduction, stress responses, and the establishment of beneficial plant-microbe interactions [1]. Ca2+ influx is the signal of the earliest responses stimulated by pathogen/microbe-associated molecular patterns (PAMPs/MAMPs) and functions in stomatal immunity [2,3]. Upon elicitor perception, the initial activation of RESPIRATORY BURST OXIDASE HOMOLOGUE (RBOH) primes the system for subsequent activation through Ca2+ signaling, leading to reactive oxygen species (ROS) production, which further enhances the immunity [2,4]. It is also a vital signal in response to abiotic stresses, such as touch, cold, drought, hypoxia, and salinity [5,6,7]. In addition, Ca2+ participates in the regulation of polarized cell growth, the elongation of pollen tubes and root hairs, and plant senescence [1,8]. Proper Ca2+ concentration could inhibit chlorophyll and protein degradation to inhibit leaf senescence, while an excessive Ca2+ influx can induce the production of reactive oxygen species (ROS) and promote aging [9].
Ca2+-ATPases belong to the superfamily of P-type ATPases forming phosphor-aspartate enzyme intermediates during the reaction cycle [10]. They play an important role in maintaining Ca2+ homeostasis in plant cells by controlling Ca2+ efflux from the cytosol to organelles and/or to the apoplast [11]. Ca2+-ATPases have four highly conserved protein domains, including a cation transporter/ATPase N-terminal (PF00690), a cation transporter/ATPase C-terminal (PF00689), an E1–E2 ATPase (PF00122), and a haloacid dehalogenase-like hydrolase (PF00702) [12]. According to their protein sequences and protein locations in plants, plant Ca2+-ATPases can be grouped into two phylogenetic subgroups, P-type IIA ER-type Ca2+-ATPases (ECAs) and P-type IIB autoinhibited Ca2+-ATPases (ACAs) [13]. ACAs are characterized by an N-terminally situated calmodulin (CaM)-binding domain (CMBD). Autoinhibition can be relieved by Ca2+/CaM-binding, thus allowing ACAs to be directly controlled by free Ca2+ levels and giving rapid feedback regulation of Ca2+ signals [14,15].
To date, ACA gene family members have been identified in multiple species. Ten AtACA genes were first identified in Arabidopsis thaliana [16]. Subsequently, ACA gene families have been reported in Solanum lycopersicum [17], Nicotiana tabacum [17], Capsicum annuum [17], Solanum tuberosum [17], Brassica rapa [18], Ricinus communis [19], Zea mays [20], and Triticum aestivum [21]. ACA genes have been proven to function in the development and growth of plants [11]. For instance, AtACA1 regulates stomatal aperture, cytosolic distribution of chloroplasts in response to light, and root gravitropic curvature [22]. AtACA2, AtACA7, AtACA9, and AtACA13 contribute to pollen germination and fitness [23,24,25]. AtACA10 is involved in vegetative growth and inflorescence structure [26]. Further, ACA genes also play important roles in biotic stress responses in plants [27]. AtACA2 and AtACA4 can alleviate hypersensitivity to salt [28]. GmSCA1 is highly and rapidly induced by salt [29]. GsACA1 positively regulates plant tolerance to both salt and carbonate-alkaline stresses in Glycine soja [30]. In moss (Physcomitrella patens), the knockout of the ACA gene PCA1 results in an increased sensitivity to salt stress [31]. TaACAs are suggested to have diverse functions in various biological processes and stress responses [21]. Recently, more specific evidence suggested that ACA genes function in defense responses in plants. GmSCA1 is induced by a fungal elicitor [29]. AtACA11 is a genetic suppressor of the programmed cell death pathway in Arabidopsis [32]. AtACA4, AtACA8, and AtACA11 mediate the pathogen effector flg22-dependent Ca2+ signaling [33,34]. AtACA8 interacts with FLAGELLIN SENSITIVE2 (FLS2), the receptor for flg22, along with AtACA10 and AtACA12 [35].
In rice, OsACA gene family members have been identified by bioinformatics methods [36,37]. OsACA1 and OsACA8 are reported to be early-glutamate-responsive in rice roots [38]. OsACA1 is positively regulated by OsTF1L, a rice HD-Zip transcription factor which promotes lignin biosynthesis and stomatal closure that improves drought tolerance [39]. OsACA6 can efficiently modulate the ROS metabolism and proline biosynthesis in response to the abiotic stresses of salt, drought, and heat in tobacco [40,41]. OsACA7 is downregulated upon the supplementation of Silicon (Si) under abiotic and biotic stresses [42,43]. Currently, the knowledge of OsACAs is fragmented, and is mainly focused on stress responses and the regulation of plant growth and development; solid evidence on the role of OsACAs in biotic responses is still limited.
Here, we focused on an especially distinctive clade of OsACA genes, based on a phylogenetic tree between Arabidopsis (A. thaliana), Tomato (S. lycopersicum), maize (Z. mays), and rice (Oryza sativa), which was reported to be related to defense response. We constructed two knockout lines of the OsACA9 gene by CRISPR/Cas9 and characterized the phenotype and resistance spectrum of the knockout mutant. We carried out transcriptome analysis by RNA-Seq to investigate the molecular mechanism of defense response mediated by Ca2+-ATPase OsACA9 in rice. Our results suggest a pleiotropic effect of OsACA9 on disease resistance and leaf senescence and provide new evidence for the role of Ca2+-ATPases in plant–pathogen interactions.

2. Results

2.1. Evolution Analysis and Expression Profile Suggest the Potential Role of OsACA9 in Various Biological Processes

Previously, eleven OsACA genes had been initially reported from rice [36,37]. However, one of the family members (OsACA10) was discarded due to the lack of the cation transporter/ATPase N-terminal domain and the E1–E2 ATPase domain [37]. To analyze the evolution of OsACA genes, a phylogenetic tree of ACA family proteins from rice (10 OsACAs), Arabidopsis (10 AtACAs), tomato (8 SlACAs), and maize (16 ZmACAs) was constructed and classified into four clades (Figure 1A). Interestingly, OsACA9 and OsACA11 were grouped with AtACA8 and AtACA10 in the same branch (Clade IV), which is especially distinctive for its role in defense response [35,44]. Moreover, the gene structure of OsACA9 and OsACA11 were similar, with multiple introns and exons, which was obviously different from other OsACA family members (Figure 1B). These results indicated that OsACA9 and OsACA11 were likely to be paralogous in evolution and might be redundant in biological function. Therefore, the expression level of OsACA9 and OsACA11 in different tissues of rice plants were analyzed by reverse transcription-quantitative real-time PCR (RT-qPCR). The results showed that OsACA9 is highly expressed in various tissues in rice, especially in the leaf and stem, while OsACA11 was not or lowly expressed in all tissues tested (Figure 1C). These results suggest that OsACA9, as the homologue to AtACA8 and AtACA10, might biologically play a more important role in defense response as well as the growth and development in rice.
To investigate the function of OsACA9 in rice, two sgRNAs (sgRNA-1, sgRNA-2) were designed (Figure 2A) and constructed into a CRISPR-Cas9 vector, respectively. The result vectors were genetically transformed into a wild-type rice cultivar Zhonghua 11 (ZH11). Two different transgenic lines, aca9#1 and aca9#2, were identified, in which a thymine (T) base was inserted by the sgRNA-1, and four bases/one base were deleted by the sgRNA-2, respectively (Figure 2B). The insertion and deletion both resulted in reading frameshift mutations (knockout) and the early termination of translation of the OsACA9 gene, and the two mutants of OsACA9 were subsequently used for phenotype identification.

2.2. Knockout of OsACA9 Broadly Enhanced the Bacterial Disease Resistance of Rice

To identify the disease resistance of the Osaca9 mutants, different strains of bacterial leaf blight (BB) and bacterial leaf streak (BLS) were used for analysis. Six Xanthomonas oryzae pv. oryzae (Xoo) strains of BB (PXO86, PXO71, POX99A, Zhe173, C2 and C5) and two Xanthomonas oryzae pv. oryzicola (Xoc) strains of BLS (BLS256 and RS105) were inoculated in the leaves of the Osaca9 mutants and the wild-type control ZH11 (susceptible cultivar), respectively. The results showed that the lesion lengths of all the Xoo and Xoc strains formed in the leaves of the Osaca9 mutants were significantly decreased compared to the wild-type control (Figure 3A–C). This result suggested that the knockout of OsACA9 could broadly enhance the resistance of rice plants to bacterial disease. Moreover, it was found that the major ROS, hydrogen peroxide (H2O2), was excessively accumulated in the leaves of the Osaca9 mutants. The histological staining using 3,3′-diaminobenzidine (DAB) showed abundant brown spots in the leaves of the Osaca9 mutants, especially at leaf margins (Figure 3D), and the content of H2O2 was higher than that of the wild-type control (Figure 3E). The result suggested that the mutation of OsACA9 caused accumulation of ROS in leaves, which may contribute to its enhanced resistance to pathogens.

2.3. OsACA9 Regulated the Leaf Senescence and Development of Rice

Based on current knowledge, ROS has been proved a driving force for the accelerated senescence of leaves [45,46]. At the filling stage, the leaves of Osaca9 mutants appeared yellowish, and the plant height significantly decreased (Figure 4A,B and Figure 5F). The net photosynthetic rate (Pn) of flag leaves in Osaca9 mutants was significantly lower than that of the wild-type control (Figure 4C). The chlorophyll a (Chla), chlorophyll b (Chlb), and carotenoid (Car) contents of the flag leaves, second top leaves, and third top leaves from a single tiller of Osaca9 mutants were significantly reduced (Figure 4D). All these results indicated that the Osaca9 mutants exhibited an early-leaf-senescence phenotype and reduced photosynthesis efficiency in leaves.
Additionally, the agronomic traits of the Osaca9 mutants and the wild-type control were assayed at the mature stage. For the Osaca9 mutants, the flag leaf length (Figure 5D) and width (Figure 5E), the panicle length (Figure 5A,G), the grain number per panicle (Figure 5B,H), and the grain thickness (Figure 5M) were significantly decreased; however, the effective tiller number per plant (Figure 5I) and the grain length (Figure 5C,K) were significantly increased. There was no significant difference in seed setting rate (Figure 5J), the grain width (Figure 5L), and 1000-grain weight (Figure 5N) between the Osaca9 mutants and the wild-type control. As a result, the yield per plant of the Osaca9 mutants was decreased by about 15% of that of the wild-type control (Figure 5O).

2.4. Transcriptomic Analysis Suggested OsACA9 Involved in MAPK Cascades in PAMP-Triggered Defense Response

To investigate the potential molecular mechanism of the OsACA9 gene involved in regulating disease resistance and early leaf senescence, total RNAs from wild-type control and the Osaca9 knockout lines were respectively extracted and used for the deep sequencing by RNA-Seq technology. A total of 291 differentially expressed genes (DEGs; Log2 radio ≥1 or ≤−1; p < 0.05; Table S1) were identified, among which 234 genes were upregulated and 57 genes were downregulated. GO (gene ontology) enrichment analysis showed that the DEGs were significantly enriched in molecular functions related to binding to organic compounds, ATP and ions, and kinase activity. Further, DEGs were also enriched in biological processes that respond to external abiotic stimuli (Supplementary Table S2).
Moreover, the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis of the DEGs revealed that these genes were mainly involved in disease resistance pathways including MAPK signaling and plant–pathogen interactions, both of whose function in plant immunity have been extensively studied and supported by solid evidence [47]. Meanwhile, multiple pathways related to basal metabolism and immune metabolites were also significantly enriched, such as the alpha-Linolenic acid metabolism pathway; the glycine, serine and threonine metabolism pathway; and the biosynthesis of terpenoids and fatty acid as well (Figure 6A). Among those DEGs, five receptor-like cytoplasmic kinases (RLCKs) genes and two cell wall-associated kinases (WAKs) genes were upregulated, which are central players in basal immunity [47]. Moreover, several functionally known genes related to rice disease resistance were also upregulated, including rice sheath blight (RSB), disease resistance gene OsRSR (disease resistance protein RPM1), MDPK (malectin domain protein kinase), and rice blast disease resistance associated gene OsRBBI2 (Bowman–Birk trypsin inhibitor). Interestingly, OsSULTR3;6 encoding a sulphate transporter in rice, which is considered as the key susceptibility gene to BLS [48], was downregulated (Supplementary Table S3). At the same time, the amounts of genes involved in the synthesis of disease resistant metabolites were also significantly upregulated, such as the 4-coumarate-CoA ligase OsACS6, the aldehyde oxidase OsAO2, the 1-aminocyclopropane-1-carboxylate oxidases Os2ODD8 and Os2ODD23, and the beta-sesquiphellandrene synthases OsTPS30 and OsTPS31. In addition, several transcription factors (TFs) involved in disease resistance were also upregulated (Figure 6B, Table 1, Supplementary Table S4). Gene expression assessed by RT-qPCR in the wild type and two lines of Osaca9 was consistent with that measured by transcriptome sequencing (Figure 6C). Thus, OsACA9 might enhance disease resistance by activating the basal immune pathways and antibacterial metabolite synthesis pathways simultaneously.

3. Discussion

The ACA gene family is identified to be involved in multiple biological process, especially in stress responses. Based on the expression profile analysis, the expression level of ACA genes varies under abiotic stresses of salt, drought, and heat in multiple species [21,27,28,40,41]. However, recently, mounting evidence shows that ACAs have an additional role in defense response in Arabidopsis [34,44], but none of the studies have reported ACA’s role in disease resistance in rice yet. According to the results of phylogenetic analysis between O. sativa, A. thaliana, S. lycopersicum, and Z. mays (Figure 1A), we focused on an especially distinctive clade (Clade IV) related to biotic response including AtACA8 and AtACA10 [35,44,68], while other clades are mainly involved in the regulation of plant growth, development, and abiotic stress responses [11,17]. We constructed two knockout lines of the OsACA9 gene by CRISPR/Cas9 and characterized the phenotype and resistance spectrum of the knockout mutants. We carried out transcriptome analysis by RNA-Seq to investigate the molecular mechanism of defense response mediated by Ca2+-ATPase OsACA9 in rice. Our results suggest a pleiotropic effect of OsACA9 on disease resistance and leaf senescence and provide new evidence for the role of Ca2+-ATPases in plant–pathogen signaling.
Based on current knowledge, AtACA8 and AtACA10 were reported to be involved in PAMP-mediated immune response and interact with the receptor kinase FLS2 [44,68]. Naturally, we expected Osaca9 to show the similar sensitive phenotype to pathogens as aca8 and aca10 mutants. However, the resistance evaluation results turned out to be absolutely the opposite, showing that two knockout lines of Osaca9 were more resistant to multiple strains of Xoo including PXO86, PXO71, POX99A, Zhe173, C2, and C5, and Xoc strains including BLS256 and RS105 (Figure 3A–C). The transcriptome analysis showed that DEGs were enriched in pathways including MAPK signaling and plant–pathogen interactions (Figure 6A). Many of the DEGs encoded disease resistance proteins, growth/development-related proteins, signaling components, and transcription factors, as well as proteins involved in protein phosphorylation (Figure 6B; Table 1). Among them were multiple genes encoding RLCKs, which have emerged as a major class of signaling proteins that regulate plant cellular activities in response to biotic stresses and endogenous extracellular signaling molecules. By associating with immune receptor kinases (RKs), RLCKs regulate multiple downstream signaling nodes to orchestrate a complex array of defense responses against microbial pathogens [47]. Interestingly, the receptor for flg22, FLS2, was upregulated, while the brassinosteroid insensitive 1–associated kinase 1 (BAK1), an important factor in BRI1-mediated BR signaling [69,70], was downregulated. These results suggest that OsACA9 may have a similar signal transduction mechanism to AtACA8 and AtACA10 by affecting RKs and RLCKs in activating immune responses. Moreover, except for BB and BLS, multiple genes positively regulating rice resistance to sheath blight (SHB) and rice blast including MDPK, OsWRKYs, OsbHLH057, OsRSR, and OsRBBI2 were upregulated (Table 1), which suggests the potential function of OsACA9 in broad-spectrum resistance to pathogens.
Meanwhile, the accumulation of ROS, an essential role in signaling plant immunity, was detected in the Osaca9 mutants (Figure 3D,E). Interestingly, among all DEGs, OsPIP2;7 encoding an aquaporin and OsPRX7 encoding a peroxidase precursor, which function in ROS scavenging [35,40,71], were both upregulated (Figure 6B, Table 1). These results suggested an enhanced ROS scavenging capability of Osaca9, which helps to prevent the excessive accumulation of ROS from causing damage to plant cells while improving disease resistance. Previous studies have also reported some of the ACAs’ roles in ROS production and scavenging [35,40]. OsACA6 was reported to efficiently modulate ROS metabolism and proline biosynthesis in response to the abiotic stresses of salt, drought, cold, and heat in tobacco [40,41]. Thus, we speculated that OsACA9 might also have a role in stress response, which must be investigated through further studies. However, the excessive accumulation of ROS has also proved harmful to organisms as it triggers protein oxidation and enzyme inactivation and leads to reduced photosynthesis efficiency and impaired chloroplast development [45,46].
In addition, Osaca9 mutants displayed an early-leaf-senescence phenotype, especially during the filling stage (Figure 3D and Figure 4A–C). We speculated that the early-leaf-senescence phenotype was partly affected by ROS accumulation, as an expense of its enhanced resistance. Further, changes in endogenous phytohormone levels influence signaling networks that function in senescence processes. Plant hormones such as ethylene, jasmonic acid (JA), salicylic acid (SA), abscisic acid (ABA), and brassinosteroid (BR) promote leaf senescence and are extensively involved in responses to various abiotic and biotic stresses [72]. According to the transcription analysis, the upregulation of EIN3 (Figure 6C), a key transcription factor in the ethylene response pathway which upregulates the transcription level of ORE1 by inhibiting the transcription of miR164 [73], and the downregulation of BAK1 (Figure 6C), an important factor in BRI1-mediated BR signaling, may have contributed to the early senescence of Osaca9. Further, master TFs-mediated transcriptional regulation plays a crucial role in the regulation of leaf senescence [74]. Members of WRKY TFs and the basic helix-loop-helix (bHLH) family TFs, which coordinate with endogenous hormones to finely regulate the leaf senescence process, were upregulated in Osaca9 (Table 1, Supplementary Table S4). Moreover, RKs and MAPKs, regulating leaf senescence by affecting the phosphorylation status of target proteins [75], were also detected in DEGs (Figure 6B,C, Table 1). These results finely explain the early-leaf-senescence phenotype of Osaca9.
In addition, the increased resistance might also cause dwarfism, through shortened leaves and panicles [76]. According to the statistical analysis of agronomic traits of the Osaca9 mutant, unsurprisingly, the length and width of the flag leaves of knockout lines decreased significantly. The length of the main panicle became shorter, leading to a decreased grain number (Figure 5A,B,D,E,G,H). However, there was no significant difference in seed setting rate and 1000-grain weight between Osaca9 and wild-type (Figure 5J,N), suggesting that Osaca9 did not affect grain filling and seed production. Interestingly, the reduced yield per plant of Osaca9 was partly compensated by increased available tiller number (Figure 5I), suggesting that OsACA9 might have an additional role in regulating the rhizomes of rice.
In plants, Ca2+ is an important nutrient and cellular secondary signaling molecule that is essential for responses to biotic and abiotic stimuli as well as growth and development [11,77]. Owing to its cytotoxicity, cytosolic Ca2+ levels must be maintained at low (10−8 to 10−7 M) levels in living cells [78]. Ca2+ homeostasis is maintained by an array of Ca2+ transport elements during unfavorable situations including Ca2+-ATPases, cyclic nucleotide-gated channels (CNGCs), and glutamate receptor homologs (GLRs), etc. [79]. ACAs are considered one major type of Ca2+-ATPases that mediate active Ca2+ transport out of the cytosol. In OsACA9 knockout mutants, the efflux of Ca2+ might be hindered, resulting in the continuous accumulation of free cytosolic Ca2+. The cytosol Ca2+ burst has been proven to promote the phosphorylation of RBOH and trigger ROS production [80], as well as the upstream activation of many responses, including the activation of MAPK kinases, as well as the induction of pathogenesis-related gene expression [81], leading to enhanced resistance in Osaca9 mutants. It can be implied that OsACA9 might balance plant growth and development with immunity response in such model (Figure 6D). In addition, a deeper understanding of the growth–defense trade-offs mediated by OsACA9 will be a powerful tool to reveal plant–bacterial interaction and can be useful for optimizing crop breeding [79]. We are currently conducting a detailed study to determine how mutations in OsACA9 activate pivotal immune pathways, especially the MAPK signaling pathway, and further explore the molecular mechanisms of how Ca2+-ATPase orchestrates external stimulus signaling through the Ca2+ signature.

4. Materials and Methods

4.1. Plant Materials

A Japonica rice cultivar named ZH11 was used for knockout of OsACA9 gene by CRISPR/Cas9. The different knockout lines of OsACA9 gene, aca9#1 and aca9#2, and the wild-type control ZH11 were planted in the field of an experimental farm in summer in South China. For the plot yield test, 200 individuals were planted in an 8 m2 (2 × 4) plot with a planting density of 20 cm × 20 cm. The agronomic traits of rice plants were measured on individual plants with at least 10 replicates at the mature stage.

4.2. Phylogenetic Analysis of OsACA Gene Family

The OsACA gene family has been reported in rice by Singh et al. (2014). For phylogenetic analysis, protein sequences of ACA family members of rice (O. sativa), Arabidopsis (A. thaliana) [37], tomato (S. lycopersicum) [17], and maize (Z. mays) [20] were used, and a phylogenetic tree was constructed by the neighbor joining method with 1000 bootstrap replicates by MEGA7.0 [82]. OsACA genes’ structure and clutter were performed using TBtools [83].

4.3. Construction of OsACA9 Knockout Rice Plants by CRISPR/Cas9

The coding sequence (CDS) of the OsACA9 gene (Os02g0176700) was analyzed by the program CRISPR Design: http://crispr.mit.edu (accessed on 28 June 2022). Two sgRNAs, sgRNA-1 (5′-GCTTCATATGAAATTTGCGG-3′) and sgRNA-2 (5′-ATACAGTGACGCTGAAAGGG-3′), were selected for OsACA9 gene editing and constructed into a CRISPR/Cas9 vector using the BGK03 kit (Biogle, Changzhou, China) according to the manufacturer’s instructions. The vector was transformed into the wild-type rice cultivar ZH11 by the Agrobacterium-mediated method [84]. The genomic DNA from the leaves of transgenic plants was extracted using the CTAB method and used for PCR amplification of OsACA9 gene editing target sites. The primer pairs used for PCR were sgRNA-1 F: 5′-CATGCAGGGGAAGCGTTT-3′ and sgRNA-1 R: 5′-CCTTTTGCCAATCAACCA-3′, and sgRNA-2 F: 5′-GCTCCAGTGAACATAAGATA-3′ and sgRNA-2 R: 5′-CTCGAAAGTATGTACTCCAGAT-3′.

4.4. Assay of Chlorophyll Content and Photosynthesis

At the filling stage, 0.1 g of fresh leaves of rice plants were sampled, and the chlorophyll pigment was extracted with equal volumes of 80% acetone. All samples were kept in a dark place at room temperature for 24 h with frequent shaking until all leaves turned white. The absorbance of the supernatant was measured by the spectrophotometer NanoDrop2000C (Thermo Fisher Scientific, Wilmington, DE, USA) at 665 nm, 649 nm, and 470 nm. Each sample was tested with at least three biological repeats. The concentrations of chlorophyll a, chlorophyll b and carotenoids were calculated according to Lichtenthaler’s method [85]. For the assay of photosynthesis, the net photosynthetic rate of the flag leaves of rice plants was measured using the LI-6400 Portable Photosynthesis System (LI-COR, Lincoln, NE, USA) following the manufacturer’s instructions.
Chlorophyll a concentration = 13.95A665 − 6.88A649
Chlorophyll b concentration = 24.96A649 − 7.32A665
Carotenoid concentration = (1000A470 − 2.05Chla − 114Chlb)/245

4.5. Disease Resistance Identification of Rice Plants

Six Xoo strains (PXO86, PXO71, POX99A, Zhe173, C2 and C5) and two Xoc strains (BLS256 and RS105) were used for disease resistance identification. Strains were cultured on a potato dextrose agar (PDA) medium that contained 20 g sucrose, 5 g peptone, 0.5 g Ca(NO3)2, 0.43 g Na2HPO4, and 0.05 g FeSO4 per liter and grew at 28 °C for 3 days. The bacterial colony was suspended in sterile distilled water at an optical density of OD600 = 1.0 and immediately used for inoculation. At the tillering stage, the leaf tip clipping method [86] was used for inoculation of Xoo, and the non-needle syringe injection method [87] was used for inoculation of Xoc. Lesions on the inoculated leaves were measured for the evaluation of resistance 2 weeks after inoculation.

4.6. DAB Staining and H2O2 Assay

Leaves of the rice plant were stained by DAB according to Thordal-Christensen’s method [88]. In brief, leaves were placed in 1 mg/mL DAB–HCI and 10 mM ascorbic acid, pH 3.8, incubated in the growth chamber for 8 h, and then cleared in boiling ethanol (96%) for 10 min. H2O2 is visualized as a reddish brown coloration. Peroxide (H2O2) Content Assay kits (Solarbio, Beijing, China) were used to measure endogenous H2O2 content. Each measurement was performed in at least three biological replicates.

4.7. RNA-Seq and Data Analysis

At the tillering stage, leaves of rice plants were sampled and sent to the Beijing Genomics Institute (Shenzhen, China) for the construction of the cDNA library and next-generation sequencing (NGS). Each sample had three biological replicates which were sequenced by the BGISEQ-500 platform. Differentially expressed genes (DEGs; log2 ratio ≥1 or ≤−1) were detected with a significance p < 0.05 between samples. Gene functional enrichments were analyzed using GO: http://www.geneontology.org/ (accessed on 11 August 2023) and KEGG: http://www.genome.jp/kegg/ (accessed on 11 August 2023).

4.8. RNA Extraction and RT-qPCR

Total RNA of rice was extracted using RNeasy Plant Mini Kit (QIAGEN, Valencia, CA, USA) and RNase-Free DNase Set (QIAGEN, Valencia, CA, USA) following the manufacturer’s instruction. Synthesis of first-strand cDNAs from RNA was carried out by the M-MLV Reverse Transcriptase (Promega, Beijing, China) according to the manufacturer’s instructions. qPCR was conducted using TB Green™ Premix Ex Taq™ (TaKaRa, Beijing, China) with a final volume of 10 μL per reaction. Each reaction mixture consisted of 5.0 μL TB Green Premix Ex Taq (TaKaRa, Beijing, China), 1.0 μL cDNA template, 0.5 μL of each primer, and 3.3 μL RNase-Free double distilled water. The qPCR was performed on the Step OnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) using the following program: 95 °C for 5 min, followed by 40 cycles at 95 °C for 5 s and 60 °C for 30 s. Relative expression level of the interest gene was calculated by the 2−ΔΔCT method using the housekeeping genes OsEF1 and OsActin as the standardization controls [89]. Each measurement was performed in at least three biological replicates. The primers are listed in Table 2.

Supplementary Materials

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

Author Contributions

Designed and supervised the experiments, X.C. and B.M.; performed the experiments, X.W., Z.W., Y.L., J.H., Z.H., J.L. and X.F.; contributed to data analysis, X.W., Z.G. and P.L.; wrote the manuscript, X.W. and P.L.; provided guidance on the whole manuscript, X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (32071987, 32101697, 32272096), the Natural Science Foundation of Zhejiang Province (LQ22C130005, LZ23C130004), and the College Students’ Science and Technology Innovation Program of Zhejiang Province (2023R404038).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data of transcriptomic sequencing can be downloaded in the NCBI Sequence Read Archive under accession number PRJNA1056136.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tian, W.; Wang, C.; Gao, Q.; Li, L.; Luan, S. Calcium Spikes, Waves and Oscillations in Plant Development and Biotic Interactions. Nat. Plants 2020, 6, 750–759. [Google Scholar] [CrossRef]
  2. Köster, P.; DeFalco, T.A.; Zipfel, C. Ca2+ Signals in Plant Immunity. EMBO J. 2022, 41, e110741. [Google Scholar] [CrossRef]
  3. Upadhyay, S.K. Calcium Channels, OST1 and Stomatal Defence: Current Status and Beyond. Cells 2022, 12, 127. [Google Scholar] [CrossRef]
  4. Dong, X.; Gifford, M.L.; Su, C. Ca2+ Signatures in Symbiosis: Another Level of Dynamism for This Key Messenger. J. Exp. Bot. 2024, 75, 508–510. [Google Scholar] [CrossRef] [PubMed]
  5. Knight, M.R.; Campbell, A.K.; Smith, S.M.; Trewavas, A.J. Transgenic Plant Aequorin Reports the Effects of Touch and Cold-Shock and Elicitors on Cytoplasmic Calcium. Nature 1991, 352, 524–526. [Google Scholar] [CrossRef] [PubMed]
  6. Knight, H.; Trewavas, A.J.; Knight, M.R. Calcium Signalling in Arabidopsis thaliana Responding to Drought and Salinity. Plant J. 1997, 12, 1067–1078. [Google Scholar] [CrossRef] [PubMed]
  7. Upadhyay, S.K. CPK12 and Ca2+-Mediated Hypoxia Signaling. Plant Signal Behav. 2023, 18, 2273593. [Google Scholar] [CrossRef] [PubMed]
  8. Lim, P.O.; Kim, H.J.; Nam, H.G. Leaf Senescence. Annu. Rev. Plant Biol. 2007, 58, 115–136. [Google Scholar] [CrossRef] [PubMed]
  9. Knight, H.; Brandt, S.; Knight, M.R. A History of Stress Alters Drought Calcium Signalling Pathways in Arabidopsis. Plant J. 1998, 16, 681–687. [Google Scholar] [CrossRef] [PubMed]
  10. Møller, J.V.; Juul, B.; le Maire, M. Structural Organization, Ion Transport, and Energy Transduction of P-Type ATPases. Biochim. Biophys. Acta (BBA)-Rev. Biomembr. 1996, 1286, 1–51. [Google Scholar] [CrossRef]
  11. García Bossi, J.; Kumar, K.; Barberini, M.L.; Domínguez, G.D.; Rondón Guerrero, Y.D.C.; Marino-Buslje, C.; Obertello, M.; Muschietti, J.P.; Estevez, J.M. The Role of P-Type IIA and P-Type IIB Ca2+-ATPases in Plant Development and Growth. J. Exp. Bot. 2020, 71, 1239–1248. [Google Scholar] [CrossRef]
  12. Wdowikowska, A.; Kłobus, G. Plant P-Type ATPases. Postepy Biochem. 2011, 57, 85–91. [Google Scholar]
  13. Pedersen, C.N.S.; Axelsen, K.B.; Harper, J.F.; Palmgren, M.G. Evolution of Plant P-Type ATPases. Front. Plant Sci. 2012, 3, 31. [Google Scholar] [CrossRef] [PubMed]
  14. Curran, A.C.; Hwang, I.; Corbin, J.; Martinez, S.; Rayle, D.; Sze, H.; Harper, J.F. Autoinhibition of a Calmodulin-Dependent Calcium Pump Involves a Structure in the Stalk That Connects the Transmembrane Domain to the ATPase Catalytic Domain. J. Biol. Chem. 2000, 275, 30301–30308. [Google Scholar] [CrossRef] [PubMed]
  15. Hwang, I.; Harper, J.F.; Liang, F.; Sze, H. Calmodulin Activation of an Endoplasmic Reticulum-Located Calcium Pump Involves an Interaction with the N-Terminal Autoinhibitory Domain. Plant Physiol. 2000, 122, 157–168. [Google Scholar] [CrossRef] [PubMed]
  16. Axelsen, K.B.; Palmgren, M.G. Inventory of the Superfamily of P-Type Ion Pumps in Arabidopsis. Plant Physiol. 2001, 126, 696–706. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, Y.X.; Shu, Y.; Chen, X.L.; Wang, A.X. Identification and Analysis of Ca2+-ATPase Gene Family in Solanaceae. Mol. Plant Breed. 2021, 4268–4277. [Google Scholar]
  18. Jiang, Q.; Wu, X.; Zhang, X.; Ji, Z.; Cao, Y.; Duan, Q.; Huang, J. Genome-Wide Identification and Expression Analysis of AS2 Genes in Brassica Rapa Reveal Their Potential Roles in Abiotic Stress. Int. J. Mol. Sci. 2023, 24, 10534. [Google Scholar] [CrossRef] [PubMed]
  19. Li, Y.; Zhang, C.; Geng, L.; Chen, Y.; Zhang, L.; Xiang, D.; Liu, P. Identification and Expression Pattern Analysis of RcACA Gene Family in Castor under Abiotic Stresses. Chin. J. Biotechnol. 2023, 39, 2861–2873. [Google Scholar]
  20. Chen, C.; Zhang, M.; Zhang, M.; Yang, M.; Dai, S.; Meng, Q.; Lv, W.; Zhuang, K. Ethylene-Insensitive 3-LIKE 2 Regulates β-Carotene and Ascorbic Acid Accumulation in Tomatoes during Ripening. Plant Physiol. 2023, 192, 2067–2080. [Google Scholar] [CrossRef]
  21. Taneja, M.; Upadhyay, S.K. Molecular Characterization and Differential Expression Suggested Diverse Functions of P-Type II Ca2+ ATPases in Triticum aestivum L. BMC Genom. 2018, 19, 389. [Google Scholar] [CrossRef]
  22. Huang, L.; Berkelman, T.; Franklin, A.E.; Hoffman, N.E. Characterization of a Gene Encoding a Ca(2+)-ATPase-like Protein in the Plastid Envelope. Proc. Natl. Acad. Sci. USA 1993, 90, 10066–10070. [Google Scholar] [CrossRef]
  23. Rahmati Ishka, M.; Brown, E.; Rosenberg, A.; Romanowsky, S.; Davis, J.A.; Choi, W.-G.; Harper, J.F. Arabidopsis Ca2+-ATPases 1, 2, and 7 in the Endoplasmic Reticulum Contribute to Growth and Pollen Fitness. Plant Physiol. 2021, 185, 1966–1985. [Google Scholar] [CrossRef]
  24. Schiøtt, M.; Romanowsky, S.M.; Baekgaard, L.; Jakobsen, M.K.; Palmgren, M.G.; Harper, J.F. A Plant Plasma Membrane Ca2+ Pump Is Required for Normal Pollen Tube Growth and Fertilization. Proc. Natl. Acad. Sci. USA 2004, 101, 9502–9507. [Google Scholar] [CrossRef] [PubMed]
  25. Iwano, M.; Igarashi, M.; Tarutani, Y.; Kaothien-Nakayama, P.; Nakayama, H.; Moriyama, H.; Yakabe, R.; Entani, T.; Shimosato-Asano, H.; Ueki, M.; et al. A Pollen Coat-Inducible Autoinhibited Ca2+-ATPase Expressed in Stigmatic Papilla Cells Is Required for Compatible Pollination in the Brassicaceae. Plant Cell 2014, 26, 636–649. [Google Scholar] [CrossRef] [PubMed]
  26. George, L.; Romanowsky, S.M.; Harper, J.F.; Sharrock, R.A. The ACA10 Ca2+-ATPase Regulates Adult Vegetative Development and Inflorescence Architecture in Arabidopsis. Plant Physiol. 2008, 146, 716–728. [Google Scholar] [CrossRef] [PubMed]
  27. Kamrul Huda, K.M.; Yadav, S.; Akhter Banu, M.S.; Trivedi, D.K.; Tuteja, N. Genome-Wide Analysis of Plant-Type II Ca(2+)ATPases Gene Family from Rice and Arabidopsis: Potential Role in Abiotic Stresses. Plant Physiol. Biochem. 2013, 65, 32–47. [Google Scholar] [CrossRef] [PubMed]
  28. Anil, V.S.; Rajkumar, P.; Kumar, P.; Mathew, M.K. A Plant Ca2+ Pump, ACA2, Relieves Salt Hypersensitivity in Yeast. Modulation of Cytosolic Calcium Signature and Activation of Adaptive Na+ Homeostasis. J. Biol. Chem. 2008, 283, 3497–3506. [Google Scholar] [CrossRef] [PubMed]
  29. Chung, W.S.; Lee, S.H.; Kim, J.C.; Heo, W.D.; Kim, M.C.; Park, C.Y.; Park, H.C.; Lim, C.O.; Kim, W.B.; Harper, J.F.; et al. Identification of a Calmodulin-Regulated Soybean Ca(2+)-ATPase (SCA1) That Is Located in the Plasma Membrane. Plant Cell 2000, 12, 1393–1407. [Google Scholar] [CrossRef] [PubMed]
  30. Sun, M.; Jia, B.; Cui, N.; Wen, Y.; Duanmu, H.; Yu, Q.; Xiao, J.; Sun, X.; Zhu, Y. Functional Characterization of a Glycine Soja Ca(2+)ATPase in Salt-Alkaline Stress Responses. Plant Mol. Biol. 2016, 90, 419–434. [Google Scholar] [CrossRef]
  31. Qudeimat, E.; Faltusz, A.M.C.; Wheeler, G.; Lang, D.; Holtorf, H.; Brownlee, C.; Reski, R.; Frank, W. A PIIB-Type Ca2+-ATPase Is Essential for Stress Adaptation in Physcomitrella Patens. Proc. Natl. Acad. Sci. USA 2008, 105, 19555–19560. [Google Scholar] [CrossRef]
  32. Boursiac, Y.; Lee, S.M.; Romanowsky, S.; Blank, R.; Sladek, C.; Chung, W.S.; Harper, J.F. Disruption of the Vacuolar Calcium-ATPases in Arabidopsis Results in the Activation of a Salicylic Acid-Dependent Programmed Cell Death Pathway. Plant Physiol. 2010, 154, 1158–1171. [Google Scholar] [CrossRef] [PubMed]
  33. Benschop, J.J.; Mohammed, S.; O’Flaherty, M.; Heck, A.J.R.; Slijper, M.; Menke, F.L.H. Quantitative Phosphoproteomics of Early Elicitor Signaling in Arabidopsis. Mol. Cell Proteom. 2007, 6, 1198–1214. [Google Scholar] [CrossRef]
  34. Hilleary, R.; Paez-Valencia, J.; Vens, C.S.; Toyota, M.; Palmgren, M.; Gilroy, S. Tonoplast-Localized Ca2+ Pumps Regulate Ca2+ Signals during Pattern-Triggered Immunity in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2020, 117, 18849–18857. [Google Scholar] [CrossRef] [PubMed]
  35. Frei dit Frey, N.; Mbengue, M.; Kwaaitaal, M.; Nitsch, L.; Altenbach, D.; Häweker, H.; Lozano-Duran, R.; Njo, M.F.; Beeckman, T.; Huettel, B.; et al. Plasma Membrane Calcium ATPases Are Important Components of Receptor-Mediated Signaling in Plant Immune Responses and Development. Plant Physiol. 2012, 159, 798–809. [Google Scholar] [CrossRef] [PubMed]
  36. Baxter, I.; Tchieu, J.; Sussman, M.R.; Boutry, M.; Palmgren, M.G.; Gribskov, M.; Harper, J.F.; Axelsen, K.B. Genomic Comparison of P-Type ATPase Ion Pumps in Arabidopsis and Rice. Plant Physiol. 2003, 132, 618–628. [Google Scholar] [CrossRef]
  37. Singh, A.; Kanwar, P.; Yadav, A.K.; Mishra, M.; Jha, S.K.; Baranwal, V.; Pandey, A.; Kapoor, S.; Tyagi, A.K.; Pandey, G.K. Genome-Wide Expressional and Functional Analysis of Calcium Transport Elements during Abiotic Stress and Development in Rice. FEBS J. 2014, 281, 894–915. [Google Scholar] [CrossRef]
  38. Kan, C.-C.; Chung, T.-Y.; Wu, H.-Y.; Juo, Y.-A.; Hsieh, M.-H. Exogenous Glutamate Rapidly Induces the Expression of Genes Involved in Metabolism and Defense Responses in Rice Roots. BMC Genom. 2017, 18, 186. [Google Scholar] [CrossRef]
  39. Bang, S.W.; Lee, D.-K.; Jung, H.; Chung, P.J.; Kim, Y.S.; Do Choi, Y.; Suh, J.-W.; Kim, J.-K. Overexpression of OsTF1L, a Rice HD-Zip Transcription Factor, Promotes Lignin Biosynthesis and Stomatal Closure That Improves Drought Tolerance. Plant Biotechnol. J. 2019, 17, 118–131. [Google Scholar] [CrossRef] [PubMed]
  40. Huda, K.M.K.; Banu, M.S.A.; Garg, B.; Tula, S.; Tuteja, R.; Tuteja, N. OsACA6, a P-Type IIB Ca2+ ATPase Promotes Salinity and Drought Stress Tolerance in Tobacco by ROS Scavenging and Enhancing the Expression of Stress-Responsive Genes. Plant J. 2013, 76, 997–1015. [Google Scholar] [CrossRef]
  41. Um, T.; Park, T.; Shim, J.S.; Kim, Y.S.; Lee, G.-S.; Choi, I.-Y.; Kim, J.-K.; Seo, J.S.; Park, S.C. Application of Upstream Open Reading Frames (UORFs) Editing for the Development of Stress-Tolerant Crops. Int. J. Mol. Sci. 2021, 22, 3743. [Google Scholar] [CrossRef]
  42. Manivannan, A.; Ahn, Y.-K. Silicon Regulates Potential Genes Involved in Major Physiological Processes in Plants to Combat Stress. Front. Plant Sci. 2017, 8, 1346. [Google Scholar] [CrossRef]
  43. Amirbakhtiar, N.; Ismaili, A.; Ghaffari, M.-R.; Mirdar Mansuri, R.; Sanjari, S.; Shobbar, Z.-S. Transcriptome Analysis of Bread Wheat Leaves in Response to Salt Stress. PLoS ONE 2021, 16, e0254189. [Google Scholar]
  44. Kadota, Y.; Liebrand, T.W.H.; Goto, Y.; Sklenar, J.; Derbyshire, P.; Menke, F.L.H.; Torres, M.-A.; Molina, A.; Zipfel, C.; Coaker, G.; et al. Quantitative Phosphoproteomic Analysis Reveals Common Regulatory Mechanisms between Effector- and PAMP-Triggered Immunity in Plants. New Phytol. 2019, 221, 2160–2175. [Google Scholar] [CrossRef]
  45. van Dongen, J.T.; Licausi, F. Oxygen Sensing and Signaling. Annu. Rev. Plant Biol. 2015, 66, 345–367. [Google Scholar] [CrossRef]
  46. Song, G.; Kwon, C.-T.; Kim, S.-H.; Shim, Y.; Lim, C.; Koh, H.-J.; An, G.; Kang, K.; Paek, N.-C. The Rice SPOTTED LEAF4 (SPL4) Encodes a Plant Spastin That Inhibits ROS Accumulation in Leaf Development and Functions in Leaf Senescence. Front. Plant Sci. 2019, 9, 1925. [Google Scholar] [CrossRef] [PubMed]
  47. Liang, X.; Zhou, J.-M. Receptor-Like Cytoplasmic Kinases: Central Players in Plant Receptor Kinase–Mediated Signaling. Annu. Rev. Plant Biol. 2018, 69, 267–299. [Google Scholar] [CrossRef] [PubMed]
  48. Cernadas, R.A.; Doyle, E.L.; Niño-Liu, D.O.; Wilkins, K.E.; Bancroft, T.; Wang, L.; Schmidt, C.L.; Caldo, R.; Yang, B.; White, F.F.; et al. Code-Assisted Discovery of TAL Effector Targets in Bacterial Leaf Streak of Rice Reveals Contrast with Bacterial Blight and a Novel Susceptibility Gene. PLoS Pathog. 2014, 10, e1003972. [Google Scholar] [CrossRef] [PubMed]
  49. Vij, S.; Giri, J.; Dansana, P.K.; Kapoor, S.; Tyagi, A.K. The Receptor-Like Cytoplasmic Kinase (OsRLCK) Gene Family in Rice: Organization, Phylogenetic Relationship, and Expression during Development and Stress. Mol. Plant 2008, 1, 732–750. [Google Scholar] [CrossRef] [PubMed]
  50. Wang, F.; Han, T.; Song, Q.; Ye, W.; Song, X.; Chu, J.; Li, J.; Chen, Z.J. The Rice Circadian Clock Regulates Tiller Growth and Panicle Development Through Strigolactone Signaling and Sugar Sensing. Plant Cell 2020, 32, 3124–3138. [Google Scholar] [CrossRef]
  51. Wang, A.; Shu, X.; Jing, X.; Jiao, C.; Chen, L.; Zhang, J.; Ma, L.; Jiang, Y.; Yamamoto, N.; Li, S.; et al. Identification of Rice (Oryza Sativa L.) Genes Involved in Sheath Blight Resistance via a Genome-wide Association Study. Plant Biotechnol. J. 2021, 19, 1553–1566. [Google Scholar] [CrossRef]
  52. Cui, Z.; Xue, C.; Mei, Q.; Xuan, Y. Malectin Domain Protein Kinase (MDPK) Promotes Rice Resistance to Sheath Blight via IDD12, IDD13, and IDD14. Int. J. Mol. Sci. 2022, 23, 8214. [Google Scholar] [CrossRef]
  53. Xu, X.; Xu, Z.; Li, Z.; Zakria, M.; Zou, L.; Chen, G. Increasing Resistance to Bacterial Leaf Streak in Rice by Editing the Promoter of Susceptibility Gene OsSULRT3;6. Plant Biotechnol. J. 2021, 19, 1101–1103. [Google Scholar] [CrossRef]
  54. Kumar, K.; Mosa, K.A.; Chhikara, S.; Musante, C.; White, J.C.; Dhankher, O.P. Two Rice Plasma Membrane Intrinsic Proteins, OsPIP2;4 and OsPIP2;7, Are Involved in Transport and Providing Tolerance to Boron Toxicity. Planta 2014, 239, 187–198. [Google Scholar] [CrossRef] [PubMed]
  55. Byeon, Y.; Back, K. Molecular Cloning of Melatonin 2-hydroxylase Responsible for 2-hydroxymelatonin Production in Rice (ORyza Sativa). J. Pineal Res. 2015, 58, 343–351. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, J.; Zhang, C.; Li, Y. Genome-Wide Identification and Expression Profiles of 13 Key Structural Gene Families Involved in the Biosynthesis of Rice Flavonoid Scaffolds. Genes 2022, 13, 410. [Google Scholar] [CrossRef]
  57. Yang, L.; Wang, C.C.; Guo, W.D.; Li, X.B.; Lu, M.; Yu, C.L. Differential Expression of Cell Wall Related Genes in the Elongation Zone of Rice Roots under Water Deficit. Russ. J. Plant Physiol. 2006, 53, 390–395. [Google Scholar] [CrossRef]
  58. Sun, Y.; Zhang, Y.-J.; Cao, G.-C.; Gu, S.-H.; Wu, K.-M.; Gao, X.-W.; Guo, Y.-Y. Rice Gene Expression Profiles Responding to Larval Feeding of the Striped Stem Borer at the 1st to 2nd Instar Stage. Insect Sci. 2011, 18, 273–281. [Google Scholar] [CrossRef]
  59. Sun, Y.; Huang, X.; Ning, Y.; Jing, W.; Bruce, T.J.A.; Qi, F.; Xu, Q.; Wu, K.; Zhang, Y.; Guo, Y. TPS46, a Rice Terpene Synthase Conferring Natural Resistance to Bird Cherry-Oat Aphid, Rhopalosiphum Padi (Linnaeus). Front. Plant Sci. 2017, 8, 110. [Google Scholar] [CrossRef] [PubMed]
  60. Zhang, L.; Xiang, Z.; Li, J.; Wang, S.; Chen, Y.; Liu, Y.; Mao, D.; Luan, S.; Chen, L. BHLH57 Confers Chilling Tolerance and Grain Yield Improvement in Rice. Plant Cell Environ. 2023, 46, 1402–1418. [Google Scholar] [CrossRef] [PubMed]
  61. Jang, S.; Cho, K.; Shibato, J.; Han, O.; Iwahashi, H.; Tamogami, S.; Majeed Zargar, S.; Kubo, A.; Masuo, Y.; Agrawal, G.K.; et al. Rice OsOPRs: Transcriptional Profiling Responses to Diverse Environmental Stimuli and Biochemical Analysis of OsOPR1. J. Plant Biol. 2009, 52, 229–243. [Google Scholar] [CrossRef]
  62. Kuroha, T.; Nagai, K.; Gamuyao, R.; Wang, D.R.; Furuta, T.; Nakamori, M.; Kitaoka, T.; Adachi, K.; Minami, A.; Mori, Y.; et al. Ethylene-Gibberellin Signaling Underlies Adaptation of Rice to Periodic Flooding. Science 2018, 361, 181–186. [Google Scholar] [CrossRef]
  63. Wang, H.; Hao, J.; Chen, X.; Hao, Z.; Wang, X.; Lou, Y.; Peng, Y.; Guo, Z. Overexpression of Rice WRKY89 Enhances Ultraviolet B Tolerance and Disease Resistance in Rice Plants. Plant Mol. Biol. 2007, 65, 799–815. [Google Scholar] [CrossRef] [PubMed]
  64. Takagi, K.; Tasaki, K.; Komori, H.; Katou, S. Hypersensitivity-Related Genes HSR201 and HSR203J Are Regulated by Calmodulin-Binding Protein 60-Type Transcription Factors and Required for Pathogen Signal-Induced Salicylic Acid Synthesis. Plant Cell Physiol. 2022, 63, 1008–1022. [Google Scholar] [CrossRef] [PubMed]
  65. Song, Y.; Niu, R.; Yu, H.; Guo, J.; Du, C.; Zhang, Z.; Wei, Y.; Li, J.; Zhang, S. OsSLA1Functions in Leaf Angle Regulation by Enhancing the Interaction between OsBRI1and OsBAK1in Rice. Plant J. 2022, 110, 1111–1127. [Google Scholar] [CrossRef] [PubMed]
  66. Shi, X.; Tian, Q.; Deng, P.; Zhang, W.; Jing, W. The Rice Aldehyde Oxidase OsAO3 Gene Regulates Plant Growth, Grain Yield, and Drought Tolerance by Participating in ABA Biosynthesis. Biochem. Biophys. Res. Commun. 2021, 548, 189–195. [Google Scholar] [CrossRef]
  67. Takai, R.; Isogai, A.; Takayama, S.; Che, F.-S. Analysis of Flagellin Perception Mediated by Flg22 Receptor OsFLS2 in Rice. Mol. Plant-Microbe Interact. 2008, 21, 1635–1642. [Google Scholar] [CrossRef]
  68. Yu, H.; Yan, J.; Du, X.; Hua, J. Overlapping and Differential Roles of Plasma Membrane Calcium ATPases in Arabidopsis Growth and Environmental Responses. J. Exp. Bot. 2018, 69, 2693–2703. [Google Scholar] [CrossRef] [PubMed]
  69. Schwessinger, B.; Roux, M.; Kadota, Y.; Ntoukakis, V.; Sklenar, J.; Jones, A.; Zipfel, C. Phosphorylation-Dependent Differential Regulation of Plant Growth, Cell Death, and Innate Immunity by the Regulatory Receptor-like Kinase BAK1. PLoS Genet. 2011, 7, e1002046. [Google Scholar] [CrossRef]
  70. Albrecht, C.; Boutrot, F.; Segonzac, C.; Schwessinger, B.; Gimenez-Ibanez, S.; Chinchilla, D.; Rathjen, J.P.; de Vries, S.C.; Zipfel, C. Brassinosteroids Inhibit Pathogen-Associated Molecular Pattern-Triggered Immune Signaling Independent of the Receptor Kinase BAK1. Proc. Natl. Acad. Sci. USA 2012, 109, 303–308. [Google Scholar] [CrossRef]
  71. Meng, X.; Zhang, S. MAPK Cascades in Plant Disease Resistance Signaling. Annu. Rev. Phytopathol. 2013, 51, 245–266. [Google Scholar] [CrossRef]
  72. Zhang, Y.-M.; Guo, P.; Xia, X.; Guo, H.; Li, Z. Multiple Layers of Regulation on Leaf Senescence: New Advances and Perspectives. Front. Plant Sci. 2021, 12, 788996. [Google Scholar] [CrossRef]
  73. Li, Z.; Peng, J.; Wen, X.; Guo, H. Ethylene-Insensitive3 Is a Senescence-Associated Gene That Accelerates Age-Dependent Leaf Senescence by Directly Repressing MiR164 Transcription in Arabidopsis. Plant Cell 2013, 25, 3311–3328. [Google Scholar] [CrossRef]
  74. Breeze, E.; Harrison, E.; McHattie, S.; Hughes, L.; Hickman, R.; Hill, C.; Kiddle, S.; Kim, Y.; Penfold, C.A.; Jenkins, D.; et al. High-Resolution Temporal Profiling of Transcripts during Arabidopsis Leaf Senescence Reveals a Distinct Chronology of Processes and Regulation. Plant Cell 2011, 23, 873–894. [Google Scholar] [CrossRef]
  75. Ahmad, S.; Guo, Y. Signal Transduction in Leaf Senescence: Progress and Perspective. Plants 2019, 8, 405. [Google Scholar] [CrossRef]
  76. Liu, M.; Shi, Z.; Zhang, X.; Wang, M.; Zhang, L.; Zheng, K.; Liu, J.; Hu, X.; Di, C.; Qian, Q.; et al. Inducible Overexpression of Ideal Plant Architecture1 Improves Both Yield and Disease Resistance in Rice. Nat. Plants 2019, 5, 389–400. [Google Scholar] [CrossRef]
  77. Lecourieux, D.; Ranjeva, R.; Pugin, A. Calcium in Plant Defence-Signalling Pathways. New Phytologist 2006, 171, 249–269. [Google Scholar] [CrossRef]
  78. Costa, A.; Navazio, L.; Szabo, I. The Contribution of Organelles to Plant Intracellular Calcium Signalling. J. Exp. Bot. 2018, 69, 4175–4193. [Google Scholar] [CrossRef]
  79. Kaur, A.; Madhu; Taneja, M.; Upadhyay, S.K. Applications of Calcium Transport Elements in Plant Improvement. In Calcium Transport Elements in Plants; Elsevier: Amsterdam, The Netherlands, 2021; pp. 427–445. [Google Scholar]
  80. Chand, V.; Gupta, V. Interaction between Ca2+ and ROS Signaling in Plants. In Calcium Transport Elements in Plants; Elsevier: Amsterdam, The Netherlands, 2021; pp. 387–410. [Google Scholar]
  81. Ranf, S.; Eschen-Lippold, L.; Pecher, P.; Lee, J.; Scheel, D. Interplay between Calcium Signalling and Early Signalling Elements during Defence Responses to Microbe- or Damage-Associated Molecular Patterns. Plant J. 2011, 68, 100–113. [Google Scholar] [CrossRef] [PubMed]
  82. 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]
  83. 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] [PubMed]
  84. Nishimura, A.; Aichi, I.; Matsuoka, M. A Protocol for Agrobacterium-Mediated Transformation in Rice. Nat. Protoc. 2006, 1, 2796–2802. [Google Scholar] [CrossRef] [PubMed]
  85. Hartmut, K. Lichtenthaler Chlorophylls and Carotenoids: Pigments of Photosynthetic Biomembranes. Methods Enzymol. 1987, 350–382. [Google Scholar]
  86. Kauffman, H.E.; Reddy, A.P.; Hsieh, S.P.; Merca, S.D. An Improved Technique for Evaluating Resistance of Rice Varieties to Xanthomonas Oryzae. Plant Dis. Report. 1973, 57, 537–541. [Google Scholar]
  87. Wu, T.; Bi, Y.; Yu, Y.; Zhou, Z.; Yuan, B.; Ding, X.; Zhang, Q.; Chen, X.; Yang, H.; Liu, H.; et al. Activated Expression of Rice DMR6-like Gene OsS3H Partially Explores the Susceptibility to Bacterial Leaf Streak Mediated by Knock-Out OsF3H04g. Int. J. Mol. Sci. 2023, 24, 13263. [Google Scholar] [CrossRef]
  88. Thordal-Christensen, H.; Zhang, Z.; Wei, Y.; Collinge, D.B. Subcellular Localization of H2O2 in Plants. H2O2 Accumulation in Papillae and Hypersensitive Response during the Barley—Powdery Mildew Interaction. Plant J. 1997, 11, 1187–1194. [Google Scholar] [CrossRef]
  89. 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]
Figure 1. Phylogenetic tree of the OsACA family members and their expression in rice tissues. (A) Phylogenetic relationship between ACA proteins from O. sativa (blue), A. thaliana (green), S. lycopersicum (red), and Z. mays (black). (B) Cluster and gene structure of OsACAs. The neighbor-joining tree was constructed using 1000 bootstrap replicates. (C) Expression pattern of OsACA9 and OsACA11 in rice different tissues by RT-qPCR.
Figure 1. Phylogenetic tree of the OsACA family members and their expression in rice tissues. (A) Phylogenetic relationship between ACA proteins from O. sativa (blue), A. thaliana (green), S. lycopersicum (red), and Z. mays (black). (B) Cluster and gene structure of OsACAs. The neighbor-joining tree was constructed using 1000 bootstrap replicates. (C) Expression pattern of OsACA9 and OsACA11 in rice different tissues by RT-qPCR.
Ijms 25 01874 g001
Figure 2. The sgRNAs and gene editing target sites of OsACA9 by CRISPR-Cas9 in rice. (A) The sequences and target sites of the two different sgRNAs, sgRNA-1 and sgRNA-2, in the OsACA9 gene. The sgRNA sequences of the target are shown in black and the sequences of PAM (protospacer-adjacent motif) are shown in blue. (B) The sequencing results of the two sgRNA targets in the OsACA9 gene from different mutation lines.
Figure 2. The sgRNAs and gene editing target sites of OsACA9 by CRISPR-Cas9 in rice. (A) The sequences and target sites of the two different sgRNAs, sgRNA-1 and sgRNA-2, in the OsACA9 gene. The sgRNA sequences of the target are shown in black and the sequences of PAM (protospacer-adjacent motif) are shown in blue. (B) The sequencing results of the two sgRNA targets in the OsACA9 gene from different mutation lines.
Ijms 25 01874 g002
Figure 3. Resistance evaluation of Osaca9 to Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas oryzae pv. oryzicola (Xoc) strains and ROS detection in the Osaca9 mutants. (AC) Defense reactions of wild-type (ZH11) and the Osaca9 mutants to Xoo strains (PXO86, PXO71, POX99A, Zhe173, C2 and C5) (A,B) and Xoc strains (BLS256 and RS105) (C). (D) DAB staining of leaves of wild-type and the Osaca9 mutants. (E) Statistical analysis of H2O2 content. WT refers to the wild-type control, and aca9#1 and aca9#2 refer to two distinct knockout lines of OsACA9. ** indicates the significant difference at p < 0.01 by t-test.
Figure 3. Resistance evaluation of Osaca9 to Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas oryzae pv. oryzicola (Xoc) strains and ROS detection in the Osaca9 mutants. (AC) Defense reactions of wild-type (ZH11) and the Osaca9 mutants to Xoo strains (PXO86, PXO71, POX99A, Zhe173, C2 and C5) (A,B) and Xoc strains (BLS256 and RS105) (C). (D) DAB staining of leaves of wild-type and the Osaca9 mutants. (E) Statistical analysis of H2O2 content. WT refers to the wild-type control, and aca9#1 and aca9#2 refer to two distinct knockout lines of OsACA9. ** indicates the significant difference at p < 0.01 by t-test.
Ijms 25 01874 g003
Figure 4. Early-leaf-senescence phenotype of the Osaca9 mutants. (A) Plant phenotype at filling stage (bar = 10 cm). (B) The flag leaf (L1), second top leaf (L2), and third top leaf (L3) at filling stage. (C) Light-response curve of net photosynthesis rate (Pn) of the flag leaves at filling stage. (D) Chlorophyll a (Chla), Chlorophyll b (Chlb), and Carotenoid (Car) content of the leaves shown in (B). WT refers to the wild-type control, and aca9#1 and aca9#2 are two distinct knockout lines of the OsACA9 gene. ** indicates the significant difference at p < 0.01 by t-test.
Figure 4. Early-leaf-senescence phenotype of the Osaca9 mutants. (A) Plant phenotype at filling stage (bar = 10 cm). (B) The flag leaf (L1), second top leaf (L2), and third top leaf (L3) at filling stage. (C) Light-response curve of net photosynthesis rate (Pn) of the flag leaves at filling stage. (D) Chlorophyll a (Chla), Chlorophyll b (Chlb), and Carotenoid (Car) content of the leaves shown in (B). WT refers to the wild-type control, and aca9#1 and aca9#2 are two distinct knockout lines of the OsACA9 gene. ** indicates the significant difference at p < 0.01 by t-test.
Ijms 25 01874 g004
Figure 5. Agronomic characters of the Osaca9 mutants. (A) Pictures of main panicles (bar = 5 cm). (B) All grains per panicle (bar = 2 cm). (C) Length of ten grains (bar = 2 cm). (D) Flag leaf length. (E) Flag leaf width. (F) Plant height. (G) Main panicle length. (H) Grain number per main panicle. (I) Effective tiller number per plant. (J) Seed setting ratio. (K) Grain length. (L) Grain width. (M) Grain thickness. (N) 1000-grain weight. (O) Yield per plant. The number on the histogram indicates the p value of the significance between aca9#1 or aca9#2 and the wild-type control (WT) by t-test.
Figure 5. Agronomic characters of the Osaca9 mutants. (A) Pictures of main panicles (bar = 5 cm). (B) All grains per panicle (bar = 2 cm). (C) Length of ten grains (bar = 2 cm). (D) Flag leaf length. (E) Flag leaf width. (F) Plant height. (G) Main panicle length. (H) Grain number per main panicle. (I) Effective tiller number per plant. (J) Seed setting ratio. (K) Grain length. (L) Grain width. (M) Grain thickness. (N) 1000-grain weight. (O) Yield per plant. The number on the histogram indicates the p value of the significance between aca9#1 or aca9#2 and the wild-type control (WT) by t-test.
Ijms 25 01874 g005
Figure 6. Transcriptome analysis of the Osaca9 mutants and the schematic diagram of the predicted mechanism mediated by OsACA9. (A) KEGG pathway enrichment analysis of the differentially expressed genes (DEGs) between the Osaca9 mutants and the wild-type (WT) control. (B) Heatmap of the relative expression levels of representative DEGs involved in plant immunity. (C) Expression analysis of disease resistance-associated genes between WT and the Osaca9 mutation lines, using OsEF1 and OsActin as the standardization controls, respectively. (D) The schematic diagram of the predicted disease resistance mechanism mediated by OsACA9. Aca9#1 and aca9#2 are two distinct knockout lines of the OsACA9 gene. * and ** indicate the significance at p < 0.05 and p < 0.01 by t-test, respectively.
Figure 6. Transcriptome analysis of the Osaca9 mutants and the schematic diagram of the predicted mechanism mediated by OsACA9. (A) KEGG pathway enrichment analysis of the differentially expressed genes (DEGs) between the Osaca9 mutants and the wild-type (WT) control. (B) Heatmap of the relative expression levels of representative DEGs involved in plant immunity. (C) Expression analysis of disease resistance-associated genes between WT and the Osaca9 mutation lines, using OsEF1 and OsActin as the standardization controls, respectively. (D) The schematic diagram of the predicted disease resistance mechanism mediated by OsACA9. Aca9#1 and aca9#2 are two distinct knockout lines of the OsACA9 gene. * and ** indicate the significance at p < 0.05 and p < 0.01 by t-test, respectively.
Ijms 25 01874 g006
Table 1. Differentially expressed genes (DEGs) related to disease resistance between the Osaca9 mutants and the wild-type (WT) control.
Table 1. Differentially expressed genes (DEGs) related to disease resistance between the Osaca9 mutants and the wild-type (WT) control.
Gene IDGene NameDefinitionFold (log2)FunctionReference
Os01g0113500OsRLCK9receptor-like kinase2.642
Os01g0116900OsRLCK17LRK-type protein; protein kinase domain containing protein3.050Resistance to Xanthomonas oryzae[49]
Os01g0116800OsRLCK18TAK144.914Resistance to abiotic stress[49]
Os11g0695000OsRLCK354leucine-rich repeat receptor protein kinase EXS precursor6.877
Os11g0691800OsRLCK355receptor-like protein kinase precursor4.635
Os04g0127500OsWAK29OsWAK receptor-like protein kinase6.259
Os04g0370100OsWAK46OsWAK receptor-like protein kinase7.100
Os02g0618200OsRSR1two-component response regulator3.034Regulate tiller growth and panicle development[50]
Os11g0229500OsRSRdisease resistance protein9.878Resistance to SHB[51]
Os01g0124000OsRBBI2BBTI2-Bowman-Birk type bran trypsin inhibitor precursor6.995
Os09g0355400MDPKmalectin domain protein kinase9.393Resistance to SHB[52]
Os01g0719300OsSULTR3;6secondary active sulfate transmembrane transporter activity−4.383Resistance to bacterial leaf streak[53]
Os06g0306300OsPRX78peroxidase precursor3.878
Os09g0541000OsPIP2;7aquaporin, plasma membrane intrinsic protein4.199Tolerance to boron toxicity[54]
Os01g0901600OsACS6AMP-binding domain containing protein2.457Phenylpropanoid biosynthesis
Os04g0182200Os2ODD231-aminocyclopropane-1-carboxylate oxidase homolog 22.366Benzoxazinoid biosynthesis [55]
Os03g0289800Os2ODD8leucoanthocyanidin dioxygenase2.387Response to low-temperature stress and salt stress[56]
Os04g0604300OsXTH9glycosyl hydrolases family 163.210Cell wall modification[57]
Os08g0167800TPS30;46terpene synthase4.295Key gene for biosynthesis of limonene, methyl salicylate[58]
Os08g0168000OsTPS31terpene synthase activity4.549Resistance to insect[59]
Os07g0543000OsbHLH057bHLH transcription factor1.557Stress tolerance and yield[60]
Os06g0215900OsOPR412-oxophytodienoate reductase8.469Plant hormone response related; Jasmonic Acid biosynthesis related[61]
Os03g0324200OsEIN3ethylene signaling regulation factors6.398Involved in environmental humidity regulation of rice blast resistance in rice[62]
Os11g0117400OsWRKY104WRKY transcription factor2.929Response to biotic and abiotic stresses[63]
Os11g0685700OsWRKY61DNA-binding transcription factor activity9.098
Os12g0458100HSR201transferase family protein, putative, expressed5.811Hypersensitivity-related; required for pathogen signal-induced salicylic acid synthesis[64]
Os11g051450OsBAK1brassinosteroid insensitive 1-associated receptor kinase 1 precursor−6.408plant growth and development; PAMP triggered immune response[65]
Os10g0138100OsAO2, AAO2aldehyde oxidase12.345Regulates plant growth, grain yield, and drought tolerance[66]
Os04g0618700OsFLS2leucine-rich repeat receptor protein kinase EXS precursor; flg22 Receptor2.101PAMP triggered immune response[67]
Table 2. Primers used for qPCR.
Table 2. Primers used for qPCR.
GeneForward Primer (5′–3′)Reverse Primer (5′–3′)
OsActinGTTACTCATTCACCACAACGGCCCTTTCAGGAGGGGCGACC
OsEF1AAGAGGAAGTCAGCGGCTAAGCAGAATGGGCAGGAAAATACA
OsACA9CGGGACGACGACGATGGCCGCGCTTGGCGGGGATGT
OsACA11TGGCGGCAAGCAGCTCTAGTCGGAATGCAGCCCTGACAAC
OsBAK1ATCAACTGGGTGGAGAGTGAATCTCCCAAGGTGTGTGGTA
OsEIN3CCTCAAGAAGGCCTGGAAGGCTTGGCGGTCATCTTGTCCT
OsRLCK354TCCGTCCTTAACCTCACCATGGGATTTGTCCAGAGAGCTG
OsTPS30TACCGCTGGCTATACTCACACTCATGGCTTCTAGGCTTGG
OsAO2TGCTGGGAGAAAGCCGAAATTGGCACTCTGCACTTGACTT
OsFLS2AGAGATTGTCGCTCCATGCCAGAGCGAGTTGTTCTGGACG
OsRSR1AACACACTCCCAATCCGCTTAGGGTGGTTTTCCCCAATCC
MDPKGGACGCAACACCTGGATACTCGTGTATGGAGTGGTGCCAT
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

Wang, X.; Wang, Z.; Lu, Y.; Huang, J.; Hu, Z.; Lou, J.; Fan, X.; Gu, Z.; Liu, P.; Ma, B.; et al. OsACA9, an Autoinhibited Ca2+-ATPase, Synergically Regulates Disease Resistance and Leaf Senescence in Rice. Int. J. Mol. Sci. 2024, 25, 1874. https://doi.org/10.3390/ijms25031874

AMA Style

Wang X, Wang Z, Lu Y, Huang J, Hu Z, Lou J, Fan X, Gu Z, Liu P, Ma B, et al. OsACA9, an Autoinhibited Ca2+-ATPase, Synergically Regulates Disease Resistance and Leaf Senescence in Rice. International Journal of Molecular Sciences. 2024; 25(3):1874. https://doi.org/10.3390/ijms25031874

Chicago/Turabian Style

Wang, Xinyu, Ziyao Wang, Yiduo Lu, Jiani Huang, Zhuoer Hu, Junlei Lou, Xinyue Fan, Zhimin Gu, Pengcheng Liu, Bojun Ma, and et al. 2024. "OsACA9, an Autoinhibited Ca2+-ATPase, Synergically Regulates Disease Resistance and Leaf Senescence in Rice" International Journal of Molecular Sciences 25, no. 3: 1874. https://doi.org/10.3390/ijms25031874

APA Style

Wang, X., Wang, Z., Lu, Y., Huang, J., Hu, Z., Lou, J., Fan, X., Gu, Z., Liu, P., Ma, B., & Chen, X. (2024). OsACA9, an Autoinhibited Ca2+-ATPase, Synergically Regulates Disease Resistance and Leaf Senescence in Rice. International Journal of Molecular Sciences, 25(3), 1874. https://doi.org/10.3390/ijms25031874

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