Comparative Genome-Wide Analysis of the Malate Dehydrogenase Gene Families in Cotton

Role of phosphorus on the biogeochemical behavior of cadmium in the contaminated soil under leaching and pot experiments

Phosphorus (P) is involved in various biochemical reactions in plant growth, so it is beneficial to plants growing in soils contaminated by metals, including cadmium (Cd). However, few studies have reported on the mechanistic roles of P in mitigating Cd toxicity to ryegrass root, and especially in alleviating the disruption of the mitochondrial function of living cells. In this study, the physiological and biochemical mechanisms associated with ryegrass growth under various Cd and P treatments were investigated using leaching and pot systems. The concentration of Cd in soil leachates showed a significantly positive relationship with redox potential (P < 0.05), but negative relationship (P < 0.05) with leachate pH values and dissolved organic carbon (DOC), indicating that exogenous P addition (as H2PO4-) may decrease Cd leaching from contaminated soil. Compared to the control (without P addition), the cumulative Cd content was reduced by 53.3% and 64.5% in the soil leachate with exogenous P application (20 mg/L and 80 mg/L), respectively. Notably, application of P decreased the Cd concentrations in the symplastic fractions and increased the Cd concentrations in the apoplastic fractions in root tips, which may help to alleviate Cd stress to the protoplast. Moreover, exogenous P was found to play a positive role in mitochondrial function and Ca2+ variation in root cells under Cd stress, which provides novel insights into the mechanisms of exogenous P in alleviating plant Cd injury.

Genome-Wide Identification of AhMDHs and Analysis of Gene Expression under Manganese Toxicity Stress in Arachis hypogaea

Malate dehydrogenase (MDH) is one kind of oxidation-reduction enzyme that catalyzes the reversible conversion of oxaloacetic acid to malic acid. It has vital functions in plant development, photosynthesis, abiotic stress responses, and so on. However, there are no reports on the genome-wide identification and gene expression of the MDH gene family in Arachis hypogaea. In this study, the MDH gene family of A. hypogaea was comprehensively analyzed for the first time, and 15 AhMDH sequences were identified. According to the phylogenetic tree analysis, AhMDHs are mainly separated into three subfamilies with similar gene structures. Based on previously reported transcriptome sequencing results, the AhMDH expression quantity of roots and leaves exposed to manganese (Mn) toxicity were explored in A. hypogaea. Results revealed that many AhMDHs were upregulated when exposed to Mn toxicity, suggesting that those AhMDHs might play an important regulatory role in A. hypogaea's response to Mn toxicity stress. This study lays foundations for the functional study of AhMDHs and further reveals the mechanism of the A. hypogaea signaling pathway responding to high Mn stress.

Genome-Wide Identification and Expression Analysis of Malate Dehydrogenase Gene Family in Sweet Potato and Its Two Diploid Relatives

Malate dehydrogenase (MDH; EC 1.1.1.37) plays a vital role in plant growth and development as well as abiotic stress responses, and it is widely present in plants. However, the MDH family genes have not been explored in sweet potato. In this study, nine, ten, and ten MDH genes in sweet potato (Ipomoea batatas) and its two diploid wild relatives, Ipomoea trifida and Ipomoea triloba, respectively, were identified. These MDH genes were unevenly distributed on seven different chromosomes among the three species. The gene duplications and nucleotide substitution analysis (Ka/Ks) revealed that the MDH genes went through segmental duplications during their evolution under purifying selection. A phylogenetic and conserved structure divided these MDH genes into five subgroups. An expression analysis indicated that the MDH genes were omni-presently expressed in distinct tissues and responded to various abiotic stresses. A transcription factor prediction analysis proved that Dof, MADS-box, and MYB were the main transcription factors of sweet potato MDH genes. These findings provide molecular features of the MDH family in sweet potato and its two diploid wild relatives, which further supports functional characterizations.

OsMDH12: A Peroxisomal Malate Dehydrogenase Regulating Tiller Number and Salt Tolerance in Rice

Salinity is an important environmental factor influencing crop growth and yield. Malate dehydrogenase (MDH) catalyses the reversible conversion of oxaloacetate (OAA) to malate. While many MDHs have been identified in various plants, the biochemical function of MDH in rice remains uncharacterised, and its role in growth and salt stress response is largely unexplored. In this study, the biochemical function of OsMDH12 was determined, revealing its involvement in regulating tiller number and salt tolerance in rice. OsMDH12 localises in the peroxisome and is expressed across various organs. In vitro analysis confirmed that OsMDH12 converts OAA to malate. Seedlings of OsMDH12-overexpressing (OE) plants had shorter shoot lengths and lower fresh weights than wild-type (WT) plants, while osmdh12 mutants displayed the opposite. At maturity, OsMDH12-OE plants had fewer tillers than WT, whereas osmdh12 mutants had more, suggesting OsMDH12's role in tiller number regulation. Moreover, OsMDH12-OE plants were sensitive to salt stress, but osmdh12 mutants showed enhanced salt tolerance. The Na+/K+ content ratio increased in OsMDH12-OE plants and decreased in osmdh12 mutants, suggesting that OsMDH12 might negatively affect salt tolerance through influencing the Na+/K+ balance. These findings hint at OsMDH12's potential as a genetic tool to enhance rice growth and salt tolerance.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Lactate Dehydrogenase Superfamily in Rice and Arabidopsis: Understanding the Molecular Evolution and Structural Diversity

Lactate/malate dehydrogenases (Ldh/Maldh) are ubiquitous enzymes involved in the central metabolic pathway of plants and animals. The role of malate dehydrogenases in the plant system is very well documented. However, the role of its homolog L-lactate dehydrogenases still remains elusive. Though its occurrence is experimentally proven in a few plant species, not much is known about its role in rice. Therefore, a comprehensive genome-wide in silico investigation was carried out to identify all Ldh genes in model plants, rice and Arabidopsis, which revealed Ldh to be a multigene family encoding multiple proteins. Publicly available data suggest its role in a wide range of abiotic stresses such as anoxia, salinity, heat, submergence, cold and heavy metal stress, as also confirmed by our qRT-PCR analysis, especially in salinity and heavy metal mediated stresses. A detailed protein modelling and docking analysis using Schrodinger Suite reveals the presence of three putatively functional L-lactate dehydrogenases in rice, namely OsLdh3, OsLdh7 and OsLdh9. The analysis also highlights the important role of Ser-219, Gly-220 and His-251 in the active site geometry of OsLdh3, OsLdh7 and OsLdh9, respectively. In fact, these three genes have also been found to be highly upregulated under salinity, hypoxia and heavy metal mediated stresses in rice.

Genome-Wide Characterization and Gene Expression Analyses of Malate Dehydrogenase (MDH) Genes in Low-Phosphorus Stress Tolerance of Chinese Fir (Cunninghamia lanceolata)

Malate dehydrogenase (MDH) genes play vital roles in developmental control and environmental stress tolerance in sessile plants by modulating the organic acid-malic acid level. However, MDH genes have not yet been characterized in gymnosperm, and their roles in nutrient deficiency are largely unexplored. In this study, 12 MDH genes were identified in Chinese fir (Cunninghamia lanceolata), namely, ClMDH-1, -2, -3, …, and -12. Chinese fir is one of the most abundant commercial timber trees in China, and low phosphorus has limited its growth and production due to the acidic soil of southern China. According to the phylogenetic analysis, MDH genes were classified into five groups, and Group 2 genes (ClMDH-7, -8, -9, and 10) were only found to be present in Chinese fir but not in Arabidopsis thaliana and Populus trichocarpa. In particular, the Group 2 MDHs also had specific functional domains-Ldh_1_N (malidase NAD-binding functional domain) and Ldh_1_C (malate enzyme C-terminal functional domain)-indicating a specific function of ClMDHs in the accumulation of malate. All ClMDH genes contained the conserved MDH gene characteristic functional domains Ldh_1_N and Ldh_1_C, and all ClMDH proteins exhibited similar structures. Twelve ClMDH genes were identified from eight chromosomes, involving fifteen ClMDH homologous gene pairs, each with a Ka/Ks ratio of <1. The analysis of cis-elements, protein interactions, and transcription factor interactions of MDHs showed that the ClMDH gene might play a role in plant growth and development, and in response to stress mechanisms. The results of transcriptome data and qRT-PCR validation based on low-phosphorus stress showed that ClMDH1, ClMDH6, ClMDH7, ClMDH2, ClMDH4, ClMDH5, ClMDH10 and ClMDH11 were upregulated under low-phosphorus stress and played a role in the response of fir to low-phosphorus stress. In conclusion, these findings lay a foundation for further improving the genetic mechanism of the ClMDH gene family in response to low-phosphorus stress, exploring the potential function of this gene, promoting the improvement of fir genetics and breeding, and improving production efficiency.

Identification and Characterization of Malate Dehydrogenases in Tomato (Solanum lycopersicum L.)

Malate dehydrogenase, which facilitates the reversible conversion of malate to oxaloacetate, is essential for energy balance, plant growth, and cold and salt tolerance. However, the genome-wide study of the MDH family has not yet been carried out in tomato (Solanum lycopersicum L.). In this study, 12 MDH genes were identified from the S. lycopersicum genome and renamed according to their chromosomal location. The tomato MDH genes were split into five groups based on phylogenetic analysis and the genes that clustered together showed similar lengths, and structures, and conserved motifs in the encoded proteins. From the 12 tomato MDH genes on the chromosomes, three pairs of segmental duplication events involving four genes were found. Each pair of genes had a Ka/Ks ratio < 1, indicating that the MDH gene family of tomato was purified during evolution. Gene expression analysis exhibited that tomato MDHs were differentially expressed in different tissues, at various stages of fruit development, and differentially regulated in response to abiotic stresses. Molecular docking of four highly expressed MDHs revealed their substrate and co-factor specificity in the reversible conversion process of malate to oxaloacetate. Further, co-localization of tomato MDH genes with quantitative trait loci (QTL) of salt stress-related phenotypes revealed their broader functions in salt stress tolerance. This study lays the foundation for functional analysis of MDH genes and genetic improvement in tomato.

Genome-Wide Identification of MDH Family Genes and Their Association with Salt Tolerance in Rice

Malate dehydrogenase (MDH) is widely present in nature and regulates plant growth and development, as well as playing essential roles, especially in abiotic stress responses. Nevertheless, there is no comprehensive knowledge to date on MDH family members in rice. In this study, a total of 12 MDH members in rice were identified through genome-wide analysis and divided into three groups on the basis of their phylogenetic relationship and protein-conserved motifs. Evolutionary analysis showed that MDH proteins from rice, maize and wheat shared a close phylogenetic relationship, and the MDH family was conserved in the long-term process of domestication. We identified two segmental duplication events involving four genes, which could be the major force driving the expansion of the OsMDH family. The expression profile, cis-regulatory elements and qRT-PCR results of these genes revealed that a few OsMDH showed high tissue specificity, almost all of which had stress response elements in the promoter region, and ten MDH members were significantly induced by salt stress. Through gene-based association analysis, we found a significant correlation between salt tolerance at the seedling stage and the genetic variation of OsMDH8.1 and OsMDH12.1. Additionally, we found that the polymorphism in the promoter region of OsMDH8.1 might be related to the salt tolerance of rice. This study aimed to provide valuable information on the functional study of the rice MDH gene family related to salt stress response and revealed that OsMDH8.1 might be an important gene for the cultivar improvement of salt tolerance in rice.

Differential expressions and enzymatic properties of malate dehydrogenases in response to nutrient and metal stresses in Stylosanthes guianensis

Malate dehydrogenase (MDH, EC 1.1.1.37) is a key enzyme that catalyzes a reversible NAD-dependent dehydrogenase reaction from oxaloacetate (OAA) to malate. Although MDH has been documented to participate in cellular metabolism and redox homeostasis in plants, the roles of MDH members in the tropical legume Stylosanthes guianensis (stylo) remain less definitive. In this study, except SgMDH1 that had been previously characterized, six novel MDH genes were isolated from stylo and were then designated as SgMDH2 to SgMDH7. All of the SgMDH proteins possessed the common features of NAD binding, dimerization interface and substrate binding sites. Expression analysis showed that three SgMDHs exhibited preferential expressions in leaves, and one SgMDH was mainly expressed in roots. Furthermore, SgMDHs were regulated by nutrient deficiencies in stylo roots, especially for phosphorus (-P) and potassium (-K) deficiencies. Differential responses of SgMDHs to trace metal stress and heavy metal toxicity were observed in stylo roots, suggesting the involvement of SgMDHs in the response of stylo to metal stresses. The six novel SgMDHs were subsequently expressed and purified from Escherichia coli to analyze their biochemical properties. Although SgMDHs exhibited variations in subcellular localizations, each SgMDH protein displayed a high level of catalytic efficiency towards OAA and NADH but a low level of catalytic efficiency towards malate and NAD⁺. In addition, the activities of recombinant SgMDH proteins were pH-dependent and temperature-sensitive, and exhibited differential regulations by various metal ions. These results together suggest the potential roles of SgMDHs in stylo coping with nutrient and metal stresses.

Protein Analysis of Pollen Tubes after the Treatments of Membrane Trafficking Inhibitors Gains Insights on Molecular Mechanism Underlying Pollen Tube Polar Growth

Pollen tube elongation is characterized by a highly-polarized tip growth process dependent on an efficient vesicular transport system and largely mobilized by actin cytoskeleton. Pollen tubes are an ideal model system to study exocytosis, endocytosis, membrane recycling, and signaling network coordinating cellular processes, structural organization and vesicular trafficking activities required for tip growth. Proteomic analysis was applied to identify Nicotiana tabacum Differentially Abundant Proteins (DAPs) after in vitro pollen tube treatment with membrane trafficking inhibitors Brefeldin A, Ikarugamycin and Wortmannin. Among roughly 360 proteins separated in two-dimensional gel electrophoresis, a total of 40 spots visibly changing between treated and control samples were identified by MALDI-TOF MS and LC-ESI-MS/MS analysis. The identified proteins were classified according to biological processes, and most proteins were related to pollen tube energy metabolism, including ammino acid synthesis and lipid metabolism, structural features of pollen tube growth as well modification and actin cytoskeleton organization, stress response, and protein degradation. In-depth analysis of proteins corresponding to energy-related pathways revealed the male gametophyte to be a reliable model of energy reservoir and dynamics.

A comprehensive phenotypic and genomic characterization of Ethiopian sorghum germplasm defines core collection and reveals rich genetic potential in adaptive traits

Understanding population genetic structure and diversity of a crop is essential in designing selection strategies in plant breeding. About 2010 Ethiopian sorghum accessions were phenotyped for different traits at multiple locations. A subset of the collection, 1628 accessions, predominantly landraces, some improved varieties, and inbred lines were genotyped by sequencing. Phenotypic data revealed association of important traits with different sorghum growing agro-climatic regions, high genetic diversity and the presence of rare natural variation in the Ethiopian sorghum germplasm. Subsequent genotypic analysis determined optimum number of sub-populations, distinct cluster groups and ancestries of each sorghum accessions. To improve utilization of germplasm, a core subset of 387 lines were selected following posteriori grouping of genotypes based on cluster groups obtained through GBS analysis followed by stratified random sampling using quantitative traits. In order to evaluate how well this new sorghum and millet innovation lab (SMIL) collection from Ethiopia is represented within the largest world sorghum collection at United States Department of Agriculture - National Plant Germplasm System (USDA-NPGS) and the sorghum association panel (SAP), comparisons were conducted based on SNP data. The SMIL collection displayed high genetic diversity with some redundancy with the USDA-NPGS germplasm but SAP showed clear distinction. Furthermore, genome-environment association analysis identified candidate genes associated with adaptation to abiotic factors, that will be important for exploitation of adaptive potential to different environments. In summary, our results described the diversity and relationship of sorghum collections, representativeness of developed core and provide novel insights into candidate genes associated to abiotic stress tolerance.

Cloning and characterization of the Bambusa oldhamii BoMDH-encoded malate dehydrogenase

Malate dehydrogenase (MDH), which is ubiquitously occurred in nature, catalyzes the interconversion of malate and oxaloacetate. Higher plants contain multiple forms of MDH that differ in coenzyme specificity, subcellular localization and physiological function. A putative Bambusa oldhamii BoMDH cDNA was screened with the specific probe from the bamboo cDNA library. Sequence alignment shows that there's a high homology between the deduced amino acid sequence of BoMDH and MDH protein in Oryza sativa glyoxysome (92%). A 57 kDa fusion protein was expressed by IPTG induction in Escherichia coli BL21 (DE3), and an obvious MDH activity was detected in the recombinant protein. The molecular mass of recombinant BoMDH was estimated to be 120 kDa, and the subunit form was 57 kDa by denatured SDS-PAGE, indicating that BoMDH presents as a homodimer. The optimum temperature and pH for BoMDH activity were 40 °C and 9.5, respectively. The K_m values of BoMDH for malate and NAD⁺ were 5.2 mM and 0.52 mM. The k_cat/K_m values of BoMDH for malate and NAD⁺ were 163 min^-1 mM^-1 and 3060 min^-1 mM^-1.

Histone Deacetylase (HDAC) Gene Family in Allotetraploid Cotton and Its Diploid Progenitors: In Silico Identification, Molecular Characterization, and Gene Expression Analysis under Multiple Abiotic Stresses, DNA Damage and Phytohormone Treatments

Histone deacetylases (HDACs) play a significant role in a plant’s development and response to various environmental stimuli by regulating the gene transcription. However, HDACs remain unidentified in cotton. In this study, a total of 29 HDACs were identified in allotetraploid Gossypium hirsutum, while 15 and 13 HDACs were identified in Gossypium arboretum and Gossypium raimondii, respectively. Gossypium HDACs were classified into three groups (reduced potassium dependency 3 (RPD3)/HDA1, HD2-like, and Sir2-like (SRT) based on their sequences, and Gossypium HDACs within each subgroup shared a similar gene structure, conserved catalytic domains and motifs. Further analysis revealed that Gossypium HDACs were under a strong purifying selection and were unevenly distributed on their chromosomes. Gene expression data revealed that G. hirsutumHDACs were differentially expressed in various vegetative and reproductive tissues, as well as at different developmental stages of cotton fiber. Furthermore, some G. hirsutum HDACs were co-localized with quantitative trait loci (QTLs) and single-nucleotide polymorphism (SNPs) of fiber-related traits, indicating their function in fiber-related traits. We also showed that G. hirsutum HDACs were differentially regulated in response to plant hormones (abscisic acid (ABA) and auxin), DNA damage agent (methyl methanesulfonate (MMS)), and abiotic stresses (cold, salt, heavy metals and drought), indicating the functional diversity and specification of HDACs in response to developmental and environmental cues. In brief, our results provide fundamental information regarding G.hirsutumHDACs and highlight their potential functions in cotton growth, fiber development and stress adaptations, which will be helpful for devising innovative strategies for the improvement of cotton fiber and stress tolerance.

Comparative Genome-wide Analysis and Expression Profiling of Histone Acetyltransferase (HAT) Gene Family in Response to Hormonal Applications, Metal and Abiotic Stresses in Cotton

Post-translational modifications are involved in regulating diverse developmental processes. Histone acetyltransferases (HATs) play vital roles in the regulation of chromation structure and activate the gene transcription implicated in various cellular processes. However, HATs in cotton, as well as their regulation in response to developmental and environmental cues, remain unidentified. In this study, 9 HATs were identified from Gossypium raimondi and Gossypium arboretum, while 18 HATs were identified from Gossypium hirsutum. Based on their amino acid sequences, Gossypium HATs were divided into three groups: CPB, GNAT, and TAF_II250. Almost all the HATs within each subgroup share similar gene structure and conserved motifs. Gossypium HATs are unevenly distributed on the chromosomes, and duplication analysis suggests that Gossypium HATs are under strong purifying selection. Gene expression analysis showed that Gossypium HATs were differentially expressed in various vegetative tissues and at different stages of fiber development. Furthermore, all the HATs were differentially regulated in response to various stresses (salt, drought, cold, heavy metal and DNA damage) and hormones (abscisic acid (ABA) and auxin (NAA)). Finally, co-localization of HAT genes with reported quantitative trait loci (QTL) of fiber development were reported. Altogether, these results highlight the functional diversification of HATs in cotton growth and fiber development, as well as in response to different environmental cues. This study enhances our understanding of function of histone acetylation in cotton growth, fiber development, and stress adaptation, which will eventually lead to the long-term improvement of stress tolerance and fiber quality in cotton.