Next Article in Journal
Quantitative Analysis and Molecular Docking Simulation of Flavonols from Eruca sativa Mill. and Their Effect on Skin Barrier Function
Next Article in Special Issue
Phenotyping and Exploitation of Kompetitive Allele-Specific PCR Assays for Genes Underpinning Leaf Rust Resistance in New Spring Wheat Mutant Lines
Previous Article in Journal
Maternal Influenza and Offspring Neurodevelopment
Previous Article in Special Issue
Wheat Water-Soluble Carbohydrate Remobilisation under Water Deficit by 1-FEH w3
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physiological and Transcriptomic Analyses Reveal Commonalities and Specificities in Wheat in Response to Aluminum and Manganese

Guangxi Key Laboratory of Agro-Environment and Agric-Products Safety, College of Agriculture, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2024, 46(1), 367-397; https://doi.org/10.3390/cimb46010024
Submission received: 29 November 2023 / Revised: 22 December 2023 / Accepted: 28 December 2023 / Published: 2 January 2024
(This article belongs to the Special Issue Advanced Research in Wheat Genome and Breeding)

Abstract

:
Aluminum (Al) and manganese (Mn) toxicity are the top two constraints of crop production in acid soil. Crops have evolved common and specific mechanisms to tolerate the two stresses. In the present study, the responses (toxicity and tolerance) of near-isogenic wheat lines (ET8 and ES8) and their parents (Carazinho and Egret) to Al and Mn were compared by determining the physiological parameters and conducting transcriptome profiling of the roots. The results showed the following: (1) Carazinho and ET8 exhibited dual tolerance to Al and Mn compared to Egret and ES8, indicated by higher relative root elongation and SPAD. (2) After entering the roots, Al was mainly distributed in the roots and fixed in the cell wall, while Mn was mainly distributed in the cell sap and then transported to the leaves. Both Al and Mn stresses decreased the contents of Ca, Mg, and Zn; Mn stress also inhibited the accumulation of Fe, while Al showed an opposite effect. (3) A transcriptomic analysis identified 5581 differentially expressed genes (DEGs) under Al stress and 4165 DEGs under Mn stress. Among these, 2774 DEGs were regulated by both Al and Mn stresses, while 2280 and 1957 DEGs were exclusively regulated by Al stress and Mn stress, respectively. GO and KEGG analyses indicated that cell wall metabolism responds exclusively to Al, while nicotianamine synthesis exclusively responds to Mn. Pathways such as signaling, phenylpropanoid metabolism, and metal ion transport showed commonality and specificity to Al and Mn. Transcription factors (TFs), such as MYB, WRKY, and AP2 families, were also regulated by Al and Mn, and a weighted gene co-expression network analysis (WGCNA) identified PODP7, VATB2, and ABCC3 as the hub genes for Al tolerance and NAS for Mn tolerance. The identified genes and pathways can be used as targets for pyramiding genes and breeding multi-tolerant varieties.

1. Introduction

Acid soils are found throughout the world. It is estimated that approximately 30~40% of the world’s total land area and more than 50% of the world’s potential arable land have a pH below 5.5; moreover, soil acidification is increasing worldwide [1]. The primary limitations of acid soils are toxic levels of aluminum (Al) and manganese (Mn), as well as suboptimal levels of phosphorous (P) [2].
Aluminum (Al) makes up 7% of the Earth’s crust and is the most abundant metal element. It is fixed in minerals or bound to mineral surfaces in non-acid soils, and thus benign to plant growth. When soil acidifies and the pH drops below 5, Al3+ is solubilized into the soil solution. The activity of Al3+ in the soil solution increases with the decreasing pH and peaks at pH 4.1 [3]. The primary symptom of Al toxicity is a rapid (beginning within minutes) inhibition of root growth, resulting in a reduced and damaged root system and limited water and mineral nutrient uptake [2].
Manganese is an essential micronutrient and a toxic element for higher plants, with a narrow concentration window between deficiency and toxicity. With a decreasing pH, the amount of exchangeable manganese, mainly in Mn2+ form, increases in the soil solution. Mn2+ can be readily transported into the root cells and translocated to the shoots, where it is finally accumulated [4]. Mn toxicity symptoms are localized to the shoot, characterized by stunted growth, chlorosis, and necrotic lesions in the leaves, accompanied by the inhibition of photosynthesis and the accumulation of excessive reactive oxygen species (ROS) [5,6].
As the co-existing metal becomes stressed, Al and Mn exert common toxicity effects on the morphology and physiology, such as inhibiting root growth, disturbing ion balance, and overproducing reactive oxygen species, subsequently resulting in oxidative damage, decreasing the chlorophyll content, and deteriorating the photosynthetic performance [2,7,8]. On the other hand, plants have evolved external exclusion and internal tolerance mechanisms to cope with both active Al and excess Mn. External exclusion mechanisms include the chelation of Al and Mn in the rhizosphere by organic acids and second metabolites [9] and the induction of a pH barrier in the rhizosphere [6,10,11]. By contrast, internal tolerance mechanisms enable plants to deal with Al3+ and Mn2+ once they enter plant cells, either by forming harmless complexes with organic acids, proteins, phenols, and other organic compound ligands, by sequestering them to organelles (e.g., vacuoles, endoplasmic reticula, and the Golgi apparatus), or by rapidly repairing any lesions incurred, including those resulting from oxidative stress [5,6,10].
Due to differences between Al and Mn in the binding site, distribution, and biological functions in plants, some plants’ responses (toxicity and resistance) to Al and Mn are specific [5,9,12] or similar but play different roles and are dominated by various genes. For example, plants transport Al and Mn from the roots to leaves or internally across the cytosol into various organelles, apoplasts, and cell walls through the activities of series transporters, representing a diverse set of metal transporter families [5,13,14]. In addition, the complex gene regulatory network of plant hormones, signal transduction, and transcription factors endows plants with specific and common tolerance mechanisms under Al and Mn stress, such as regulating root growth, cell wall modification, and organic acid secretion [5,15,16].
Traditional strategies used to maintain crop production in acid soils include applying lime to raise the soil pH and planting tolerant crop species and cultivars. Although the liming of soil can ameliorate some Al and Mn toxicity and increase calcium and magnesium concentrations in many acid soils worldwide, it is expensive and ineffective in the subsoil. In some cases, heavy lime application may cause a deleterious effect on the soil structure and reduce the availability of other cationic micronutrients that are essential for plant growth [17,18]. Plant species vary in their tolerance to stresses in acid soils. Native plants usually have more substantial tolerances to adverse factors than crops [19]. Moreover, there is a wide range of tolerance to either Al [20] or Mn [21] among genotypes within the same species. Ideal crop cultivars growing in acid soil should tolerate both Al and Mn toxicity simultaneously. However, plants may employ some specific mechanisms to tolerate Al or Mn despite their coexistence. Many cultivars tolerate only one toxicity [4,21]. Two strategies could be used to obtain duel-tolerant varieties. One approach is to select and breed manganese-tolerant varieties from aluminum-tolerant germplasms, or vice versa, but this is generally restricted by a genetic bottleneck. Another strategy is to design smart plants and pyramid genes to fortify the common ones and compensate for the differential tolerances. Uncovering the commonalities and the specificities of Al and Mn toxicity and tolerant mechanisms is a prerequisite for identifying and pyramiding hub genes and breeding multi-tolerant genotypes. Although the physiological mechanisms by which plants coordinately adapt to multiple stresses in acid soils have been extensively studied [4,21,22], the molecular mechanisms are still largely unknown. Many genes may have co-evolutionary mechanisms in their regulatory pathways and function in the context of multiple stresses.
In the face of increasing environmental challenges, understanding the molecular and physiological responses to abiotic stresses is paramount for ensuring global food security. Aluminum (Al) and manganese (Mn) stresses pose significant threats to crop productivity [4], particularly in wheat, a staple food for a large portion of the world’s population [22]. High-throughput RNA sequencing will be employed to capture the dynamic changes in gene expression profiles, and the transcriptome analysis will provide a genome-wide perspective on the molecular events occurring in wheat plants under Al and Mn stresses [23,24]. This study integrated physiological and transcriptome analyses to reveal the specificity and commonality in the toxicity and tolerance mechanisms of wheat under Al and Mn stresses, providing a new perspective on the molecular network of plants responding to multiple stresses in acid soil and further discovering candidate genes for wheat Al and Mn tolerance, providing genetic resources.

2. Materials and Methods

2.1. Plant Materials and Treatment

Near-isogenic (over 95%) wheat lines differing in Al tolerance at the ALMT1 locus, ET8 (homozygous Al-tolerant), and ES8 (homozygous Al-sensitive) [25], as well as their parents Carazinho (Al tolerant) and Egret (Al sensitive), were used as plant materials. The ET8 and ES8 lines were derived from a cross between Carazinho and Egret, with the resulting progeny being backcrossed eight times to Egret or derivatives of Egret. Grain yield increased in Egret-derived lines when Al tolerance from Carazinho was introduced. The seeds were provided by Peter Ryan from CSIRO Plant Industry, Canberra, Australia.
Wheat seeds were surface-sterilized with 0.5% NaClO (v/v) solution for 20 min, rinsed with sterile water, and allowed to germinate in wet paper towel for 4~5 days at 4 °C before the seedlings were transplanted into continuously aerated Magnavaca hydroponic solution (pH 4.2) for a 2-day acclimation to low pH in a walk-in growth chamber with a stable temperature of 23 °C, a 14 h photoperiod of 400 μmol m−2 s−1 illumination level, and a relative humidity of 60–80%. Then, the seedlings were stressed by replacing fresh Magnavaca solutions with an addition of 100 µM KAl(SO)4 (as Al toxicity), 500 µM MnCl2 (as Mn toxicity), or no addition (as control). Free Al3+ and Mn2+ concentrations in the solution were 25 µM and 267 µM, respectively, as predicted using the WinIAP program (https://www.sequentix.de/index.php#Winiap accessed on 10 October 2023). The experiment was arranged in a completely randomized block design with at least three replicates per batch.

2.2. Relative Root Length

Before and four days after stress, the fresh roots of different treatment groups were collected and completely expanded onto the scanner platform (Epson Expression 1000XL, Seiko Epson Corporation, Nagano, Japan). The total length per plant was evaluated using image analysis with WinRhizo Ver. 2.0 (Regent Instruments Inc., Quebec, QC, Canada). The RRL, the ratio of net root growth under Al or Mn treatment compared to the control, was calculated.

2.3. SPAD and Brown Spots

After four days of stress, the relative content of chlorophyll in fully developed leaves, as dictated by SPAD (Soil and Plant Analyzer Development) value, was determined with a portable Konica SPAD-502 Plus instrument (Konica Minolta Holdings Inc., Tokyo, Japan). Brown spots on leaves under Mn stress were checked and photographed.

2.4. Metal Ions Accumulation and Subcellular Distribution

Al and Mn concentrations in roots and shoots and Ca, Mg, Zn, and Fe concentrations in whole plant were determined after four days of stress. The plant materials were dried in an oven at 105 °C for 1 h, followed by drying at 80 °C for 3 h until a constant weight was achieved. Then, they were ground using a JC-FW100 grinder (Juchuan Co., Qingdao, China). After being dried and pulverized, the plant materials were conventionally digested with mixed acid (concentrated nitric acid/perchloric acid = 4:1) in a microwave digestion instrument (MARS, CEM company, Charlotte, NC, USA) with a digestion procedure of 150 °C for 8 min, 180 °C for 5 min, and 200 °C for 5 min, and then determined using inductively coupled plasma atomic emission spectrometry (ICP-AES, Juguang Co., Beijing, China).
Cell sap in root apices preparation and the extraction of ions (Al and Mn) from the residual cell wall were conducted according to [26]. Al and Mn concentrations were determined by ICP-AES. The mixture of standards was prepared using certified reference materials (CRM) purchased from Sigma-Aldrich (St. Louis, MO, USA) for Al, Mn, Ca, Mg, Zn, and Fe determination.

2.5. RNA Extraction, Library Preparation, Sequencing, and Read Mapping

Total RNA was extracted from 1 cm root apices using TRIzol Reagent (Invitrogen Co., Carlsbad, CA, USA) according to the manufacturer’s instructions. The RNA quality was determined using 5300 Bioanalyser (Agilent Co., Santa Clara, CA, USA) and quantified using the ND-2000 (Thermo Fisher NanoDrop, Waltham, MA, USA). Only high-quality RNA sample (OD260/280 = 1.8~2.2, OD260/230 ≥ 2.0, RIN ≥ 6.5, 28S:18S ≥ 1.0, > 1 μg) was used to construct sequencing library and verified by qRT-PCR.
RNA purification, reverse transcription, library construction, and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Thirty-six RNA-seq transcriptome libraries (4 varieties × 3 treatments × 3 biological replicates) were prepared and sequenced.
The raw paired-end reads were trimmed and quality controlled using Fastp [27] with default parameters. The clean reads were separately aligned to the wheat genome with orientation mode using HISAT2 software Ver. 2.2.1 [28]. The mapped reads of each sample were assembled using StringTie Ver. 1.3.6 [29] in a reference-based approach. All assembled genes were annotated against Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-Prot, Protein family (Pfam), Gene Ontology (GO), Clusters of Orthologous Groups of Proteins (COG), and NCBI non-redundant protein sequences (NR) libraries.

2.6. Differential Expression and Functional Enrichment Analysis

The expression level of each transcript was calculated according to the fragments per kilobase of exons model per million mapped reads (FPKM) method. RSEM was used to quantify gene abundances [30]. Differential expression analysis was performed using the DESeq2 [31]. DEGs with |log2FC| ≥ 1 and FDR ≤ 0.001 were considered significantly differently expressed genes. Cluster heat diagrams were drawn using Toolkit for Biologists (TB) tools Ver. 2. 030 with default settings [32]. Those data were added to the Comprehensive Gene Expression Database with the accession number PRJNA1031207. In addition, GO annotation and KEGG pathway analysis were carried out using Goatools and KOBAS [33] to identify which DEGs were significantly enriched in GO terms and metabolic pathways at Bonferroni-corrected p-value ≤ 0.05 compared with the whole-transcriptome background.

2.7. Quantitative Real-Time PCR (qRT-PCR) Validation

Using the identical RNA/cDNAs for RNA-seq as templates, qRT-PCR was performed on a Bio-Rad CFX96 (Bio-Rad Laboratories, Hercules, CA, USA) to verify the authenticity of transcriptomic profile expression patterns. The 10 μL reaction system contained 5 μL of TB Green®Premix Ex Taq™ II (Takara standard Co., Osaka, Japan); 10 μM of primers, each 0.2 μL; 1 μL of cDNA template, and 3.6 μL of ddH2O. The amplification procedure was initially 95 °C for 30 s, followed by 40 cycles of 95 °C for 10 s, and 60 °C for 30 s (two-step thermal cycling). Twenty-one candidate DEGs involved in various processes were randomly picked up as target genes, and the housekeeping gene TaActin was used as an internal control. The primers used for qRT-PCR are listed in Table S1. The relative level of expression was calculated using the 2−ΔΔCT formula [34] (Livak and Schmittgen, 2001).

2.8. Weighted Gene Co-Expression Network Analysis (WGCNA)

Co-expression networks were constructed based on pairwise correlations of gene expression across all samples to find clusters (modules) of highly correlated genes and explore the connection of the modules with physiological traits and the hub genes [35]. After filtering out the genes with a low expression (FPKM < 2; CV > 0.15), 9924 genes were retained for the WGCNA. We introduced the appropriate soft threshold power (β = 7) in this study according to scale independence and mean connectivity curves produced using the software. Then, the adjacency was transformed into a topological overlap matrix (TOM), and average linkage hierarchical clustering was conducted according to the TOM-based dissimilarity measure, along with merged genes with a minimum size of 30 and merge cut height of 0.25 as a module. Modules were defined as clusters of highly interconnected genes, and genes within the same cluster had high correlation coefficients among them. The total connectivity and intramodular connectivity were calculated with weighted and correlation coefficients in each module. Biologically meaningful co-expression modules were subsequently analyzed through the heatmap of the correlation coefficients between trait and module. Co-expression network analysis uses the top 30 gene nodes with connectivity, in which each node represents a gene, and the connecting lines (edges) between genes represent co-expression correlations. Interaction network visualization for each module was performed using Cytoscape version 3.9.1 [36]. Hub genes have the highest Kme (Eigengene connectivity) value, indicating the most connections in the network.

2.9. Statistical Analysis

Statistical analysis was carried out using SPSS version 19.0 for analysis of variance (ANOVA), followed by Duncan’s multiple range test (DMRT) to compare means among treatments if the ANOVA result was significant (p ≤ 0.05). Before ANOVA, the data were tested for normality using Chi-squared analysis and transformed where necessary.
Since RRL is a ratio of means, while the measured net root growth values (x) were not paired between control and stress groups, the mean value and standard derivation of RRG were calculated according to the following Equation [37]:
R R L ¯ = x ¯ A l   o r   M n x ¯ C K ;   S D R R L = R R L ¯ × ( S D A l   o r   M n x ¯ A l   o r   M n ) 2 + ( S D C K x ¯ C K ) 2

3. Results

3.1. Tolerance of Wheat Varieties to Al and Mn Toxicity

The inhibition of root growth, as indicated by the relative root length (RRL), is the initial and primary symptom of early Al toxicity. Aluminum tolerance was precisely the breeding target for creating NILs ES8 and ET8. Not surprisingly, the RRLs of the aluminum-tolerant parent, Carazinho, and the near-isogenic line, ET8, were significantly higher than those of Egret and ES8. Interestingly, the change order of the root length (Carazinho > ET8 > Egret > ES8) under Mn stress was the same as that under Al stress, and the change amplitude was similar (Figure 1).
Al and Mn toxicity decreased the chlorophyll level, as indicated by the SPAD value, and the decreasing amplitude values in Carazinho and ET8 were lower than those of Egret and ES8. There was no significant difference in the inhibitory effect of Mn and Al in Carazinho and ET8, while Mn showed a more substantial effect than Al in the Egret and ES8 varieties. Furthermore, Mn toxicity symptoms like brown spots on the old leaves or crinkles on the young leaves were found only on the leaves of Egret and ES8 after 4 d of Mn stress (Figure 2).
Collectively, nearly equal damage to the wheat was caused by 25 µM of Al and 267 µM of Mn. The concentrations of active Al and Mn stresses screened from the preliminary experiments were suitable for a further comparison of the phenotypes and mechanisms of the two stresses. The consistency of tolerance to Al and Mn in near-isogenic wheat lines (ET8 and ES8) and their parents (Carazinho and Egret) confirmed the possibility of co-adaptation and breeding multi-resistant varieties.

3.2. Al and Mn Subcellular Distribution in Root Cell and Contents in Plants

After an influx from the rhizosphere and entry into root cells, Al was mainly fixed in the cell wall, with a proportion variation of 82% to 93% among the genotypes, and only a small part was distributed in the cell saps. The proportions of Al fixed by the cell wall in the tolerant varieties were lower than those of the sensitive varieties, possibly resulting in a low degree of damage to the cell wall rigidification. By contrast, most of the Mn in the root cell was distributed in the cell sap rather than in the cell wall. There was a slight genotypic variation, but there was no correlation to Mn tolerance (Figure 3).
Most of the Al absorbed by the wheat accumulated in the root, and only a small part of Al (<20%) was transported to the shoot, while Mn was transported to the shoot, and only a small amount of Mn (<10%) stayed in the root (Figure 3). The root-to-shoot ratio (R/T) of the Al and Mn contents did not differ significantly among the cultivars. Of note, due to the small biomass, the Mn and Al contents in the root, shoot, and whole plant of ES8 and ET8 were significantly lower than those of their parents, Carazinho and Egret. There was no significant difference in concentration based on dry weight among cultivars.

3.3. Al and Mn Effectiveness on Four Metal Ions Contents

The variance analysis results revealed that for every element (Ca, Mg, Fe, and Zn), the contents in wheat were significantly different among stress treatments (F = 27.336~83.02, F0.01(2, 24) = 5.61), varieties (F = 10.9~116.62, F0.01(3, 24) = 4.72), and their combinations (F = 19.69~52.90, F0.01(11, 24) = 3.09). The Ca, Mg, and Zn levels in the wheat were dramatically decreased by both Al and Mn stresses. The decreasing effects of Al and Mn toxicity on the Ca content were roughly equivalent. However, the decreasing effects of Al stress on the Mg and Zn contents were more substantial than those of Mn toxicity. Al toxicity significantly increased the Fe contents, while Mn toxicity had the opposite effect. The significant differences in the Ca, Mg, Fe, and Zn contents among the varieties were not correlated to Al or Mn tolerance. The significance of multiple differences among stress combinations with varieties was labeled in Figure 4.

3.4. Wheat Root Transcriptome Profiling in Response to Al and Mn

Whole genome transcriptome sequencing analysis was applied to investigate the molecular responses of wheat roots to Al and Mn toxicity. Thirty-six cDNA libraries were constructed, combining three stresses, four varieties, and three replicates. These libraries yielded 1.8 G clean reads. The Q30 values of all samples exceeded 93.5% with a confidence of 99.9%, and the minimum mapping rate was over 90.55%, indicating that the sequencing was of high quality and qualification for the subsequent gene expression analysis.
A total of 115,534 expressed genes (91,609 known and 23,925 new) and 193,406 transcripts (113,693 known and 79,713 new) were assembled. The functions of 91,313 genes and 113,446 transcripts were finally annotated while referencing the KEGG, Swiss-Prot, Pfam, GO, COG, and NR libraries.
To validate the RNA-seq’s accuracy of gene expression, qRT-PCR was performed using 21 randomly selected genes (Table S1) and identical cDNA templates. The results of the qRT-PCR showed a very significant positive correlation with the transcriptome data, indicating that the transcriptome data were reliable (Figure 5).

3.5. Identification of Differentially Expressed Genes in Response to Al and Mn Toxicity

The FPKMs of genes expressed under Al and Mn toxicity were compared with the control using the criteria |log2FC| ≥ 1 and p < 0.001. The DEGs from various varieties were pooled within the stress treatments. A total of 5461 DEGs were identified under Al stress, comprising 3227 upregulated and 2234 downregulated DEGs (Figure 6A). Similarly, 4165 DEGs were identified under Mn stress, comprising 2617 upregulated and 1548 downregulated genes (Figure 6B).
The Venn diagram categorizes the 7699 genes responded by Al and Mn into eight sets, classified into four types based on expression patterns (Figure 6C,D, Table S2). Among these, 1950 genes expressed in the same pattern under Al and Mn toxicity (1283 genes upregulated and 667 genes downregulated by both Al and Mn; Type I) were identified as possible candidate genes contributing to common toxicity or tolerant mechanisms. In contrast, seven genes were differentially expressed under Al and Mn toxicity but in the opposite regulatory direction (Type Il), among which one gene was annotated as a solute carrier family 40 member I, and it was downregulated by Al but upregulated by Mn, and six genes were upregulated by Al but downregulated by Mn, including three genes encoding a bidirectional sugar transporter, one gene encoding ricin B-like lectin, one gene encoding cytochrome P450, and one gene encoding cysteine synthase. In addition, 3504 genes (1938 up, 1566 down) were exclusively regulated by Al (Type Ill), and 1950 genes (1333 up, 875 down) were exclusively regulated by Mn (Type IV). The latter three types of genes were studied to analyze the potential specific toxicity or tolerant mechanisms (Figure 6C,D).

3.6. GO Functional Annotations and KEGG Pathway Analysis

The GO annotation results showed that binding (GO: 0005488) and catalytic activities (GO: 0003824) in the molecular function (MF) categories; the membrane part (GO: 0044425), cell part (GO: 0044464), and organelle (GO: 0043226) in the cellular component (CC) categories; and the metabolic process (GO: 0008152), cellular process (GO: 0009987), and biological regulation (GO: 0065007) in the biological process (BP) categories were the most abundant DEGs of the GO annotation terms. In the transporter activity (GO: 0005215) term, 193 genes (88 up, 105 down) responded exclusively to Al, 161 genes (95 up, 66 down) exclusively responded to Mn, and 94 genes (49 up, 45 down) responded to both Al and Mn, expressed in different trends. In the transcription regulator activity (GO: 0140110) term, 145 genes (145 up, 35 down) responded exclusively to Al, 161 genes (109 up, 22 down) responded exclusively to Mn, and 94 genes (115 up, 14 down) responded to both Al and Mn, indicating that most transcription factors were upregulated in response to Al and Mn. In the antioxidant activity (GO: 0016209) term, 90 genes (32 up, 58 down) exclusively responded to Al, 59 genes (15 up, 44 down) exclusively responded to Mn, and 94 genes (5 up, 38 down) responded to both Al and Mn (Figure 7).
The top 20 KEGG enrichment pathways with the smallest FDR values are displayed in Figure 8. The results showed that the genes upregulated by Al or Mn were enriched into five categories, including metabolism (M,14 pathways), genetic information processing (GIP, 4 pathways), environmental information processing (EIP, 2 pathways), organismal systems (OS, 1 pathway), and cellular process (CP, 1 pathway). In contrast, the genes downregulated by Al or Mn were mainly grouped into the metabolism (M) category, which contains 22 pathways. The genes that were upregulated by both Al and Mn were enriched in protein processing in the endoplasmic reticulum (GIP, 137 genes), spliceosome (32 GIP), plant–pathogen interaction (OS, 31 genes), endocytosis (CP, 16 genes), and phenylalanine metabolism (M, 15 genes). The genes that were downregulated by both Al and Mn were enriched in seven pathways, such as phenylpropanoid biosynthesis (M, 46 genes), starch and sucrose metabolism (M, 10 genes), and base excision repair (GIPs, 5 genes). Genes that were upregulated by Al were enriched in diterpenoid biosynthesis (M, 16 genes) and aminoacyl-tRNA biosynthesis (GIP, 2 genes). The genes that were downregulated by Al were enriched into 12 pathways, such as fatty acid degradation (M, 10 genes), terpenoid backbone biosynthesis (M, 16 genes), and zeatin biosynthesis (M, 10 genes). Genes that were upregulated by Mn were enriched into plant hormone signal transduction (EIP, 28 genes); alpha-linolenic acid metabolism (M, 8 genes); phenylalanine, tyrosine, and tryptophan metabolism (M, 7 genes); and ABC transporter (6 genes) pathways. Genes that were downregulated by Mn were enriched into cutin, suberine, and wax biosynthesis (M, 14 genes); benzoxazinoid biosynthesis (M, 6 genes); and glycerolipid metabolism (M, 9 genes). Genes that were upregulated by Al and Mn were enriched in the MAPK signaling pathway (EIP, 23 genes upregulated by Al and 50 genes upregulated by Mn) and brassinosteroid biosynthesis (M, 5 genes upregulated by Al and 4 genes upregulated by Mn), and genes that were downregulated by Al and Mn were enriched in galactose metabolism (M, 9 genes downregulated by Al and 9 genes downregulated by Mn). Genes that were enriched in fatty acid (M, 8 genes upregulated by Al and 10 genes downregulated by Mn) and cysteine and methionine metabolism (M, 17 genes downregulated by Al and 23 genes upregulated by Mn) were oppositely regulated by Al and Mn.

3.7. Correlations between Physiological Traits and Expressed Module Eigengenes

The correlation analysis showed that RRL and SPAD were significantly positively correlated under Mn stress but not under Al stress. The Al contents in the root, stem, and total plant were significantly positively correlated, while they were not correlated with the distribution ratio between the cell wall and cell sap. Similarly, the Mn contents in the root, stem, and total plant were also significantly positively correlated. Besides that, they all were negatively regulated with the distribution in the cell wall and thus positively correlated with the distribution in the cell sap. There was no significant correlation between the Ca, Mg, Fe, and Zn variation amplitudes under Mn stress. However, the variation of Mg was positively correlated with Ca and negatively correlated with Fe and Zn under Al stress (Figure 9).
A total of 9924 filtered genes in 36 samples were clustered and divided into 18 module eigengenes (ME) decorated with diacritical colors, respectively, based on a scale-free topological model, β = 7 (Figure 10 and Figure S1). The correlation coefficient and correlation between the module and trait were calculated to explore the genes and expression patterns involved in the physiological trait and tolerance and to find specific expression genes in the samples. Five among eighteen MEs were notably correlated with at least one trait. MEmidnightblue (95 genes) was positively correlated with Mn traits (contents and distribution). MEblack (409 genes) was positively correlated with Al traits (contents and distribution). MElightcyan (93 genes) was positively correlated with Al traits and Fe content, while it was negatively correlated with Mn traits. MEred (642 genes) was negatively correlated with RRL. MEturquoise (2075 genes) was positively correlated with RRL, while it was negatively correlated with Al traits (Figure 10). The expression and annotation of genes in the above five modules are listed in Table S3-1–S3-5. The module is a highly interconnected set of genes that co-express to function and involve a biological process. However, the genes within a module are not necessarily differentially expressed genes. For instance, 2075 genes were clustered into MEturquoise; none were upregulated, and only a small portion was weakly downregulated by Al and Mn. RRL is the primary indicator for Al and Mn toxicity and tolerance screening. Forty-two genes in MEturquoise were directly related to root development, such as putative auxin transporter-like protein (LAX3), which was related to lateral root formation; bifunctional L-3-cyanoalanine synthase (PALP), which was related to root hair cell development; germin-like protein (GLP), which was related to the regulation of root development; protein root hair defective 3 (RHD3); root phototropism protein 3 (NPH3); and root-specific lectin (ChtBD) (Table S3-5). These genes did not respond to Al and Mn, but provided additive information that was not unfolded by the DEGs.
A connectivity analysis was conducted with the top 30 genes with high connectivity in each ME (MEmidnightblue, MEblack, MElightcyan, MEred, and MEturquoise) using Cytoscape to visualize the gene interaction regulation network (Figure 11).
The top five genes with a co-expression weight greater than 0.15, more connectivity, and higher Kme value in the interaction networks were screened. The candidate hub genes in the module were designated based on the annotated gene function and the correlation between the ME and traits. The connectivity, expressions, and annotations of the candidate hub genes are shown in Table 1.
In MEmidnightblue, all of the top five genes encoding nicotianamine synthase were upregulated exclusively by Mn. In MEblack, two genes that functioned as zinc finger transcription factors and three genes that were involved in intracellular trafficking and secretion were upregulated strongly by Al, while they were upregulated marginally by Mn. In MElightcyan, two genes functioned as inorganic ion transporters and were involved in glutathione metabolism; one gene encoding the TCP family transcription factor and one gene that was involved in protein processing in the endoplasmic reticulum presented more connectivities. However, their expressions were constitutive, and not strongly regulated by Al nor Mn. This result provides new insights into the mechanisms underlying Al and Mn toxicity and tolerance beyond differentially expressed genes (DEGs). In MEred, three genes involved in genetic information processing and two genes that encoded HSP20 were significantly upregulated by Al and Mn. In MEturquoise, two genes involved in proton and auxin transport and two genes involved in cell wall biogenesis and phenylpropanoid biosynthesis were slightly downregulated by Al and Mn. These genes may be a hub in their module and play a vital role in regulating the corresponding physiological trait.

4. Discussion

4.1. Cell Wall Biogenesis and Macromolecule Metabolism Respond Exclusively to Al Stress

The root tip is the first site for aluminum absorption, and aluminum toxicity inhibits the expansion and elongation of root cells and subsequent cell division [2]. Cell walls are the first barrier to Al uptake at the root apex. More than 80% of the total accumulated Al in plant roots is tightly bound to the cell walls, with only a tiny fraction entering the cell sap (cytoplasm), while Mn is mainly distributed in the cell sap (Figure 3). The primary cell wall components, pectin, and hemicellulose can bind Al; the binding of Al to cell walls decreases extensibility and cell elongation and increases rigidity [5]. Al stress significantly increases the contents of cell wall polysaccharides (pectin, hemicellulose I, and hemicellulose II) and the activities of pectin methylesterase and pectin demethylation at the root tips, increasing the fixation of Al by cell walls. The pectin methylation degree is determined by the pectin methylesterase (PME) activity [15]. Xyloglucan is the most active Hemicellulose component in the cell wall; it is biosynthesized in the Golgi apparatus by a series of glycan synthases and glycosyltransferases (XTHs) before being exported to the wall [38]. XTHs are considered important factors in controlling the strength and ductility of the cell wall and are extremely sensitive to Al [39]. Becnel et al. (2006) found that at least twenty-one XTH members had strong expressions in the root, especially XTH14, XTH15, and XTH31 [40]. The root-specific expressions in XTH5, XTH12, XTH13, XTH14, XTH17, XTH18, XTH19, XTH20, XTH26, and XTH31 were differentially expressed in response to Al [41].
This study found that genes regulating cell wall modification were specifically induced by Al stress but not significantly by Mn. Four of the five genes encoding PME were upregulated, while one was downregulated, which may explain the increase in cell wall Al binding (Figure 12). Similarly, XTHs were differentially expressed with a divergent pattern: XTH22, XTH23, XTH25, XTH27, and XTH28 were upregulated, while XTH8, XTH16, XTH24, XTH26, XTH30, and XTH31 were downregulated. In addition, the NAC transcription factor, which directly binds to the XTH31 promoter region in Arabidopsis [42], was also highly upregulated in our study. In vivo, the localization of XTH activity showed that Al greatly inhibited this enzyme activity within 30 min of exposure, which was concomitant with Al-induced callose deposition in the roots [39]. Under Al stress, plants may increase resistance by reducing the production of xyloglucans through the downregulation of the XTH genes, wherein XTH31 plays a significant role; the expression of other XTH genes may be involved in synergistic XTH31 expression or provide energy [39].
The WGCNA results showed that MEturquoise and MEblack were strongly correlated with Al, while they were weakly correlated with Mn. Numerous genes in the two MEs were involved in synthesizing and transporting cell wall components (Table S3-2 and S3-5). For instance, peroxidase (POD) P7, encoded by the hub gene TraesCS7D02G369700, is mainly involved in clearing root cell ROS and synthesizing lignin. Another hub gene, TraesCS3B02G280000, annotated as V-type proton ATPase bundle B2, regulates the synthesis and transport of lignin monomers (Table 1) [43]. Lignin synthesis is a typical defense mechanism in plant stress response. Lignin accumulates in the cell wall of the root endothelial layer, inhibiting the entry of heavy metals into the xylem or their transport from vascular bundles to the aboveground. OsSTAR1 and OsSTAR2, encoding a nucleotide-binding domain and a transmembrane domain of an ABC transporter, improve Al tolerance in rice through cell wall modification [44,45]. The TraesCS3A02G129000 gene encoding the ABC transporter C family member 3-like was exactly the hub gene in MEblack (Table 1).
The hardening of cell walls has a dual effect: on one hand, it inhibits the elongation of roots, and on the other hand, it prevents harmful metals from entering the cytoplasm. The modification of root cell wall properties can be reversible, so it is worth discussing whether cell wall stiffening is the cause or consequence of root growth inhibition.

4.2. Nicotianamine Synthesis Responds Exclusively to Mn Stress

Nicotianamine (NA) is a low-molecular-weight metal-bound ligand that maintains the homeostasis of Fe, Zn, and Mn. It protects cells from metal-induced damage, such as oxidative stress [46]. L-methionine is the precursor for NA biosynthesis, which binds with ATP to synthesize S-adenosylmethionine (SAM). Subsequently, three S-adenosylmethionine molecules are condensed to NA by nicotianamine synthase (NAS). NA is finally converted into 2′-deoxymugineic acids (DMA) and other mugineic acids (MAs) under the catalysis of nicotianamine aminotransferase (NAAT) [47]. The overexpression of MxNAS2 and MxNAS3 genes from M. xiaojinensis increases the concentration of NA-Mn in the flowers and leaves of tobacco and Arabidopsis, respectively [46]. Yellow strip-like protein (YSL) is involved in the uptake, transport, and relocation of metal complexes such as Mn-NA, Mn-HMA, and Mn-DMA to maintain Mn ion balance in plants [6,48].
The results included nine genes encoding NAS, four genes encoding NAAT, and two genes encoding YSL, which were specifically and strongly upregulated by Mn stress but not significantly regulated by Al (Figure 13). For instance, the expressions of NAS4, NAS1, NAS2, and NAATA increased by 3040, 420, 342, and 577-fold under Mn stress. In addition, the YSL2 gene (TraesCS6D02G223000), which transports Mn-NA and Mn-DMA [48,49,50], was also upregulated exclusively by Mn (Figure 14).
WGCNA showed that the genes encoding NAS, NAAT, and YSL2 were all significantly upregulated by Mn in MEmidnightblue; additionally, three genes encoding heavy-metal-associated domain (HMA) and one gene encoding ZIP were identified (Table S3-1). Most importantly, five genes encoding NAS were all included in the top 30 genes with the highest connectivity (Table 1), playing a hub role in the regulatory network to Mn stress.
Increasing the synthesis of NA and its derivatives and enhancing the long-distance transport of Mn-MAS is a specific strategy to tolerate Mn.

4.3. Metal Ion Transportation and Accumulation under Al and Mn Stresses

The Al absorption site is located at the root tip, and Al induces a series of transporter genes to coordinate the absorption, transport, and redistribution of Al. NRAMP aluminum transporter 1 (NRAT1) plays a vital role in rice Al tolerance by reducing the level of toxic Al in the root cell wall and transporting Al ions (not Al-citrate complex) into the root cell, where they are ultimately sequestered in the vacuole [26,51]. The ALS protein family is responsible for transporting excess Al to the vacuole for isolation and thus enhancing tolerance to Al in rice [52], tea [53], and buckwheat [54]. In this study, under Al stress, two nodulin 26-like intrinsic protein (NIP) genes, TaNIP1;1 (TraesCS7B02G122600) and TaNIP1;2 (TraesCS7A02G215700), were significantly upregulated by Al stress. This may promote the transport of the malate–Al complex (not Al ion) from the cell wall to the symplasm and further transport to the leaf vacuole for segregation [14,45,55]. WGCNA showed that MEblack and MElightcyan were significantly correlated with the content and distribution of Al. Some genes upregulated in MEblack may be involved in Al absorption, including those annotated as ABC transporters and one gene encoding NRAMP (Table S3-2).
Manganese is essential for plant growth, with a narrow concentration window between deficiency and toxicity. Therefore, Mn transport systems should be tightly regulated in plants. Plants must deliberately control the uptake, transportation, and distribution of manganese to prevent toxicity caused by excessive amounts of this element. Transporters play a vital role in these processes and thus help plants to tolerate Mn toxicity. Mn absorption occurs in the mature root zone and can be easily transported through the xylem to the shoot [56]. The Arabidopsis root plasma membrane-localized transporter natural resistance-associated macrophage protein (AtNRAMP1) mediates Mn2+ absorption. Two members of the zinc/iron-regulated transporter protein (ZIP) family, AtZIP1 and AtZIP2, are involved in the transport of Mn from the root to stem [8,57]; OsYSL2 is responsible for the long-distance transport and distribution of Mn [50]. In this study, the genes of TaYSL2 (TraesCS6D02G223000), TaZIP1 (TraesCS2A02G424200), and TaZIP2 (TraesCS6A02G158700) were significantly upregulated under Mn stress, and thus increased Mn absorption (Figure 14). Mn shares some ion channel proteins with other metal ions. For instance, OsNRAMP5, a plasma membrane protein, is involved in constitutive Fe and Mn uptake and inducive transport during flowering and seed development [58]. In the present study, the expression of reported multifunctional transporter genes, such as CAX2, CCX3, MTP9, MTP10, MTP11, MTP11, NRAMP3, NRAMP5, and YSL6 [8], were not regulated by Mn. The diversity of genes involved in Mn uptake and transport and the versatility of genes lead to the complexity of the regulation of Mn ion balance.
Al and Mn inhibit root growth, subsequently leading to the deprivation of essential nutrient elements. Due to its higher charge, the Al ion inhibits the absorption of bivalent cations (Ca2+, Mg2+, Fe2+, and Zn2+) more strongly than other toxic ions [1]. It removes bivalent cations from the plasma membrane and cell wall [59]. High concentrations of Mn2+ compete with other nutrient cations, thus hindering their uptake [8]. On the other hand, roots can improve tolerance to Al and Mn toxicity stresses by controlling the uptake of mineral nutrients and maintaining ion homeostasis, including calcium, magnesium, zinc, and iron [6,8,10].
Our results showed that Al and Mn stresses decreased the Ca, Mg, and Zn contents, and decreased the absorption and accumulation of Fe by Al, while Mn exerted the opposite effects. Al and Mn significantly downregulated one gene (TraesCS5D02G030700) encoding the cation/calcium exchanger (CCX1), contributing to the decrease in the Ca content. Plants utilize calmodulin-activated Ca2+ pumps (ACA) at the plasma membrane, endoplasmic reticulum, and vacuole [60]. The ACA7 (TraesCSU02G003600) and ACA5 (TraesCS5D02G284100) might play essential roles in altering the production of cellular Ca2+ signals and thus mediate adaptive responses (Figure 14).
Al and Mn shared hydrating ionic radii with Mg, competing for binding sites on apoplasm and transporters in the plasma membrane [61]. Aquaporin (PIP) regulates water and Mg2+ transport [62]. Two PIPs (TraesCS5A02G100000 and TraesCS7B02G002000) were significantly downregulated, which was related to a decrease in Mg2+ absorption. The decrease in the Mg content under Al and Mn stress may be the reason for the decrease in the chlorophyll content (SPAD). The overexpression of MGT1 can enhance Al tolerance [63]. The expression of one MGT gene (TraesCS3A02G380600) was upregulated under Al stress but not Mn due to an adaptive response related to the regulation of membrane potential balance by H+-ATPase activity under Al stress (Figure 14).
Fe is an essential trace element for plant growth and development. However, free Fe in cells produces toxic ROS through the Fenton reaction [64], necessitating strict control of iron homeostasis through various strategies. Al stress increased Fe absorption in tea plants by upregulating Fe ion transporter genes, including FRO, YSL, and IRT, maintaining Fe homeostasis [53]. Mn and Fe have antagonistic solid effects, and excessive Mn2+ takes advantage of Fe on the binding of Fe/Mn co-transporters [8,65]. Most of the CDF/MTP and VIT families can transport Fe2+ and Mn2+, and excess Mn2+ blocks the absorption of Fe2+ [8]. In this study, three VIT family genes and one YSL5 gene were significantly downregulated under Al and Mn stresses, reducing Fe absorption. One IRT1 gene (TraesCS4D02G017600) and one MGP gene (TraesCS3D02G257900) were upregulated by Al stress but downregulated by Mn stress (Figure 14). These results indicate that Al stress increased Fe2+ absorption by upregulating IRT and MRP. MTP4 (TraesCS4A02G077100) and MTP7 (TraesCS3D02G038000) were downregulated and upregulated by Mn stress, respectively, which affected the absorption and transport of Mn and Fe.
High active Al3+ competes with Zn2+ for binding sites on the root plasma membrane, and excess Mn2+ shares a similar transporter with Zn2+, resulting in the plant inhibition of zinc absorption [1,6]. Heavy-metal-associated domain (HMA) and ZIP proteins are involved in the absorption and transport of Zn [8]. Our study indicated that seven genes from the HMA family were significantly downregulated under Al and Mn stresses. The expression of genes from the ZIP family was diverse, with the significant upregulation of ZIP1 (TraesCS2A02G424200) and the downregulation of ZIP8 (TraesCS2B02G533800) under Al and Mn stresses. In addition, ZIP2 (TraesCS6A02G158700) was upregulated by Mn stress but downregulated by Al, which is favorable for the absorption of Mn and Zn (Figure 14).
In MEturquoise (Table S3-5), the genes encoding the calcium load-activated calcium channel (CAEEL); the Fe transporters, VIT and YSL; and the Zn transporter, HMA, were downregulated by Al and Mn. In contrast, Zn transporter genes (ZIP) in the MEred were upregulated (Table S3-4). MEturquoise was positively correlated with RRL, while MEred was negatively correlated with RRL. Reducing the relative root length may decrease the absorption area and gene expression levels, ultimately inhibiting nutrient uptake.

4.4. Phenylpropanoid Biosynthesis under Al and Mn Stresses

The KEGG results indicated that genes involved in phenylalanine metabolism commonly responded to Al and Mn, and most genes involved in phenylpropanoid metabolism were downregulated by both Al and Mn; some genes involved in flavonoid metabolism specifically responded to Al (Figure 8). There is the same expression trend of genes involved in the metabolism pathway from phenylalanine to sinapoy-choline and scopolin under Al and Mn, with the downregulation of seven genes encoding phenylalanine ammonia-lyase (PAL), three genes encoding putrescine hydroxycinnamoyltransferase (PHT), one gene encoding caffeoylshikimate esterase (CSE), two genes encoding tricetin 3′,4′,5′-O-trimethyltransferase (COMT), three genes encoding Cytochrome P450 (CYP), and three genes encoding Anthocyanin 3′-O-beta-glucosyltransferase (TOGT), and the upregulation of one gene (TraesCS3A02G329200) encoding serine carboxypeptidase-like (SCPL); however, one gene (TraesCS7A02G398900) encoding CYP45073A and one gene (TraesCS5B02G522400) encoding TOGT1 were exclusively upregulated by Al. Two genes (TraesCS6D02G136600 and TraesCS6A02G147300) encoding cinnamoyl-CoA reductase (CCR) were upregulated by Al and Mn. Among seven genes encoding beta-glucosidase (BGLG), three were downregulated by Al and Mn, BGLG31 (TraesCS5D02G302700) and BGLG6 (TraesCS4A02G056400) showed significant upregulation by Al, and BGLB12 (TraesCS2A02G329000) and BGLB3 (TraesCS5D02G404700) were upregulated by Mn but downregulated by Al. One gene (TraesCS3D02G388100) encoding 4-coumarate-CoA ligase (4CL) and one gene (TraesCS7D02G328300) encoding cinnamyl alcohol dehydrogenase (CAD) were exclusively up- and downregulated by Al, respectively (Figure 15).
Phenylpropanoid biosynthesis is an important pathway that produces plant secondary metabolites, such as phenylpropanoids, flavonoids, and diterpenoids, which are secreted in response to metal ions [66,67]. PAL catalyzes phenylalanine cinnamic acid, the first rate-limiting enzyme in the phenylalanine biosynthesis pathway [68]. Our results showed that genes involved in the pathway from cinnamic acid to sinapoy-choline and scopoline were downregulated, while PAL and CCR were upregulated by both Al and Mn. Similar to our study, genes encoding PAL and CCR were upregulated in phenylpropanoid metabolism and are activated to counteract Mn stress in Arabis paniculate [69]. In Populus euphratica, the content of phenylalanine and the antioxidant enzyme activity increases under Mn stress [70]. Phenylpropanoid metabolism exhibits different expression patterns in the roots and leaves under Mn stress in Stylosanthes. The expressions of PAL1 and chalcone-flavanone isomerase (CFI) were upregulated in the leaves, while PAL2 and chalcone synthase (CHS) were downregulated in the roots [71]. It was also found that phenylpropanoid biosynthesis and phenylalanine metabolism played a critical role in the defense against aluminum [72]. Lettuce protected roots from oxidative damage by upregulating phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase, and 4-coumarate coenzyme A ligase [73]. These results suggest that phenylpropanoid biosynthesis enhanced the antioxidant capacity for the tolerance against both Al and Mn.
In addition, Al specifically upregulated 4CL and CYP, which are involved in the synthetic precursors of flavonoid metabolism, p-coumaroyl-CoA, and cinnamoyl-CoA. The KEGG results showed that Al exclusively downregulated some genes involved in flavonoid metabolism, and two genes (TraesCS2D02G582000 and TraesCS4A02G003100) encoding tryptamine benzoyltransferase 1 (TBT1) and one gene (TraesCS2A02G482200) encoding anthocyanidin reductase (ANR) were exclusively upregulated (Table S2-3). The flavonoid biosynthesis site is located at the upper end of the surface cells in the plants’ root growth area and participates in internal and external Al detoxification by forming solid complexes with toxic Al ions [74]. Alfalfa increased the resistance to Al stress by increasing flavonoid accumulation in roots and secretion from the root tips [75]. Our results suggest that wheat has a tolerance to Al by promoting flavonoids to chelate Al ions.

4.5. Signaling in Roots under Al and Mn Stresses

According to the KEGG results, the majority of genes involved in signaling, including those related to the MAPK signaling pathway, plant–pathogen interaction, and plant hormones like IAA, ETH, and JA, exhibited similar expression patterns. However, the specificity of ZT and BR in response to Al and Mn was noticeable.
Genes that are commonly expressed in plant hormones (Figure 8), including two genes (TraesCS2B02G210600 and TraesCS2D02G191800) encoding indole-3-acetic acid-amido synthetase GH3, one gene (TraesCS3D02G411300) encoding auxin-responsive factor (ARF), two genes (TraesCS2D02G391400 and TraesCS1A02G328800) encoding ethylene-responsive transcription factor (ERF), and six genes encoding jasmonate ZIM domain protein (JAZ), were upregulated. Two genes (TraesCS3A02G348400 and TraesCS3D02G342000) encoding abscisic acid receptor PYL were downregulated by both Al and Mn, with a pronounced downregulation by Al than Mn. While some genes exclusively responded to Al or Mn, three genes involved in brassinosteroid (BR) were exclusively upregulated by Al, while one gene (TraesCS5A02G214800) encoding CYP92A6 was exclusively upregulated by Mn. Furthermore, three genes involved in zeatin (ZT) were exclusively downregulated by Al, including two genes (TraesCS3A02G263300 and TraesCS3D02G263000) encoding adenylate isopentenyl transferase 1 (IPT) and one gene (TraesCS2B02G013000) encoding protein NRT1/ PTR family 8.5 (PTR) (Figure 16).
IAA and JA affected the antioxidant defense mechanisms. They altered the cellular redox homeostasis for Al and Mn, showing the commonality of plant hormones in response to Al and Mn stresses [76]. In addition, plants respond to both stresses in some specific mechanism. The primary lesion of Al toxicity is the rapid inhibition of root elongation, with plant hormones primarily reacting to Al stress in the roots [10]. Genes involved in BR and ZT were specifically reposed to Al, which might have contributed to root growth and organic acid. The activated ethylene signal by Al promoted IAA or cytokinin (CTK) accumulation and led to a modification in the cell wall structure through auxin-responsive factors (ARFs). It was reported that BZ signaling is involved in the expressions of ALMT1 and STOP. It also stimulated the synthesis of the jasmonic acid signaling to modulate microtubule (MT) polymerization to regulate root growth [15]. Auxin-regulated GmMATE and H+-ATPase enhance Al tolerance by modifying the cell walls and promoting citrate exudation. Abscisic acid (ABA) was found to mediate the Al-induced citrate exudation in soybean and ethylene-mediated JA signal-regulated ALMT1-controlled malate exudation [15]. These results suggest that plants participating in Al tolerance mainly focus on the roots. While the changes in plant hormone content were influenced by Mn accumulation, moderate Mn increased the IAA, ZT, and 6-BA (6-benzylaminopurine) contents but decreased the ABA content in the leaves. A high Mn accumulation significantly increased the leaf JA and ABA contents, which inhibited plant growth and stimulated stress tolerance [77]. Mn enhanced the resistance to Mn toxicity based on IAA and ABA synthesis by regulating growth and scavenging ROS for antioxidants and Mn transporter [78].
Most genes involved in the MAPK signaling pathway and plant–pathogen interaction expressed the same trend under Al and Mn stresses. Six genes encoding calmodulin-like proteins (CML) and one gene (TraesCS4B02G327800) encoding respiratory burst oxidase homolog protein (RBOH) were significantly upregulated by Al and Mn. Among three genes encoding mitogen-activated protein kinase (MAPK), M2K9 and M3K17 showed significant upregulation, while Al and Mn stresses downregulated M3K18. Moreover, one gene (TraesCS7D02G161200) encoding a pathogenesis-related protein (PRP) and MEKK1 was exclusively upregulated by Mn (Figure 16).
Signaling, such as Ca2+, ROS, and NO, is a late response to Al stress. Al promoted the binding of Ca2+ to CML24, leading to malate exudation by regulating the expressions of STOP1, WRYK46, and ALMT1 [79]. The MAPK pathway positively regulates STOP1, thus conferring Al resistance [80]. Similarly, the response of these signals is vital for Mn. Mn2+ increased the catalase activity by regulating mitogen-activated protein kinases (MAPKs) and calmodulins in the presence of calcium to reduce Mn toxicity [81]. The production of one of the ROS, hydrogen peroxide (H2O2), is often generated contemporaneously with NO, and is regulated by CML and RBOH, respectively. NO was shown to play a familiar role in reducing Al toxicity by regulating the root growth and antioxidant capacity in wheat roots [82,83] and alleviating Mn toxicity by preventing oxidative stress in rice [84].

4.6. Transcription Factors Regulate Wheat Tolerance to Al and Mn Stresses

The GO annotation analysis showed that transcription factors (TFs) were predominantly upregulated by Al and Mn (Figure 7). Transcription factors are the terminal points of stress signal transduction and molecular switches for downstream gene expression [85]. Numerous TFs, such as WRKY, MYB, GATA, bZIP, bHLH, and ERF, have been reported to regulate Al and Mn stress responses [86,87,88]. For instance, sensitive to proton toxicity (STOP1) and its rice homolog Al resistance transcription factor (ART1) are both members of the C2H2 zinc finger family, which enhances Al tolerance by regulating ALMT1 [89,90]. The overexpression of GsERF1 enhances aluminum tolerance in Arabidopsis thaliana through ethylene-mediated pathways and ABA signaling pathways [91]. Similarly, ZmbHLH105 improves tolerance to manganese stress by regulating antioxidant mechanism-mediated ROS clearance and Mn/Fe-related transporter expression in plants [92]. Some transcription factors that regulate genes associated with aluminum (Al) tolerance may also have similar or opposite roles in manganese (Mn) tolerance. For instance, C2H2 transcription factors activate Al-tolerant genes, including glutamate dehydrogenase 1 (GDH1), pectin methylesterase inhibitor (PMI), malic enzyme (ME), aluminum-activated malate transporter (ALMT), multidrug and toxic compound exclusion protein (MATE), and tonoplast dicarboxylate transporter (TDT) [93,94]. These homologous genes were upregulated by Mn toxicity [95]. WRKY46 is a negative regulator of ALMT1 in Arabidopsis thaliana, and the destruction of WRKY46 leads to increased malate secretion, resulting in Al tolerance [96]. However, WRKY gene transcripts in stylo are enhanced by Mn stress [95], suggesting that there are conserved and divergent regulatory networks mediated by WRKY-mediated gene expression under Al and Mn toxicity.
An expression analysis showed that 640 TFs were significantly regulated by at least one stress, covering 29 TFs families (Table 2). The MYB family accounted for the most, followed by the WRKY, AP2, bHLH, and NAC families. The differentially expressed TFs from the BES1, GRF, HD-ZIP, RAV, TCP, and NF-YA families were exclusively regulated by Al. In contrast, two TFs from calmodulin-binding transcription activators (CAMTAs) and ethylene-insensitive 3-like (EIL) families were exclusively regulated by Mn, indicating the specificity of these TFs in sensing stress and regulating downstream genes. A total of 455 differentially expressed TFs were found under Al stress, comprising 291 upregulated and 82 downregulated TFs. Similarly, 373 TFs were differentially expressed under Mn stress, including 414 upregulated and 212 downregulated TFs. There were far more upregulated TFs than downregulated TFs by either Al or Mn, suggesting that wheat tends to activate TFs to adapt to toxicity. Numerous differentially expressed TFs were found under only one stress, with 267 TFs exclusively regulated by Al and 185 by Mn. A total of 188 differentially expressed TFs, including 156 upregulated and 32 downregulated TFs, were identified under two stresses in the same trend, demonstrating that the most co-expressed TFs function similarly under Al and Mn stresses.

5. Conclusions

The current study analyzed the physiology and transcriptomics of roots from wheat seedlings that were tolerant to Al and Mn stresses. The results showed that Al was mainly fixed in the root cell walls, regulated by PEMs and XTHs. Mn was distributed primarily on the cell saps and then transferred to the shoots by NA and YSL. The expression of metal cation transporter genes was altered by Al and Mn, leading to a reduced absorption of Ca, Mg, and Zn. Additionally, Mn decreased Fe absorption, while Al had the opposite effect. The identified DEGs associated with commonality and specific pathways based on the KEGG enrichment analysis might represent an adaptive strategy of wheat to cope with Al and Mn. Numerous genes related to TFs, such as MYB, WRKY, and AP2 families, were activated to tolerate Al and Mn stresses. In addition, WGCNA identified the hub genes PODP7, VATB2, and ABCC3 for Al tolerance and NAS for Mn tolerance. For these reasons, a further study is recommended to determine the ability of pyramid genes to breed multi-tolerant varieties, to understand the interaction of TFs with other proteins, and to regulate the target genes.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cimb46010024/s1, Figure S1: Network topology of different soft threshold powers; Table S1: Primers used for qRT-PCR; Table S2: Expression and functional annotation of DEGs in response to Al and Mn (6 sheets contained); Table S3: Expression and functional annotation of genes in 5 MEs (5 sheets contained).

Author Contributions

Conceptualization, D.D. and Y.X.; methodology, C.X. and D.L.; investigation, C.X. and D.L.; formal analysis, W.Z.; validation, Y.Q., Q.L. and S.S. writing—original draft preparation, D.D.; writing—review and editing, D.D. and M.U.; funding acquisition, D.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 31960633, and the Guangxi Natural Science Foundation, grant number 2023GXNSFAA026445.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in this study are available in the article or Supplementary Materials. The RNA-seq raw data can be found on the NCBI repository, accession number PRJNA1031207.

Acknowledgments

We are highly grateful to Peter Ryan form the CSIRO Plant Industry for providing the wheat seeds.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Moustaka, J.; Ouzounidou, G.; Baycu, G.; Moustakas, M. Aluminum resistance in wheat involves maintenance of leaf Ca2+ and Mg2+ content, decreased lipid peroxidation and Al accumulation, and low photosystem II excitation pressure. Biometals 2016, 29, 611–623. [Google Scholar] [CrossRef] [PubMed]
  2. Kochian, L.V.; Hoekenga, O.A.; Pineros, M.A. How do crop plants tolerate acid soils?—Mechanisms of aluminum tolerance and phosphorous efficiency. Annu. Rev. Plant Biol. 2004, 55, 459–493. [Google Scholar] [CrossRef] [PubMed]
  3. Kinraide, T.B. Toxicity factors in acidic forest soils: Attempts to evaluate separately the toxic effects of excessive Al3+ and H+ and insufficient Ca2+ and Mg2+ upon root elongation. Eur. J. Soil Sci. 2003, 54, 323–333. [Google Scholar] [CrossRef]
  4. Muhammad, N.; Cai, S.; Shah, J.M.; Zhang, G. The combined treatment of Mn and Al alleviates the toxicity of Al or Mn stress alone in barley. Acta Physiol. Plant. 2016, 38, 1–7. [Google Scholar] [CrossRef]
  5. Wang, X.; Ai, S.; Liao, H. Deciphering interactions between phosphorus status and toxic metal exposure in plants and rhizospheres to improve crops reared on acid soil. Cells 2023, 12, 441. [Google Scholar] [CrossRef] [PubMed]
  6. Millaleo, R.; Reye-Diaz, M.; Ivanov, A.G.; Mora, M.L.; Alberdi, M. Manganese as essential and toxic element for plants: Transport, accumulation and resistance mechanisms. J. Soil Sci. Plant Nutr. 2010, 10, 470–481. [Google Scholar] [CrossRef]
  7. Nahar, K.; Hasanuzzaman, M.; Suzuki, T.; Fujita, M. Polyamines-induced aluminum tolerance in mung bean: A study on antioxidant defense and methylglyoxal detoxification systems. Ecotoxicology 2017, 26, 58–73. [Google Scholar] [CrossRef]
  8. Alejandro, S.; Hoeller, S.; Meier, B.; Peiter, E. Manganese in plants: From acquisition to subcellular allocation. Front. Plant Sci. 2020, 11, 300. [Google Scholar] [CrossRef]
  9. Morita, A.; Yanagisawa, O.; Maeda, S.; Takatsu, S.; Ikka, T. Tea plant (Camellia sinensis L.) roots secrete oxalic acid and caffeine into medium containing aluminum. Soil Sci. Plant Nutr. 2011, 57, 796–802. [Google Scholar] [CrossRef]
  10. Kochian, L.V.; Pineros, M.A.; Liu, J.; Magalhaes, J.V. Plant adaptation to acid soils: The molecular basis for crop aluminum resistance. Annu. Rev. Plant Biol. 2015, 66, 571–598. [Google Scholar] [CrossRef]
  11. Liu, J.; Pineros, M.A.; Kochian, L.V. The role of aluminum sensing and signaling in plant aluminum resistance. J. Integr. Plant Biol. 2014, 56, 221–230. [Google Scholar] [CrossRef] [PubMed]
  12. Lytle, C.M.; Lytle, F.W.; Smith, B.N. Use of XAS to determine the chemical speciation of bioaccumulated manganese in Potamogeton pectinatus. J. Environ. Qual. 1996, 25, 311–316. [Google Scholar] [CrossRef]
  13. Che, J.; Zhao, X.Q.; Shen, R.F. Molecular mechanisms of plant adaptation to acid soils: A review. Pedosphere 2023, 33, 14–22. [Google Scholar] [CrossRef]
  14. Wang, Y.; Li, R.; Li, D.; Jia, X.; Zhou, D.; Li, J.; Lyi, S.M.; Hou, S.; Huang, Y.; Kochian, L.V.; et al. NIP1; 2 is a plasma membrane-localized transporter mediating aluminum uptake, translocation, and tolerance in Arabidopsis. Proc. Natl. Acad. Sci. USA 2017, 114, 5047–5052. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, P.; Wan, N.; Horst, W.J.; Yang, Z.B. From stress to responses: Aluminium-induced signalling in the root apex. J. Exp. Bot. 2023, 74, 1358–1371. [Google Scholar] [CrossRef]
  16. Leskova, A.; Javot, H.; Giehl, R.F.H. Metal crossroads in plants: Modulation of nutrient acquisition and root development by essential trace metals. J. Exp. Bot. 2022, 73, 1751–1765. [Google Scholar] [CrossRef]
  17. Moreira, A.; Moraes, L.A.C.; Navroski, D. Lime and micronutrients interaction in soybean genotypes adapted to tropical and subtropical conditions. Commun. Soil Sci. Plant Anal. 2017, 48, 792–800. [Google Scholar] [CrossRef]
  18. Rao, K.S.; Rao, B.S.; Vishnuvardhan, D.; Prasad, K.V. Alteration of superhelical state of DNA by aluminium (Al). BBA-Biomembr. 1993, 1172, 17–20. [Google Scholar] [CrossRef]
  19. Watanabe, T.; Osaki, M. Mechanisms of adaptation to high aluminum condition in native plant species growing in acid soils: A review. Commun. Soil Sci. Plant 2002, 33, 1247–1260. [Google Scholar] [CrossRef]
  20. Zishiri, R.M.; Mutengwa, C.S.; Tandzi, L.N.; Manyevere, A. Growth response and dry matter partitioning of quality protein maize (Zea mays L.) genotypes under aluminum toxicity. Agronomy 2022, 12, 1262. [Google Scholar] [CrossRef]
  21. Khabaz-Saberi, H.; Rengel, Z.; Wilson, R.; Setter, T.L. Variation of tolerance to manganese toxicity in Australian hexaploid wheat. J. Plant Nutr. Soil Sci. 2010, 173, 103–112. [Google Scholar] [CrossRef]
  22. Khabaz-Saberi, H.; Rengel, Z. Aluminum, manganese, and iron tolerance improves performance of wheat genotypes in waterlogged acidic soils. J. Plant Nutr. Soil Sci. 2010, 173, 461–468. [Google Scholar] [CrossRef]
  23. Salvador-Moreno, N.; Ryan, P.R.; Holguin, I.; Delhaize, E.; Benito, C.; Gallego, F.J. Transcriptional profiling of wheat and wheat-rye addition lines to identify candidate genes for aluminum tolerance. Biol. Plant. 2018, 62, 741–749. [Google Scholar] [CrossRef]
  24. Campos, C.; Nobre, T.; Goss, M.J.; Faria, J.; Barrulas, P.; Carvalho, M. Transcriptome analysis of wheat roots reveals a differential regulation of stress responses related to arbuscular mycorrhizal fungi and soil disturbance. Biology 2019, 8, 93. [Google Scholar] [CrossRef] [PubMed]
  25. Sasaki, T.; Yamamoto, Y.; Ezaki, B.; Katsuhara, M.; Ahn, S.J.; Ryan, P.R.; Delhaize, E.; Matsumoto, H. A wheat gene encoding an aluminum-activated malate transporter. Plant J. 2004, 37, 645–653. [Google Scholar] [CrossRef] [PubMed]
  26. Li, J.Y.; Liu, J.P.; Dong, D.K.; Jia, X.M.; McCouch, S.R.; Kochian, L.V. Natural variation underlies alterations in Nramp aluminum transporter (NRAT1) expression and function that play a key role in rice aluminum tolerance. Proc. Natl. Acad. Sci. USA 2014, 111, 6503–6508. [Google Scholar] [CrossRef]
  27. Chen, S.F.; Zhou, Y.Q.; Chen, Y.R.; Gu, J. FASTP: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, 884–890. [Google Scholar] [CrossRef] [PubMed]
  28. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef]
  29. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef]
  30. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef]
  31. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  33. Xie, C.; Mao, X.; Huang, J.; Ding, Y.; Wu, J.; Dong, S.; Kong, L.; Gao, G.; Li, C.-Y.; Wei, L. KOBAS 2.0: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011, 39, W316–W322. [Google Scholar] [CrossRef]
  34. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using realtime quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  35. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [PubMed]
  36. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape stringapp: Network analysis and visualization of proteomics data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef] [PubMed]
  37. Ryan, P.R.; Dong, D.; Teuber, F.; Wendler, N.; Muehling, K.H.; Liu, J.; Xu, M.; Salvador Moreno, N.; You, J.; Maurer, H.-P.; et al. Assessing how the aluminum-resistance traits in wheat and rye transfer to hexaploid and octoploid triticale. Front. Plant Sci. 2018, 9, 1334. [Google Scholar] [CrossRef] [PubMed]
  38. Rocha, J.; Ciceron, F.; de Sanctis, D.; Lelimousin, M.; Chazalet, V.; Lerouxel, O.; Breton, C. Structure of Arabidopsis thaliana FUT1 reveals a variant of the GT-B class fold and provides insight into xyloglucan fucosylation. Plant Cell 2016, 28, 2352–2364. [Google Scholar] [CrossRef]
  39. Yang, J.L.; Zhu, X.F.; Peng, Y.X.; Zheng, C.; Li, G.X.; Liu, Y.; Shi, Y.Z.; Zheng, S.J. Cell wall hemicellulose contributes significantly to aluminum adsorption and root growth in Arabidopsis. Plant Physiol. 2011, 155, 1885–1892. [Google Scholar] [CrossRef]
  40. Becnel, J.; Natarajan, M.; Kipp, A.; Braam, J. Developmental expression patterns of Arabidopsis XTH genes reported by transgenes and genevestigator. Plant Mol. Biol. 2006, 61, 451–467. [Google Scholar] [CrossRef]
  41. Yokoyama, R.; Nishitani, K. A comprehensive expression analysis of all members of a gene family encoding cell-wall enzymes allowed us to predict cis-regulatory regions involved in cell-wall construction in specific organs of Arabidopsis. Plant Cell Physiol. 2001, 42, 1025–1033. [Google Scholar] [CrossRef] [PubMed]
  42. Tao, Y.; Wan, J.X.; Liu, Y.S.; Yang, X.Z.; Shen, R.F.; Zhu, X.F. The NAC transcription factor ANAC017 regulates aluminum tolerance by regulating the cell wall-modifying genes. Plant Physiol. 2022, 189, 2517–2534. [Google Scholar] [CrossRef]
  43. Tsuyama, T.; Kawai, R.; Shitan, N.; Matoh, T.; Sugiyama, J.; Yoshinaga, A.; Takabe, K.; Fujita, M.; Yazaki, K. Proton-dependent coniferin transport, a common major transport event in differentiating xylem tissue of woody plants. Plant Physiol. 2013, 162, 918–926. [Google Scholar] [CrossRef] [PubMed]
  44. Huang, C.F.; Yamaji, N.; Mitani, N.; Yano, M.; Nagamura, Y.; Ma, J.F. A bacterial-type ABC transporter is involved in aluminum tolerance in rice. Plant Cell 2009, 21, 655–667. [Google Scholar] [CrossRef] [PubMed]
  45. Larsen, P.B.; Geisler, M.J.; Jones, C.A.; Williams, K.M.; Cancel, J.D. ALS3 encodes a phloem-localized ABC transporter-like protein that is required for aluminum tolerance in Arabidopsis. Plant J. 2005, 41, 353–363. [Google Scholar] [CrossRef] [PubMed]
  46. Seregin, I.V.; Kozhevnikova, A.D. Nicotianamine: A key player in metal homeostasis and hyperaccumulation in plants. Int. J. Mol. Sci. 2023, 24, 10822. [Google Scholar] [CrossRef] [PubMed]
  47. Gautam, T.; Jan, I.; Batra, R.; Singh, K.; Pandey, R.; Sharma, P.K.; Balyan, H.S.; Gupta, P.K. Further studies on nicotianamine aminotransferase (NAAT) genes involved in biofortification in bread wheat (Triticum aestivum L.). Plant Gene 2023, 33, 100389. [Google Scholar] [CrossRef]
  48. Araki, R.; Murata, J.; Murata, Y. A novel barley yellow stripe 1-like transporter (HvYSL2) localized to the root endodermis transports metal-phytosiderophore complexes. Plant. Cell Physiol. 2011, 52, 1931–1940. [Google Scholar] [CrossRef]
  49. Sasaki, A.; Yamaji, N.; Xia, J.; Ma, J.F. OsYSL6 is involved in the detoxification of excess manganese in rice. Plant Physiol. 2011, 157, 1832–1840. [Google Scholar] [CrossRef]
  50. Ishimaru, Y.; Masuda, H.; Bashir, K.; Inoue, H.; Tsukamoto, T.; Takahashi, M.; Nakanishi, H.; Aoki, N.; Hirose, T.; Ohsugi, R.; et al. Rice metal-nicotianamine transporter, OsYSL2, is required for the long-distance transport of iron and manganese. Plant J. 2010, 62, 379–390. [Google Scholar] [CrossRef]
  51. Xia, J.; Yamaji, N.; Kasai, T.; Ma, J.F. Plasma membrane-localized transporter for aluminum in rice. Proc. Natl. Acad. Sci. USA 2010, 107, 18381–18385. [Google Scholar] [CrossRef] [PubMed]
  52. Huang, C.F.; Yamaji, N.; Chen, Z.; Ma, J.F. A tonoplast-localized half-size ABC transporter is required for internal detoxification of aluminum in rice. Plant J. 2012, 69, 857–867. [Google Scholar] [CrossRef] [PubMed]
  53. Hao, J.; Peng, A.; Li, Y.; Zuo, H.; Li, P.; Wang, J.; Yu, K.; Liu, C.; Zhao, S.; Wan, X.; et al. Tea plant roots respond to aluminum-induced mineral nutrient imbalances by transcriptional regulation of multiple cation and anion transporters. BMC Plant Biol. 2022, 22, 203. [Google Scholar] [CrossRef] [PubMed]
  54. Lei, G.J.; Yokosho, K.; Yamaji, N.; Fujii-Kashino, M.; Ma, J.F. Functional characterization of two half-size ABC transporter genes in aluminium-accumulating buckwheat. New Phytol. 2017, 215, 1080–1089. [Google Scholar] [CrossRef] [PubMed]
  55. Larsen, P.B.; Cancel, J.; Rounds, M.; Ochoa, V. Arabidopsis ALS1 encodes a root tip and stele localized half type ABC transporter required for root growth in an aluminum toxic environment. Planta 2007, 225, 1447–1458. [Google Scholar] [CrossRef] [PubMed]
  56. Shao, J.F.; Yamaji, N.; Shen, R.F.; Ma, J.F. The key to vn homeostasis in plants: Regulation of Mn transporters. Trends. Plant Sci. 2017, 22, 215–224. [Google Scholar] [CrossRef] [PubMed]
  57. Milner, M.J.; Seamon, J.; Craft, E.; Kochian, L.V. Transport properties of members of the ZIP family in plants and their role in Zn and Mn homeostasis. J. Exp. Bot. 2013, 64, 369–381. [Google Scholar] [CrossRef]
  58. Ishimaru, Y.; Bashir, K.; Nakanishi, H.; Nishizawa, N.K. OsNRAMP5, a major player for constitutive iron and manganese uptake in rice. Plant Signal. Behav. 2012, 7, 763–766. [Google Scholar] [CrossRef]
  59. Quinteiro Ribeiro, M.A.; de Almeida, A.A.F.; Mielke, M.S.; Gomes, F.P.; Pires, M.V.; Baligar, V.C. Aluminum effects on growth, photosynthesis, and mineral nutrition of cacao genotypes. J. Plant Nutr. 2013, 36, 1161–1179. [Google Scholar] [CrossRef]
  60. Boursiac, Y.; Lee, S.M.; Romanowsky, S.; Blank, R.; Sladek, C.; Chung, W.S.; Harper, J.F. Disruption of the vacuolar calcium-ATPases in Arabidopsis results in the activation of a salicylic acid-dependent programmed cell death pathway. Plant Physiol. 2010, 154, 1158–1171. [Google Scholar] [CrossRef]
  61. Bose, J.; Babourina, O.; Rengel, Z. Role of magnesium in alleviation of aluminium toxicity in plants. J. Exp. Bot. 2011, 62, 2251–2264. [Google Scholar] [CrossRef] [PubMed]
  62. Ma, Q.; Feng, Y.; Luo, S.; Cheng, L.; Tong, W.; Lu, X.; Li, Y.; Zhang, P. The aquaporin MePIP2;7 improves MeMGT9-mediated Mg2+ acquisition in cassava. J. Integr. Plant Biol. 2023, 65, 2349–2367. [Google Scholar] [CrossRef] [PubMed]
  63. Chen, Z.C.; Yamaji, N.; Motoyama, R.; Nagamura, Y.; Ma, J.F. Up-regulation of a magnesium transporter gene OsMGT1 is required for conferring aluminum tolerance in rice. Plant Physiol. 2012, 159, 1624–1633. [Google Scholar] [CrossRef] [PubMed]
  64. Kim, S.A.; Punshon, T.; Lanzirotti, A.; Li, L.; Alonso, J.M.; Ecker, J.R.; Kaplan, J.; Guerinot, M.L. Localization of iron in Arabidopsis seed requires the vacuolar membrane transporter VIT1. Science 2006, 314, 1295–1298. [Google Scholar] [CrossRef] [PubMed]
  65. Pittman, J.K. Managing the manganese: Molecular mechanisms of manganese transport and homeostasis. New Phytol. 2005, 167, 733–742. [Google Scholar] [CrossRef] [PubMed]
  66. Sharma, A.; Shahzad, B.; Rehman, A.; Bhardwaj, R.; Landi, M.; Zheng, B. Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress. Molecules 2019, 24, 2452. [Google Scholar] [CrossRef]
  67. Izbianska, K.; Arasimowicz-Jelonek, M.; Deckert, J. Phenylpropanoid pathway metabolites promote tolerance response of lupine roots to lead stress. Ecotoxicol. Environ. Saf. 2014, 110, 61–67. [Google Scholar] [CrossRef]
  68. Feduraev, P.; Skrypnik, L.; Riabova, A.; Pungin, A.; Tokupova, E.; Maslennikov, P.; Chupakhina, G. Phenylalanine and tyrosine as exogenous precursors of wheat (Triticum aestivum L.) secondary metabolism through PAL-associated pathways. Plants 2020, 9, 476. [Google Scholar] [CrossRef]
  69. Tang, T.; Tao, F.; Li, W. Characterisation of manganese toxicity tolerance in Arabis paniculata. Plant Divers. 2021, 43, 163–172. [Google Scholar] [CrossRef]
  70. Lei, Y.; Korpelainen, H.; Li, C. Physiological and biochemical responses to high Mn concentrations in two contrasting Populus cathayana populations. Chemosphere 2007, 68, 686–694. [Google Scholar] [CrossRef]
  71. Liu, P.; Huang, R.; Hu, X.; Jia, Y.; Li, J.; Luo, J.; Liu, Q.; Luo, L.; Liu, G.; Chen, Z. Physiological responses and proteomic changes reveal insights into Stylosanthes response to manganese toxicity. BMC Plant Biol. 2019, 19, 212. [Google Scholar] [CrossRef] [PubMed]
  72. Zhou, P.; Su, L.; Lv, A.; Wang, S.; Huang, B.; An, Y. Gene expression analysis of alfalfa seedlings response to acid-aluminum. Int. J. Genom. 2016, 2016, 2095195. [Google Scholar] [CrossRef] [PubMed]
  73. Chen, Y.; Huang, L.; Liang, X.; Dai, P.; Zhang, Y.; Li, B.; Lin, X.; Sun, C. Enhancement of polyphenolic metabolism as an adaptive response of lettuce (Lactuca sativa) roots to aluminum stress. Environ. Pollut. 2020, 261, 114230. [Google Scholar] [CrossRef] [PubMed]
  74. Ezaki, B.; Higashi, A.; Nanba, N.; Nishiuchi, T. An S-adenosyl methionine synthetase (SAMS) gene from Andropogon virginicus L. confers aluminum stress tolerance and facilitates epigenetic gene regulation in Arabidopsis thaliana. Front. Plant Sci. 2016, 7, 1627. [Google Scholar] [CrossRef] [PubMed]
  75. Su, L.; Lv, A.; Wen, W.; Fan, N.; Li, J.; Gao, L.; Zhou, P.; An, Y. MsMYB741 is involved in alfalfa resistance to aluminum stress by regulating flavonoid biosynthesis. Plant J. 2022, 112, 756–771. [Google Scholar] [CrossRef] [PubMed]
  76. Choudhury, S.; Mazumder, M.K.; Moulick, D.; Sharma, P.; Tata, S.K.; Ghosh, D.; Ali, H.M.; Siddiqui, M.H.; Brestic, M.; Skalicky, M.; et al. A computational study of the role of secondary metabolites for mitigation of acid soil stress in cereals using dehydroascorbate and mono-dehydroascorbate reductases. Antioxidants 2022, 11, 458. [Google Scholar] [CrossRef] [PubMed]
  77. Jiang, P.; Liu, J.; Yu, G.; Lei, L.; Jiang, X. Moderate Mn accumulation enhances growth and alters leaf hormone contents in the hyperaccumulator Celosia argentea Linn. Environ. Exp. Bot. 2021, 191, 104603. [Google Scholar] [CrossRef]
  78. Wu, Q.; Lin, X.; Li, S.; Liang, Z.; Wang, H.; Tang, T. Endophytic Bacillus sp. AP10 harboured in Arabis paniculata mediates plant growth promotion and manganese detoxification. Ecotoxicol. Environ. Saf. 2023, 262, 115170. [Google Scholar] [CrossRef]
  79. Zhu, X.; Wang, P.; Bai, Z.; Herde, M.; Ma, Y.; Li, N.; Liu, S.; Huang, C.F.; Cui, R.; Ma, H.; et al. Calmodulin-like protein CML24 interacts with CAMTA2 and WRKY46 to regulate ALMT1-dependent Al resistance in Arabidopsis thaliana. New Phytol. 2022, 233, 2471–2487. [Google Scholar] [CrossRef]
  80. Zhou, F.; Singh, S.; Zhang, J.; Fang, Q.; Li, C.; Wang, J.; Zhao, C.; Wang, P.; Huang, C.F. The MEKK1-MKK1/2-MPK4 cascade phosphorylates and stabilizes STOP1 to confer aluminum resistance in Arabidopsis. Mol. Plant 2023, 16, 337–353. [Google Scholar] [CrossRef]
  81. Bataller, S.; Villacastin, A.J.; Shen, Q.J.; Bergman, C. Comparative nutritional assessment and metabolomics of a WRKY rice mutant with enhanced germination rates. Agronomy 2023, 13, 1149. [Google Scholar] [CrossRef]
  82. Faria-Lopes, R.O.; Muniz, D.R.; Chaves, I.S.; Franca, M.G.C.; Modolo, L.V. Nitric oxide precursors prevent Al-triggered auxin flow inhibition in Triticum aestivum roots. J. Adv. Res. 2019, 15, 27–36. [Google Scholar] [CrossRef] [PubMed]
  83. Sun, C.; Liu, L.; Yu, Y.; Liu, W.; Lu, L.; Jin, C.; Lin, X. Nitric oxide alleviates aluminum-induced oxidative damage through regulating the ascorbate-glutathione cycle in roots of wheat. J. Integr. Plant Biol. 2015, 57, 550–561. [Google Scholar] [CrossRef] [PubMed]
  84. Srivastava, S.; Dubey, R.S. Nitric oxide alleviates manganese toxicity by preventing oxidative stress in excised rice leaves. Acta Physiol. Plant. 2012, 34, 819–825. [Google Scholar] [CrossRef]
  85. Thakro, V.; Malik, N.; Basu, U.; Srivastava, R.; Narnoliya, L.; Daware, A.; Varshney, N.; Mohanty, J.K.; Bajaj, D.; Dwivedi, V.; et al. A superior gene allele involved in abscisic acid signaling enhances drought tolerance and yield in chickpea. Plant Physiol. 2023, 191, 1884–1912. [Google Scholar] [CrossRef] [PubMed]
  86. Melo, J.O.; Martins, L.G.C.; Barros, B.A.; Pimenta, M.R.; Lana, U.G.P.; Duarte, C.E.M.; Pastina, M.M.; Guimaraes, C.T.; Schaffert, R.E.; Kochian, L.V.; et al. Repeat variants for the SbMATE transporter protect sorghum roots from aluminum toxicity by transcriptional interplay in cis and trans. Proc. Natl. Acad. Sci. USA 2019, 116, 313–318. [Google Scholar] [CrossRef] [PubMed]
  87. Barros, V.A.; Chandnani, R.; de Sousa, S.M.; Maciel, L.S.; Tokizawa, M.; Guimaraes, C.T.; Magalhaes, J.V.; Kochian, L.V. Root adaptation via common genetic factors conditioning tolerance to multiple stresses for crops cultivated on acidic tropical soils. Front. Plant Sci. 2020, 11, 565339. [Google Scholar] [CrossRef]
  88. Ding, S.; Zhang, H.; Zhou, C.; Bao, Y.; Xu, X.; Chen, Y.; Shen, Z.; Chen, C. Transcriptomic, epigenomic and physiological comparisons reveal key factors for different manganese tolerances in three Chenopodium ambrosioides L. populations. Plant Physiol. Biochem. 2023, 201, 107883. [Google Scholar] [CrossRef]
  89. Iuchi, S.; Koyama, H.; Iuchi, A.; Kobayashi, Y.; Kitabayashi, S.; Kobayashi, Y.; Ikka, T.; Hirayama, T.; Shinozaki, K.; Kobayashi, M. Zinc finger protein STOP1 is critical for proton tolerance in Arabidopsis and coregulates a key gene in aluminum tolerance. Proc. Natl. Acad. Sci. USA 2007, 104, 9900–9905. [Google Scholar] [CrossRef]
  90. Yamaji, N.; Huang, C.F.; Nagao, S.; Yano, M.; Sato, Y.; Nagamura, Y.; Ma, J.F. A zinc finger transcription factor ART1 regulates multiple genes implicated in aluminum tolerance in rice. Plant Cell 2009, 21, 3339–3349. [Google Scholar] [CrossRef]
  91. Li, L.; Li, X.; Yang, C.; Cheng, Y.; Cai, Z.; Nian, H.; Ma, Q. GsERF1 enhances Arabidopsis thaliana aluminum tolerance through an ethylene-mediated pathway. BMC Plant Biol. 2022, 22, 258. [Google Scholar] [CrossRef] [PubMed]
  92. Sun, K.; Wang, H.; Xia, Z. The maize bHLH transcription factor bHLH105 confers manganese tolerance in transgenic tobacco. Plant Sci. 2019, 280, 97–109. [Google Scholar] [CrossRef] [PubMed]
  93. Kobayashi, Y.; Ohyama, Y.; Kobayashi, Y.; Ito, H.; Iuchi, S.; Fujita, M.; Zhao, C.R.; Tanveer, T.; Ganesan, M.; Kobayashi, M.; et al. STOP2 activates transcription of several genes for Al- and low pH-tolerance that are regulated by STOP1 in Arabidopsis. Mol. Plant 2014, 7, 311–322. [Google Scholar] [CrossRef] [PubMed]
  94. Sawaki, Y.; Iuchi, S.; Kobayashi, Y.; Kobayashi, Y.; Ikka, T.; Sakurai, N.; Fujita, M.; Shinozaki, K.; Shibata, D.; Kobayashi, M.; et al. STOP1 regulates multiple genes that protect Arabidopsis from proton and aluminum toxicities. Plant Physiol. 2009, 150, 281–294. [Google Scholar] [CrossRef]
  95. Jia, Y.; Li, X.; Liu, Q.; Hu, X.; Li, J.; Dong, R.; Liu, P.; Liu, G.; Luo, L.; Chen, Z. Physiological and transcriptomic analyses reveal the roles of secondary metabolism in the adaptive responses of Stylosanthes to manganese toxicity. BMC Genom. 2020, 21, 861. [Google Scholar] [CrossRef]
  96. Ding, Z.J.; Yan, J.Y.; Xu, X.Y.; Li, G.X.; Zheng, S.J. WRKY46 functions as a transcriptional repressor of ALMT1, regulating aluminum-induced malate secretion in Arabidopsis. Plant J. 2013, 76, 825–835. [Google Scholar] [CrossRef]
Figure 1. Relative root length of wheat varieties under Al and Mn stresses. Data are expressed as the mean ± standard deviation (SD). Histogram columns within a cluster labeled with different lowercases were significantly different (p < 0.05).
Figure 1. Relative root length of wheat varieties under Al and Mn stresses. Data are expressed as the mean ± standard deviation (SD). Histogram columns within a cluster labeled with different lowercases were significantly different (p < 0.05).
Cimb 46 00024 g001
Figure 2. (A) Chlorophyll level under Al and Mn stresses. (B) Leaf brown spots (labeled in red circles) under Mn stress. Data are expressed as the mean ± standard deviation (SD). Histogram columns within a cluster labeled with different lowercases were significantly different (p < 0.05).
Figure 2. (A) Chlorophyll level under Al and Mn stresses. (B) Leaf brown spots (labeled in red circles) under Mn stress. Data are expressed as the mean ± standard deviation (SD). Histogram columns within a cluster labeled with different lowercases were significantly different (p < 0.05).
Cimb 46 00024 g002
Figure 3. Al and Mn contents and subcellular distribution in root cell. Data are expressed as the mean ± standard deviation (SD). Histogram columns within a cluster labeled with different lowercases were significantly different (p < 0.05).
Figure 3. Al and Mn contents and subcellular distribution in root cell. Data are expressed as the mean ± standard deviation (SD). Histogram columns within a cluster labeled with different lowercases were significantly different (p < 0.05).
Cimb 46 00024 g003
Figure 4. Nutrient element contents in wheat plants under Al and Mn toxicity. Data are expressed as the mean ± standard deviation (SD). Histogram columns labeled with different lowercase letters are significantly different (p < 0.05).
Figure 4. Nutrient element contents in wheat plants under Al and Mn toxicity. Data are expressed as the mean ± standard deviation (SD). Histogram columns labeled with different lowercase letters are significantly different (p < 0.05).
Cimb 46 00024 g004
Figure 5. Correlation analysis of differentially expressed genes’ (DEGs) expressions between transcriptome data and qRT-PCR results. *** indicates significance at p < 0.001, every black dot represents one gene expressed in the same sample determined by RNA-seq and qRT-PCR methods.
Figure 5. Correlation analysis of differentially expressed genes’ (DEGs) expressions between transcriptome data and qRT-PCR results. *** indicates significance at p < 0.001, every black dot represents one gene expressed in the same sample determined by RNA-seq and qRT-PCR methods.
Cimb 46 00024 g005
Figure 6. Al and Mn responsive gene expression profiles. The volcanic map of DEGs under Al and Mn stresses, respectively (A,B). Venn diagrams converge the DEGs under Al and Mn stresses (C,D).
Figure 6. Al and Mn responsive gene expression profiles. The volcanic map of DEGs under Al and Mn stresses, respectively (A,B). Venn diagrams converge the DEGs under Al and Mn stresses (C,D).
Cimb 46 00024 g006
Figure 7. GO annotation analysis of four types of DEGs under Al and Mn stress. Type I_DEGs: genes expressed in the same direction under Al and Mn stresses; Type II_DEGs: 7 genes were expressed in opposite directions but were not annotated; Type III_DEGs: genes exclusively responded to Al; Type IV_DEGs: genes exclusively responded to Mn; MF: molecular function; CC: catalytic activity; BP: biological process.
Figure 7. GO annotation analysis of four types of DEGs under Al and Mn stress. Type I_DEGs: genes expressed in the same direction under Al and Mn stresses; Type II_DEGs: 7 genes were expressed in opposite directions but were not annotated; Type III_DEGs: genes exclusively responded to Al; Type IV_DEGs: genes exclusively responded to Mn; MF: molecular function; CC: catalytic activity; BP: biological process.
Cimb 46 00024 g007
Figure 8. KEGG enrichment for eight DEGs set under Al and Mn stresses. The bubble size indicates the number of genes involved in the terms and pathways. DEGs set under Al_up/Mn_down and Al_down/Mn_up were not enriched into any pathways. Note: Al (Mn)_up: upregulated by Al (Mn); Al (Mn)_down: downregulated by Al (Mn); Al (Mn)_ns: did not respond to Al (Mn); M: metabolism; GIP: genetic information processing; EIP: environmental information processing; OS: organismal systems; CP: cellular process.
Figure 8. KEGG enrichment for eight DEGs set under Al and Mn stresses. The bubble size indicates the number of genes involved in the terms and pathways. DEGs set under Al_up/Mn_down and Al_down/Mn_up were not enriched into any pathways. Note: Al (Mn)_up: upregulated by Al (Mn); Al (Mn)_down: downregulated by Al (Mn); Al (Mn)_ns: did not respond to Al (Mn); M: metabolism; GIP: genetic information processing; EIP: environmental information processing; OS: organismal systems; CP: cellular process.
Cimb 46 00024 g008
Figure 9. Correlation between tolerance, toxicity, and nutrient ion contents under Al (left bottom) and Mn (right upper) stresses. The parameters in yellow line (from left upper to right bottom) are as follow: relative root length, SPAD, Al or Mn contents in root, shoot, total plant, ratio of root to shoot, cell wall and cell sap, and nutrients (Ca, Mg, Fe, and Zn) content in total plant, respectively.
Figure 9. Correlation between tolerance, toxicity, and nutrient ion contents under Al (left bottom) and Mn (right upper) stresses. The parameters in yellow line (from left upper to right bottom) are as follow: relative root length, SPAD, Al or Mn contents in root, shoot, total plant, ratio of root to shoot, cell wall and cell sap, and nutrients (Ca, Mg, Fe, and Zn) content in total plant, respectively.
Cimb 46 00024 g009
Figure 10. Correlation between gene expression module and physiological trait. (A) Hierarchical clustering tree showing co-expression modules. Each leaf in the tree represents one gene. The major tree branches constitute the modules labeled with different colors. (B) Matrix of module–trait correlation. A total of 18 modules are shown on the left, and the relevance color scale from −1 to 1 is shown on the right. The high correlation between a specific module and the sample is shown in dark red or dark green. The numbers in parentheses represent the significance (p), and the numbers above represent the correlation coefficient (r).
Figure 10. Correlation between gene expression module and physiological trait. (A) Hierarchical clustering tree showing co-expression modules. Each leaf in the tree represents one gene. The major tree branches constitute the modules labeled with different colors. (B) Matrix of module–trait correlation. A total of 18 modules are shown on the left, and the relevance color scale from −1 to 1 is shown on the right. The high correlation between a specific module and the sample is shown in dark red or dark green. The numbers in parentheses represent the significance (p), and the numbers above represent the correlation coefficient (r).
Cimb 46 00024 g010
Figure 11. Visualization of gene co-expression network in five modules. Each node represents a gene, and the nodes are connected by lines. The genes at both ends of the lines have similar expression patterns. The top five genes with the highest connectivity were identified as candidate hub genes in each module.
Figure 11. Visualization of gene co-expression network in five modules. Each node represents a gene, and the nodes are connected by lines. The genes at both ends of the lines have similar expression patterns. The top five genes with the highest connectivity were identified as candidate hub genes in each module.
Cimb 46 00024 g011
Figure 12. The expression levels of DEGs involved in cell wall biogenesis and macromolecule metabolism under Al and Mn stresses. PAE, PME: Pectin methylesterase; XTH: Glycan synthases and glycosyltransferases; NAC: NAM, ATAF1/2, and cup-shaped cotyledon-2.
Figure 12. The expression levels of DEGs involved in cell wall biogenesis and macromolecule metabolism under Al and Mn stresses. PAE, PME: Pectin methylesterase; XTH: Glycan synthases and glycosyltransferases; NAC: NAM, ATAF1/2, and cup-shaped cotyledon-2.
Cimb 46 00024 g012
Figure 13. The expression levels of DEGs involved in nicotianamine synthesis under Al and Mn stresses. NAS: nicotianamine synthase; NAAT: nicotianamine aminotransferase.
Figure 13. The expression levels of DEGs involved in nicotianamine synthesis under Al and Mn stresses. NAS: nicotianamine synthase; NAAT: nicotianamine aminotransferase.
Cimb 46 00024 g013
Figure 14. The expression levels of DEGs involved in metal ion transport under Al and Mn stresses. ACA: calmodulin-activated Ca2+ pumps; CCX: cation/calcium exchanger; MGT: magnesium transporter; PIP: aquaporin; IRT: iron-regulated transporter; VIT: vacuole iron transporter; YSL: yellow stripe-like protein; MRP: multidrug resistance-associated protein; MTP: metal tolerance protein; ZIP: zinc/iron-regulated transporter protein; HMA: heavy-metal-associated domain.
Figure 14. The expression levels of DEGs involved in metal ion transport under Al and Mn stresses. ACA: calmodulin-activated Ca2+ pumps; CCX: cation/calcium exchanger; MGT: magnesium transporter; PIP: aquaporin; IRT: iron-regulated transporter; VIT: vacuole iron transporter; YSL: yellow stripe-like protein; MRP: multidrug resistance-associated protein; MTP: metal tolerance protein; ZIP: zinc/iron-regulated transporter protein; HMA: heavy-metal-associated domain.
Cimb 46 00024 g014
Figure 15. The expression levels of DEGs involved in phenylpropanoid metabolism under Al and Mn stresses. The gene expression values of different samples or experimental conditions were normalized and scale-mapped to [−1,1] from minimum (green) to maximum (red). PAL: Phenylalanine ammonia-lyase; 4CL: 4-coumarate--CoA ligase; CCR: Cinnamoyl-CoA reductase; BGLG: Beta-glucosidase; CYP: Cytochrome P450; CAD: Cinnamyl alcohol dehydrogenase; TOGT: Anthocyanin 3′-O-beta-glucosyltransferase; PHT: Putrescine hydroxycinnamoyl transferase; CSE: Caffeoylshikimate esterase; COMT: Tricetin 3′,4′,5′-O-trimethyltransferase; SCPL: Serine carboxypeptidase-like protein.
Figure 15. The expression levels of DEGs involved in phenylpropanoid metabolism under Al and Mn stresses. The gene expression values of different samples or experimental conditions were normalized and scale-mapped to [−1,1] from minimum (green) to maximum (red). PAL: Phenylalanine ammonia-lyase; 4CL: 4-coumarate--CoA ligase; CCR: Cinnamoyl-CoA reductase; BGLG: Beta-glucosidase; CYP: Cytochrome P450; CAD: Cinnamyl alcohol dehydrogenase; TOGT: Anthocyanin 3′-O-beta-glucosyltransferase; PHT: Putrescine hydroxycinnamoyl transferase; CSE: Caffeoylshikimate esterase; COMT: Tricetin 3′,4′,5′-O-trimethyltransferase; SCPL: Serine carboxypeptidase-like protein.
Cimb 46 00024 g015
Figure 16. The expression levels of DEGs involved in signaling under Al and Mn stresses. ARF: auxin-responsive factor; GH: indole-3-acetic acid-amido synthetase GH; ERF: ethylene-responsive transcription factor; PYL: abscisic acid receptor PYL; JAZ: jasmonate ZIM domain protein; BR: brassinosteroid; IPT: adenylate isopentenyltransferase; PTR: protein NRT1/PTR family; CML: calmodulin-like protein; M3K, M2K, MEKK: mitogen-activated protein kinase; PR: pathogenesis-related protein; RBOH: respiratory burst oxidase homolog protein.
Figure 16. The expression levels of DEGs involved in signaling under Al and Mn stresses. ARF: auxin-responsive factor; GH: indole-3-acetic acid-amido synthetase GH; ERF: ethylene-responsive transcription factor; PYL: abscisic acid receptor PYL; JAZ: jasmonate ZIM domain protein; BR: brassinosteroid; IPT: adenylate isopentenyltransferase; PTR: protein NRT1/PTR family; CML: calmodulin-like protein; M3K, M2K, MEKK: mitogen-activated protein kinase; PR: pathogenesis-related protein; RBOH: respiratory burst oxidase homolog protein.
Cimb 46 00024 g016
Table 1. Connectivity, expressions, and annotations of candidate hub genes’ correlation between ME and traits.
Table 1. Connectivity, expressions, and annotations of candidate hub genes’ correlation between ME and traits.
ModuleGene_idConnectivityKme
Value
Fold ChangeGene Description
Al/CKMn/CK
midnightblueTraesCS5A02G55200050.1460.9860.892124.101Nicotianamine synthase
TraesCS3B02G06850050.0640.9860.830329.475Nicotianamine synthase
TraesCS4D02G18490050.0320.9940.69453.203Nicotianamine synthase
TraesCS4B02G18390049.7290.9870.429178.986Nicotianamine synthase
TraesCS2A02G04990049.6450.9850.486343.721Nicotianamine synthase
blackTraesCS3A02G129000109.3350.9752.3301.080ABC transporter C family member 3-like
TraesCS2B02G323600105.7430.9693.9241.598RING-type E3 ubiquitin transferase; zinc finger of C3HC4-type
TraesCS3A02G231600105.4720.9706.4591.459Homeobox-leucine zipper protein HOX3
TraesCS3D02G129900104.6210.9662.7111.007ABC subfamily C transporter
TraesCS4B02G307300104.4570.9697.4541.341Intracellular trafficking, secretion, and vesicular transport
lightcyanTraesCS2D02G07050033.8910.9671.5970.592Inorganic ion transport; ascorbate peroxidase; glutathione metabolism
TraesCS2B02G08740032.2810.9551.5080.604Inorganic ion transport; ascorbate peroxidase; glutathione metabolism
TraesCS3A02G49300030.6100.9441.4970.660Function unknown hypothetical protein TRIUR3
TraesCS7D02G25830029.4170.9431.9210.560Protein processing in endoplasmic reticulum; hydrolase
TraesCS7A02G25750029.3760.9411.7970.594TCP family transcription factor; hydrolase
redTraesCS4D02G214100226.1240.9825.4404.174Positive regulation of mRNA splicing via spliceosome
TraesCS5A02G548000225.2900.98040.42235.575BAG family molecular chaperone regulator 6; IQ calmodulin-binding motif
TraesCS4B02G208900225.0950.97514.93212.665Nucleotide exchange factor Fes1B
TraesCS4D02G145500223.5120.965403.199310.521Hsp20/alpha crystallin family
TraesCS6B02G374100222.5310.97081.00457.046Hsp20/alpha crystallin family
turquoiseTraesCS7D02G369700748.1960.9890.3010.412Peroxidase P7; phenylpropanoid biosynthesis
TraesCS3B02G280000728.8920.9850.4340.626Proton-transporting V-type ATPase
TraesCS1D02G404600727.5110.9840.3460.525Plant-type secondary cell wall biogenesis; Fasciclin-like arabinogalactan protein
TraesCS7A02G453100727.0240.9800.2600.448Function unknown
TraesCS7A02G492400722.9600.9860.3070.517Auxin efflux carrier component; intracellular trafficking, secretion, and vesicular transport
Table 2. Statistics of differentially expressed transcription factors under Al and Mn stresses.
Table 2. Statistics of differentially expressed transcription factors under Al and Mn stresses.
TFs FamilyAl_upAl_downAl_nrTotal
Mn_upMn_nrMn_downMn_nrMn_upMn_down
MYB30309241413120
WRKY352103541114
AP2163401220385
bHLH7611821457
NAC819434745
HSF255022034
bZIP017511933
HB216120324
GRAS57038023
B336061117
LBD33033517
Dof24022111
MIKC34002110
C2H24302009
SRS6300009
DBB2021016
GATA0202015
M_type0201104
ARF1300004
ZF-HD3000003
GRF0002002
BES10100001
HD-ZIP0100001
RAV0100001
EIL0000101
TCP0100001
CAMTA0000101
FAR11000001
NF-YA0100001
Total156190327713550640
Notes: Al (Mn)_up: upregulated by Al (Mn); Al (Mn)_down: downregulated by Al (Mn); Al (Mn)_ns: did not respond to Al (Mn).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, D.; Xian, C.; Zhang, W.; Qin, Y.; Li, Q.; Usman, M.; Sun, S.; Xing, Y.; Dong, D. Physiological and Transcriptomic Analyses Reveal Commonalities and Specificities in Wheat in Response to Aluminum and Manganese. Curr. Issues Mol. Biol. 2024, 46, 367-397. https://doi.org/10.3390/cimb46010024

AMA Style

Luo D, Xian C, Zhang W, Qin Y, Li Q, Usman M, Sun S, Xing Y, Dong D. Physiological and Transcriptomic Analyses Reveal Commonalities and Specificities in Wheat in Response to Aluminum and Manganese. Current Issues in Molecular Biology. 2024; 46(1):367-397. https://doi.org/10.3390/cimb46010024

Chicago/Turabian Style

Luo, Daozhen, Chunnuan Xian, Wenjie Zhang, Ying Qin, Qing Li, Muhammad Usman, Shiheng Sun, Yongxiu Xing, and Dengfeng Dong. 2024. "Physiological and Transcriptomic Analyses Reveal Commonalities and Specificities in Wheat in Response to Aluminum and Manganese" Current Issues in Molecular Biology 46, no. 1: 367-397. https://doi.org/10.3390/cimb46010024

APA Style

Luo, D., Xian, C., Zhang, W., Qin, Y., Li, Q., Usman, M., Sun, S., Xing, Y., & Dong, D. (2024). Physiological and Transcriptomic Analyses Reveal Commonalities and Specificities in Wheat in Response to Aluminum and Manganese. Current Issues in Molecular Biology, 46(1), 367-397. https://doi.org/10.3390/cimb46010024

Article Metrics

Back to TopTop