- Research
- Open access
- Published:
Pan-transcriptomic Profiling Demarcates Serendipita Indica-Phosphorus Mediated Tolerance Mechanisms in Rice Exposed to Arsenic Toxicity
Rice volume 16, Article number: 28 (2023)
Abstract
Inadvertent accumulation of arsenic (As) in rice (Oryza sativa L.) is a concern for people depending on it for their subsistence, as it verily causes epigenetic alterations across the genome as well as in specific cells. To ensure food safety, certain attempts have been made to nullify this highest health hazard encompassing physiological, chemical and biological methods. Albeit, the use of mycorrhizal association along with nutrient reinforcement strategy has not been explored yet. Mechanisms of response and resistance of two rice genotypes to As with or without phosphorus (P) nutrition and Serendipita indica (S. indica; S.i) colonization were explored by root transcriptome profiling in the present study. Results revealed that the resistant genotype had higher auxin content and root plasticity, which helped in keeping the As accumulation and P starvation response to a minimum under alone As stress. However, sufficient P supply and symbiotic relationship switched the energy resources towards plant’s developmental aspects rather than excessive root proliferation. Higher As accumulating genotype (GD-6) displayed upregulation of ethylene signaling/biosynthesis, root stunting and senescence related genes under As toxicity. Antioxidant defense system and cytokinin biosynthesis/signaling of both genotypes were strengthened under As + S.i + P, while the upregulation of potassium (K) and zinc (Zn) transporters depicted underlying cross-talk with iron (Fe) and P. Differential expression of phosphate transporters, peroxidases and GSTs, metal detoxification/transport proteins, as well as phytohormonal metabolism were responsible for As detoxification. Taken together, S. indica symbiosis fortified with adequate P-fertilizer can prove to be effective in minimizing As acquisition and accumulation in rice plants.
Background
Arsenic (As) becomes a constituent of food chain via its ingestion through contaminated water or food laden with As, which is either geogenically and/or anthropogenically released in the environment (Norton et al. 2012). Rice (Oryza sativa L.), due to its heavy irrigation requirement and cultivation in anaerobic soils, gets exposed to polluted water, making its grains tainted with inorganic arsenic species ‘Asi’ (i.e., arsenate—AsV and arsenite— AsIII), and the organic forms dimethylarsinic acid, with trace amount of monomethylarsinic acid and tetramethylarsonium (Hansen et al. 2011). In China alone, almost 50% agricultural soil samples have verified As concentration of more than 40 mg kg− 1, which exceeds the threshold value as listed in “China Soil Environmental Quality (GB15618–2018)” (Lin et al. 2015). The related acute/chronic health damages for people on a rice-subsistence diet and the cumulative ravaging impacts on living cells, due to its competitive interference in P-reliant biological processes, make arsenic’s exclusion much more critical (Sehar et al. 2022). Because plant roots are the first organs coming in contact with xenobiotic compounds/stressors, and function at the interface, it is imminent to study the avoidance or tolerance mechanisms initiated at this end to cope with As stress. Having a strong inter-ionic competition between them (Boorboori & Zhang 2022), the chemical similarity of AsV with P leads to its uptake, rather indistinguishably, through phosphate transporters. In Arabidopsis root, disarray of the P transporters, e.g., Pho1;1 and Pho1;4, provided tolerance to AsV (Shin et al. 2004). Kamiya et al. (2013) described PO43− deficiency signaling pathway as a target of AsV and the involvement of OsPT1 in shoot As accumulation in rice. Moreover, knockout of OsPT1 decreased AsV uptake and Asi accumulation in rice grains (Cao et al. 2017). As the P/As molar ratio is integral for the regulation of As mobility and accumulation, controlling the interactions between P and As in soil, however, is very complicated.
The results obtained from transcriptomic studies involving different cell types help us gain a deeper understanding regarding the constituents of a specific cell type, normal cellular functions, and the aberrations that may result from deviated gene activity (Huang et al. 2012). In order to understand the As toxicity mechanism, several microarray-based gene expression studies have been conducted (Dangleben et al. 2013). Arsenic exposure leads to adaptational changes, which encompass chromatin structure, transcriptome and the production of specific gene isoforms that, even upon withdrawal, do not lead to the complete recovery of chromatin structures and/or gene expression patterns (Riedmann et al. 2015). Norton et al. (2008) studied whole genome of rice roots exposed to arsenate in hydroponic conditions and denoted glutathione conjugation and AsV methylation as the most crucial biochemical responses to AsV stress. A separate analysis of AsV and AsIII stresses at transcriptome level revealed differential expression of transporters, defense and stress-responsive genes, metallothioneins, sulfate-metabolizing and heat-shock proteins as well as regulatory genes, while the induction of cytochrome P450 genes occurred only under AsV stress (Chakrabarty et al. 2009). Furthermore, ethylene induced lignification occurred along with the upregulation of ATP-binding cassette superfamily transporters, ROS response network genes, mitogen-activated protein kinases (MAPKs) and calcium-dependent protein kinases (CDPKs), MYB and zinc-finger proteins, and senescence related (GARP-G2-like) transcription factors (Huang et al. 2012). Arsenate disrupts IAA biosynthesis, transport, and localization by negatively affecting the expression of IAA biosynthesis, i.e., OsASA2 and OsYUCCA2, and transporter genes, i.e., AUX1 and PIN5b (Ronzan et al. 2019). Arsenic notoriously impose iron (Fe) deficiency in plants, which in turn interferes with chlorophyll biosynthesis (Das and Sarkar 2018).
Comparatively, technology of bioremediation is taking precedence in attenuating As absorption in rice. According to prior reports, microbes could meliorate As accumulation by regulating As speciation, its bioavailability, and plant’s metabolism (Shri et al. 2019). Besides, the transcript levels of OsLsi1 and OsLsi2 were found to be repressed by arbuscular mycorrhizal fungi (AMF) in rice plants, thereby reducing As uptake/translocation (Shri et al. 2019). Transcriptional studies involving mycorrhizal fungus Rhizophagus irregularis colonization in rice roots resulted in a decrease in secondary cell-wall (SCW) related transcripts correlated with modified SCW composition; also, the hormonal profile and transport-related transcripts across different root types demonstrated a potential impact of symbiosis on root functioning and architecture (Gutjahr et al. 2015). In recent years, the focus of research has advanced to a specific group of AMFs, namely ‘Sebacinales’, including two phylogenetic-subgroups, i.e., Sebacinaceae and Serendipitaceae (Weiß et al. 2011); the latter contains several species, including Serendipita indica (formerly known as Piriformospora indica) found in Thar deserts of India (Varma et al. 1999), and others associated with plant family Orchidaceae (Oktalira et al. 2021). This endophyte is capable of establishing mutualistic association with many monocotyledonous and dicotyledonous plants (Franken 2012). Colonization of soybean roots with S. indica have been reported to instigate an upregulatory expression of abscisic acid-, gibberellin- and auxin- synthesis, signaling, or transport genes; where the auxin responsive genes displayed potential of reshaping root development and architecture (Bajaj et al. 2018). It has proven to adroitly tolerate, detoxify as well as immobilize As by decreasing its bioavailability via accumulation, adsorption and precipitation (Mohd et al. 2017; Ghorbani et al. 2021). Also, the remarkable P solubilization and utilization capabilities of this fungus assure its survival in the As-spiked growth conditions (Kushwaha et al. 2022). Pertinently, the employment of S. indica facilitates the uptake of available nutrients but cannot compensate the deficit thereof. Though several studies have delivered single puzzle pieces to the whole picture, the S. indica colonization-mediated arsenic tolerance in addition to P nutrition and related molecular mechanisms are quite understudied. In this work, the transcriptional changes of the two rice genotypes differing in their responses to As accumulation accompanied by P and S. indica colonization as nutrient management strategy and myco-remediative tool have been focused on, respectively.
Materials and Methods
Plant Material Serendipita Indica Culture Preparation and Inoculation
Healthy seeds of two rice genotypes, Zhongzhe You-1 (ZZY-1; G1; low As accumulator) and Guodao-6 (GD-6; G2; high As accumulator), were disinfected using 3% perhydroxic acid for about 15 min, then rinsed 3 times with ddH2O. Thereafter, the seeds of each genotype were soaked in a labeled germination box overnight at room temperature (25 °C). Afterward, the well moist seeds were placed for germination in a rice growth chamber at 29 °C ± 1 °C during the day (12 h) and 26 °C ± 1 °C at night. A hydroponic setup was arranged and 10 days old, uniform sized seedlings were transferred to half strength nutrient solution (Zeng et al. 2008), for 3 days followed by full strength solution for 4 days before fungal inoculation. Detailed composition of nutrient solution is mentioned in Additional file 1: Table S1. Serendipita indica was first grown on agar Petri plates with Aspergillus medium; circular agar disks (5 mm in diameter) inclosing active spores of S. indica (Hill and Käfer 2001), were placed on solidified medium and kept in absolute dark for 7 days at 30 °C ± 1 °C. For inoculation, S. indica was grown in a liquid Käfer medium from previously inoculated Petri plates; a small circular section of about 1 mm diameter of the inoculated fungus was introduced into a 500 mL glass flask and then incubated for 12 days at 30 °C ± 1 °C in a shaking incubator (200 rpm). Thereafter, the liquid containing active spores were filtered through sterile muslin cloth then used for plant inoculation by root soaking method (Chadha et al. 2015).
Confirmation of Serendipita Indica Root Colonization and Treatment Application
Seven-day post-inoculation, successful colonization was confirmed by microscopic visualization of S. indica spores in the root cortical tissues of the host plant. Inoculated plant roots were first gathered and rinsed thoroughly with water before cutting them into 1 cm long pieces and then immersed in KOH solution (10%) overnight at room temperature (25 °C; Chadha et al. 2015). The roots were thereafter, washed three times with ddH2O and then immersed again, this time with 1% HCl for about 3 min prior to trypan blue (0.05%) staining. Followed by successful colonization, plants were kept either as positive control (S.i) or treated with alone 10 µM As L− 1 (As + S. indica), 50 µM P L− 1 (P + S. indica) and their combination (As + P + S. indica). A set of un-inoculated plants were also treated (As and As + P), meanwhile seedlings without any treatment were designated as control. With the aim of clarifying the transcriptomic profile of two genotypes, ZZY-1 (G1) and GD-6 (G2), under As (T1), As + P (T2), As + S.i (T3) and As + S.i + P (T4) treatments along with control (T5), cDNA from total 30 root samples (5 × 3 × 2, treatments×replicates×genotypes, respectively) was prepared and clean paired-end reads were subjected to sequencing. Arsenic (10 µM) was provided as Na2HAsO4.7H2O (0.0078 g L− 1), whereas for 50 µM phosphorus, 0.25 mM NaH2PO4.2H2O (0.039 g L− 1) was added (stock solutions).
RNA Extraction, Library Preparation and Sequencing
Extraction of total RNA from the root tissues was carried out using TRIzol® reagent according to the instructions provided by manufacturer (Magen). Quantitative examination of RNA samples involved absorbance ratio analysis at A260/A280 through a Nanodrop ND-2000 system (Thermo Scientific, USA), and the integrity of RNA was determined by an Agilent Bioanalyzer 4150 system (Agilent Technologies, CA, USA). Paired-end library construction required ABclonal mRNA-seq Lib Prep Kit (ABclonal, China). Briefly, the purification of mRNA was done using l µg total RNA along with oligo (dT) magnetic beads, and subsequent fragmentation using divalent cations in ABclonal First Strand Synthesis Reaction Buffer at elevated temperatures. Accordingly, random hexamers were needed for the synthesis of first-strand cDNAs, while mRNA fragments served as templates for Reverse Transcriptase (RNase H), followed by synthesis of second-strand cDNA, which utilized DNA polymerase I, RNAseH, buffer and dNTPs. Post-synthesis, the double stranded cDNA fragments were adapter-ligated for paired-end library preparation, and later on used for PCR amplification. The resultant products were purified (AMPure XP system) and Agilent Bioanalyzer 4150 system assisted in assessing their quality. Libraries were sequenced using Illumina Novaseq 6000 (or MGISEQ-T7), which helped generate 150 bp paired-end reads. Afterwards, the bioinformatics analyses were performed using an in-house pipeline from Shanghai Applied Protein Technology Co., Ltd. The key softwares and procedures are as follows:
Quality Control and Mapping
Firstly, the raw reads (in fastq format) were processed to obtain clean ones through perl scripts, which involved removal of the adapter sequence and filtering of low-quality reads (the number of lines with a string quality value ≤ 25 accounts for more than 60% of the entire reading), keeping N (undetermined base information) ratio > 5% of reads. The acquired reads were then aligned to reference genome separately, using HISAT2 software (http://daehwankim1ab.eithub.io/hisat2/), to obtain mapped reads.
Quantitative Analysis of Differential Gene Expression Levels
Numbers of reads mapped to each gene were counted via FeatureCounts (http://subread.sourceforue.net/). In due course, the FPKM value for each gene was determined predicated on the length of the gene and corresponding mapped reads’ count. Analysis of differential expression was carried out using DESeq2 (http://bioconductor.ore/packages/release/bioc/html/DESeq2.html), where the DEGs with | log2FC | > 1 and Padj < 0.05 were designated as significantly different.
Enrichment and Transcription Factor Analysis
Differential genes were evaluated and functionally annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to explain the differences between samples. To serve this purpose, clusterProfiler R software package was employed, where the values displaying P < 0.05 were considered significantly enriched for the GO or KEGG function in question. Correspondingly, TF analysis of DEGs involved direct data extraction from PlantTFDB database (http://planttfdb.cbi.pku.edu.cn/), whereas InterProScan tool (https://www.ebi.ac.uk/interpro/) assisted in annotating gene functions, and the predictive screening of putative transcription factors was carried out using DBD database (https://db1p.uni-trier.de/rec/journals/nar/Kummerfe1dT06.html).
Verification of RNA-Seq Results with qRT-PCR
The RNA-Seq results were verified by performing qRT-PCR analysis following the methods discussed in detail by Adil et al. (2020). For internal reference, the rice actin gene (accession number XM_015774830) was used with the primer sequences: F— 5´TCCTCCGTGGAGAAGAGCTA3´ and R—5´GCAATGCCAGGGAACATAGT3´. Oligonucleotide primers for the nine genes related to sugar and metal transport, detoxification and cell-wall components were used, and the results determined by the 2(−ΔΔCt) method (Schmittgen and Livak 2008), agreed with the whole transcriptome findings (displaying a correlation coefficient = 0.893) (Additional file 2: Fig. S1 and Additional file 1: Table S2).
Statistical Analysis
Data were expressed as mean ± standard error of three biological replicates for each sample. Statistical analyses were performed using the Graphpad Prism (Graphpad Software, San Diego, CA, USA). Statistical significance was assessed at P < 0.05 significance level using Duncan’s multiple comparison.
Results
With the aim of clarifying the transcriptomic profile of two genotypes, ZZY-1 (G1) and GD-6 (G2), under As (T1), As + P (T2), As + S.i (T3) and As + S.i + P (T4) treatments along with control (T5), cDNA from total 30 root samples was prepared and clean paired-end reads were subjected to sequencing. About 8.9 Gb reads sample− 1 were generated with a GC content of 46.9–51.66%, Q30 ˃ 92.8%, and an alignment efficiency of 69.5–90.1%, from which higher than 67% reads were uniquely mapped (Additional file 1: Table S3). Cumulatively, 59,833,545 clean reads were generated from each library. Most importantly, the sequencing quality was quite satisfactory to proceed analyzing further (Additional file 2: Fig. S2). Novel transcripts were merged with the reference sequence of Oryza sativa spp. japonica, thereafter clean reads were mapped and gene expression levels were calculated for each sample under control and four treatments exploiting RSEM tool (RNA-Seq by Expectation-Maximization). As shown in Additional file 1: Table S4, a total of 3860 (11.2%) novel genes were registered encoding proteins with unknown functions, including some that matched the description of putative CorA-like Mg2+ transporter protein [Oryza sativa Japonica Group] protein kinase domain (Novel.18,755), probable aquaporin TIP3-1 (Novel.30,957), NADH dehydrogenase subunit 6 (mitochondrion; Novel.10,674), mitochondrial outer membrane protein porin 2 (Novel.24,621), chloroplastic NADH-glutamate synthase 2 precursor (Novel.25,513).
Gene Expression Analysis and Identification of DEGs
Out of 34,575 total genes, the number of differentially expressing genes was found highest in G2T5 vs. G2T1 (29.8%), followed by G1T1 vs. G2T1 (Fig. 1A). Pairwise comparisons displayed more downregulated genes than upregulated ones overall. Succinctly, the comparison, i.e., G1T1 vs. G1T2, G1T1 vs. G1T4, G1T1 vs. G2T1, G1T3 vs. G2T3, G2T1 vs. G2T2, G2T1 vs. G2T3 and G2T1 vs. G2T4, brought about a high number of upregulated DEGs, i.e., 65.2%, 64.5%, 50.2%, 51.1%, 50.2%, 50% and 51.5%, respectively. There were 2577 commonly expressed, while 4860 and 1907 uniquely expressed genes under As (G1T1 vs. G2T1) and As + S.i + P (G1T4 vs. G2T4), respectively, indicating that a huge number of genes expressed upon arsenic toxicity stimulus (Fig. 1B and C). Among the genotypes across the treatments, 839, 1281, 1821 and 1884 genes were uniquely expressed for G1T4 vs. G2T4, G1T3 vs. G2T3, G1T2 vs. G2T2 and G1T1 vs. G2T1, respectively, while 182 DEGs were expressed mutually.
Comparative stats of differentially expressed transcripts and their counts (A) across the treatments between ZZY-1 and GD-6. Following treatments were employed: As10µM (T1), As10µM+P50µM (T2), As10µM+S.i (T3), As10µM+S.i + P50 µM (T4) and control (T5). Venn diagram representing number of differentially expressed genes classified into differential interactions of As/S.i/P; (B) Overlaps among DEGs within group comparisons between ZZY-1 vs. GD-6 (G1T1 vs. G2T1, G1T2 vs. G2T2, G1T3 vs. G2T3 and G1T4 vs. G2T4); and (C) G1T1 vs. G2T1 and G1T4 vs. G2T4.
Genotypic comparisons showed that ZZY-1 genotype had 353, 107, 828 and 2978 unique genes expressed differentially for G1T1 vs. G1T2, G1T1 vs. G1T3, G1T1 vs. G1T4 and G1T5 vs. G1T1, respectively, while 175 genes expressed mutually. Guodao-6 displayed 571, 156, 780 and 3464 DEGs when compared for G2T1 vs. G2T2, G2T1 vs. G2T3, G2T1 vs. G2T4 and G2T5 vs. G2T1, respectively, with 956 co-expressing genes (Additional file 2: Fig. S3). Collectively, the control group comparison of ZZY-1 and GD-6 (G1T5 vs. G2T5) presented 5486 DEGs, which setting against G1T1 vs. G2T1, G1T2 vs. G2T2, G1T3 vs. G2T3 and G1T4 vs. G2T4, shared 2854, 3190, 2730 and 2366 DEGs, respectively. On the one hand, when ZZY-1 control against As10μM group (G1T5 vs. G1T1; 4343 DEGs) was compared with its own treatment groups, i.e., G1T1 vs. G1T2, G1T1 vs. G1T3 and G1T1 vs. G1T4, it shared 490, 107 and 1121 DEGs, respectively. On the other hand, GD-6 control against As10μM group (G2T5 vs. G2T1; 10,303 DEGs) shared 4685, 818 and 5751 DEGs with G2T1 vs. G2T2, G2T1 vs. G2T3 and G2T1 vs. G2T4 groups, respectively.
Correspondingly, principle component analysis (PCA) revealed that samples of the two genotypes were clearly separated by PC1 and PC2, which could explain for 64.7% of the total variation, signifying the differential genotypic behavior (Additional file 2: Fig. S4); moreover, the controls were clearly separated from treatments, indicating that arsenic toxicity posed a substantial impact on the transcript levels. Relevantly, the volcano plots for the pairwise treatment comparisons depicted that the highest number of down-regulated DEGs belonged to G1T2 vs. G2T2 (4114; 57.6%); whereas G1T4 vs. G2T4 reportedly exhibited the highest percentage (71.7%) of down-regulated genes (Fig. 2A-D). Among the genotype based pairwise comparisons, G1T1 vs. G1T3 presented 62.6% downregulated DEGs as opposed to the G1T1 vs. G1T2, G1T1 vs. G1T4 and G1T5 vs. G1T1 (34.8%, 35.4% and 57.6%, respectively). In case of GD-6, G2T5 vs. G2T1 comparison had higher percentage of down-regulated DEGs (54.9%) compared to its other three counterparts. Contingent on expression levels, all genes could be placed into 6 sub clusters (Additional file 2: Fig. S5) and 7 color modules (Additional file 2: Fig. S6).
Volcanogram representing differentially expressed transcripts analyzed by DEseq2 in roots of two rice genotypes (G1—ZZY-1 and G2—GD-6) in response to As, P and S. indica. (A), G1T1 vs. G2T1; (B), G1T2 vs. G2T2; (C), G1T3 vs. G2T3; and (D), G1T4 vs. G2T4. X axis represents log2 transformed fold change. Y axis represents -log10 (adjusted P value). T1 = As10µM, T2 = As10µM+P50µM, T3 = As10µM+S.i and T4 = As10µM+S.i + P50 µM. Maroon points represent upregulated DEGs, blue points represent downregulated ones and gray dots represent non-DEGs.
Pathway Analyses and Functional Enrichment of DEGs via KEGG
Information related to the expressed genes and their putative functions were determined through KEGG pathway classification/functional enrichment analysis (Fig. 3A and S7). During sorting, the smallest P-values (< 0.05) were given weightage and the number of genes enriched in the first 20 entries were counted. Among the top 20 KEGG pathways, ‘linoleic acid metabolism’, ‘sesquiterpenoid and triterpenoid biosynthesis’ and ‘phenylpropanoid biosynthesis’ were enriched mutually for all treatment wise comparisons. ‘Starch/ sucrose metabolism’ and ‘α-linolenic acid metabolism’ were enriched mutually for G1T1 vs. G2T1, G1T3 vs. G2T3 and G1T4 vs. G2T4, while ‘flavone and flavonol biosynthesis’, plant hormone signal transduction’ ‘MAPK signaling’, ‘biosynthesis of secondary metabolites’, ‘diterpenoid biosynthesis’, ‘isoquinoline alkaloid biosynthesis’ and ‘plant-pathogen interaction’ were inclusively and significantly enriched KEGG pathways for G1T2 vs. G2T2, G1T3 vs. G2T3 and G1T4 vs. G2T4 comparisons. G1T1 vs. G2T1 and G1T2 vs. G2T2 (Additional file 2: Fig. S7) had common and significant enrichment of ‘ribosome’ and ‘pyruvate metabolism’; G1T3 vs. G2T3 and G1T4 vs. G2T4 were enriched for ‘metabolic pathways’ and ‘cysteine and methionine metabolism’; for the comparisons of G1T2 vs. G2T2 and G1T4 vs. G2T4, it was related to ‘glutathione metabolism’ and ‘photosynthesis’; G1T2 vs. G2T2 and G1T3 vs. G2T3 displayed significant enrichment of ‘cutin, suberin and wax biosynthesis’; additionally, G1T1 vs. G2T1 and G1T4 vs. G2T4 shared ‘terpenoid backbone biosynthesis’, whereas for G1T1 vs. G2T1 and G1T3 vs. G2T3, ‘fatty acid biosynthesis’, ‘fatty acid metabolism’, ‘DNA replication’ and ‘biotin metabolism’ were expressed as commonly enriched pathways. In particular, the G1 vs. G2 arsenic stress (T1) group showed ‘purine metabolism’, ‘fatty acid degradation’, ‘phagosome’, ‘cyanoamino acid metabolism’, ‘benzoxazinoid biosynthesis’, ‘steroid biosynthesis’, ‘valine, leucine and isoleucine degradation’ and ‘glycosphingolipid biosynthesis - lacto and neolacto series’, to be significantly enriched. Furthermore, ‘flavonoid biosynthesis’, ‘pentose and glucuronate interconversions’, ‘carotenoid biosynthesis’, ‘tropane, piperidine and pyridine alkaloid biosynthesis’, and ‘tyrosine metabolism’ enriched for G1T2 vs. G2T2; whereas, the uniquely enriched pathway for G1T3 vs. G2T3 was ‘nitrogen metabolism’. ‘Zeatin biosynthesis’, ‘betalain biosynthesis’, and ‘glucosinolate biosynthesis’ were exclusively enriched for G1T4 vs. G2T4 (Fig. 3A and Additional file 2: Fig. S7).
Gene Ontological Classification of DEGs
In pursuit of functionally annotating the DEGs further, gene ontological categories were investigated including the domains: biological processes (BPs), cellular components (CCs), and molecular functions (MFs) (Fig. 3B and Additional file 2: Fig. S8). The enrichment results implied that ‘ADP binding’ was commonly enriched across the treatment wise paired comparisons among the molecular functions; however, ‘structural constituent of ribosome’, ‘structural molecule activity’, ‘rRNA binding’, ‘large ribosomal subunit rRNA binding’, ‘nutrient reservoir activity’, ‘O-acetyltransferase activity’, ‘small ribosomal subunit rRNA binding’, ‘D-threo-aldose 1-dehydrogenase activity’, ‘oxidoreductase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor’, were significantly and exclusively enriched for G1T1 vs. G2T1. The enrichment of gene ontological terms, i.e., ‘heme binding’, ‘tetrapyrrole binding’, ‘oxidoreductase activity’, ‘oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen’ and ‘monooxygenase activity’, were found to be mutual for G1T2 vs. G2T2, G1T3 vs. G2T3 and G1T4 vs. G2T4. The GO terms: ‘oxidoreductase activity, acting on peroxide as acceptor’ and ‘peroxidase activity’ were highlighted for G1T2 vs. G2T2 and G1T3 vs. G2T3, inclusively. While, G1T2 vs. G2T2 had an exclusive enrichment of ‘hydroquinone: oxygen oxidoreductase activity’ and ‘glucosyltransferase activity’; G1T3 vs. G2T3 had ‘microtubule binding’, ‘tubulin binding’; and G1T4 vs. G2T4 had ‘iron ion binding’, ‘pattern binding’, ‘polysaccharide binding’, as well as ‘carbohydrate binding’.
KEGG analysis of DEGs presented as scatter plot showing 20 most significantly enriched pathways (A), and relative comparative gene ontology (GO) in ZZY-1 under As10µM + S.i + P50µM in comparison to GD-6 (B). Each circle in the dot plot represents a pathway and the size of the circle represents the number of genes enriched in the pathway. The ordinate represents the pathway name and the abscissa represents the rich factor, which compares the ratio of genes annotated to a pathway among the DEGs to the ratio of genes annotated to that pathway among all genes. The color of the circle represents the q-value, where smaller values equate to more reliable enrichment significances of the DEGs in the pathway. GO terms were sorted based on q-values (< 0.05) and belonged to biological processes, cellular components, and molecular functions were shown in green, red and blue color, respectively
The G1T1 vs. G2T1 did not share any term in relation to cellular components or biological functions among the top 10 enriched GO terms with the other three comparisons, whereas ‘cell periphery’, ‘extracellular region’, ‘intrinsic component of membrane’, ‘integral component of membrane’ and ‘plasma membrane’ were the commonly CCs enriched for G1T2 vs. G2T2, G1T3 vs. G2T3, G1T4 vs. G2T4. For G1T2 vs. G2T2 and G1T3 vs. G2T3 comparisons, the shared terms were ‘apoplast’, ‘external encapsulating structure’, ‘cell-wall’, ‘anchored component of membrane’ and ‘anchored component of plasma membrane’. Exclusively, ‘membrane part’, ‘extracellular region part’, ‘extracellular space’, ‘lipid droplet’, and ‘membrane’ were enriched for G1T4 vs. G2T4. Among the most enriched and mutually expressed biological process term for G1T2 vs. G2T2, G1T3 vs. G2T3 and G1T4 vs. G2T4, was ‘defense response’. Additionally, genes governing ‘response to stress’ were commonly highlighted for G1T2 vs. G2T2 and G1T4 vs. G2T4. In case of G1T2 vs. G2T2 and G1T3 vs. G2T3, ‘cell-wall organization or biogenesis’ expressed mutually. Interestingly, the BPs up-regulated for the G1T4 vs. G2T4 involved ‘telomere maintenance via telomerase’, ‘telomere maintenance via telomere lengthening’, ‘inflorescence development’, ‘reproductive shoot system development’, ‘regulation of salicylic acid biosynthetic process’, ‘salicylic acid biosynthetic process’, ‘regulation of salicylic acid metabolic process’, along with ‘RNA-dependent DNA biosynthetic process’, ‘antibiotic biosynthetic process’ and ‘response to hydrogen peroxide’ terms, which are unique terms for this treatment-wise comparison.
Other groups displayed biological processes: ‘cell cycle’, ‘mitotic cell cycle’, ‘mitotic cell cycle process’, ‘cell cycle process’, ‘microtubule-based process’, ‘cell division’ exclusive for G1T3 vs. G2T3; ‘lignin catabolic process’, ‘phenylpropanoid catabolic process’, ‘cell-wall biogenesis’, ‘plant type cell-wall biogenesis’, ‘plant-type secondary cell-wall biogenesis’ as well as ‘external encapsulating structure organization’ for G1T2 vs. G2T2; however for G1T1 vs. G2T1, most BPs were unique and associated with ‘peptide biosynthetic process’, ‘amide biosynthetic process’, ‘peptide metabolic process’, ‘organonitrogen compound biosynthetic process’, ‘cellular nitrogen compound biosynthetic process’, ‘ribosome assembly’, ‘cellular amide metabolic process’, ‘cellular macromolecule biosynthetic process’ and ‘macromolecule biosynthetic process’ (Additional file 2: Fig. S8).
Transcription Factors
A number of transcription families were exclusively/ inclusively expressed, including WRKY, FAR1, NAC, bZIP, bHLH, B3, C2H2, ERF, M-type-MADS, MIKC-MADS, MYB- and -related. It is worth pointing out that a number of genes related to FAR1 TF (57) were identified for G1T1 vs. G2T1 (out of 452 genes), while out of 598 genes associated with TFs under G1T2 vs. G2T2, WRKY and FAR1 were on a par (59 genes each). Moreover, from 548 DEGs, 58 were identified belonging to FAR1 followed by B3 TF with 53 genes for G1T3 vs. G2T3. A comparatively smaller number of TF related gene expression was observed (Fig. 4) for G1T4 vs. G2T4 as out of 332 DEGs, most were associated with FAR1 (39) and B3 (33).
Heat Shock and Stress Responsive DEGs
A total of 215 stress-related DEGs were found being inclusively and exclusively expressed, including heat shock proteins and heat stress TFs, LEA proteins, zinc finger stress associated proteins, drought-, salt- and cold-tolerance related genes, PIPs and others. Regulation of stress related DEGs occurred in a treatment and genotype dependent manner (Additional file 2: Fig. S9). Heat shock protein 17.9 kDa (Os11g0244200) and OsHsp16.9B (Os01g0136200) were not expressed in the controls of both genotypes but the expression was observed across all treatments, more so in ZZY1; whereas, a heat shock protein 70 (Os11g0703900) did not express under the control (T5) and As10µM + S.i + P50µM treatment (T4) for both the genotypes. A remarkable increase in the expressions of OsHSP82A (Os08g0500700), OsHsp71.1 (Os03g0276500), OsHsp17.3 (Os03g0266900), OsctHSP70-1 (Os05g0460000), and OsMed37_1 (Os01g0840100) was witnessed across all treatments in both genotypes against their respective controls (Fig. 5). Exclusively under GD-6, OsHsp18.0 (Os01g0184100) was expressed, moreover a protein involved in negative salt stress tolerance (OCPI2; Os01g0615100) had higher expression under arsenic treatment (T1) in GD-6 than ZZY-1. Stress response related lipid transport protein OsACBP1 (Os08g0162800) expressions were found much higher in ZZY-1 genotype compared to GD-6 under all treatments. Meanwhile, the expression of OsCCD1 (Os06g0683400) was higher in controls of both genotypes, albeit the exposure of arsenic undoubtedly reduced the levels regardless of the amelioratory variables.
Detoxification Related DEGs
Inclusively and exclusively, around 167 genes were identified to be expressed differentially under this category. It was observed that a higher number of detoxification related genes were expressed under S.i colonization (T3) as compared to other treatments, reaching more than 6 log2-FC, with higher percentage of upregulated expression levels. Numerous peroxidases (POD), glutathione S-transferases, glutathione synthases, glutaredoxins were regulated by the treatments applied. The log2-FPKM values showed higher expression levels of Os1-CysPrxA (Os07g0638300) in ZZY-1 than GD-6 across all treatments against controls, which plainly depicts the genotypic difference in response to As exposure. Compared to ZZY-1, GD-6 displayed higher expression of OsAPx5 (Os12g0178200). Furthermore, OsGSTU35 (Os10g0368100) and OsGSTU50 (Os10g0530900) were also expressed highly under all treatments for GD-6 in comparison with its own control and ZZY-1 (Fig. 6 and Additional file 1: Table S5).
Transporters and Heavy Metal Related DEGs
A number of transport related DEGs (52) were identified in this study, including ABC (ATP-binding cassette) transporter related DEGs that belonged to Multidrug Resistance (MDR) and members of A, B, C, E, G and I families, phosphate (12), Zn (6), potassium (12), sulfate (9), copper (2), ammonium (2), nitrate (22), silicon (2), aluminium-induced malate transporters (7) and heavy metal related genes (112). Differences were observed in the expression of the genes among treatments and between genotypes (Fig. 7 and Additional file 1: Table S6). The expression of ABC transporters varied under different treatments for both genotypes; for genotype GD-6, these transporters had elevated expression levels under alone As10 stress than in combined treatments. However, OsABCG16 (Os06g0731200) and OsABCG21 (Os09g0401100) expressed highly under T4, OsABCG13 (Os05g0384600), OsABCG5 (Os03g0281900), OsABCG3 (Os01g0836600) under T3, whereas OsABCG2 (Os01g0615500) and OsABCB25 (Os03g0755900) under T2 (Additional file 2: Fig. S10 and Fig. S11).
Upon As addition, the expression of OsSPX-MFS1 (Os04g0573000) was increased comparatively more under As alone treatment in GD-6 than ZZY-1. Another phosphate transporter OsPT1 (Os03g0150600) was altered noticeably, experiencing downregulation for ZZY-1 under all treatments against its control, albeit an increased expression was observed under all treatments including control in GD-6. The transcript levels of OsPT14, which function in Pi transport in subcellular organelles, increased under T1, T3 and T4 treatments for ZZY-1 but in case of GD-6, this increase was observed mainly under T1 and T2. Ammonium transporter 1 (Os02g0620600) seemed to have a higher expression in ZZY-1 under T1 in comparison with its T2, T3 and T4 counterparts as well as GD-6 genotype. However, OsAMT2 (Os01g0831300) and OsDUR3 (Os10g0580400) expressions were positively highlighted under As stress for GD-6. Copper transporter 5 (Os05g042470) displayed expression under T2 and T4 for both genotypes, more prominently in GD-6, while OsCOPT6 (Os08g0455900) had higher expressions under T1 and T3 for GD-6, and T1 and T2 for ZZY-1. Aluminum-activated malate transporters were also found to express differentially, exhibiting high expressions generally in GD-6. Interestingly, a couple of differentially expressed genes associated with potassium transport were registered, with an elevated expression for both genotypes but under different treatments; SKC1 (Os01g0307500) and OsHAK4 (Os08g0466200) portrayed a higher expression in the controls of ZZY-1 than GD-6. Although, OsHAK5 (Os01g0930400), OsHAK7 (Os07g0669700), OsHAK13 (Os06g0671000) and OsHAK21 (Os03g0576200) were expressed more under T1-T3 for GD-6, OsHAK8 (Os03g0337500) and OsHAK22 (Os07g0102100) expressed for ZZY-1 under the same three treatments. Furthermore, OsAKT1 (Os01g0648000) was expressed the most under T4 for ZZY-1. Among the NO3– transporters, similar to nitrate-induced NOI protein (Os11g0226400), expressed only for GD-6, whereas silicon transporters LSI1 (Os02g0745100) and LSI2 (Os03g0107300), being highly regulated under control conditions, mostly downregulated under T1 and T3, though T4 in case of both genotypes and T2 in case of GD-6 only brought about higher expressions of these transporters. On the contrary, LSI6 (Os06g0228200) significantly upregulated under all treatments for both genotypes, more profoundly for ZZY-1 (Additional file 1: Table S6).
DEGs Associated with Phytohormones
Owing to their crucial roles in abiotic stress regulation, a mutual comparison of DEGs related to auxin (Additional file 2: Fig. S12), ethylene (Additional file 2: Fig. S13), abscisic acid, jasmonic acid, salicylic acid (Fig. 8), cytokinin and gibberellic acid (Fig. 9), was performed. There were 255 inclusively and exclusively expressed DEGs, and a very large number of auxin related genes were upregulated for ZZY-1 regardless of the treatment applied, mostly being under T2, but DEGs such as OsSAUR37 (Os09g0437100), similar to SAR1 (suppressor of auxin resistance1; Os02g0142950) and OsSAUR57 (Os12g0609600), were prominently expressed under T4 for GD-6 genotype. Moreover, similar to auxin-binding protein, OsABP20 (Os04g0288100) expressed exclusively under GD-6. In terms of DEGs involved in ethylene biosynthesis or signaling, GD-6 exhibited innately higher expression as compared to ZZY-1, though an increase in expression was found under T1 and T3 for ZZY-1 as well, mostly under T2 and T3. However, ACO4 (Os11g0186900) and ACO5 (Os05g0149400) were highly expressed for this genotype under T4. Also, similar to ethylene-responsive transcriptional coactivator (Os06g0592500) gene was dramatically increased in its expression with the onset of As added treatments, where T1 and T3 enforced higher expressions in GD-6 compared to its T2 and T4 counterparts but also against ZZY-1. The expressions of brassinosteroids associated DEGs were comparatively higher in ZZY-1, specially under T1, T2 and T3; moreover, brassinosteroid insensitive 1-associated receptor kinase 1 (Os08g0255500), expressed exclusively for GD-6. Differentially expressed genes involved in gibberellin regulation were significantly expressed for ZZY-1 genotype, while GD-6 had a fair share of upregulated cytokinin related DEGs, quite noticeable under T2, T3 and T4. A relatively higher number of ABA and SA related DEGs were expressed significantly and positively for ZZY-1 compared to GD-6 genotype, which presented an affinity towards jasmonic acid associated DEGs’ upregulation under its As stressed treatments as opposed to ZZY-1. The comparisons revealed that G1T3 vs. G2T3 had a very strong auxin signaling response (Figs. 8 and 9, Additional file 1: Table S7).
Hierarchical heat map showing the expression profile of inclusively expressed genes related to the biosynthesis and signaling of JA, SA and ABA in response to As alone (T1), As + P (T2), As + S.i (T3) and As + S.i + P (T4) treatments along with control (T5), in the two rice genotypes ZZY-1 (G1) and GD-6 (G2)
Hierarchical heat map illustrating the expression profile of inclusively expressed genes related to gibberellin, cytokinin and brassinosteroid biosynthesis and signaling in response to As alone (T1), As + P (T2), As + S.i (T3) and As + S.i + P (T4) treatments along with control (T5), in the two rice genotypes ZZY-1 (G1) and GD-6 (G2)
Additionally, GDSL esterase/lipase (GELPs) associated DEGs, such as OsGELP93 (Os07g0642200), OsGELP56 (Os04g0507700), OsGELP67 (Os05g0419800) and OsGELP50 (Os03g0365900) were highly expressed under T4 for GD-6. In regard to ZZY-1, the expressions were positively aggregated under T2 followed by T3, T1 and then T4. Certain cell-wall biosynthesis related genes (expansins) and different channels were expressed differentially and their expression profiles have been presented as heat map (Fig. 10 and Additional file 2: Fig. S14).
Discussion
The Efficient Detoxification and Exclusion of Heavy Metals
Arsenic toxicity can adversely affect the growth of a living organism by the accumulation of lipid peroxidation products and damaging the membrane integrity, however in the presence of such xenobiotics, plants minimize this damage by increasing the cellular activity of detoxification enzymes (Schroder et al. 2009). Brinke et al. (2018), in their transcriptome study on rice roots after exposure to AsIII-tainted sediments, found that a GO term “glutathione biosynthesis” was prominent. According to Verbruggen et al. (2009), glutathione biosynthesis is induced under heavy metal stress, providing a substrate for phytochelatins’ formation that constitutively participate in detoxifying several metal(loid)s. Transcripts related to ROS-detoxification were distinctly upregulated under As + P and As + S.i of ZZY-1 as well as T1 and T3 of GD-6. The conversion of AsV to AsIII involves oxidation of glutathione (GSH) to GSSG, therefore it is not uncommon to observe glutathione-S-transferase, glutathione reductase (GR), glutaredoxin (Grx) and thioredoxin (Trx) associated transcripts being upregulated in rice roots upon AsV exposure (Huang et al. 2012). Arsenic stress induces lignification (Huang et al. 2012), which solidifies the cell-wall and restricts root growth and the production of ROS contribute to this AsV-induced rigidity of the cell-wall. Genotype ZZY-1 presented a phenotype that, under the adverse As-laden growth conditions, managed to produce new root growth and adopted better avoidance strategies than GD-6. The addition of P with As led to suberization of roots internal tissues. In Arabidopsis roots, ABCG2, ABCG6 and ABCG20 transporters partake in endodermal suberin lamellae formation (Yadav et al. 2014). Moreover, OsABCG5 in rice as well as StABCG1 in potato tuber and its periderm have also been identified and designated as suberin monomer transporters (Chen et al. 2022). Moreover, the exudation of the HM-chelating organic acids, ‘malate’ and ‘citrate’ into the rhizosphere, by Al-activated malate transporter (ALMT) and multidrug and toxin extrusion (MATE) transporters, respectively, serve as a HM-resistance mechanism in plants (Liang et al. 2013; Liu & Li 2022). The transporter-facilitated secretion of organic acids not only detoxifies/extrudes toxic metal ions from roots but also improves P and Fe availability (Delhaize et al. 2012; Kobayashi and Nishizawa 2012). Transmembrane proteins, like P-type ATPases, function as pumps across the plasma or intracellular membranes, coupling ATP hydrolysis to the efflux of a cation out of the cytosol (Nuria 2021). In this study, the expression of P-type ATPase, K/Mg/Cd/Cu/Zn/Na/Ca/Na/H-transporter family protein (Os11g0485200), was found to be upregulated under all treatments of ZZY-1, showing the capability of this genotype as a strong effluxer.
Preferential Growth of Lateral Roots and Potential Modification at Chromosome Level
Roots provide anchorage to the plants and the finer lateral roots in particular, explore the surrounding medium and forage for nutrients and water (Gutjahr et al. 2015). Given experimental conditions where AsV was abundant, the As tolerant plant exhibited a trait of cellular proliferation, producing lateral roots, which helped them to avoid taking up AsV, as reported in Pterris vittata (Wang et al. 2011), where new roots tended to accumulate and transport less AsV than primary roots, which on the contrary, uptake more AsIII. Generally, AMF colonization is distributed asymmetrically within the root system (Gutjahr and Paszkowski 2013). This becomes evident in rice roots that consist of three major root types, i.e., crown roots that are weakly colonized, large lateral roots that get strongly colonized, and fine lateral roots (FLRs), which are preferably not colonized (Gutjahr et al. 2009; Rebouillat et al. 2009). In current study, high number of lateral root production by the As tolerant genotype facilitated the colonization of S. indica, which corroborates a separate study conducted on rice employing another AMF Rhizophagus irregularis (Gutjahr et al. 2015). The study conducted by Choi et al. (2021), investigated the hereditary foundation of the naturally varying telomere length in A. thaliana, maize and rice and found a negative correlation between telomere length and flowering time. Both phosphorus nutrition and S. indica are reported to facilitate the transitional phase of vegetative growth towards maturity. It is not surprising that the combination of both is actually helping rice plants to reduce the flowering time by regulating the telomere lengthening and at the same time avoiding As-induced programmed cell death, as in plants the longer the length of telomere, the shorter the time taken up by plant to reach maturity, and vice versa. This knowledge of the two tested ameliorative treatments could be used in scientific development integrating different perspectives. Another aspect that differentiated the tolerant genotype from the susceptible one was the ribosomal damage repair under As alone treatment (G1T1 vs. G2T1), which if not done promptly and appropriately, would cause cellular functions to cease and programmed cell death.
Phosphorus and S. indica Co-Treatment Induced Reprogramming
In plants, GDSL-type esterase/lipase enzymes (GELP) play key functions encompassing anther and pollen development, along with responses to biotic and abiotic stresses (Chepyshko et al. 2012; Cenci et al. 2022). A study carried out by Riemann et al. (2007), explained their response in association with light and jasmonate signaling; moreover, OsGELP110/OsGELP115 were characterized as contributors to the formation of rice exine, whereas the double knockout (osgelp110/osgelp115) experienced male sterility (Zhang et al. 2020). Yuan et al. (2020) worked with an ER-localized GELP ‘OsGELP34’, and observed its essentiality in male reproduction and pollen development. Although, control plants of both genotypes displayed a certain degree of expression, arsenic alone and in combination with other treatments brought about a surge in genes related to GELP category (Additional file 2: Fig. S14), which represents the strategy of producing reproductive organs for the continuation of species.
The genotypes differed in their phytohormonal response as a large number of auxin responsive and biosynthesis related genes were upregulated under all the treatments of ZZY-1 on one hand (Additional file 2: Fig. S12), while on the other hand GD-6 predominantly experienced an upregulation of genes related to ethylene (Additional file 2: Fig. S13). We know from the results that ZZY-1 overall performed better than GD-6, which could be explained in a way that this genotype presented more plasticity due to higher auxin levels that inhibit elongation of primary root and promote formation of lateral root (Alarcón et al. 2019). It is interesting that ethylene and auxin work in synchronicity (Vanstraelen and Benková 2012). Ethylene acts as a central negative regulator of auxin levels in cells, restricting elongation, and this all happens in the root epidermis (Vaseva et al. 2018), which being in continuous contact with the environment surrounding it, encounters external cues, such as heavy metal toxicity.
In terms of nutrient transporters, phosphate-starvation related gene OsSPX-MFS1 (Wang et al. 2012; He et al. 2021), was highly upregulated in GD-6 under As stress, which suggests that the competition of the As vs. P had much severe impact on this genotype; additionally, the control (T5), As + S.i + P (T4) and As + P (T2) treatments of the same genotype utilized OsPT1, whereas under As + S.i (T3) another P transporter, OsPT2 was employed, which shows that the symbiotic association modified the phosphate uptake mechanism. Genotype ZZY-1, on the contrary, made use of only OsPT14 under T1, T3 and T4. It is worthy to mention here that Shi et al. (2013), proposed OsPHT2;1 (OsPT14) as a marker gene for screening high phosphorus efficiency genotypes, whose overexpression could increase foliar P concentrations and the biomass of rice plants. Ammonium is a preferred form of nitrogen for uptake in rice plants (Konishi and Ma 2021), and arsenic stress is known to adversely affect nitrogen metabolism, therefore it was not surprising to observe an upregulation of ammonium transporters (OsAMT1 in case of ZZY-1 and OsAMT2 for GD-6 along with a urea transporter OsDUR3), whose expression mellowed down under T4 and become somewhat closer to control values, showing that the pressure of altered ammonium nutrition was adjusted accordingly. Also, for ZZY-1 plants, the expression of the high-affinity nitrate transporter OsNAR2.1 was upregulated under As + S.i and As + S.i + P, which is known to share a role in nitrate absorption and translocation in rice, providing plant with drought tolerance under water deficit conditions (Chen et al. 2019). Pertaining to As stress tolerance, GD-6 also had a few tricks up its sleeves, as it displayed a higher expression of OsUGE1 (UDP-glucose 4-epimerase 1; Os05g0595100) under As alone stress, which is putatively involved in cell-wall carbohydrate partitioning amidst conditions where nitrogen is deficient (Guevara et al. 2014). Also, the gene expression of OsNiR, which encodes nitrite reductase and is indispensable for N assimilation in rice (Yu et al. 2021), was upregulated in GD-6 under all treatments, more prominently under As + P. Sulphur pool is required for the synthesis of many enzymes and antioxidant compounds, such as GSH and phytochelatins (Clemens and Ma 2016); differential upregulation of OsSultr1 and OsSultr3 represent the genotypic variation in regard to dealing with the As-induced disruption in S-metabolism. Abercrombie et al. (2008) demonstrated increased transcript levels of genes encoding chloroplastic and a cytosolic Cu/ZnSOD upon AsV exposure, which could explain the higher activity of copper and zinc transporters in this study, particularly for GD-6. Both genotypes responded in a certain way when it came to K+ nutrition, as numerous related genes were upregulated in both genotypes to cope with As-toxicity. Specific to ZZY-1, the expression of OsAKT1 was induced higher under T4 than other treatments, which indicates the role of K+ in close connection with P against As toxicity (Ródenas et al. 2019). Having a closer look at the differential expression of metal transporters, involvement of YSL family genes, IRT2 and FRDL1 suggest the involvement of Fe nutrition; also, the two NRAMPs (natural resistance-associated macrophage proteins), i.e., OsNramp1 and OsNramp5 show transport activity for Fe2+, albeit their total Fe2+ uptake contribution is comparatively small (Takahashi et al. 2011; Sasaki et al. 2012). Maintenance of nutritional homeostasis is critical to plant’s growth as disruption/deficiency in one could lead to excess in the other or vice versa (Xie et al. 2019). Expressions of OsENA1 under all As stressed treatments for ZZY-1 indicated the employment of salt stress response, as this gene in barely expresses under normal growth conditions but gets intensified under elevated salt concentrations and alkaline pH (Ruiz et al. 2003).
Construing Stress Response and Phytohormonal Tuning
To produce functional proteins along with removing misfolded/aggregated harmful proteins, a system of protein quality control operates inside the cells that involves certain molecular chaperones, folding catalysts and proteases, which only assist and do not become a component of final products (Ellis 1993). These chaperons are categorized by their molecular-sizes, i.e., Hsp100s, Hsp90s, Hsp70s, Hsp60s, Hsp40s, and small heat-shock proteins, where each member assists in folding, refolding, oligomeric assembly, translocation, and degradation of a specific substrate proteins (Hartl et al. 2011). A large number of HSPs and related TFs were upregulated for ZZY-1 and GD-6 (under T1, T2, T3), however the excessive protein repair was reduced under T4, much more in ZZY-1. It is worth noting that Hsp16.9 A was upregulated mainly for ZZY-1, and higher expression of OsHsp18 were witnessed in GD-6. These small heat shock proteins actually contribute positively to certain biotic/abiotic defense responses in rice plants (Kuang et al. 2017). According to a study conducted by Huang et al. (2012) almost 40% of the OsHsf TFs were observed to be induced in rice during AsV stress. In line with the stress responsive transcripts, gene expression of PP2C1 and NAC58 reduced in ZZY-1 compared to its control counterparts. Noteworthily, PP2Cs function as negative regulators of the ABA signaling network, which mediate stomatal closure (Jung et al. 2020), therefore reducing the expression of PP2C1 could trigger ABA related signaling cascade. Additionally, the NAC58 is a transcription activator that acts as a positive regulator of the jasmonate (JA) pathway, mediating leaf senescence and may directly regulate JA biosynthesis genes, such as LOX2, AOC, AOS2, AOC1 and OPR7 (Zhou et al. 2013). The obtained results show that most genes related to JA biosynthesis were not induced much in ZZY-1 than GD-6, rather transcript levels of salicylic acid and gibberellin biosynthesis genes were upregulated in ZZY-1 (Fig. 8). Under chilling stress, the combined application of the phytohormone salicylic acid and H2O2 raised the GA/ABA ratio, accelerating seed germination of maize (Li et al. 2017). Moreover, both genotypes under T4 showed lesser inclination towards ethylene responsive and biosynthesis genes expression (Additional file 2: Fig. S13), which suggests the differential fine tuning of phytohormonal operational mechanism, as an increase in ethylene response factors is a characteristic of arsenic toxicity in rice (Huang et al. 2012). Arsenic exposure tends to downregulate the expression of cytokinin signaling genes, such as OsCRL4 and OsRR20 (Huang et al. 2012), but genotype ZZY-1 displayed an upregulatory trend in the expression of OsCRL1 under all treatments, as well as of RR9 and RR10 specifically under As + S.i + P (Fig. 9). Cytokinins are enriched in the shoot apical meristem and immature leaves, regulating nutrient sink and source activity and promote crop productivity by the activation of inflorescence meristems (Choi et al. 2011). Moreover, CK and BR interaction reportedly promotes ovule initiation and increases seed number per silique in case of Arabidopsis (Zu et al. 2022). Another gene worthy of mention among the stress responsive genes is OsCIA, which is present under normal growth conditions in panicles, leaves, and seedling tissues, however cold stress also induces its expression in anthers (Imin et al. 2006); the expressions were prominent under T1, T2 and T3 of ZZY-1, which implies that arsenic stress could also be sensed by the plants as chilling stress. ATP-binding cassette (ABC) transporters transform the energy gained from hydrolysis of ATP into trans-bilayer substrate movement either into or outside the cytoplasm (Locher 2009). Genotype GD-6, particularly under As alone stress expressed a number of ABC transporters, which indicates that in an attempt to exclude or sequester the accumulated As into the vacuoles, this genotype spent a lot of energy.
Shift in Number of Genes Under Certain Transcriptional Factors
Among a number of TFs, FAR1, bHLH, WRKY, B3 and ERF were quite prominent, though the number of genes presented under them changed remarkably with added ameliorative treatments (Fig. 4). FAR-RED IMPAIRED RESPONSE1 (FAR1) partakes in plant growth and development processes, such as regulating phytochrome-signaling processes, branching/tillering, shade avoidance, and flowering time in Arabidopsis thaliana and rice (Genoud et al. 2008); in addition, it is reported by Dai et al. (2022), to positively regulate the responses to salt and temperature stress in Eucalyptus grandis. Similarly, ‘B3-domain’ family of transcription factors have crucial and diverse functions, such as seed growth, development and stress. The mutations of B3’s master regulators, i.e., abi3, lec2, and fus3, revealed the association of this TF in embryonic Fe distribution and FERRITIN gene expression (Grant-Grant et al. 2022).
Schematic illustration of differentially expressed transcripts associated with major KEGG and GO annotated networks expressed during arsenic stress along with phosphorus fortification and S. indica symbiosis. These transcripts have been significantly detected when ZZY-1 under the above-mentioned treatment was compared with GD-6.
Furthermore, a huge number of expansins were detected to be upregulated, mainly in ZZY-1 under T1, T2 and T3. The expansins mediate cell-wall loosening, and thereby, are responsible for apical growth enhancement (Choi et al. 2003), which might explain the plasticity of ZZY-1 genotype. The strategies adopted by plants against As uptake, its translocation and accumulation vary phenomenally, however prophylactic measures must be taken to ensure the safety of humans and the collective environment. Stress responses, like senescence and PCD, could be reoriented towards inflorescence/early maturity, which gives the plant a chance to reproduce under unpropitious conditions. Schematic diagram of differentially expressed transcripts associated with major KEGG and GO annotated networks expressed during arsenic stress along with phosphorus fortification and S. indica symbiosis is illustrated in Fig. 10.
Conclusions
Arsenic toxicity to the plants can be reduced by using self-renewing and cost-effective methods, such as AMF incorporation, which is less intrusive and potentially delays phosphate starvation response with adequate P supply. This study showed that the As + S.i + P interactive treatment mellowed down the chaotic increase in heat shock proteins, and potentially induced early flowering genes by lengthening the telomere. Keeping the above information into consideration, genes stimulating endogenous BR (CycD3) and their cross-talk with CK (Os11g0143300, Os12g0139400, Os01g0952500 and Os01g0197700), as well as SA (Os11g0663100, Os03g0778000 and Os11g0669100), metal effluxer (P-type ATPase, K/Mg/Cd/Cu/Zn/Na/Ca/Na/H-transporter family protein) and detoxification related genes (glutathione-S-transferases), elevation of K+ transporter AKT1, Fe+ transporters (OsYSL2 and OsYSL15), and phosphate transporter PHT14, greatly contributed to As tolerance acquisition. Supply of P remarkably changed the response of plants colonized with S. indica, which raises questions regarding the efficiency of this AMF and the untapped potential that could be exploited by combining this tool with other HM remediation techniques.
Data Availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].
References
Abercrombie JM, Halfhill MD, Ranjan P, Rao MR, Saxton AM, Yuan JS, Stewart CN (2008) Transcriptional responses of Arabidopsis thaliana plants to As (V) stress. BMC Plant Biol 8:87. https://doi.org/10.1186/1471-2229-8-87
Adil MF, Sehar S, Chen G, Jilani G, Chen ZH, Chaudhry AN, Shamsi IH (2020) Cadmium-zinc cross- talk delineates toxicity tolerance in rice via differential genes expression and physiological/ultrastructural adjustments. Ecotoxicol Environ Saf 190:110076. https://doi.org/10.1016/j.ecoenv.2019.110076
Alarcón MV, Salguero J, Lloret PG (2019) Auxin modulated initiation of lateral roots is linked to pericycle cell length in maize. Front Plant Sci 10:11. https://doi.org/10.3389/fpls.2019.00011
Bajaj R, Huang Y, Gebrechristos S, Mikolajczyk B, Brown H, Prasad R, Varma A, Bushley KE (2018) Transcriptional responses of soybean roots to colonization with the root endophytic fungus Piriformospora indica reveals altered phenylpropanoid and secondary metabolism. Sci Rep 8:10227. https://doi.org/10.1038/s41598-018-26809-3
Boorboori MR, Zhang HY (2022) The role of Serendipita indica (Piriformospora indica) in improving plant resistance to drought and salinity stresses. Biology 11(7):952. https://doi.org/10.3390/biology11070952
Brinke A, Reifferscheid G, Klein R, Feiler U, Buchinger S (2018) Transcriptional changes measured in rice roots after exposure to arsenite-contaminated sediments. Environ Sci Pollut Res Int 25(3):2707–2717. https://doi.org/10.1007/s11356-017-0515-z
Cao Y, Sun D, Ai H, Mei H, Liu X, Sun S, Xu G, Liu Y, Ma LQ (2017) Knocking out OsPT4 gene decreases arsenate uptake by rice plants and inorganic arsenic accumulation in rice grains. Environ Sci Technol 51:12131–12138. https://doi.org/10.1021/acs.est.7b03028
Cenci A, Concepción-Hernández M, Guignon V, Angenon G, Rouard M (2022) Genome-wide classification and phylogenetic analyses of the GDSL-type esterase/lipase (GELP) family in flowering plants. Int J Mol Sci 23:12114. https://doi.org/10.3390/ijms232012114
Chadha N, Mishra M, Rajpal K, Bajaj R (2015) An ecological role of fungal endophytes to ameliorate plants under biotic stress. Arch Microbiol 197:869–881. https://doi.org/10.1007/s00203-015-1130-3
Chakrabarty D, Trivedi PK, Misra P, Tiwari M, Shri M, Shukla D, Kumar S, Rai A, Pandey A, Nigam D, Tripathi RD, Tuli R (2009) Comparative transcriptome analysis of arsenate and arsenite stresses in rice seedlings. Chemosphere 74:688–702. https://doi.org/10.1016/j.chemosphere.2008.09.082
Chen JG, Qi TT, Hu Z, Fan XR, Zhu LL, Iqbal MF, Yin XM, Xu GH, Fan XR (2019) OsNAR2.1 positively regulates drought tolerance and grain yield under drought stress conditions in rice. Front Plant Sci 10:197. https://doi.org/10.3389/fpls.2019.00197
Chen AL, Liu T, Wang Z, Chen XP (2022) Plant root suberin: a layer of defence against biotic and abiotic stresses. Front Plant Sci 13:1056008. https://doi.org/10.3389/fpls.2022.1056008
Chepyshko H, Lai CP, Huang LM, Liu JH, Shaw JF (2012) Multifunctionality and diversity of GDSL esterase/lipase gene family in rice (Oryza sativa L. japonica) genome: new insights from bioinformatics analysis. BMC Genom 13:309. https://doi.org/10.1186/1471-2164-13-309
Choi D, Lee Y, Cho HT, Kende H (2003) Regulation of expansin gene expression affects growth and development in transgenic rice plants. Plant Cell 15(6):1386–1398. https://doi.org/10.1105/tpc.011965
Choi JY, Choi D, Lee S, Ryu CM, Hwang I (2011) Cytokinins and plant immunity: old foes or new friends? Trends Plant Sci 16(7):388–394. https://doi.org/10.1016/j.tplants.2011.03.003
Choi JY, Abdulkina LR, Yin, Chastukhina IB, Lovell JT, Agabekian IA, Young PG, Razzaque S, Shippen DE, Juenger TE, Shakirov EV, Purugganan MD (2021) Natural variation in plant telomere length is associated with flowering time. Plant Cell 33(4):1118–1134. https://doi.org/10.1093/plcell/koab022
Clemens S, Ma JF (2016) Toxic heavy metal and metalloid accumulation in crop plants and foods. Annu Rev Plant Biol 67:489–512. https://doi.org/10.1146/annurev-arplant-043015-112301
Dai J, Sun J, Peng W, Liao W, Zhou Y, Zhou X-R, Qin Y, Cheng Y, Cao S (2022) FAR1/FHY3 transcription factors positively regulate the salt and temperature stress responses in Eucalyptus grandis. Front Plant Sci 13:883654. https://doi.org/10.3389/fpls.2022.883654
Dangleben NL, Skibola CF, Smith MT (2013) Arsenic immunotoxicity: a review. Environ Health 12(1):73. https://doi.org/10.1186/1476-069X-12-73
Das J, Sarkar P (2018) Remediation of arsenic in mung bean (Vigna radiata) with growth enhancement by unique arsenic-resistant bacterium Acinetobacter lwofi. Sci Total Environ 624:1106–1118. https://doi.org/10.1016/j.scitotenv.2017.12.157
Delhaize E, Ma J, Ryan PR (2012) Transcriptional regulation of aluminium tolerance genes. Trends Plant Sci 17(6):341–348. https://doi.org/10.1016/j.tplants.2012.02.008
Ellis RJ (1993) The general concept of molecular chaperones. Philos Trans R Soc Lond B Biol Sci 339:257–261. https://doi.org/10.1098/rstb.1993.0023
Franken P (2012) The plant strengthening root endophyte Piriformospora indica: potential application and the biology behind. Appl Microbiol Biotechnol 96:1455–1464. https://doi.org/10.1007/s00253-012-4506-1
Genoud T, Schweizer F, Tscheuschler A, Debrieux D, Casal JJ, Schafer E, Hiltbrunner A, Fankhauser C (2008) FHY1 mediates nuclear import of the light-activated phytochrome A photoreceptor. PLoS Genet 4(8):e1000143. https://doi.org/10.1371/journal.pgen.1000143
Ghorbani A, Tafteh M, Roudbari N, Pishkar L, Zhang W, Wu C (2021) Piriformospora indica augments arsenic tolerance in rice (Oryza sativa) by immobilizing arsenic in roots and improving iron translocation to shoots. Ecotoxicol Environ Saf 209:111793. https://doi.org/10.1016/j.ecoenv.2020.111793
Grant-Grant S, Schaffhauser M, Baeza-Gonzalez P, Gao F, Conéjéro G, Vidal EA, Gaymard F, Dubos C, Curie C, Roschzttardtz H (2022) B3 transcription factors determine iron distribution and FERRITIN gene expression in embryo but do not control total seed iron content. Front Plant Sci 13:870078. https://doi.org/10.3389/fpls.2022.870078
Guevara DR, El-Kereamy A, Yaish MW, Yong MB, Rothstein SJ (2014) Functional characterization of the rice UDP-glucose 4-epimerase 1, OsUGE1: a potential role in cell wall carbohydrate partitioning during limiting nitrogen conditions. PLoS ONE 9(5):e96158. https://doi.org/10.1371/journal.pone.0096158
Gutjahr C, Paszkowski U (2013) Multiple control levels of root system remodeling in arbuscular mycorrhizal symbiosis. Front Plant Sci 4:204. https://doi.org/10.3389/fpls.2013.00204
Gutjahr C, Casieri L, Paszkowski U (2009) Glomus intraradices induces changes in root system architecture of rice independently of common symbiosis signaling. New Phytol 182(4):829–837. https://doi.org/10.1111/j.1469-8137.2009.02839.x
Gutjahr C, Sawers RJH, Marti G, Andrés-Hernández L, Yang SY, Casieri L, Angliker H, Oakeley EJ, Wolfender JL, Abreu-Goodger C et al (2015) Transcriptome diversity among rice root types during asymbiosis and interaction with arbuscular mycorrhizal fungi. Proc Natl Acad Sci USA112:6754–6759. https://doi.org/10.1073/pnas.1504142112
Hansen HR, Raab A, Price AH, Duan G, Zhu YG, Norton GJ, Feldmann J, Meharg AA (2011) Identification of tetramethylarsonium in rice grains with elevated arsenic content. J Environ Monit 13(1):32–34. https://doi.org/10.1039/C0EM00460J
Hartl FU, Bracher A, Hayer-Hartl M (2011) Molecular chaperones in protein folding and proteostasis. Nature 475:324–332. https://doi.org/10.1038/nature10317
He QJ, Lu H, Guo HX, Wang Y, Zhao P, Li Y, Wang F, Xu JM, Mo XR, Mao CZ (2021) OsbHLH6 interacts with OsSPX4 and regulates the phosphate starvation response in rice. Plant J 105:649–667. https://doi.org/10.1111/tpj.15061
Hill TW, Käfer E (2001) Improved protocols for aspergillus minimal medium: trace element and minimal medium salt stock solutions. Fungal Genet Rep 48:20–21. https://doi.org/10.4148/1941-4765.1173
Huang TL, Nguyen QTT, Fu SF, Lin CY, Chen YC, Huang HJ (2012) Transcriptomic changes and signalling pathways induced by arsenic stress in rice roots. Plant Mol Biol 80(6):587–608. https://doi.org/10.1007/s11103-012-9969-z
Imin N, Kerim T, Weinman JJ, Rolfe BG (2006) Low temperature treatment at the young microspore stage induces protein changes in rice anthers. Mol Cell Proteomics 5(2):274–292. https://doi.org/10.1074/mcp.M500242-MCP200
Jung C, Nguyen NH, Cheong JJ (2020) Transcriptional regulation of protein phosphatase 2 C genes to modulate abscisic acid signaling. Int J Mol Sci 21(24):9517. https://doi.org/10.3390/ijms21249517
Kamiya T, Islam R, Duan G, Uraguchi S, Fujiwara T (2013) Phosphate deficiency signaling pathway is a target of arsenate and phosphate transporter, OsPT1, is involved in as accumulation in shoots of rice. Soil Sci Plant Nutr 59:580–590. https://doi.org/10.1080/00380768.2013.804390
Kobayashi T, Nishizawa NK (2012) Iron uptake, translocation, and regulation in higher plants. Annu Rev Plant Biol 63:131–152. https://doi.org/10.1146/annurev-arplant-042811-105522
Konishi N, Ma JF (2021) Three polarly localized ammonium transporter 1 members are cooperatively responsible for ammonium uptake in rice under low ammonium condition. New Phytol 232:1778–1792. https://doi.org/10.1111/nph.17679
Kuang J, Liu JZ, Mei J, Wang CC, Hu HT, Zhang YJ, Sun MH, Ning X, Xiao LT, Yang L (2017) A class II small heat shock protein OsHsp18.0 plays positive roles in both biotic and abiotic defense responses in rice. Sci Rep 7:11333. https://doi.org/10.1038/s41598-017-11882-x
Kushwaha AS, Thakur RS, Patel DK, Kumar M (2022) Impact of arsenic on phosphate solubilization, acquisition and poly-phosphate accumulation in endophytic fungus Serendipita indica. Microbiol Res 259:127014. https://doi.org/10.1016/j.micres.2022.127014
Li Z, Xu J, Gao Y, Wang C, Guo G, Luo Y, Huang YT, Hu WM, Sheteiwy MS, Guan YJ, Hu J (2017) The synergistic priming effect of exogenous salicylic acid and H2O2 on chilling tolerance enhancement during maize (Zea mays L.) seed germination. Front Plant Sci 8:1153. https://doi.org/10.3389/fpls.2017.01153
Liang CY, Piñeros MA, Tian J, Yao ZF, Sun LL, Liu JP, Shaff J, Coluccio A, Kochian LV, Liao H (2013) Low pH, aluminum, and phosphorus coordinately regulate malate exudation through GmALMT1 to improve soybean adaptation to acid soils. Plant Physiol 161(3):1347–1361. https://doi.org/10.1104/pp.112.208934
Lin CY, Wang J, Cheng HG, Ouyang W (2015) Arsenic profile distribution of the wetland argialbolls in the Sanjiang plain of northeastern China. Sci Rep 5:10766. https://doi.org/10.1038/srep10766
Liu Q, Li Y (2022) Nutrient and non-nutrient factors associated with the arsenic uptake and buildup in rice: a review. J Soil Sci Plant Nutr 22:4798–4815. https://doi.org/10.1007/s42729-022-00961-2
Locher KP (2009) Structure and mechanism of ATP-binding cassette transporters. Philos Trans R Soc Lond B Biol Sci 364(1514):239–245. https://doi.org/10.1098/rstb.2008.0125
Mohd S, Shukla J, Kushwaha AS, Mandrah K (2017) Endophytic fungi Piriformospora indica mediated protection of host from arsenic toxicity. Front Microbiol 8:1–14. https://doi.org/10.3389/fmicb.2017.00754
Norton GJ, Lou-Hing DE, Meharg AA, Price AH (2008) Rice-arsenate interactions in hydroponics: whole genome transcriptional analysis. J Exp Bot 59:2267–2276. https://doi.org/10.1093/jxb/ern097
Norton GJ, Pinson SRM, Alexander J, McKay S, Hansen H, Duan GL, Islam MR et al (2012) Variation in grain arsenic assessed in a diverse panel of rice (Oryza sativa) grown in multiple sites. New Phytol 193:650–664. https://doi.org/10.1111/j.1469-8137.2011.03983.x
Nuria F (2021) Membrane transporters, an overview of the arbuscular mycorrhizal fungal transportome. Encyclopedia of Mycology, Elsevier, 1:44–53. https://doi.org/10.1016/B978-0-12-819990-9.00021-4
Oktalira FT, May TW, Dearnaley JDW, Linde CC (2021) Seven new Serendipita species associated with australian terrestrial orchids. Mycologia 113(5):968–987. https://doi.org/10.1080/00275514.2021.1919848
Rebouillat J, Dievart A, Verdeil JL, Escoute J, Giese G, Breitler JC, Gantet P, Espeout S, Guiderdoni E, Perin C (2009) Molecular genetics of rice root development. Rice 2:15–34. https://doi.org/10.1007/s12284-008-9016-5
Riedmann C, Ma Y, Melikishvili M, Godfrey SG, Zhang Z, Chen KC, Rouchka EC, Fondufe-Mittendorf YN (2015) Inorganic Arsenic-induced cellular transformation is coupled with genome wide changes in chromatin structure, transcriptome and splicing patterns. BMC Genom 16:212. https://doi.org/10.1186/s12864-015-1295-9
Riemann M, Gutjahr C, Korte A, Riemann M, Danger B, Muramatsu T, Bayer U, Waller F, Furuya M, Nick P (2007) GER1, a GDSL motif-encoding gene from rice is a novel early light- and jasmonate-induced gene. Plant Biol 9(1):32–40. https://doi.org/10.1055/s-2006-924561
Ródenas R, Martínez V, Nieves-Cordones M, Rubio F (2019) High external K+ concentrations impair pi nutrition, induce the phosphate starvation response, and reduce arsenic toxicity in Arabidopsis plants. Int J Mol Sci 20:2237. https://doi.org/10.3390/ijms20092237
Ronzan M, Piacentini D, Fattorini L, Caboni E, Eiche E, Ziegler J, Hause B, Riemann B, Betti C, Altamura MM, Falasca G (2019) Auxin-jasmonate crosstalk in Oryza sativa L. root system formation after cadmium and/or arsenic exposure. Environ Exp Bot 165:59–69. https://doi.org/10.1016/j.envexpbot.2019.05.013
Ruiz A, Yenush L, Ariño J (2003) Regulation of ENA1 Na+-ATPase gene expression by the Ppz1 protein phosphatase is mediated by the calcineurin pathway. Eukaryot Cell 2(5):937–948. https://doi.org/10.1128/EC.2.5.937-948.2003
Sasaki A, Yamaji N, Yokosho K, Ma JF (2012) Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice. Plant Cell 24(5):2155–2167. https://doi.org/10.1105/tpc.112.096925
Schmittgen TD, Livak KJ (2008) Analyzing real-time PCR data by the comparative CT method. Nat Protoc 3:1101–1108. https://doi.org/10.1038/nprot.2008.73
Schroder P, Lyubenova L, Huber C (2009) Do heavy metals and metalloids influence the detoxification of organic xenobiotics in plants? Environ Sci Pollut Res 16(7):795–804. https://doi.org/10.1007/s11356-009-0168-7
Sehar S, Feng Q, Adil MF, Sahito FS, Ibrahim Z, Baloch DM, Ullah N, Ouyang Y, Guo Y, Shamsi IH (2022) Tandem application of endophytic fungus Serendipita indica and phosphorus synergistically recuperate arsenic induced stress in rice. Front Plant Sci 13:982668. https://doi.org/10.3389/fpls.2022.982668
Shi SL, Wang DF, Yan Y, Zhang F, Wang HD, Gu M et al (2013) Function of phosphate transporter OsPHT2;1 in improving phosphate utilization in rice. China J Rice Sci 27:457–465. https://doi.org/10.3969/j.issn.1001-7216.2013.05.002
Shin H, Shin HS, Dewbre GR, Harrison MJ (2004) Phosphate transport in Arabidopsis: Pht1;1 and Pht1;4 play a major role in phosphate acquisition from both low- and high-phosphate environments. Plant J 39:629–642. https://doi.org/10.1111/j.1365-313X.2004.02161.x
Shri M, Singh PK, Kidwai M, Gautam N, Dubey S, Verma G, Chakrabarty D (2019) Recent advances in arsenic metabolism in plants: current status, challenges and highlighted biotechnological intervention to reduce grain arsenic in rice. Metallomics 11:519–532. https://doi.org/10.1039/C8MT00320C
Takahashi R, Ishimaru Y, Senoura T, Shimo H, Ishikawa S, Arao T, Nakanishi H, Nishizawa NK (2011) The OsNRAMP1 iron transporter is involved in cd accumulation in rice. J Exp Bot 62(14):4843–4850. https://doi.org/10.1093/jxb/err136
Vanstraelen M, E Benková (2012) Hormonal interactions in the regulation of plant development. Annu Rev Cell Dev Biol 28:463–487. https://doi.org/10.1146/annurev-cellbio-101011-155741
Varma A, Verma S, Sahay NS, Butehorn B, Franken P (1999) Piriformospora indica, a cultivable plant- growth-promoting root endophyte. Appl Environ Microbiol 65:2741–2744. https://doi.org/10.1128/AEM.65.6.2741-2744.1999
Vaseva II, Qudeimat E, Potuschak T, Du YL, Genschik P, Vandenbussche F, Van Der Straeten D (2018) The plant hormone ethylene restricts Arabidopsis growth via the epidermis. Proc Natl Acad Sci USA 115(17):E4130–E4139. https://doi.org/10.1073/pnas.1717649115
Verbruggen N, Hermans C, Schat H (2009) Mechanisms to cope with arsenic or cadmium excess in plants. Curr Opin Plant Biol 12(3):364–372. https://doi.org/10.1016/j.pbi.2009.05.001
Wang X, Ma LQ, Rathinasabapathi B, Cai Y, Liu YG, Zeng GM (2011) Mechanisms of efficient arsenite uptake by arsenic hyperaccumulator Pteris vittata. Environ Sci Technol 45(22):9719–9725. https://doi.org/10.1021/es2018048
Wang C, Huang W, Ying YH, Li S, Secco D, Tyerman S, Whelan J, Shou HX (2012) Functional characterization of the rice SPX-MFS family reveals a key role of OsSPX-MFS1 in controlling phosphate homeostasis in leaves. New Phytol 196(1):139–148. https://doi.org/10.1111/j.1469-8137.2012.04227.x
Weiß M, Sýkorová Z, Garnica S, Riess K, Martos F, Krause C, Oberwinkler F, Bauer R, Redecker D (2011) Sebacinales everywhere: previously overlooked ubiquitous fungal endophytes. PLoS ONE 6:e16793. https://doi.org/10.1371/journal.pone.0016793
Xie XN, Hu WT, Fan XN, Chen H, Tang M (2019) Interactions between phosphorus, zinc, and iron homeostasis in nonmycorrhizal and mycorrhizal plants. Front Plant Sci 10:1172. https://doi.org/10.3389/fpls.2019.01172
Yadav V, Molina I, Ranathunge K, Castillo IQ, Rothstein SJ, Reed JW (2014) ABCG transporters are required for suberin and pollen wall extracellular barriers in Arabidopsis. Plant Cell 26(9):3569–3588. https://doi.org/10.1105/tpc.114.129049
Yu J, Xuan W, Tian Y, Fan L, Sun J, Tang W, Chen G, Wang B, Liu Y, Wu W, Liu X, Jiang X, Zhou C, Dai Z, Xu D, Wang C, Wan J (2021) Enhanced OsNLP4-OsNiR cascade confers nitrogen use efficiency by promoting tiller number in rice. Plant Biotechnol J 19(1):167–176. https://doi.org/10.1111/pbi.13450
Yuan GQ, Zou T, Zhang X, Liu MM, Luo TT, He ZY, Tao Y, Zhou D, Zhao HF, Liu SJ, Liu R, Zhou ML, Zhu J, Liang YY, Deng QM, Wang SQ, Zheng AP, Wang AJ, Liu HN, Wang LX, Li P, Li SC (2020) A rice GDSL esterase/lipase protein (GELP) is required for anther and pollen development. Mol Breed 40(9):90. https://doi.org/10.1007/s11032-020-01170-4
Zeng FR, Chen S, Miao Y, Wu FB, Zhang GP (2008) Changes of organic acid exudation and rhizosphere pH in the rice plants under chromium stress. Environ Pollut 155:284–289. https://doi.org/10.1016/j.envpol.2007.11.019
Zhang HH, Wang ML, Li YQ, Yan W, Chang ZY, Ni HL, Chen ZF, Wu JX, Xu CJ, Deng XW, Tang XY (2020) GDSL esterase/lipases OsGELP34 and OsGELP110/OsGELP115 are essential for rice pollen development. J Integr Plant Biol 62(10):1574–1593. https://doi.org/10.1111/jipb.12919
Zhou Y, Huang WF, Liu L, Chen TY, Zhou F, Lin YJ (2013) Identification and functional characterization of a rice NAC gene involved in the regulation of leaf senescence. BMC Plant Biol 13:132. https://doi.org/10.1186/1471-2229-13-132
Zu SH, Jiang YT, Chang JH, Zhang YJ, Xue HW, Lin WH (2022) Interaction of brassinosteroid and cytokinin promotes ovule initiation and increases seed number per silique in Arabidopsis. J Integr Plant Biol 64(3):702–716. https://doi.org/10.1111/jipb.13197
Acknowledgements
This research work was financially supported by the Sino-Pakistan Project NSFC Grant No. 31961143008, National Natural Science Foundation of China, International (Regional) Cooperation and Exchange Program, Research fund for International young scientists Grant Nos. 31750110462, 32250410280 and Jiangsu Collaborative Innovation Centre for Modern Crop Production (JCIC-MCP) China.
Author information
Authors and Affiliations
Contributions
Shafaque Sehar and Muhammad Faheem Adil: Conceptualization, Methodology, Experimental work, Data analysis and Writing—original draft. Qidong Feng and Dongming Wei: Methodology and Experimental work. Syed Muhammad Hassan Askri, and Falak Sehar Sahito: Writing - review & editing. Imran Haider Shamsi: Conceptualization, Supervision, Funding acquisition, Correspondence and Writing—review & editing. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic Supplementary Material
Below is the link to the electronic supplementary material.
Additional file 1:
Table S1. List of salts used in making nutrient solution. Table S2. Sequences of oligonucleotide primers used for validating RNA-Seq data. Table S3. Summary statistics of rice transcriptome generated by Illumina NovaSeq 6000. Table S4. List of novel genes detected during As/S. indica/P interaction in ZZY-1 and GD-6. Table S5. Inclusively expressed detoxification-related differential transcripts (log2 fold change) observed in response to As/S. indica/P interaction in ZZY-1 and GD-6. Table S6. Inclusively expressed nutrient and metal transport-related differential transcripts (log2 fold change) observed in response to As/S. indica/P interaction in ZZY-1 and GD-6. Table S7. Inclusively expressed phytohormone-related differential transcripts (log2 fold change) observed in response to As/S. indica/P interaction in ZZY-1 and GD-6.
Additional file 2:
Fig. S1. Validation of data obtained from RNA-Seq with real-time quantitative reverse transcription PCR. Fig. S2. Summary of gene mapping ratio as observed across the samples including all treatments between ZZY-1 (G1) and GD-6 (G2). Fig. S3. Venn diagram representing exclusive and inclusive DEGs expressed for G1 = ZZY1 (A), and G2 = GD-6 (B), in response to As alone, As + P, As + S.i and As + S.i + P treatments along with control. Fig. S4. Principal component analysis (PCA) of As, P and S. indica-induced transcripts across the treatments in two rice genotypes. Fig. S5. Time series analysis illustrating Mfuzz results. T1 = As10µM, T2 = As10µM + P50µM, T3 = As10µM + S.i, T4 = As10µM + S.i + P50µM and T5 = control. Fig. S6. Distribution of differentially expressed transcripts using color module (A), and cluster dendrogram (B). Fig. S7. KEGG analysis of DEGs presented as scatter plot showing 20 most significantly enriched pathways for: (A) G1T1 vs. G2T1, (B) G1T2 vs. G2T2, and (C) G1T3 vs. G2T3 comparisons. Fig. S8. Relative comparative gene ontology (GO) in ZZY-1 under As10 (A), As10µM + P50µM (B), and As10µM + S.i (C) in comparison to GD-6. Fig. S9. Heat map presenting the expression profile of mutually expressed stress responsive DEGs. Fig. S10. Heat map showing the expression profile of inclusively expressed heavy metal associated transporters. Fig. S11. Heat map representing the expression profile of inclusively expressed ABC transporters. Fig. S12. Hierarchical heat map depicting the expression profile of inclusively expressed genes related to auxin biosynthesis and signaling. Fig. S13. Hierarchical heat map illustrating the expression profile of inclusively expressed genes related to the biosynthesis and signaling of ethylene. Fig. S14. Hierarchical heatmap representation of differentially expressed GELP transcripts (inclusively/exclusively) observed in response to As alone, As + P, As + S.i and As + S.i + P along with control, in two rice genotypes.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sehar, S., Adil, M.F., Askri, S.M.H. et al. Pan-transcriptomic Profiling Demarcates Serendipita Indica-Phosphorus Mediated Tolerance Mechanisms in Rice Exposed to Arsenic Toxicity. Rice 16, 28 (2023). https://doi.org/10.1186/s12284-023-00645-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12284-023-00645-0