Chlorophyll decomposition is accelerated in banana leaves after the long-term magnesium deficiency according to transcriptome analysis

Abstract

Magnesium (Mg) is an essential macronutrient for plant growth and development. Physiological and transcriptome analyses were conducted to elucidate the adaptive mechanisms to long-term Mg deficiency (MD) in banana seedlings at the 6-leaf stage. Banana seedlings were irrigated with a Mg-free nutrient solution for 42 days, and a mock control was treated with an optimum Mg supply. Leaf edge chlorosis was observed on the 9^th leaf, which gradually turned yellow from the edge to the interior region. Accordingly, the total chlorophyll content was reduced by 47.1%, 47.4%, and 53.8% in the interior, center and edge regions, respectively, and the net photosynthetic rate was significantly decreased in the 9^th leaf. Transcriptome analysis revealed that MD induced 9,314, 7,425 and 5,716 differentially expressed genes (DEGs) in the interior, center and edge regions, respectively. Of these, the chlorophyll metabolism pathway was preferentially enriched according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The expression levels of the five candidate genes in leaves were consistent with what is expected during chlorophyll metabolism. Our results suggest that changes in the expression of genes related to chlorophyll synthesis and decomposition result in the yellowing of banana seedling leaves, and these results are helpful for understanding the banana response mechanism to long-term MD.

Introduction

Crops absorb magnesium ions (Mg²⁺) from the soil mainly through their roots, which are usually affected by various factors, such as soil texture, cation exchangeable capacity and climatic factors [1]. The unique chemical property of Mg²⁺ results in a bond with negatively charged root cell walls, promoting ion exchange with the soil [2]. However, soil saturated with cations would lead to Mg deficiency (MD), particularly in acidic soils in high rainfall areas [3]. MD has a major negative consequence on crop production, leading to a detrimental effect on yields and poor food and feed quality [4]. Thus, maintaining Mg contents for agricultural production is very important. In plants, approximately 5–35% of Mg is detected in chloroplasts, and Mg is the central element of the tetrapyrrole ring in chlorophyll [5, 6]. Mg usually accounts for ~0.15%-0.35% of the dry composition of vegetative organs but reaches ~6%-20% in chlorophyll [7]. MD has become a serious problem as a result of heavy rainfall and the inappropriate use of nitrogen, phosphorus, and potassium fertilizers [8, 9].

A substantial pool of total cellular Mg is required to synthesize chlorophyll in photosynthetic tissues, while the rest is used for ribosome bridging during translation and for chelation with nucleotides, nucleic acids and other phosphate-containing compounds [10]. As the central atom of chlorophyll, Mg is essential for photosynthesis in green plants. Chlorophyll acts in pigment-protein complexes to capture light and transfer the electrons in PSI and PSII [11, 12]. The first step for chlorophyll synthesis is the insertion of Mg²⁺ into protoporphyrin IX catalyzed by Mg-chelatase (MgCh) [13]. Mg-chelatase is a key enzyme at the branch point in the synthetic pathway of heme and chlorophyll, which is composed of the ChlI, ChlH, and ChlD submits in plants [14–16]. A series of genes encoding the enzymes involved in the Mg branch of the chlorophyll biosynthesis pathway, viz., CHLH, CHL1, CHLD, CHLM, CRD1, DVR, POR A/B/C and CHLG, have been identified using in vitro and in vivo approaches [17].

The typical chlorosis symptom of MD appears in the leaf intervein and usually appears first on the lower and older leaves owing to the mobile nature of Mg [18, 19]. At the early stage of MD in plants, photoassimilates accumulate in source leaves before photosynthesis is suppressed, resulting in an excessive accumulation of carbohydrates and enhanced production of reactive oxygen species (ROS) [6]. Later, excessive sucrose regulates the chlorophyll A/B binding protein 2 (CAB2) gene in a feedback mode, leading to a decreased chlorophyll concentration [20]. In Arabidopsis, AtSGR1 and AtSGR2, which encode the enzyme Mg-dechelatase, are responsible for the breakdown of chlorophyll [21]. Under MD conditions, the expression of OsSGR is negatively regulated by ROS, affecting chlorophyll degradation in rice leaves [22]. miRNAs are also involved in plant MD tolerance [23]. As a result, MD affects dry matter production and carbohydrate partitioning in sugar bean, Arabidopsis, barley and Citrus sinensis [20, 24–26].

Banana is one of the world’s major food crops in the tropical and subtropical regions, where soil acidification has been a serious growing problem [27, 28]. As a result, MD has been a major risk factor for banana production. The impact of MD on the allocation of carbohydrates to different organs of the banana plant has been reported [29]. In contrast to model species such as Arabidopsis, rice, maize and soybean, the molecular mechanism of banana leaf chlorosis is unclear [21, 22, 30, 31]. Therefore, in the present study, an attempt was made to decipher the molecular mechanism of MD-induced banana leaf chlorosis through physiological and transcriptomic characterizations of banana (BaXi Jiao, Musa acuminata, AAA, cultivar Cavendish) seedling leaves exposed to long-term MD.

Materials and methods

Plant material and growth conditions

Banana (BaXi Jiao, M. acuminata, AAA, cultivar Cavendish) seedlings were purchased from the Danzhou Rapid Propagation of Banana Breeding Base (Hainan, China), cultivated in potting soil (Pindstrup, Denmark) until the 6-leaf stage and then used for MD treatment.

For MD treatment, all seedlings were transplanted to sand after rinsing the roots with double distilled water and then seedlings of uniform growth were randomly divided into two groups, viz., the MD and mock control groups. The modified Hoagland’s nutrient solution contained 4 mM Ca(NO₃)₂·4H₂O, 2 mM NH₄H₂PO₄, 4 mM KCl, 60 μM Fe-EDTA, 25 μM H₃BO₃, 2 μM MnSO₄·H₂O, 2 μM ZnSO₄·7H₂O, 0.5 μM CuSO₄·5H₂O, and 0.05 μM H₂MoO₄. For the mock control group, normal Hoagland’s nutrient solution containing 1 mM MgSO₄·7H₂O and 6 mM KNO₃ was used, while 1 mM K₂SO₄, 4 mM KNO₃ and 1 mM NH₄NO₃ were used to maintain the same K, S and N supply in the MD group. The two kinds of Hoagland’s nutrient solutions were supplied every 2 days. For each group, at least 15 plantlets were used, and 3 biological replicates were set. All seedlings were grown in the greenhouse of Hainan University. After 42 days of treatment, the 9^th leaves and the roots were collected, immediately frozen in liquid nitrogen and then stored at -80°C for RNA extraction. The 9^th-leaf samples were divided into the edge (LM_E), center (LM_C) and interior (LM_I) regions.

Photosynthesis and chlorophyll content measurement

After 42 days of treatment, photosynthetic parameters were measured in the 9^th leaf of the banana seedlings at 9:00–11:00 a.m. on a sunny day. All measurements were carried out using a CIRAS-3 portable photosynthesis system (PP Systems, USA). The photosynthetic photon flux density provided by a red/blue LED light source was amounted to 1200 μmol m^-2s^-1, the ambient CO₂ concentration was adjusted to 390 μmol mol^-1 by CO₂ injection and the leaf temperature was maintained at 27°C.

The relative chlorophyll content of each leaf was measured using a chlorophyll meter (SPAD-502 Plus; Konica Minolta), and the SPAD value of the 9^th leaf in each plant was recorded from 31 to 60 days after treatment. The measurements were taken from at least three biological replicates per treatment, and the values were averaged.

Biomass and Mg content analysis in leaves after long-term MD

After harvest, the plants were washed with deionized water to remove any residual ions. The roots, stems and leaves were heated at 105°C for 15–20 min, and then dried at 80°C until they reached a constant weight. The dry weight was measured and recorded.

For Mg content measurement, leaves were dried at 60°C for 3 days and then digested with nitric acid using a microwave digestion system (MILESTONE Ethos UP). After dilution in deionized water, the metal content in the samples was determined by ICP-MS (Agilent 7000 series).

Transcriptome library construction and analysis of differentially expressed genes (DEGs)

Fresh samples were ground into powder in liquid nitrogen, and total RNA (at least 1 μg) was extracted from LM_I, LM_C, LM_E, and control samples using the RNA-prep Pure Plant Plus Kit (#DP441, TIANGEN, Beijing, China) according to the manufacturer’s instructions. After pre-processing and rRNA removal, a cDNA library was constructed using the Illumina NovaSeq 6000 platform (Majorbio Biology Company, Shanghai, China).

Clean reads were mapped based on the Musa acuminate genome by using the HISAT2 program, and the read count for each gene was obtained from mapping results [32]. Gene expression levels were estimated using the RSEM program [33]. DEGs between the two groups (MD and Mg sufficiency) were analyzed using the DESeq2 R package (1.24.0). Genes with an adjusted P value < 0.05 and |log2FC| ≥ 1 filtering condition following DESeq2 analysis were considered to be differentially expressed between groups. Gene Ontology (GO) enrichment analysis of the DEGs was conducted using GOATOOLS to determine over- and underrepresented terms [34]. KOBAS was used to test the significance of the enriched DEGs in special Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [35].

Real-time quantitative polymerase chain reaction (qRT-PCR)

Five candidate genes related to chlorophyll metabolism were selected, and the primers were designed by Premier Primer 5.0 software (Premier Biosoft International, Palo Alto, CA, USA). The primer sequences and reference gene are listed in a supplemental table (S1 Table). cDNA synthesis was performed with a HiScript II 1st Strand cDNA Synthesis Kit (+gDNA wiper) (Vazyme, R212-01). Gene expression levels were analyzed by qRT-PCR on an ABI Q7 (Applied Biosystems, USA) using the ChamQ Universal SYBR qPCR Master Mix (Vazyme, Q711–02). The reaction mixtures, with a final volume of 5 μL, included 2× SYBR Mix, cDNA (1 μL), 0.2 μL each of forward and reverse primers, and PCR-grade water (3.6 μL). qPCR was performed under the following conditions: 95°C for 5 min, followed by 40 cycles at 95°C for 30s, 55°C for 30s, and 72°C for 40s. The qPCR was repeated three times for each gene in each sample, and experiments were performed in three biological replicates. The expression levels were calculated using the previous method [36]. The MaRPS2 gene was used as an internal control.

Results

Long-term MD affected the growth of banana

Mg is essential to plant growth and development, especially for photosynthesis. To investigate the physiological changes in banana leaves in response to long-term MD stress, banana (BaXi Jiao, M. acuminata, AAA, cultivar Cavendish) seedlings at the 6-leaf-stage were planted in quartz sand without a Mg supply for 42 days (Fig 1A). The 9^th leaf (L9) was visible etiolated and showed symptoms of chlorosis (Fig 1A and 1B), which gradually decreased in decreased in degree from the edge to the inside regions (Fig 1B).

[Figure omitted. See PDF.]

(A) The phenotype of banana seedlings after MD for 42 days. For the control groups, 1 mM Mg²⁺ Hoagland nutrient solutions were used, and 0 mM Mg²⁺ solutions were used for the long-term MD treatment (-Mg). Scale bars = 10 cm. (B) The phenotype of the 7^th– 11^th banana seedling leaves between the control and the MD treatment. Scale bars = 10 cm. (C) Mg content in different parts of the 9^th leaves (L9) (n = 3). (D) Fresh weight and dry weight of 9^th leaf. FW: fresh weight, DW: dry weight. Values represent the mean ± standard error (SE) (n = 3). (E) The different parts of the 9^th leaf were collected for RNA-seq. LM_I represents leaf interior, LM_E represents leaf edge and LM_C represents leaf center. Data analysis was performed by two-tailed t test, * means significantly different.

Compared to the control group, the Mg content of L9 decreased significantly (p < 0.01) under long-term MD in the LM_I, LM_C and LM_E regions (Fig 1C and 1E). Moreover, the Mg content gradually increased from the 8^th leaf (L8) to the 11^th leaf (L11) under the MD treatment, but showed the opposite trend in the control group (S1 Fig). However, MD stress had no significant (p = 0.99, 0.79) effect on the fresh and dried biomass of the 9^th leaf (Fig 1D). In conclusion, long-term MD led to chlorosis symptoms and lower Mg contents in banana leaves.

Chlorophyll content and photosynthesis rate under long-term MD

To explore the physiological process of banana leaves in response to long-term MD, the relative chlorophyll concentration (SPAD) was measured. The SPAD value of L9 showed a significant (p < 0.01) decrease after 60 days of MD treatment, with the differences starting at 40 days (Fig 2A). Moreover, the SPAD values were recorded from L11 to L6, and decreased significantly (p < 0.01) from L9 to L6 in the MD group after 42 days of Mg starvation treatment (Fig 2B). As Mg is a movable element in plants, it can be transferred from older leaves to younger leaves. The results showed that L9 was the first young leaf to respond to MD among all the leaves. Therefore, we chose L9 for subsequent analysis.

[Figure omitted. See PDF.]

(A) SPAD values of L9 recorded from 31 to 60 days in the control and MD treatments (n = 5). (B) SPAD value variation from L6 to L11 on the 42nd day in the control and MD treatments (n = 24). (C) Internal CO₂ concentrations of the MD treatment in L9 (n = 10). (D) Net photosynthetic rate (Pn) of the MD treatment in L9 (n = 10). (E-G) Chlb, Chla and Car contents in different regions of L9 (n = 3). Chla: chlorophyll a; Chlb: chlorophyll b; Car: carotenoid. Data analysis was performed by two-tailed t test, * means significantly different.

The intercellular CO₂ concentration was significantly (p = 0.02) increased, and the net photosynthetic rate of L9 was significantly (p < 0.01) reduced. Moreover, the contents of chlorophyll a, chlorophyll b and carotenoids in different regions of L9 leaves were significantly (p = 0.02, 0.02, 0.01) lower under long-term MD treatment in the interior region than in the control group (Fig 2E–2G). Similarly, the same trend was observed in L10 leaves (S2A–S2C Fig).

RNA sequencing (RNA-seq) and DEG analysis

To further reveal the molecular response mechanism of banana leaves to MD stress, banana L9 under long-term MD for 42 days were subjected to RNA-seq. We constructed four cDNA libraries, namely LM_E, LM_C, LM_I and Control, which represent the edge, center and interior regions of L9 in the MD and control samples. After removing the sequencing adaptor reads and low-quality data, we obtained 41,524,877 reads in the LM_E library, 44,306,550 reads in the LM_C library, 44,682,441 reads in the LM_I library, 36,332,613 reads in the Control library, and 84.64 Gb of RNA-seq data. More than 93.07% of the reads had a quality score of Q30 (Table 1).

[Figure omitted. See PDF.]

Principal component analysis (PCA) was used to evaluate the reliability of samples to explore differences within or between groups. There was a significant difference of 74.54% between LM_E and the control group, while the intragroup difference was 10.02% (Fig 3A). Replicate samples are generally considered reliable and suitable for further analysis. Samples from different groups were clustered together according to the treatment conditions. The control samples were clearly separated from the others, while the LM_I and LM_C samples were clustered together.

[Figure omitted. See PDF.]

(A) Principal component analysis (PCA) of gene expression. (B) Venn diagram of different samples from the MD and control treatments. (C) Venn diagram of the number of shared DEGs between different groups. (D) Upregulated and Downregulated DEGs for LM_I vs. Control, LM_C vs. Control and LM_E vs. Control.

A total of 20,114 transcripts were detected in all samples. In total, 80, 148, 278 and 595 transcripts were specifically expressed in the LM_E, LM_C, LM_I and control, respectively (Fig 3B). After using the filtering criteria to determine the significance of the differences in gene expression levels, there were 2,586 coexpressed DEGs in different regions of the L9, while 2,765, 871 and 1,395 DEGs were specifically expressed in LM_I, LM_C and LM_E, respectively (Fig 3C). Compared with the control group, there were 4,724, 3,894 and 2,135 DEGs were upregulated, and 4,592, 3,533 and 3,583 DEGs downregulated in LM_I, LM_C and LM_E, respectively (Fig 3D). Among the 2,586 common DEGs, 1,121, 1,029 and 1,119 DEGs were upregulated, and 1,465, 1,557, 1,467 DEGs were downregulated in LM_I, LM_C and LM_E, respectively (S3 Fig). The results showed that more common DEGs were downregulated than upregulated.

GO and KEGG enrichment analyses of DEGs

All DEGs found were classified into the cellular component (CC), molecular function (MF) and biological process (BP) GO categories. For the MF category, the most significant terms included ’binding’, ’catalytic activity’, ’transporter activity’ and ’nucleic acid binding transcription factor activity’. For BP terms, the most significantly overrepresented terms were ’metabolic process’, ’cellular process’, ’single-organism process’, ’biological regulation’ and ’regulation of biological process’. In addition, ’cell’, ’membrane’, ’cell part’, ’membrane part’ and ’organelle’ were significantly overrepresented terms in the CC category (Fig 4).

[Figure omitted. See PDF.]

In this study, KEGG pathway enrichment analysis of the DEGs revealed 131 pathways that were enriched, including ’carotenoid biosynthesis’, ’carbon fixation in photosynthetic organisms’, ’glycolysis/gluconeogenesis’, ’carbon metabolism’, and ’porphyrin and chlorophyll metabolism’ (Fig 5). According to the KEGG pathway analysis, there were 33 DEGs enriched in the carbon fixation in photosynthetic organisms (ko00710) pathway, including 16 downregulated DEGs and 9 upregulated DEGs. Additionally, 16 DEGs were enriched in the pathway of porphyrin and chlorophyll metabolism pathway (ko00860), of which 10 DEGs were downregulated and 6 DEGs were upregulated.

[Figure omitted. See PDF.]

DEGs involved in photosynthesis and chlorophyll metabolism

Physiological and transcriptome analyses revealed that MD treatment inhibited the photosynthetic capacity and chlorophyll content in banana leaves. A heatmap of 64 candidate genes differentially expressed between treatments (MD and control) was generated to visualize the expression pattern. A total of 21 genes were associated with photosynthesis-antenna proteins, while 17 genes were upregulated in LM_E (Fig 6).

[Figure omitted. See PDF.]

Red represents upregulated genes, and blue represents downregulated genes.

To confirm the accuracy of the RNA-seq data, 5 genes (Ma06_g24000 (MaCHLI), Ma09_g24880 (MaCHLM), Ma10_g13280 (MaCHLG), Ma02_g05320 (MaSGR1) and Ma07_g18910 (MaSGR2)) related to the chlorophyll metabolism pathway were selected to investigate their expression profiles using qRT-PCR. The verification results of these genes showed similar expression patterns in the RNA-seq data (Fig 7B). Three genes (MaCHLI, MaCHLM and MaCHLG) were downregulated, while 2 genes were upregulated under long-term MD in different banana leaf regions.

[Figure omitted. See PDF.]

(A) DEGs identified by RNA-seq are shown in colored blocks. Blue blocks represent downregulated genes, and red blocks represent upregulated genes. (B) qRT-PCR analysis and RNA-seq data of MaCHLI, MaCHLM, MaCHLG, MaSGR1 and MaSGR2 expression.

Discussion

Mg²⁺ is transported from roots to aerial tissues by MGR transporters for photosynthesis [37, 38]. Long-term MD triggers dramatic molecular responses in plants including microtubule-based movement, signal transduction, protein phosphorylation and regulation of light harvesting, and the photosynthesis antioxidant system [39, 40]. We characterized the decrease in biomass partitioning and carbohydrate distribution of banana leaves under MD [29]. To gain insight into the molecular response mechanism involved in long-term MD in different regions of banana leaves with the gradual inward extension of yellowing symptoms, RNA-seq analysis was performed to identify key genes and pathways that respond to long-term MD.

For crop plants, the first visual MD symptom is usually interveinal chlorosis in the old leaves owing to the relatively mobile nature of Mg in plants [41]. However, recent studies have shown that Mg remobilization is more vigorous in young mature leaves than in old leaves, which has been verified in sugar beet, Arabidopsis and rice [7, 9, 20]. In our experiment, young mature leaves also suffered from marginal chlorosis (Fig 1). These results indicate that Mg transport correlates with leaf vigor. Moreover, MD is known to inhibit the photosynthetic rate and lead to growth retardation and low production [42, 43]. A reduction in the photosynthetic pigment content is always accompanied by photosynthetic inhibition, while the molecular mechanism needs further in-depth study [44]. Here, we showed that the chlorophyll content of banana seedlings under long-term MD significantly decreased, while the photosynthetic pigment contents gradually decreased from the interior regions to the edge (Fig 2). The decrease in the photosynthetic pigment content in response to MD is a common phenomenon [22, 26, 45]. This observation indicates that long-term MD causes a loss of reaction centers connected with the light harvesting complex and photosystem and results in a reduction in the net photosynthetic rate and chlorophyll content.

Photosynthetic processes within chloroplasts require substantial amounts of Mg [46]. MD inhibits photosynthetic capacity and CO₂ assimilation, which leads to slow metabolism, thus affecting plant growth and development [20, 47, 48]. In Arabidopsis, early-term MD altered the expression of circadian clock genes in roots and triggered abscisic acid-responsive genes, whereas long-term MD altered the expression of genes involved in the ethylene biosynthetic pathway, reactive oxygen species detoxification and photoprotection of the photosynthetic apparatus [39, 49]. Long-term MD altered the expression of genes involved in signal transduction, the stress response, carbohydrate and energy metabolism, cell transport, cell wall and cytoskeleton metabolism, and nucleic acid and protein metabolisms in Citrus reticulata [50]. MD also regulated genes involved in lignin biosynthesis pathways, regulation of the cell cycle and cell wall metabolism, resulting in lignification, enlargement and cracking of the veins in the lower leaves of Citrus sinensis [51]. In our study, the carotenoid biosynthesis, carbon fixation in photosynthetic organisms, glycolysis/gluconeogenesis, carbon metabolism, and porphyrin and chlorophyll metabolism pathways were the most significantly affected according to the KEGG pathway analysis (Fig 5). In contrast to previous studies, we found that the pathway of porphyrin and chlorophyll metabolism pathways were dramatically changed. These results suggested that DEGs induced by long-term MD are involved in the energy and substance metabolism, which affects the growth of banana seedlings. Moreover, different plants responded to MD over varying period time.

According to the RNA-seq data, the majority of altered DEGs were located in the interior regions, whereas the edge had the least number of altered DEGs in L9 under long-term MD (Fig 3). These results suggest that the leaves may suffer from programmed cell death from the edge to the interior regions. In general, the best recognized function of Mg in plants is the formation of chlorophyll pigments, where the red substrate Mg-protoporphyrin IX monomethyl ester is converted to the green product 3,8-divinyl protochlorophyllide a [52]. MgCh catalyzes the first step committed to the synthesis of chlorophyll [53]. MgCh consists of three subunits—CHLH, CHLI, and CHLD—all of which undergo transcriptional and posttranslational modifications in plants and algae [54, 55]. Previous studies showed that the expression of the chlorophyll synthesis genes CHLI, ChlM and CHLG was downregulated under MD [20, 56–58]. Consistent with the KEGG pathway analysis, we found that 3 genes (MaCHLI, MaCHLM and MaCHLG) were simultaneously downregulated in different regions of banana (Fig 6). These results indicate that the chlorophyll biosynthesis capacity is reduced under long-term MD. Moreover, chlorophyll degradation is an important part of nutrient recycling and redistribution during plant stress [59]. OsSGR encodes a new chloroplast protein senescence-related gene and regulates chlorophyll degradation in chloroplasts [60]. The qRT-PCR results demonstrated that the expression of the two SGR genes was upregulated, leading to chlorophyll degradation of (Fig 7A). As described, long-term MD treatment inhibited chlorophyll synthesis and promoted chlorophyll degradation, resulting in leaf chlorosis and a decrease in photosynthesis in bananas.

Conclusion

This study revealed that at the physiological level and transcription levels, MD affected plant chlorosis and senescence processes, namely, the photosynthesis rate, CO₂ concentration and chlorophyll content, chlorophyll synthesis and degradation-related pathways. The expression of genes related to chlorophyll synthesis and degradation in the three leaf parts—edge, center and interior—was complicated, and involved many plant physiological responses and growth regulatory mechanisms. The results of this study can provide references for the study the functions of Mg in other crops and lay a foundation for the study of the physiological responses and molecular mechanisms of MD in banana.

Supporting information

S1 Fig. Ratio of magnesium ion content to leaf edge (LM_E), center (LM_C) and interior (LM_I) in L7-11.

https://doi.org/10.1371/journal.pone.0270610.s001

(TIF)

S2 Fig. Chlorophyll content in L10 after 42 days of treatment.

(A) Chla content in different L10 samples. (B) Chlb content in different L10 samples. C. Car content in different L10 samples.

https://doi.org/10.1371/journal.pone.0270610.s002

(TIF)

S3 Fig. Histogram showing the common upregulated and downregulated DEGs in LM_I, LM_C and LM_E, respectively.

https://doi.org/10.1371/journal.pone.0270610.s003

(TIF)

S1 Table. List of oligonucleotide sequences used as qRT-PCR primers.

https://doi.org/10.1371/journal.pone.0270610.s004

(DOCX)

Acknowledgments

The authors thank the lab members for assistance.

Citation: Kan B, Yang Y, Du P, Li X, Lai W, Hu H (2022) Chlorophyll decomposition is accelerated in banana leaves after the long-term magnesium deficiency according to transcriptome analysis. PLoS ONE 17(6): e0270610. https://doi.org/10.1371/journal.pone.0270610

About the Authors:

Baolin Kan

Contributed equally to this work with: Baolin Kan, Yong Yang

Roles: Formal analysis, Investigation, Visualization, Writing – original draft

Affiliation: Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, College of Tropical Crops, Hainan University, HaiKou, China

Yong Yang

Contributed equally to this work with: Baolin Kan, Yong Yang

Roles: Investigation, Visualization, Writing – review & editing

Affiliation: Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, College of Tropical Crops, Hainan University, HaiKou, China

Pengmeng Du

Roles: Formal analysis, Investigation, Writing – review & editing

Affiliation: Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, College of Tropical Crops, Hainan University, HaiKou, China

Xinping Li

Roles: Investigation, Writing – review & editing

Affiliation: Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, College of Tropical Crops, Hainan University, HaiKou, China

Wenjie Lai

Roles: Investigation, Writing – review & editing

Affiliation: Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, College of Tropical Crops, Hainan University, HaiKou, China

Haiyan Hu

Roles: Conceptualization, Data curation, Formal analysis, Writing – review & editing

E-mail: [email protected]

Affiliation: Hainan Key Laboratory for Sustainable Utilization of Tropical Bioresources, College of Tropical Crops, Hainan University, HaiKou, China

https://orcid.org/0000-0001-9197-3526

References

1. Hariadi YC, Shabala S. Screening broad beans (Vicia faba) for magnesium deficiency. I. Growth characteristics, visual deficiency symptoms and plant nutritional status. Funct Plant Biol. 2004;31(5):529–37. pmid:32688924

2. White PJ, Broadley MR, El-Serehy HA, George TS, Neugebauer K. Linear relationships between shoot magnesium and calcium concentrations among angiosperm species are associated with cell wall chemistry. Annals of Botany. 2018;122(2):221–6. pmid:29722830

3. Yan B, Hou Y. Effect of Soil Magnesium on Plants: a Review. IOP Conference Series: Earth and Environmental Science. 2018;170:022168.

4. Verbruggen N, Hermans C. Physiological and molecular responses to magnesium nutritional imbalance in plants. Plant Soil. 2013;368(1):87–99.

5. Cakmak I, Yazici AM. Magnesium: a forgotten element in crop production. Better crops with plant food. 2010;94(2):23–5.

6. de Bang TC, Husted S, Laursen KH, Persson DP, Schjoerring JK. The molecular–physiological functions of mineral macronutrients and their consequences for deficiency symptoms in plants. New Phytol. 2021;229(5):2446–69. pmid:33175410

7. Tanoi K, Kobayashi NI. Leaf Senescence by Magnesium Deficiency. Plants (Basel). 2015;4(4):756–72. Epub 2015/01/01. pmid:27135350

8. Gransee A, Fuhrs H. Magnesium mobility in soils as a challenge for soil and plant analysis, magnesium fertilization and root uptake under adverse growth conditions. Plant Soil. 2013;368(1):5–21.

9. Chen ZC, Peng WT, Li J, Liao H. Functional dissection and transport mechanism of magnesium in plants. Semin Cell Dev Biol. 2018;74:142–52. Epub 2017/08/22. pmid:28822768

10. Kleczkowski LA, Igamberdiev AU. Magnesium Signaling in Plants. Int J Mol Sci. 2021;22(3). Epub 2021/01/29. pmid:33503839

11. Jansson S. The light-harvesting chlorophyll a/b-binding proteins. Biochim Biophys Acta. 1994;1184(1):1–19. pmid:8305447

12. Pietrzykowska M, Suorsa M, Semchonok DA, Tikkanen M, Boekema EJ, Aro EM, et al. The light-harvesting chlorophyll a/b binding proteins Lhcb1 and Lhcb2 play complementary roles during state transitions in Arabidopsis. Plant Cell. 2014;26(9):3646–60. Epub 2014/09/07. pmid:25194026

13. Walker CJ, Weinstein JD. The magnesium-insertion step of chlorophyll biosynthesis is a two-stage reaction. Biochem J. 1994;299(1):277–84. pmid:8166650

14. Rissler HM, Collakova E, DellaPenna D, Whelan J, Pogson BJ. Chlorophyll biosynthesis. Expression of a second chl I gene of magnesium chelatase in Arabidopsis supports only limited chlorophyll synthesis. Plant Physiol. 2002;128(2):770–9. Epub 2002/02/14. pmid:11842180

15. Wu Z, Zhang X, He B, Diao L, Sheng S, Wang J, et al. A chlorophyll-deficient rice mutant with impaired chlorophyllide esterification in chlorophyll biosynthesis. Plant Physiol. 2007;145(1):29–40. Epub 2007/05/31. pmid:17535821

16. Shen Y-Y, Wang X-F, Wu F-Q, Du S-Y, Cao Z, Shang Y, et al. The Mg-chelatase H subunit is an abscisic acid receptor. Nature. 2006;443(7113):823–6. pmid:17051210

17. Masuda T. Recent overview of the Mg branch of the tetrapyrrole biosynthesis leading to chlorophylls. Photosynth Res. 2008;96(2):121–43. pmid:18273690

18. Billard V, Maillard A, Coquet L, Jouenne T, Cruz F, Garcia-Mina J-M, et al. Mg deficiency affects leaf Mg remobilization and the proteome in Brassica napus. Plant Physiol Biochem. 2016;107:337–43. pmid:27362297

19. Farhat N, Ivanov AG, Krol M, Rabhi M, Smaoui A, Abdelly C, et al. Preferential damaging effects of limited magnesium bioavailability on photosystem I in Sulla carnosa plants. Planta. 2015;241(5):1189–206. pmid:25637102

20. Hermans C, Verbruggen N. Physiological characterization of Mg deficiency in Arabidopsis thaliana. J Exp Bot. 2005;56:2153–61. pmid:15983014

21. Shimoda Y, Ito H, Tanaka A. Arabidopsis STAY-GREEN, Mendel’s Green Cotyledon Gene, Encodes Magnesium-Dechelatase. The Plant Cell. 2016;28(9):2147–60. pmid:27604697

22. Peng YY, Liao LL, Liu S, Nie MM, Li J, Zhang LD, et al. Magnesium Deficiency Triggers SGR–Mediated Chlorophyll Degradation for Magnesium Remobilization. Plant Physiology. 2019;181(1):262–75. pmid:31289214

23. Islam W, Tauqeer A, Waheed A, Zeng F. MicroRNA Mediated Plant Responses to Nutrient Stress. Int J Mol Sci. 2022;23(5):2562. pmid:35269700

24. Fischer ES, Bremer E. Influence of magnesium deficiency on rates of leaf expansion, starch and sucrose accumulation, and net assimilation in Phaseolus vulgaris. physiologia plantarum. 1993;89(2):271–6.

25. Tränkner M, Jákli B, Tavakol E, Geilfus C-M, Cakmak I, Dittert K, et al. Magnesium deficiency decreases biomass water-use efficiency and increases leaf water-use efficiency and oxidative stress in barley plants. Plant Soil. 2016;406(1):409–23.

26. Li C-P, Qi Y-P, Zhang J, Yang L-T, Wang D-H, Ye X, et al. Magnesium-deficiency-induced alterations of gas exchange, major metabolites and key enzymes differ among roots, and lower and upper leaves of Citrus sinensis seedlings. Tree Physiol. 2017;37(11):1564–81. pmid:28633366

27. Sanchez PA, Logan TJ. Myths and Science about the Chemistry and Fertility of Soils in the Tropics. Myths and Science of Soils of the Tropics 1992. p. 35–46.

28. Amini Khoozani A, Birch J, Bekhit AEA. Production, application and health effects of banana pulp and peel flour in the food industry. J Food Sci Technol. 2019;56(2):548–59. Epub 2019/03/25. pmid:30906012

29. He H, Jin X, Ma H, Deng Y, Huang J, Yin L. Changes of plant biomass partitioning, tissue nutrients and carbohydrates status in magnesium-deficient banana seedlings and remedy potential by foliar application of magnesium. Scientia Horticulturae. 2020;268:109377.

30. Li H, Wang N, Ding J, Liu C, Du H, Huang K, et al. The maize CorA/MRS2/MGT-type Mg transporter, ZmMGT10, responses to magnesium deficiency and confers low magnesium tolerance in transgenic Arabidopsis. Plant Mol Biol. 2017;95(3):269–78. pmid:28871377

31. Yang N, Jiang J, Xie H, Bai M, Xu Q, Wang X, et al. Metabolomics Reveals Distinct Carbon and Nitrogen Metabolic Responses to Magnesium Deficiency in Leaves and Roots of Soybean [Glycine max (Linn.) Merr.]. frontiers in plant science. 2017;8:2091. pmid:29312369

32. Hubert O, Piral G, Galas C, Baurens FC, Mbeguie AMD. Changes in ethylene signaling and MADS box gene expression are associated with banana finger drop. Plant Sci. 2014;223:99–108. Epub 2014/04/29. pmid:24767119

33. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12(1):323. pmid:21816040

34. Klopfenstein DV, Zhang L, Pedersen BS, Ramírez F, Warwick Vesztrocy A, Naldi A, et al. GOATOOLS: A Python library for Gene Ontology analyses. Scientific Reports. 2018;8(1):10872. pmid:30022098

35. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(suppl_2):W316–W22. pmid:21715386

36. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. Epub 2002/02/16. pmid:11846609

37. Meng SF, Zhang B, Tang RJ, Zheng XJ, Chen R, Liu CG, et al. Four plasma membrane-localized MGR transporters mediate xylem Mg2+ loading for root-to-shoot Mg2+ translocation in Arabidopsis. Mol Plant. 2022. Epub 2022/01/23. pmid:35063662

38. Bakshi A, Gilroy S. Moving magnesium. Mol Plant. 2022. Epub 2022/04/16. pmid:35422405

39. Hermans C, Vuylsteke M, Coppens F, Cristescu SM, Harren FJM, Inzé D, et al. Systems analysis of the responses to long-term magnesium deficiency and restoration in Arabidopsis thaliana. New Phytol. 2010;187(1):132–44. Epub 2010/04/20. pmid:20412444

40. Yang L-T, Zhou Y-F, Wang Y-Y, Wu Y-M, Ye X, Guo J-X, et al. Magnesium Deficiency Induced Global Transcriptome Change in Citrus sinensis Leaves Revealed by RNA-Seq. Int J Mol Sci. 2019;20(13):3129. pmid:31248059

41. Robinson JBD. Nutrient Deficiencies and Toxicities in Crop Plants. Edited By W. F. Bennett. St. Paul, Minnesota: APS Press (1993), pp. 202, $59.00 in USA, elsewhere US$74.00. ISBN 0-89054-151-5. Exp Agric. 1995;31(3):391-. Epub 2008/10/03.

42. Cakmak I, Kirkby EA. Role of magnesium in carbon partitioning and alleviating photooxidative damage. Physiologia plantarum. 2008;133(4):692–704. pmid:18724409

43. Hermans C, Bourgis F, Faucher M, Strasser RJ, Delrot S, Verbruggen N. Magnesium deficiency in sugar beets alters sugar partitioning and phloem loading in young mature leaves. Planta. 2005;220(4):541–9. pmid:15580527

44. Huang J-H, Xu J, Ye X, Luo T-Y, Ren L-H, Fan G-C, et al. Magnesium deficiency affects secondary lignification of the vascular system in Citrus sinensis seedlings. Trees. 2019;33(1):171–82.

45. Tang N, Li Y, Chen L-S. Magnesium deficiency–induced impairment of photosynthesis in leaves of fruiting Citrus reticulata trees accompanied by up-regulation of antioxidant metabolism to avoid photo-oxidative damage. J Plant Nutr Soil Sci. 2012;175(5):784–93.

46. Jamali Jaghdani S, Jahns P, Tränkner M. Mg deficiency induces photo-oxidative stress primarily by limiting CO(2) assimilation and not by limiting photosynthetic light utilization. Plant Sci. 2021;302:110751. Epub 2020/12/09. pmid:33287999

47. Farhat N, Rabhi M, Falleh H, Lengliz K, Smaoui A, Abdelly C, et al. Interactive effects of excessive potassium and Mg deficiency on safflower. Acta Physiol Plant. 2013;35(9):2737–45.

48. Li J, Yokosho K, Liu S, Cao HR, Yamaji N, Zhu XG, et al. Diel magnesium fluctuations in chloroplasts contribute to photosynthesis in rice. Nature Plants. 2020;6(7):848–59. pmid:32541951

49. Hermans C, Vuylsteke M, Coppens F, Craciun A, Inzé D, Verbruggen N. Early transcriptomic changes induced by magnesium deficiency in Arabidopsis thaliana reveal the alteration of circadian clock gene expression in roots and the triggering of abscisic acid-responsive genes. New Phytol. 2010;187(1):119–31. pmid:20406411

50. Jin XL, Ma CL, Yang LT, Chen LS. Alterations of physiology and gene expression due to long-term magnesium-deficiency differ between leaves and roots of Citrus reticulata. J Plant Physiol. 2016;198:103–15. Epub 2016/05/11. pmid:27163764

51. Ye X, Huang H-Y, Wu F-L, Cai L-Y, Lai N-W, Deng C-L, et al. Molecular mechanisms for magnesium-deficiency-induced leaf vein lignification, enlargement and cracking in Citrus sinensis revealed by RNA-Seq. Tree Physiol. 2021;41(2):280–301. pmid:33104211

52. Chen GE, Hitchcock A, Mareš J, Gong Y, Tichý M, Pilný J, et al. Evolution of Ycf54-independent chlorophyll biosynthesis in cyanobacteria. Proceedings of the National Academy of Sciences. 2021;118(10):e2024633118. pmid:33649240

53. Jensen PE, Gibson LCD, Henningsen KW, Hunter CN. Expression of the chlI, chlD, and chlH Genes from the Cyanobacterium Synechocystis PCC6803 in Escherichia coli and Demonstration That the Three Cognate Proteins Are Required for Magnesium-protoporphyrin Chelatase Activity*. J Biol Chem. 1996;271(28):16662–7. pmid:8663186

54. Willows RD. Biosynthesis of Chlorophylls from Protoporphyrin IX. ChemInform. 2003;20(3):327–41. pmid:12828371

55. Gao Y-S, Wang Y-L, Wang X, Liu L. Hexameric structure of the ATPase motor subunit of magnesium chelatase in chlorophyll biosynthesis. Protein Sci. 2020;29(4):1026–32. pmid:31891428

56. Bhor SA, Tateda C, Mochizuki T, Sekine KT, Yaeno T, Yamaoka N, et al. Inducible expression of magnesium protoporphyrin chelatase subunit I (CHLI)-amiRNA provides insights into cucumber mosaic virus Y satellite RNA-induced chlorosis symptoms. Virusdisease. 2017;28(1):69–80. Epub 2017/05/04. pmid:28466058

57. Wang Z, Hong X, Hu K, Wang Y, Wang X, Du S, et al. Impaired Magnesium Protoporphyrin IX Methyltransferase (ChlM) Impedes Chlorophyll Synthesis and Plant Growth in Rice. Front Plant Sci. 2017;8:1694. Epub 2017/10/17. pmid:29033966

58. Shalygo N, Czarnecki O, Peter E, Grimm B. Expression of chlorophyll synthase is also involved in feedback-control of chlorophyll biosynthesis. Plant Mol Biol. 2009;71(4–5):425–36. Epub 2009/08/15. pmid:19680747

59. Bresson J, Bieker S, Riester L, Doll J, Zentgraf U. A guideline for leaf senescence analyses: from quantification to physiological and molecular investigations. J Exp Bot. 2018;69(4):769–86. Epub 2017/10/11. pmid:28992225

60. Park SY, Yu JW, Park JS, Li J, Yoo SC, Lee NY, et al. The senescence-induced staygreen protein regulates chlorophyll degradation. Plant Cell. 2007;19(5):1649–64. Epub 2007/05/22. pmid:17513504