Genome-wide exploration of oil biosynthesis genes in cultivated olive tree varieties (Olea europaea): insights into regulation of oil biosynthesis.

Genome-wide oil biosynthesis was explored by de novo sequencing two cultivated olive tree (Olea europaea) varieties (cv. Ayvalik and Picual). This is the first report of the former variety sequencing. As outgroups, raw reads of cv. Leccino and scaffold-level assembly of cv. Farga were also retrieved. Each of these four cultivars was chromosome-scale assembled into 23 pseudochromosomes, with 1.31 Gbp (Farga), 0.93 Gbp (Ayvalik), 0.7 Gbp (Picual), and 0.54 Gbp (Leccino) in size. Ab initio gene finding was performed on these assemblies, using wild olive tree (oleaster)-trained parameters. High numbers of gene models were predicted and anchored to the pseudochromosomes: 69,028 (Ayvalik), 55,073 (Picual), 63,785 (Farga), and 40,449 (Leccino). Using previously reported oil biosynthesis genes from wild olive tree genome project, the following homologous sequences were identified: 1,355 (Ayvalik), 1,269 (Farga), 812 (Leccino), and 774 (Picual). Of these, 358 sequences were commonly shared by all cultivars. Besides, some sequences were cultivar unique: Ayvalik (126), Farga (118), Leccino (46), and Picual (52). These putative sequences were assigned to various GO terms, ranging from lipid metabolism to stress tolerance, from signal transactions to development, and to many others, implicating that oil biosynthesis is synergistically regulated with involvement of various other pathways.

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.

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 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.

Genome-Wide Identification and Comparative Analysis of Copper Transporter Genes in Plants.

Copper (Cu) transporters have primary importance in maintenance of physiological limits of Cu homeostasis in plants. However, structural characterization of Cu transporters in many plant species is still limited. In this study, a total of 78 potential Cu transporter genes were identified from 18 different plant species. Study revealed that Cu transporters could be characterized with a CTR protein family (PF04145) domain, three putative transmembrane domains (TMDs), a single exon number, and a basic character. Met-rich motifs at N-terminal region, MXXXM motif in TMD-2, and GXXXG motif in TMD-3 could be essential for Cu transport since they were highly conserved in all analyzed species. In phylogeny, a clear distinction was observed between Cu transporter sequences of lower and higher plants. General topological features of Cu transporters in higher plants-monocots and dicots-were highly conserved compared to lower plants. Identification of Cu transporter homologous in various plant species and their comparative analysis at gene and protein levels will become valuable theoretical basis for future studies aiming to further characterization and molecular manipulation of Cu transporters.

Essential and Beneficial Trace Elements in Plants, and Their Transport in Roots: a Review.

The essentiality of 14 mineral elements so far have been reported in plant nutrition. Eight of these elements were known as micronutrients due to their lower concentrations in plants (usually ≤100 mg/kg/dw). However, it is still challenging to mention an exact number of plant micronutrients since some elements have not been strictly proposed yet either as essential or beneficial. Micronutrients participate in very diverse metabolic processes, including from the primary and secondary metabolism to the cell defense, and from the signal transduction to the gene regulation, energy metabolism, and hormone perception. Thus, the attempt to understand the molecular mechanism(s) behind their transport has great importance in terms of basic and applied plant sciences. Moreover, their deficiency or toxicity also caused serious disease symptoms in plants, even plant destruction if not treated, and many people around the world suffer from the plant-based dietary deficiencies or metal toxicities. In this sense, shedding some light on this issue, the 13 mineral elements (Fe, B, Cu, Mn, Mo, Si, Zn, Ni, Cl, Se, Na, Al, and Co), required by plants at trace amounts, has been reviewed with the primary focus on the transport proteins (transporters/channels) in plant roots. So, providing the compiled but extensive information about the structural and functional roles of micronutrient transport genes/proteins in plant roots.

Isolation of a transcription factor DREB1A gene from Phaseolus vulgaris and computational insights into its characterization: protein modeling, docking and mutagenesis.

A transcription factor DREB1A (dehydration-responsive element-binding 1A) gene was amplified and sequenced from the common bean (Phaseolus vulgaris). PvDREB1A had a 777 base pairs (bp) open reading frame encoding a protein of 225 residues. The protein sequence contained a conserved DNA-binding AP2 domain of about 60 residues and a nuclear localization signal (NLS) at N-terminus site. PvDREB1A demonstrated high homology with other DREB1 members only in AP2 domain and NLS site. The phylogenetic distribution of different DREB members showed three main groups as DREB1-3 and PvDREB1A was a member of DREB1 group. Homology modeling and secondary structure analyses revealed that PvDREB1A AP2 domain was packed into the three-stranded antiparallel beta sheets (ß1-3) and an alpha helix (α1) almost parallel to these beta sheets. Moreover, DNA-binding AP2 domain of PvDREB1A and GCC-box containing double helix DNA were docked. The docking analysis showed that PvDREB1A AP2 domain could bind to the major groove of the DNA by three-stranded antiparallel beta sheets, with residues Gly86 or Thr87 in ß1-sheet and Arg63 or Arg64 in ß3-sheet. The docked complex also indicated that AP2 domain has a preferential for the binding of GCC stretches in the double helix DNA. A total of 36 reliably estimated hot spots residues were identified with high mutability grade but none of these residues was essential for the protein function since they are located at outside the DNA-binding AP2 domain of PvDREB1A.

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.

Identification and Comparative Analysis of H2O2-Scavenging Enzymes (Ascorbate Peroxidase and Glutathione Peroxidase) in Selected Plants Employing Bioinformatics Approaches.

Among major reactive oxygen species (ROS), hydrogen peroxide (H2O2) exhibits dual roles in plant metabolism. Low levels of H2O2 modulate many biological/physiological processes in plants; whereas, its high level can cause damage to cell structures, having severe consequences. Thus, steady-state level of cellular H2O2 must be tightly regulated. Glutathione peroxidases (GPX) and ascorbate peroxidase (APX) are two major ROS-scavenging enzymes which catalyze the reduction of H2O2 in order to prevent potential H2O2-derived cellular damage. Employing bioinformatics approaches, this study presents a comparative evaluation of both GPX and APX in 18 different plant species, and provides valuable insights into the nature and complex regulation of these enzymes. Herein, (a) potential GPX and APX genes/proteins from 18 different plant species were identified, (b) their exon/intron organization were analyzed, (c) detailed information about their physicochemical properties were provided, (d) conserved motif signatures of GPX and APX were identified, (e) their phylogenetic trees and 3D models were constructed, (f) protein-protein interaction networks were generated, and finally (g) GPX and APX gene expression profiles were analyzed. Study outcomes enlightened GPX and APX as major H2O2-scavenging enzymes at their structural and functional levels, which could be used in future studies in the current direction.

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 of galactinol synthase (GolS) genes in Solanum lycopersicum and Brachypodium distachyon.

GolS genes stand as potential candidate genes for molecular breeding and/or engineering programs in order for improving abiotic stress tolerance in plant species. In this study, a total of six galactinol synthase (GolS) genes/proteins were retrieved for Solanum lycopersicum and Brachypodium distachyon. GolS protein sequences were identified to include glyco_transf_8 (PF01501) domain structure, and to have a close molecular weight (36.40-39.59kDa) and amino acid length (318-347 aa) with a slightly acidic pI (5.35-6.40). The sub-cellular location was mainly predicted as cytoplasmic. S. lycopersicum genes located on chr 1 and 2, and included one segmental duplication while genes of B. distachyon were only on chr 1 with one tandem duplication. GolS sequences were found to have well conserved motif structures. Cis-acting analysis was performed for three abiotic stress responsive elements, including ABA responsive element (ABRE), dehydration and cold responsive elements (DRE/CRT) and low-temperature responsive element (LTRE). ABRE elements were found in all GolS genes, except for SlGolS4; DRE/CRT was not detected in any GolS genes and LTRE element found in SlGolS1 and BdGolS1 genes. AU analysis in UTR and ORF regions indicated that SlGolS and BdGolS mRNAs may have a short half-life. SlGolS3 and SlGolS4 genes may generate more stable transcripts since they included AATTAAA motif for polyadenylation signal POLASIG2. Seconder structures of SlGolS proteins were well conserved than that of BdGolS. Some structural divergences were detected in 3D structures and predicted binding sites exhibited various patterns in GolS proteins.