Genome-Wide Identification and Expression Analysis of Metal Tolerance Protein Gene Family in Medicago truncatula Under a Broad Range of Heavy Metal Stress

Transcriptomics, proteomics, and metabolomics interventions prompt crop improvement against metal(loid) toxicity

The escalating challenges posed by metal(loid) toxicity in agricultural ecosystems, exacerbated by rapid climate change and anthropogenic pressures, demand urgent attention. Soil contamination is a critical issue because it significantly impacts crop productivity. The widespread threat of metal(loid) toxicity can jeopardize global food security due to contaminated food supplies and pose environmental risks, contributing to soil and water pollution and thus impacting the whole ecosystem. In this context, plants have evolved complex mechanisms to combat metal(loid) stress. Amid the array of innovative approaches, omics, notably transcriptomics, proteomics, and metabolomics, have emerged as transformative tools, shedding light on the genes, proteins, and key metabolites involved in metal(loid) stress responses and tolerance mechanisms. These identified candidates hold promise for developing high-yielding crops with desirable agronomic traits. Computational biology tools like bioinformatics, biological databases, and analytical pipelines support these omics approaches by harnessing diverse information and facilitating the mapping of genotype-to-phenotype relationships under stress conditions. This review explores: (1) the multifaceted strategies that plants use to adapt to metal(loid) toxicity in their environment; (2) the latest findings in metal(loid)-mediated transcriptomics, proteomics, and metabolomics studies across various plant species; (3) the integration of omics data with artificial intelligence and high-throughput phenotyping; (4) the latest bioinformatics databases, tools and pipelines for single and/or multi-omics data integration; (5) the latest insights into stress adaptations and tolerance mechanisms for future outlooks; and (6) the capacity of omics advances for creating sustainable and resilient crop plants that can thrive in metal(loid)-contaminated environments.

Genetics and metabolic responses of Artemisia annua L to the lake of phosphorus under the sparingly soluble phosphorus fertilizer: evidence from transcriptomics analysis

The medicinal herb Artemisia annua L. is prized for its capacity to generate artemisinin, which is used to cure malaria. Potentially influencing the biomass and secondary metabolite synthesis of A. annua is plant nutrition, particularly phosphorus (P). However, most soil P exist as insoluble inorganic and organic phosphates, which results to low P availability limiting plant growth and development. Although plants have developed several adaptation strategies to low P levels, genetics and metabolic responses to P status remain largely unknown. In a controlled greenhouse experiment, the sparingly soluble P form, hydroxyapatite (Ca5OH(PO4)3/CaP) was used to simulate calcareous soils with low P availability. In contrast, the soluble P form KH2PO4/KP was used as a control. A. annua's morphological traits, growth, and artemisinin concentration were determined, and RNA sequencing was used to identify the differentially expressed genes (DEGs) under two different P forms. Total biomass, plant height, leaf number, and stem diameter, as well as leaf area, decreased by 64.83%, 27.49%, 30.47%, 38.70%, and 54.64% in CaP compared to KP; however, LC-MS tests showed an outstanding 37.97% rise in artemisinin content per unit biomass in CaP contrary to KP. Transcriptome analysis showed 2015 DEGs (1084 up-regulated and 931 down-regulated) between two P forms, including 39 transcription factor (TF) families. Further analysis showed that DEGs were mainly enriched in carbohydrate metabolism, secondary metabolites biosynthesis, enzyme catalytic activity, signal transduction, and so on, such as tricarboxylic acid (TCA) cycle, glycolysis, starch and sucrose metabolism, flavonoid biosynthesis, P metabolism, and plant hormone signal transduction. Meanwhile, several artemisinin biosynthesis genes were up-regulated, including DXS, GPPS, GGPS, MVD, and ALDH, potentially increasing artemisinin accumulation. Furthermore, 21 TF families, including WRKY, MYB, bHLH, and ERF, were up-regulated in reaction to CaP, confirming their importance in P absorption, internal P cycling, and artemisinin biosynthesis regulation. Our results will enable us to comprehend how low P availability impacts the parallel transcriptional control of plant development, growth, and artemisinin production in A. annua. This study could lay the groundwork for future research into the molecular mechanisms underlying A. annua's low P adaptation.

Integrative application of silicon and/or proline improves Sweet corn (Zea mays L. saccharata) production and antioxidant defense system under salt stress condition

Silicon (Si) and/or proline (Pro) are natural supplements that are considered to induce plants' stress tolerance against various abiotic stresses. Sweet corn (Zea mays L. saccharata) production is severely afflicted by salinity stress. Therefore, two field tests were conducted to evaluate the potential effects of Si and/or Pro (6mM) used as seed soaking (SS) and/or foliar spray (FS) on Sweet corn plant growth and yield, physio-biochemical attributes, and antioxidant defense systems grown in a saline (EC = 7.14dS m-1) soil. The Si and/or Pro significantly increased growth and yield, photosynthetic pigments, free proline, total soluble sugars (TSS), K+/Na+ratios, relative water content (RWC), membrane stability index (MSI), α-Tocopherol (α-TOC), Ascorbate (AsA), glutathione (GSH), enzymatic antioxidants activities and other anatomical features as compared to controls. In contrast, electrolytes, such as SS and/or FS under salt stress compared to controls (SS and FS using tap water) were significantly decreased. The best results were obtained when SS was combined with FS via Si or Pro. These alterations are brought about by the exogenous application of Si and/or Pro rendering these elements potentially useful in aiding sweet corn plants to acclimate successfully to saline soil.

Genome-Wide Identification and Expression Analysis of MTP (Metal Ion Transport Proteins) Genes in the Common Bean

MTP/CDF carriers, called metal ion transport proteins, act as substrates for the transmission of micronutrients such as iron (Fe), zinc (Zn), and manganese (Mn) to membrane carriers in plants. In this study, genome-wide analysis of the MTP gene family in the common bean genome, expression analysis of the PvMTP4, PvMTP5, and PvMTP12 genes after Fe and Zn treatments, and the effects of Fe and Zn applications on iron and zinc content were investigated. This study used common bean genotypes assumed to have high or low Fe and Zn accumulation ability. PvMTP genes were defined as containing conserved catalytic domains with molecular weights and protein lengths ranging from 41.35 to 91.05 kDa and from 369 to 813 amino acids (aa), respectively. As a result of the phylogenetic analysis, three main clusters containing seven subgroups were formed. In this study, the first characterization of the MTP gene family of beans was performed, and the responses of three different PvMTP genes in the Zn-CDF group to Fe and Zn applications were revealed. The obtained findings are thought to constitute pioneering resources for future research on common bean biofortification studies, plant breeding related to Fe and Zn, and the functional characterization of the MTP gene family.

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses

Metal tolerance proteins (MTPs) as Me²⁺/H⁺(K⁺) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd²⁺ and Cu²⁺ exposure might be responsible for metal translocation from roots to leaves.

RNA sequencing in Artemisia annua L explored the genetic and metabolic responses to hardly soluble aluminum phosphate treatment

Artemisia annua L. is a medicinal plant valued for its ability to produce artemisinin, a molecule used to treat malaria. Plant nutrients, especially phosphorus (P), can potentially influence plant biomass and secondary metabolite production. Our work aimed to explore the genetic and metabolic response of A. annua to hardly soluble aluminum phosphate (AlPO₄, AlP), using soluble monopotassium phosphate (KH₂PO₄, KP) as a control. Liquid chromatography-mass spectrometry (LC-MS) was used to analyze artemisinin. RNA sequencing, gene ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied to analyze the differentially expressed genes (DEGs) under poor P conditions. Results showed a significant reduction in plant growth parameters, such as plant height, stem diameter, number of leaves, leaf areas, and total biomass of A. annua. Conversely, LC-MS analysis revealed a significant increase in artemisinin concentration under the AlP compared to the KP. Transcriptome analysis revealed 762 differentially expressed genes (DEGs) between the AlP and the KP. GH3, SAUR, CRE1, and PYL, all involved in plant hormone signal transduction, showed differential expression. Furthermore, despite the downregulation of HMGR in the artemisinin biosynthesis pathway, the majority of genes (ACAT, FPS, CYP71AV1, and ALDH1) were upregulated, resulting in increased artemisinin accumulation in the AlP. In addition, 12 transcription factors, including GATA and MYB, were upregulated in response to AlP, confirming their importance in regulating artemisinin biosynthesis. Overall, our findings could contribute to a better understanding the parallel transcriptional regulation of plant hormone transduction and artemisinin biosynthesis in A. annua L. in response to hardly soluble phosphorus fertilizer.

Primary plant nutrients modulate the reactive oxygen species metabolism and mitigate the impact of cold stress in overseeded perennial ryegrass

Overseeded perennial ryegrass (Lolium perenne L.) turf on dormant bermudagrass (Cynodon dactylon Pers. L) in transitional climatic zones (TCZ) experience a severe reduction in its growth due to cold stress. Primary plant nutrients play an important role in the cold stress tolerance of plants. To better understand the cold stress tolerance of overseeded perennial ryegrass under TCZ, a three-factor and five-level central composite rotatable design (CCRD) with a regression model was used to study the interactive effects of nitrogen (N), phosphorus (P), and potassium (K) fertilization on lipid peroxidation, electrolyte leakage, reactive oxygen species (ROS) production, and their detoxification by the photosynthetic pigments, enzymatic and non-enzymatic antioxidants. The study demonstrated substantial effects of N, P, and K fertilization on ROS production and their detoxification through enzymatic and non-enzymatic pathways in overseeded perennial ryegrass under cold stress. Our results demonstrated that the cold stress significantly enhanced malondialdehyde, electrolyte leakage, and hydrogen peroxide contents, while simultaneously decreasing ROS-scavenging enzymes, antioxidants, and photosynthetic pigments in overseeded perennial ryegrass. However, N, P, and K application mitigated cold stress-provoked adversities by enhancing soluble protein, superoxide dismutase, peroxide dismutase, catalase, and proline contents as compared to the control conditions. Moreover, N, P, and, K application enhanced chlorophyll a, chlorophyll b, total chlorophyll, and carotenoids in overseeded perennial ryegrass under cold stress as compared to the control treatments. Collectively, this 2-years study indicated that N, P, and K fertilization mitigated cold stress by activating enzymatic and non-enzymatic antioxidants defense systems, thereby concluding that efficient nutrient management is the key to enhanced cold stress tolerance of overseeded perennial ryegrass in a transitional climate. These findings revealed that turfgrass management will not only rely on breeding new varieties but also on the development of nutrient management strategies for coping cold stress.

Genome-wide identification and expression analysis of metal tolerance protein (MTP) gene family in soybean (Glycine max) under heavy metal stress

AIM

Plant metal tolerance proteins (MTPs) are plant membrane divalent cation transporters that specifically contribute to heavy metal stress resistance and mineral uptake. However, little is known about this family's molecular behaviors and biological activities in soybean.

METHODS AND RESULTS

A total of 20 potential MTP candidate genes were identified and studied in the soybean genome for phylogenetic relationships, chromosomal distributions, gene structures, gene ontology, cis-elements, and previous gene expression. Furthermore, the expression of MTPs has been investigated under different heavy metals treatments. All identified soybean MTPs (GmaMTPs) contain a cation efflux domain or a ZT dimer and are further divided into three primary cation diffusion facilitator (CDF) groups: Mn-CDFs, Zn-CDFs, and Fe/Zn-CDFs. The developmental analysis reveals that segmental duplication contributes to the GmaMTP family's expansion. Tissue-specific expression profiling revealed comparative expression profiling in similar groups, although gene expression differed between groups. GmaMTP genes displayed biased responses in either plant leaves or roots when treated with heavy metal. In the leaves and roots, nine and ten GmaMTPs responded to at least one metal ion treatment. Furthermore, in most heavy metal treatments, GmaMTP1.1, GmaMTP1.2, GmaMTP3.1, GmaMTP3.2, GmaMTP4.1, and GmaMTP4.3 exhibited significant expression responses.

CONCLUSION

Our findings provided insight into the evolution of MTPs in soybean. Overall, our findings shed light on the evolution of the MTP gene family in soybean and pave the path for further functional characterization of this gene family.

Interplay of silymarin and clove fruit extract effectively enhances cadmium stress tolerance in wheat (Triticum aestivum)

Introduction

Osmoprotectant supplementation can be used as a useful approach to enhance plant stress tolerance. However, the effect of silymarin and clove fruit extract (CFE) on wheat plants grown under cadmium (Cd) stress has not been studied.

Methods

Wheat seeds were planted in plastic pots filled with ions-free sand. A ½-strength Hoagland's nutrient solution was used for irrigation. Pots were treated with eight treatments thirteen days after sowing: 1) Control, 2) 0.5 mM silymarin foliar application [silymarin], 3) 2% CFE foliar application [CFE], 4) CFE enriched with silymarin (0.24 g silymarin L-1 of CFE) [CFE-silymarin], 5) Watering wheat seedlings with a nutritious solution of 2 mM Cd [Cd]. 6) Cadmium + silymarin, 7) Cadmium + CFE, and 8) Cadmium + CFE-silymarin. The experimental design was a completely randomized design with nine replicates.

Results and discussion

The Cd stress decreased grain yield, shoot dry weight, leaf area, carotenoids, chlorophylls, stomatal conductance, net photosynthetic rate, transpiration rate, membrane stability index, nitrogen, phosphorus, and potassium content by 66.9, 60.6, 56.7, 23.8, 33.5, 48.1, 41.2, 48.7, 42.5, 24.1, 39.9, and 24.1%, respectively. On the other hand, Cd has an Application of CFE, silymarin, or CEF-silymarin for wheat plants grown under Cd stress, significantly improved all investigated biochemical, morphological, and physiological variables and enhanced the antioxidant enzyme activities. Applying CFE and/or silymarin enhanced plant tolerance to Cd stress more efficiently. Our findings suggest using CFE-silymarin as a meaningful biostimulator for wheat plants to increase wheat plants' tolerance to Cd stress via enhancing various metabolic and physiological processes.

Individual and mutual effects of elevated carbon dioxide and temperature on salt and cadmium uptake and translocation by rice seedlings

Plant kingdoms are facing increasingly harsh environmental challenges marked by the coexposure of salinity and pollution in the pedosphere and elevated CO₂ and temperature in the atmosphere due to the rapid acceleration of industrialization and global climate change. In this study, we deployed a hydroponics-based experiment to explore the individual and mutual effects of different temperatures (low temperature, T1: 23°C; high temperature, T2: 27°C) and CO₂ concentrations (ambient CO₂: 360 ppm; medium CO₂: 450 ppm; high CO₂: 700 ppm) on the uptake and translocation of sodium chloride (NaCl, 0.0, 0.2, 0.6, and 1.1 g Na/L) and cadmium nitrate (Cd(NO₃)₂·4H₂O, 0.0, 0.2, 1.8, and 5.4 mg Cd/L) by rice seedlings. The results indicated that Cd and Na exposure significantly (P< 0.05) inhibited plant growth, but T2 and medium/high CO₂ alleviated the effects of Cd and Na on plant growth. Neither significant synergistic nor antagonistic effects of Cd and Na were observed, particularly not at T1 or high CO₂. At increasing temperatures, relative growth rates increased despite higher concentrations of Cd and Na in both rice roots and shoots. Similarly, higher CO₂ stimulated the growth rate but resulted in significantly lower concentrations of Na, while the Cd concentration was highest at medium CO₂. Coexposure experiments suggested that the concentration of Cd in roots slightly declined with additional Na and more at T2. Overall, our preliminary study suggested that global climate change may alter the distribution of mineral and toxic elements in rice plants as well as the tolerance of the plants.

Alternative splicing: transcriptional regulatory network in agroforestry

Alternative splicing (AS) in plants plays a key role in regulating the expression of numerous transcripts from a single gene in a regulatory pathway. Variable concentrations of growth regulatory hormones and external stimuli trigger alternative splicing to switch among different growth stages and adapt to environmental stresses. In the AS phenomenon, a spliceosome causes differential transcriptional modifications in messenger RNA (mRNAs), resulting in partial or complete retention of one or more introns as compared to fully spliced mRNA. Differentially expressed proteins translated from intron-retaining messenger RNA (mRNA^ir) perform vital functions in the feedback mechanism. At the post-transcriptional level, AS causes the remodeling of transcription factors (TFs) by the addition or deletion of binding domains to activate and/or repress transcription. In this study, we have summarized the specific role of AS in the regulation of gene expression through repression and activation of the transcriptional regulatory network under external stimuli and switch among developmental stages.

Schizosaccharomyces pombe Grx4, Fep1, and Php4: In silico analysis and expression response to different iron concentrations

Due to iron's essential role in cellular metabolism, most organisms must maintain their homeostasis. In this regard, the fission yeast Schizosaccharomyces pombe (sp) uses two transcription factors to regulate intracellular iron levels: spFep1 under iron-rich conditions and spPhp4 under iron-deficient conditions, which are controlled by spGrx4. However, bioinformatics analysis to understand the role of the spGrx4/spFep1/spPhp4 axis in maintaining iron homeostasis in S. pombe is still lacking. Our study aimed to perform bioinformatics analysis on S. pombe proteins and their sequence homologs in Aspergillus flavus (af), Saccharomyces cerevisiae (sc), and Homo sapiens (hs) to understand the role of spGrx4, spFep1, and spPhp4 in maintaining iron homeostasis. The three genes' expression patterns were also examined at various iron concentrations. A multiple sequence alignment analysis of spGrx4 and its sequence homologs revealed a conserved cysteine residue in each PF00085 domain. Blast results showed that hsGLRX3 is most similar to spGrx4. In addition, spFep1 is most closely related in sequence to scDal80, whereas scHap4 is most similar to spFep1. We also found two highly conserved motifs in spFep1 and its sequence homologs that are significant for iron transport systems because they contain residues involved in iron homeostasis. The scHap4 is most similar to spPhp4. Using STRING to analyze protein-protein interactions, we found that spGrx4 interacts strongly with spPhp4 and spFep1. Furthermore, spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that the three proteins play cooperative roles in iron homeostasis. At the highest level of Fe, spgrx4 had the highest expression, followed by spfep1, while spphp4 had the lowest expression; a contrast occurred at the lowest level of Fe, where spgrx4 expression remained constant. Our findings support the notion that organisms develop diverse strategies to maintain iron homeostasis.

Genotoxicity and Trace Elements Contents Analysis in Nile Tilapia (Oreochromis niloticus) Indicated the Levels of Aquatic Contamination at Three Egyptian Areas

The toxic waste and pollutants of heavy metals continuously pollute freshwater aquatic reservoirs, which have severe implications on aquatic life and human health. The present work aims to evaluate trace elements (Zn, Mn, Cu, Cd, and Pb) along with three sites, Mariout Lake, Abbassa, and River Nile Aswan in Egypt, using Nile tilapia (Oreochromis niloticus) as bioindicator. The quality assurance, health-risk assessment, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), quantitative polymerase chain reaction (qPCR), and micronucleus test were performed to investigate the effect of different trace elements on Hsp70 gene level and micronuclei formation. We observed the highest expression of Hsp70 protein band of 70 KD and stress-responsive Hsp70 gene in the liver followed by gills of Nile tilapia caught from Mariout and Abbassa, but the lowest expression was in Nile tilapia caught from Aswan. Obvious micronuclei were observed under the microscope in erythrocytes, and their number was gradually decreased in the following manner: Mariout > Abbassa > Aswan. Noticeably, Cu, Zn, and Mn contents were low. Still, Pb and Cd contents were higher than the toxicity level recommended by the Food and Agriculture Organization (FAO), The World Health Organization (WHO), and the European Commission (EC). These results showed that Hsp70's appearance at the two levels of mRNA and protein is an effective indicator for aquatic pollution besides the aberration at the chromosome level represented in the micronucleus test. Furthermore, these results showed that Nile tilapia of the Aswan region had comparatively low trace elements contamination and were suitable for consumption.

Elbaiomy R, Chen L, Wang Z, Jiao J, Zhu J, Zhou W, Chen B, Soaud SA, Abbas M, Lin N, El-Sappah AH. Evolution analysis of FRIZZY PANICLE (FZP) orthologs explored the mutations in DNA coding sequences in the grass family (Poaceae)

FRIZZY PANICLE (FZP), an essential gene that controls spikelet differentiation and development in the grass family (Poaceae), prevents the formation of axillary bud meristems and is closely associated with crop yields. It is unclear whether the FZP gene or its orthologs were selected during the evolutionary process of grass species, which possess diverse spike morphologies. In the present study, we adopted bioinformatics methods for the evolutionary analysis of FZP orthologs in species of the grass family. Thirty-five orthologs with protein sequences identical to that of the FZP gene were identified from 29 grass species. Analysis of conserved domains revealed that the AP2/ERF domains were highly conserved with almost no amino acid mutations. However, species of the tribe Triticeae, genus Oryza, and C4 plants exhibited more significant amino acid mutations in the acidic C-terminus region. Results of the phylogenetic analysis showed that the 29 grass species could be classified into three groups, namely, Triticeae, Oryza, and C4 plants. Within the Triticeae group, the FZP genes originating from the same genome were classified into the same sub-group. When selection pressure analysis was performed, significant positive selection sites were detected in species of the Triticeae and Oryza groups. Our results show that the FZP gene was selected during the grass family's evolutionary process, and functional divergence may have already occurred among the various species. Therefore, researchers investigating the FZP gene's functions should take note of the possible presence of various roles in other grass species.

In Silico Functional Prediction and Expression Analysis of C2H2 Zinc-Finger Family Transcription Factor Revealed Regulatory Role of ZmZFP126 in Maize Growth

The C2H2-zinc finger proteins (ZFP) comprise a large family of transcription factors with various functions in biological processes. In maize, the function regulation of C2H2- zine finger (ZF) genes are poorly understood. We conducted an evolution analysis and functional prediction of the maize C2H2-ZF gene family. Furthermore, the ZmZFP126 gene has been cloned and sequenced for further favorable allelic variation discovery. The phylogenetic analysis of the C2H2-ZF domain indicated that the position and sequence of the C2H2-ZF domain of the poly-zinc finger gene are relatively conserved during evolution, and the C2H2-ZF domain with the same position is highly conserved. The expression analysis of the C2H2-ZF gene family in 11 tissues at different growth stages of B73 inbred lines showed that genes with multiple transcripts were endowed with more functions. The expression analysis of the C2H2-ZF gene in P1 and P2 inbred lines under drought conditions showed that the C2H2-ZF genes were mainly subjected to negative regulation under drought stress. Functional prediction indicated that the maize C2H2-ZF gene is mainly involved in reproduction and development, especially concerning the formation of important agronomic traits in maize yield. Furthermore, sequencing and correlation analysis of the ZmZFP126 gene indicated that this gene was significantly associated with the SDW-NAP and TDW-NAP. The analysis of the relationship between maize C2H2-ZF genes and C2H2-ZF genes with known functions indicated that the functions of some C2H2-ZF genes are relatively conservative, and the functions of homologous genes in different species are similar.