Genes involved in mRNA surveillance are induced in Brachypodium distachyon under cadmium toxicity.

BACKGROUND: Cd accumulation in plant cells results in dramatic problems including oxidative stress and inhibition of vital enzymes. It also affects mineral uptakes by disrupting membrane permeability. Interaction among Cd and other plant nutrient elements changes the nutritional contents of crops and reduces their yield. METHODS AND RESULTS: In the present study, Cd stress in Brachypodium distachyon led to the upregulation of some heavy metal transport genes (influx or efflux) encoding cation-efflux proteins, heavy metal-associated proteins and NRAMP proteins. The Arabidopsis orthologs of the differentially expressed B. distachyon genes (DEGs) under Cd toxicity were identified, which exhibited Bradi4g26905 was an ortholog of AtALY1-2. Detailed co-expression network and gene ontology analyses found the potential involvement of the mRNA surveillance pathway in Cd tolerance in B. distachyon. These genes were shown to be downregulated by sulfur (S) deficiency. CONCLUSIONS: This is the first transcriptomic study investigating the effect of Cd toxicity in B. distachyon, a model plant for genomic studies in Poaceae (Gramineae) species. The results are expected to provide valuable information for more comprehensive research related to heavy metal toxicity in plants.

Mitochondrial iron transporter (MIT) gene in potato (Solanum tuberosum): comparative bioinformatics, physiological and expression analyses in response to drought and salinity.

Mitochondrial iron transporter (MIT) genes are essential for mitochondrial acquisition/import of iron and vital to proper functioning of mitochondria. Unlike other organisms, research on the MITs in plants is limited. The present study provides comparative bioinformatics assays for the potato MIT gene (StMIT) as well as gene expression analyses. The phylogenetic analyses revealed monocots-dicot divergence in MIT proteins and it was also found clade specific motif diversity. In addition, docking analyses indicated that Asp172 and Gly100 residues to be identified as the closest residues binding to ferrous iron. The percentage of structure overlap of the StMIT 3D protein model with Arabidopsis, maize and rice MIT proteins was found between 80.18% and 85.71%. The transcript analyses exhibited that the expression of StMIT was triggered under drought and salinity stresses. The findings of the present study would provide valuable leads for further studies targeting specifically the MIT gene and generally the plant iron metabolism.

Genome-wide identification of serine acetyltransferase (SAT) gene family in rice (Oryza sativa) and their expressions under salt stress.

BACKGROUND: Assimilation of sulfur to cysteine (Cys) occurs in presence of serine acetyltransferase (SAT). Drought and salt stresses are known to be regulated by abscisic acid, whose biosynthesis is limited by Cys. Cys is formed by cysteine synthase complex depending on SAT and OASTL enzymes. Functions of some SAT genes were identified in Arabidopsis; however, it is not known how SAT genes are regulated in rice (Oryza sativa) under salt stress. METHODS AND RESULTS: Sequence, protein domain, gene structure, nucleotide, phylogenetic, selection, gene duplication, motif, synteny, digital expression and co-expression, secondary and tertiary protein structures, and binding site analyses were conducted. The wet-lab expressions of OsSAT genes were also tested under salt stress. OsSATs have underwent purifying selection. Segmental and tandem duplications may be driving force of structural and functional divergences of OsSATs. The digital expression analyses of OsSATs showed that jasmonic acid (JA) was the only hormone inducing the expressions of OsSAT1;1, OsSAT2;1, and OsSAT2;2 whereas auxin and ABA only triggered OsSAT1;1 expression. Leaf blade is the only plant organ where all OsSATs but OsSAT1;1 were expressed. Wet-lab expressions of OsSATs indicated that OsSAT1;1, OsSAT1;2 and OsSAT1;3 genes were upregulated at different exposure times of salt stress. CONCLUSIONS: OsSAT1;1, expressed highly in rice roots, may be a hub gene regulated by cross-talk of JA, ABA and auxin hormones. The cross-talk of the mentioned hormones and the structural variations of OsSAT proteins may also explain the different responses of OsSATs to salt stress.

Ammonium transporter 1 (AMT1) gene family in tomato (Solanum lycopersicum L.): Bioinformatics, physiological and expression analyses under drought and salt stresses.

Nitrogen (N) is an essential macronutrient for plants, and mainly taken from the soil as ammonium (NH+4). It is particularly transported into the plants by AMmonium Transporters (AMTs), which are plasma membrane proteins. In the present study, genome-wide identification, physiological and expression analyses of tomato (Solanum lycopersicum L.) ammonium transporters 1 (SlAMT1) genes under drought and salt stresses were performed. Sequence analyses revealed the presence of variations in SlAMT1s at nucleotide and protein levels. While all the SlAMT1s comprise an ammonium transporter domain (PF00909), the numbers of their transmembrane helices were found to be diverse. Digital expression analyses proved that SlAMT1-3 gene had different expression patterns compared to the others, suggesting its functional diversities. The expression analyses revealed that SlAMT1 genes were 0.16 and 5.94 -fold down-regulated under drought and salt stresses, respectively. The results suggested that expression of SlAMT1 genes were adversely affected by abiotic stress conditions.

Pathogenesis related protein-1 (PR-1) genes in tomato (Solanum lycopersicum L.): Bioinformatics analyses and expression profiles in response to drought stress.

The pathogenesis-related protein 1 (PR-1) gene family play important roles in the plant metabolism in response to biotic and abiotic stresses. The present study aimed genome-wide identification and bioinformatics analyses of PR-1 genes in tomato (Solanum lycopersicum L.). The analyses resulted in the identification of 13 novel SlPR-1 genes, each of which produce a protein belonging to the CAP superfamily (PF00188). The KEGG annotation analyses revealed that the SlPR-1 proteins functioned in the environmental information processing (09130). The expression patterns of the PR-1 genes and some stress-related physiological parameters were investigated in Fusarium oxysporum sensitive and tolerant tomato varieties under drought stress. The drought stress leaded upregulation of all SlPR-1 genes, reaching up to 50 folds. The results indicate that the SlPR-1 genes play active roles in response to drought. This is the first study exhibiting the expression profiles of SlPR-1 genes under an abiotic stress, drought, in tomato.

Investigation of PIC1 (permease in chloroplasts 1) gene's role in iron homeostasis: bioinformatics and expression analyses in tomato and sorghum.

Iron (Fe) is a crucial micronutrient in plant metabolism; thus, iron homeostasis is critical for plant development. Permease in chloroplast 1 (PIC1) is the first protein determined in the chloroplast playing a role iron homeostasis. In the present study, the PIC1 gene was investigated at a genome-wide scale in four plant genomes; Arabidopsis, tomato, maize and sorghum. Based on the gene ontology database, 21 GO terms were found related to the PIC1 gene, most of which were involved in iron hemostasis and transport. The digital expression data revealed that the expression of the majority of PIC1 genes (62.5%) in Arabidopsis decreased under abiotic stress conditions. Expression profiles of tomato PIC1 (SlPIC1) and sorghum PIC1 (SbPIC1) genes were also analyzed under salt and drought stress conditions using Real Time-quantitative PCR (RT-qPCR). Our wet-lab studies demonstrated that the SbPIC1 gene was down-regulated under salt and drought stresses in all tissues, while SlPIC1 was up-regulated in all but root tissue under drought stress. Some structural variations were detected in predicted 3D structures of PIC1 proteins and the structural similarity values varied between 0.23 and 0.35. Consequently, these results may contribute to the understanding of the PIC1 gene in iron transport and homeostasis in plants.

Genome-wide analyses of ATP sulfurylase (ATPS) genes in higher plants and expression profiles in sorghum (Sorghum bicolor) under cadmium and salinity stresses.

ATP sulfurylase (ATPS, EC: 2.7.7.4) is a crucial enzyme for sulfate assimilation pathway in both plastids and cytosol in plants. In this study, genome-wide and comparative analyses of ATPSs in 11 higher plant species, including sequence and structural analyses have been performed. Expression of ATPS genes in sorghum under cadmium (Cd) and salinity (NaCl) stresses were also investigated to provide a model experimental data for the regulation of ATPS genes under stress conditions. Thirty-one ATPS genes from 11 plant species were found. It showed that ATPSs from different species have high sequence divergences, which cause structural differences among them. Phylogenetic analysis has shown that there are two major types of ATPSs evolved in dicots while monocots were evolved to have one type of ATPs. Finally, expression analysis of ATPS genes revealed tissue and stress dependent expression pattern, which indicates expressions of ATPS genes are tightly regulated.

Whirly (Why) transcription factors in tomato (Solanum lycopersicum L.): genome-wide identification and transcriptional profiling under drought and salt stresses.

Whirly (Why) transcription factor (TFs) constitute one of the important TF families which plays essential roles in plant metabolism to cope with environmental stresses. In the present study, Why genes were identified at genome-wide scale in tomato (Solanum lycopersicum), and bioinformatics analyses were implemented. Validation of Why genes expressions under drought and salt stresses were also performed using RT-qPCR. The analyses revealed the presence of two Why genes in tomato genome, SlWhy1 (Solyc05g007100.2.1) and SlWhy2 (Solyc11g044750.1.1). Both genes contained Whirly transcription factor domain structure (PF08536), and Why proteins were in basic character (pI ≥ 7). While the lengths of the proteins ranged from 268 to 236 amino acid residues for SlWhy1 and SlWhy2 respectively, exon numbers identified in both genes were seven. According to the digital expression data, SlWhy genes are expressed at medium level in different anatomical parts and developmental stages. In the promotor sequence analysis, 13 types of putative TF binding sites were identified, and the highest motif number was 46, found for GATA TF. Gene co-expression analyses revealed that complex networks for SlWhy genes, which are connected with various metabolic pathways. Based on the RT-qPCR data, both SlWhy1 and SlWhy2 genes were up-regulated under salt and drought stresses. 3D structure analyses revealed that SlWhy1 protein had a more diverged structure than SlWhy2 protein, based on their comparisons in Arabidopsis and potato. The results obtained in the present study could be a useful scientific basis for understanding Why genes in tomato and their functions under abiotic stress conditions.

Identification of O-acetylserine(thiol)lyase (OASTL) genes in sorghum (Sorghum bicolor) and gene expression analysis under cadmium stress.

Cysteine (Cys) is the first identified molecule in plant metabolism which includes both sulfur and nitrogen. It can be synthesized in three cellular compartments, containing chloroplast, cytoplasm and mitochondrion. The final step of cysteine biosynthesis is catalyzed by the O-acetylserine(thiol)lyase enzyme (OASTL, E.C. 4.2.99). In the present study, seven members of the OASTL gene family in the sorghum (Sorghum bicolor) genome were identified at a genome-wide scale and comparative bioinformatics analyses were performed between sorghum and Arabidopsis OASTLs. In all OASTL proteins, a pyridoxal-phosphate dependent domain structure (PALP, PF00291) was identified. The gene ontology annotations also revealed that all sorghum OASTL genes have KOG1252 (Cystathionine beta-synthase and related enzyme) and K01738 (cysteine synthase A) activities. In promotor sequences of OASTL genes, diverse cis-acting elements were found, including hormone and light responsiveness, abiotic stress responsiveness, and tissue-specific ones (meristem and endosperm). Sorghum OASTL genes demonstrated medium or high level expressions in anatomical parts and developmental stages based on the digital expression data. Expression of OASTL genes were also analyzed under cadmium (Cd) stress in sorghum by Real Time-quantitative PCR (RT-qPCR). The results exclusively showed that OASTL A1-2 gene was 1.12 fold up-regulated in roots, whereas cysteine synthase 26 was 2.25 fold down-regulated in leaves. The predicted 3D structure of OASTLs indicated some structural diversities as well as variations in the secondary structures.

A key gene bHLH115 in iron homeostasis: comprehensive bioinformatics analyses in Arabidopsis, tomato, rice, and maize.

bHLH115 transcription factor (TF) is a positive regulator of the Fe-deficiency and plays essential roles in the stress-related regulation network. In this study, orthologous bHLH115 genes in Arabidopsis, tomato, rice, and maize were analyzed using in silico methods. All bHLH115 proteins contained PF00010 (HLH: Helix-loop-helix DNA-binding domain) domain structure and their sub-cellular localizations were predicted as nucleus. The bHLH115 orthologues in monocots and dicots clearly diverged from each other. The expression analyses revealed that orthologous genes of bHLH115 in queried species were highly expressed in seed parts, leaf, stem, and flower parts. The bHLH115 genes were co-expressed with genes in plant defense system, and with genes involving in biotic and abiotic stress responses. In terms of protein structures, OsbHLH115 and ZmbHLH115, and AtbHLH115 and SlbHLH115 had the highest protein structure similarities. In addition, bHLH115 proteins have bZIP, bHLH and MYB transcription factor binding sites strengthens their engagement in various metabolic ways. Molecular docking analyses showed the different binding sites based on plant species, suggesting functional flexibilities of bHLH115 gene.

Genome-wide and comparative analysis of bHLH38, bHLH39, bHLH100 and bHLH101 genes in Arabidopsis, tomato, rice, soybean and maize: insights into iron (Fe) homeostasis.

Iron (Fe) is an essential element for plant life. Its deficiency impedes growth and development and excessive iron can cause the toxic effect via the Fenton reaction. Thus, plants have developed various mechanisms to acquire, distribute and utilize Fe for the maintenance of their iron homeostasis at cellular and systemic levels. A basic helix-loop-helix (bHLH) transcription factor family plays essential roles in many regulatory and development processes in plants. In this study, we aimed to understand the roles of bHLH38, bHLH39, bHLH100 and bHLH101 genes for Fe homeostasis in Arabidopsis, tomato, rice, soybean and maize species by using bioinformatics approaches. The gene/protein sequence analyses of these genes demonstrated that all bHLH proteins comprised helix-loop-helix DNA binding domain (PF00010) with varied exon numbers between 2 and 13. The phylogenetic analysis did not reveal a clear distinction between monocot and dicot plants. A total of 61 cis-elements were found in promotor sequences, including biotic and abiotic stress responsiveness, hormone responsiveness, and tissue specific expressions. The some structural divergences were identified in predicted 3D structures of bHLH proteins with different channels numbers. The co-expression network analysis demonstrated that bHLH39 and bHLH101 played more important roles in Fe regulation in Arabidopsis. The digital expression analysis showed various expression profiles of bHLH genes which were identified in developmental stages, anatomical parts, and perturbations. Particularly, bHLH39 and bHLH101 genes were found to be more active genes in Fe homeostasis. As a result, our findings can contribute to understanding of bHLH38, bHLH39, bHLH100 and bHLH101 genes in Fe homeostasis in plants.

Genome-wide identification and cadmium induced expression profiling of sulfate transporter (SULTR) genes in sorghum (Sorghum bicolor L.).

Sulfur is an essential element for all living organisms. Plants can convert inorganic sulfur into organic sulfur compounds by complex enzymatic steps. In this study, we conducted a genome-wide analysis of sulfate transporter genes (SULTRs) in the sorghum (Sorghum bicolor) genome and examined expression profiles of SbSULTR genes under 200 µM cadmium (Cd) exposure. As a result of sorghum genome analysis, 11 SULTR genes were identified, including SbSULTR1;1, SbSULTR1;2, SbSULTR1;3, SbSULTR2;1, SbSULTR2;2, SbSULTR3;1, SbSULTR3;2, SbSULTR3;3, SbSULTR3;4, SbSULTR3;5, and SbSULTR4. Given names are based on phylogeny and chromosomal locations. Except SbSULTR4, all SbSULTR proteins contained Sulfate_transp (PF00916), STAS (PF01740) domains and 12 trans-membrane domains. Phylogenetic analysis revealed that four major groups were identified such as SULTR1, 2, 3, and 4 groups and SULTR4 group was separated to other SULTR groups. In promotor sequences of SbSULTR genes, many diverse cis-acting elements were found mainly related with physiological processes such as light, stress and hormone responsiveness. The expression profiles of SbSULTR genes showed that SULTR1;2, 1;3, 3;3, and 3;5 genes up-regulated in root, while expression level of SULTR4 decreased under 200 µM Cd exposure. The predicted 3D structures of SULTR proteins showed some conformational changes, suggesting functional diversities of SbSULTRs. Finally, results of this study may contribute towards understanding SbSULTR genes and their regulations and roles in Cd stress in sorghum.

Genome-wide identification and expression profiling of EIL gene family in woody plant representative poplar (Populus trichocarpa).

This study aimed to improve current understanding on ethylene-insensitive 3-like (EIL) members, least explored in woody plants such as poplar (Populus trichocarpa Torr. & Grey). Herein, seven putative EIL members were identified in P. trichocarpa genome and were roughly annotated either as EIN3-like sequence associated with ethylene pathway or EIL3-like sequences related with sulfur (S)-pathway. Motif-distribution pattern of proteins also corroborated this annotation. They were distributed on six chromosomes (chr1, 3, 4 and 8-10), and were revealed to encode a protein of 509-662 residues with nuclear localization. The presence of ethylene insensitive 3 (EIN3; PF04873) domain (covering first 80-280 residues from N-terminus) was confirmed by Hidden Markov Model-based search. The first half of EIL proteins (∼80-280 residues including EIN3 domain) was substantially conserved. The second half (∼300-600 residues) was considerably diverged. Additionally, first half of proteins harbored acidic, proline-rich and glutamine-rich sites, and supported the essentiality of these regions in the transcriptional-activation and protein-function. Moreover, identified six segmental and one-tandem duplications demonstrated the negative or purifying selective nature of mutations. Furthermore, expression profile analysis indicated the possibility of a crosstalk between EIN3- and EIL3-like genes, and co-expression networks implicated their interactions with very diverse panels of biological molecules.

Genome-wide exploration of metal tolerance protein (MTP) genes in common wheat (Triticum aestivum): insights into metal homeostasis and biofortification.

Metal transport process in plants is a determinant of quality and quantity of the harvest. Although it is among the most important of staple crops, knowledge about genes that encode for membrane-bound metal transporters is scarce in wheat. Metal tolerance proteins (MTPs) are involved in trace metal homeostasis at the sub-cellular level, usually by providing metal efflux out of the cytosol. Here, by using various bioinformatics approaches, genes that encode for MTPs in the hexaploid wheat genome (Triticum aestivum, abbreviated as Ta) were identified and characterized. Based on the comparison with known rice MTPs, the wheat genome contained 20 MTP sequences; named as TaMTP1-8A, B and D. All TaMTPs contained a cation diffusion facilitator (CDF) family domain and most members harbored a zinc transporter dimerization domain. Based on motif, phylogeny and alignment analysis, A, B and D genomes of TaMTP3-7 sequences demonstrated higher homology compared to TaMTP1, 2 and 8. With reference to their rice orthologs, TaMTP1s and TaMTP8s belonged to Zn-CDFs, TaMTP2s to Fe/Zn-CDFs and TaMTP3-7s to Mn-CDFs. Upstream regions of TaMTP genes included diverse cis-regulatory motifs, indicating regulation by developmental stage, tissue type and stresses. A scan of the coding sequences of 20 TaMTPs against published miRNAs predicted a total of 14 potential miRNAs, mainly targeting the members of most diverged groups. Expression analysis showed that several TaMTPs were temporally and spatially regulated during the developmental time-course. In grains, MTPs were preferentially expressed in the aleurone layer, which is known as a reservoir for high concentrations of iron and zinc. The work identified and characterized metal tolerance proteins in common wheat and revealed a potential involvement of MTPs in providing a sink for trace element storage in wheat grains.

Genome-wide exploration of silicon (Si) transporter genes, Lsi1 and Lsi2 in plants; insights into Si-accumulation status/capacity of plants.

Silicon (Si) is a nonessential, beneficial micronutrient for plants. It increases the plant stress tolerance in relation to its accumulation capacity. In this work, root Si transporter genes were characterized in 17 different plants and inferred for their Si-accumulation status. A total of 62 Si transporter genes (31 Lsi1 and 31 Lsi2) were identified in studied plants. Lsi1s were 261-324 residues protein with a MIP family domain whereas Lsi2s were 472-547 residues with a citrate transporter family domain. Lsi1s possessed characteristic sequence features that can be employed as benchmark in prediction of Si-accumulation status/capacity of the plants. Silicic acid selectivity in Lsi1s was associated with two highly conserved NPA (Asn-Pro-Ala) motifs and a Gly-Ser-Gly-Arg (GSGR) ar/R filter. Two NPA regions were present in all Lsi1 members but some Ala substituted with Ser or Val. GSGR filter was only available in the proposed high and moderate Si accumulators. In phylogeny, Lsi1s formed three clusters as low, moderate and high Si accumulators based on tree topology and availability of GSGR filter. Low-accumulators contained filters WIGR, AIGR, FAAR, WVAR and AVAR, high-accumulators only with GSGR filter, and moderate-accumulators mostly with GSGR but some with A/CSGR filters. A positive correlation was also available between sequence homology and Si-accumulation status of the tested plants. Thus, availability of GSGR selectivity filter and sequence homology degree could be used as signatures in prediction of Si-accumulation status in experimentally uncharacterized plants. Moreover, interaction partner and expression profile analyses implicated the involvement of Si transporters in plant stress tolerance.

Insights into a key sulfite scavenger enzyme sulfite oxidase (SOX) gene in plants.

Sulfite oxidase (SOX) is a crucial molybdenum cofactor-containing enzyme in plants that re-oxidizes the sulfite back to sulfate in sulfite assimilation pathway. However, studies of this crucial enzyme are quite limited hence this work was attempted to understand the SOXs in four plant species namely, Arabidopsis thaliana, Solanum lycopersicum, Populus trichocarpa and Brachypodium distachyon. Herein studied SOX enzyme was characterized with both oxidoreductase molybdopterin binding and Mo-co oxidoreductase dimerization domains. The alignment and motif analyses revealed the highly conserved primary structure of SOXs. The phylogeny constructed with additional species demonstrated a clear divergence of monocots, dicots and lower plants. In addition, to further understand the phylogenetic relationship and make a functional inference, a structure-based phylogeny was constructed using normalized RMSD values in five superposed models from four modelled plant SOXs herein and one previously characterized chicken SOX structure. The plant and animal SOXs showed a clear divergence and also implicated their functional divergences. Based on tree topology, monocot B. distachyon appeared to be diverged from other dicots, pointing out a possible monocot-dicot split. The expression patterns of sulfite scavengers including SOX were differentially modulated under cold, heat, salt and high light stresses. Particularly, they tend to be up-regulated under high light and heat while being down-regulated under cold and salt stresses. The presence of cis-regulatory motifs associated with different stresses in upstream regions of SOX genes was thus justified. The protein-protein interaction network of AtSOX and network enrichment with gene ontology (GO) terms showed that most predicted proteins, including sulfite reductase, ATP sulfurylases and APS reductases were among prime enzymes involved in sulfite pathway. Finally, SOX-sulfite docked structures indicated that arginine residues particularly Arg374 is crucial for SOX-sulfite binding and additional two other residues such as Arg51 and Arg103 may be important for SOX-sulfite bindings in plants.

Genome-wide identification and expression analysis of sulfate transporter (SULTR) genes in potato (Solanum tuberosum L.).

MAIN CONCLUSION: Solanum tuberosum genome analysis revealed 12 StSULTR genes encoding 18 transcripts. Among genes annotated at group level ( StSULTR I-IV), group III members formed the largest SULTRs-cluster and were potentially involved in biotic/abiotic stress responses via various regulatory factors, and stress and signaling proteins. Employing bioinformatics tools, this study performed genome-wide identification and expression analysis of SULTR (StSULTR) genes in potato (Solanum tuberosum L.). Very strict homology search and subsequent domain verification with Hidden Markov Model revealed 12 StSULTR genes encoding 18 transcripts. StSULTR genes were mapped on seven S. tuberosum chromosomes. Annotation of StSULTR genes was also done as StSULTR I-IV at group level based mainly on the phylogenetic distribution with Arabidopsis SULTRs. Several tandem and segmental duplications were identified between StSULTR genes. Among these duplications, Ka/Ks ratios indicated neutral nature of mutations that might not be causing any selection. Two segmental and one-tandem duplications were calculated to occur around 147.69, 180.80 and 191.00 million years ago (MYA), approximately corresponding to the time of monocot/dicot divergence. Two other segmental duplications were found to occur around 61.23 and 67.83 MYA, which is very close to the origination of monocotyledons. Most cis-regulatory elements in StSULTRs were found associated with major hormones (such as abscisic acid and methyl jasmonate), and defense and stress responsiveness. The cis-element distribution in duplicated gene pairs indicated the contribution of duplication events in conferring the neofunctionalization/s in StSULTR genes. Notably, RNAseq data analyses unveiled expression profiles of StSULTR genes under different stress conditions. In particular, expression profiles of StSULTR III members suggested their involvement in plant stress responses. Additionally, gene co-expression networks of these group members included various regulatory factors, stress and signaling proteins, and housekeeping and some other proteins with unknown functions.

Molecular docking of Glycine max and Medicago truncatula ureases with urea; bioinformatics approaches.

Urease (EC 3.5.1.5) is a nickel-dependent metalloenzyme catalyzing the hydrolysis of urea into ammonia and carbon dioxide. It is present in many bacteria, fungi, yeasts and plants. Most species, with few exceptions, use nickel metalloenzyme urease to hydrolyze urea, which is one of the commonly used nitrogen fertilizer in plant growth thus its enzymatic hydrolysis possesses vital importance in agricultural practices. Considering the essentiality and importance of urea and urease activity in most plants, this study aimed to comparatively investigate the ureases of two important legume species such as Glycine max (soybean) and Medicago truncatula (barrel medic) from Fabaceae family. With additional plant species, primary and secondary structures of 37 plant ureases were comparatively analyzed using various bioinformatics tools. A structure based phylogeny was constructed using predicted 3D models of G. max and M. truncatula, whose crystallographic structures are not available, along with three additional solved urease structures from Canavalia ensiformis (PDB: 4GY7), Bacillus pasteurii (PDB: 4UBP) and Klebsiella aerogenes (PDB: 1FWJ). In addition, urease structures of these species were docked with urea to analyze the binding affinities, interacting amino acids and atom distances in urease-urea complexes. Furthermore, mutable amino acids which could potentially affect the protein active site, stability and flexibility as well as overall protein stability were analyzed in urease structures of G. max and M. truncatula. Plant ureases demonstrated similar physico-chemical properties with 833-878 amino acid residues and 89.39-90.91 kDa molecular weight with mainly acidic (5.15-6.10 pI) nature. Four protein domain structures such as urease gamma, urease beta, urease alpha and amidohydro 1 characterized the plant ureases. Secondary structure of plant ureases also demonstrated conserved protein architecture, with predominantly α-helix and random coil structures. In structure-based phylogeny, plant ureases from G. max, M. truncatula and C. ensiformis were clearly diverged from bacterial ureases of B. pasteurii and K. aerogenes. Glu, Thr, His and Gly were commonly found as interacting residues in most urease-urea docking complexes while Glu was available in all docked structures. Besides, Ala and Arg residues, which are reported in active-site architecture of plant and bacterial ureases were present in G. max urea-urease complex but not present in others. Moreover, Arg435 and Arg437 in M. truncatula and G. max, respectively were identified as highly mutable hotspot residues residing in amidohydro 1 domain of enzyme. In addition, a number of stabilizing residues were predicted upon mutation of these hotspot residues however Cys and Thr made strong implications since they were also found in codon-aligned sequences as substitutions of hotspot residues. Comparative analyses of primary sequence and secondary structure in 37 different plants demonstrated quite conserved natures of ureases in plant kingdom. Structure-based phylogeny indicated the presence of a possible prokaryote-eukaryote split and implicated the subjection of bacterial ureases to heavy selection in prokaryotic evolution compared to plants. Urea-urease docking complexes suggested that different species could share common interacting residues as well as may have some other uncommon residues at species-dependent way. In silico mutation analyses identified mutable amino acids, which were predicted to reside in catalytic site of enzyme therefore mutagenesis at these sites seemed to have adverse effects on enzyme efficiency or function. This study findings will become valuable preliminary resource for future studies to further understand the primary, secondary and tertiary structures of urease sequences in plants as well as it will provide insights about various binding features of urea-urease complexes.

In silico analysis of Mn transporters (NRAMP1) in various plant species.

Manganese (Mn) is an essential micronutrient in plant life cycle. It may be involved in photosynthesis, carbohydrate and lipid biosynthesis, and oxidative stress protection. Mn deficiency inhibits the plant growth and development, and causes the various plant symptoms such as interveinal chlorosis and tissue necrosis. Despite its importance in plant life cycle, we still have limited knowledge about Mn transporters in many plant species. Therefore, this study aimed to identify and characterize high affinity Arabidopsis Mn root transporter NRAMP1 orthologs in 17 different plant species. Various in silico methods and digital gene expression data were used in identification and characterization of NRAMP1 homologs; physico-chemical properties of sequences were calculated, putative transmembrane domains (TMDs) and conserved motif signatures were determined, phylogenetic tree was constructed, 3D models and interactome map were generated, and gene expression data was analyzed. 49 NRAMP1 homologs were identified from proteome datasets of 17 plant species using AtNRAMP1 as query. Identified sequences were characterized with a NRAMP domain structure, 10-12 putative TMDs with cytosolic N- and C-terminuses, and 10-14 exons encoding a protein of 500-588 amino acids and 53.8-64.3 kDa molecular weight with basic characteristics. Consensus transport residues, GQSSTITGTYAGQY(/F)V(/I)MQGFLD(/E/N) between TMD-8 and 9 were identified in all sequences but putative N-linked glycosylation sites were not highly conserved. In phylogeny, NRAMP1 sequences demonstrated divergence in lower and higher plants as well as in monocots and dicots. Despite divergence of lower plant Physcomitrella patens in phylogeny, it showed similarity in superposed 3D models. Phylogenetic distribution of AtNRAMP1 and 6 homologs inferred a functional relationship to NRAMP6 sequences in Mn transport, while distribution of OsNRAMP1 and 5 homologs implicated an involvement of NRAMP1 sequences in Mn transport or a cross-talk between in Fe-Mn homeostasis. Interactome analysis further confirmed this cross-talk between Mn and Fe pathways. Gene expression profile of AtNRAMP1 under Fe-, K-, P- and S-deficiencies, and cold, drought, heat and salt stresses revealed various proteins involving in transcription regulation, cofactor biosynthesis, diverse developmental roles, carbohydrate metabolism, oxidation-reduction reactions, cellular signaling and protein degradation pathways. Mn deficiency or toxicity could cause serious adverse effects in plants as well as in humans. To reduce these adversities mainly rely on understanding the molecular mechanisms underlying Mn uptake from the soil. However, we still have limited knowledge regarding the structural and functional roles of Mn transporters in many plant species. Therefore, identification and characterization of Mn root uptake transporter, NRAMP1 orthologs in various plant species will provide valuable theoretical knowledge to better understand Mn transporters as well as it may become an insight for future studies aiming to develop genetically engineered and biofortified plants.

Genome-Wide Identification, Characterization and Expression Profiling of Potato (Solanum tuberosum) Frataxin (FH) Gene.

Frataxin (FH) plays a crucial role in the biogenesis of mitochondria and the regulation of iron in the cells of various organisms. However, there has been very little research on FH in plants. In this study, the potato FH gene (StFH) was identified and characterized using a genome-wide approach, and its sequence was compared to those of FH genes from Arabidopsis, rice, and maize. The FH genes were found to have a lineage-specific distribution and were more conserved in monocots than in dicots. While multiple copies of FH genes have been reported in some species, including plants, only one isoform of FH was found in potato. The expression of StFH in leaves and roots was analyzed under two different abiotic stress conditions, and the results showed that StFH was upregulated more in leaves and that its expression levels increased with the severity of the stress. This is the first study to examine the expression of an FH gene under abiotic stress conditions.