Molecular docking of <Emphasis Type="Italic">Glycine max</Emphasis> and <Emphasis Type="Italic">Medicago truncatula</Emphasis> 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.     

Estimates of conserved microsynteny among the genomes of <Emphasis Type="Italic">Glycine max</Emphasis>, <Emphasis Type="Italic">Medicago truncatula</Emphasis> and <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

A growing body of research indicates that microsynteny is common among dicot genomes. However, most studies focus on just one or a few genomic regions, so the extent of microsynteny across entire genomes remains poorly characterized. To estimate the level of microsynteny between Medicago truncatula (Mt) and Glycine max (soybean), and also among homoeologous segments of soybean, we used a hybridization strategy involving bacterial artificial chromosome (BAC) contigs. A Mt BAC library consisting of 30,720 clones was screened with a total of 187 soybean BAC subclones and restriction fragment length polymorphism (RFLP) probes. These probes came from 50 soybean contig groups, defined as one or more related BAC contigs anchored by the same low-copy probe. In addition, 92 whole soybean BAC clones were hybridized to filters of HindIII-digested Mt BAC DNA to identify additional cases of cross-hybridization after removal of those soybean BACs found to be repetitive in Mt. Microsynteny was inferred when at least two low-copy probes from a single soybean contig hybridized to the same Mt BAC or when a soybean BAC clone hybridized to three or more low-copy fragments from a single Mt BAC. Of the 50 soybean contig groups examined, 54% showed microsynteny to Mt. The degree of conservation among 37 groups of soybean contigs was also investigated. The results indicated substantial conservation among soybean contigs in the same group, with 86.5% of the groups showing at least some level of microsynteny. One contig group was examined in detail by a combination of physical mapping and comparative sequencing of homoeologous segments. A TBLASTX similarity search was performed between 1,085 soybean sequences on the 50 BAC contig groups and the entire Arabidopsis genome. Based on a criterion of sequence homologues <100 kb apart, each with an expected value of < or =1e-07, seven of the 50 soybean contig groups (14%) exhibited microsynteny with Arabidopsis.     

An Italian functional genomic resource for <Emphasis Type="Italic">Medicago truncatula</Emphasis>

Background

Medicago truncatula is a model species for legumes. Its functional genomics have been considerably boosted in recent years due to initiatives based both in Europe and US. Collections of mutants are becoming increasingly available and this will help unravel the genetic control of important traits for many species of legumes.

Findings

Our report is on the production of three complementary mutant collections of the model species Medicago truncatula produced in Italy in the frame of a national genomic initiative. Well established strategies were used: Tnt1 mutagenesis, TILLING and activation tagging. Both forward and reverse genetics screenings proved the efficiency of the mutagenesis approaches adopted, enabling the isolation of interesting mutants which are in course of characterization. We anticipate that the reported collections will be complementary to the recently established functional genomics tools developed for Medicago truncatula both in Europe and in the United States.

     

Fine mapping of a <Emphasis Type="Italic">Phytophthora</Emphasis>-resistance gene <Emphasis Type="Italic">RpsWY</Emphasis> in soybean (<Emphasis Type="Italic">Glycine max</Emphasis> L.) by high-throughput genome-wide sequencing

Key message

Using a combination of phenotypic screening, genetic and statistical analyses, and high-throughput genome-wide sequencing, we have finely mapped a dominant Phytophthora resistance gene in soybean cultivar Wayao.

Abstract

Phytophthora root rot (PRR) caused by Phytophthora sojae is one of the most important soil-borne diseases in many soybean-production regions in the world. Identification of resistant gene(s) and incorporating them into elite varieties are an effective way for breeding to prevent soybean from being harmed by this disease. Two soybean populations of 191 F₂ individuals and 196 F_7:8 recombinant inbred lines (RILs) were developed to map Rps gene by crossing a susceptible cultivar Huachun 2 with the resistant cultivar Wayao. Genetic analysis of the F₂ population indicated that PRR resistance in Wayao was controlled by a single dominant gene, temporarily named RpsWY, which was mapped on chromosome 3. A high-density genetic linkage bin map was constructed using 3469 recombination bins of the RILs to explore the candidate genes by the high-throughput genome-wide sequencing. The results of genotypic analysis showed that the RpsWY gene was located in bin 401 between 4466230 and 4502773 bp on chromosome 3 through line 71 and 100 of the RILs. Four predicted genes (Glyma03g04350, Glyma03g04360, Glyma03g04370, and Glyma03g04380) were found at the narrowed region of 36.5 kb in bin 401. These results suggest that the high-throughput genome-wide resequencing is an effective method to fine map PRR candidate genes.

     

Comparative genomic analysis of CIPK gene family in <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Populus</Emphasis>

Calcium serves as a second messenger in various signal transduction pathways in plants. CBL-interacting protein kinases (CIPKs), which have a variety of functions, are involved in calcium signal transduction. Previous, the studies on CIPK family members focused on Arabidopsis and rice. Here, we present a comparative genomic analysis of the CIPK gene family in Arabidopsis and poplar, a model tree species. Twenty-seven potential CIPKs were identified from poplar using genome-wide analysis. Like the CIPK gene family from Arabidopsis, CIPK genes from poplar were also divided into intron-free and intron-harboring groups. In the intron-harboring group, the intron distribution of CIPKs is rather conserved during the genome evolutionary process. Many homologous gene pairs were found in the CIPK gene family, indicating duplication events might contribute to the amplification of this gene family. The phylogenetic comparison of CIPKs in combination with intron distribution analysis revealed that CIPK genes from both Arabidopsis and poplar might have an ancient origin, which formed earlier than the separation of these two eudicot species. Our genomic and bioinformatic analysis will provide an important foundation for further functional dissection of the CBL-CIPK signaling network in poplars. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetic and physical localization of an anthracnose resistance gene in <Emphasis Type="Italic">Medicago truncatula</Emphasis>

Anthracnose of alfalfa, caused by the fungal pathogen Colletotrichum trifolii, is one of the most destructive diseases of alfalfa worldwide. An improved understanding of the genetic and molecular mechanisms underlying host resistance will facilitate the development of resistant alfalfa cultivars, thus providing the most efficient and environmentally sound strategy to control alfalfa diseases. Unfortunately, cultivated alfalfa has an intractable genetic system because of its tetrasomic inheritance and out-crossing nature. Nevertheless, the model legume Medicago truncatula, a close relative of alfalfa, has the potential to serve as a surrogate to map and clone the counterparts of agronomically important genes in alfalfa—particularly, disease resistance genes against economically important pathogens. Here we describe the high-resolution genetic and physical mapping of RCT1, a host resistance gene against C. trifolii race 1 in M. truncatula. We have delimited the RCT1 locus within a physical interval spanning ∼200 kb located on the top of M. truncatula linkage group 4. RCT1 is part of a complex locus containing numerous genes homologous to previously characterized TIR-NBS-LRR type resistance genes. The result presented in this paper will facilitate the positional cloning of RCT1 in Medicago.     

<Emphasis Type="Italic">GmDim1</Emphasis> Gene Encodes Nucleolar Localized U5-Small Nuclear Ribonucleoprotein in <Emphasis Type="Italic">Glycine max</Emphasis>

The dim1⁺ gene family is essential for G2/M transition during mitosis and encodes a small nuclear ribonucleoprotein that functions in the mRNA splicing machinery of eukaryotes. However, the plant homolog of DIM1 gene has not been defined yet. Here, we identified a gene named GmDim1 positioned on chromosome 9 of soybean (Glycine max (L.) Merr.) with 80% homology to other eukaryotic dim1⁺ family genes. A domain of soybean DIM1 protein was primarily conserved with U5 snRNP protein family and secondarily aligned with mitotic DIM1 protein family. The GmDim1 gene was expressed constitutively in all soybean organs. The transgenic Arabidopsis thaliana (L.) plants overexpressing GmDim1 showed early flowering and stem elongation, produced multiple shoots and continued flowering after the post-flowering stage. DIM1 proteins transiently expressed in onion cells were localized in the nucleus with dense deposition in the nucleolus. Therefore, we propose that the soybean GmDim1 gene is a component of plant U5 snRNP involved in mRNA splicing and normal progress of plant growth.     

IS<Emphasis Type="Italic">4</Emphasis> family goes genomic

Background  

Insertion sequences (ISs) are small, mobile DNA entities able to expand in prokaryotic genomes and trigger important rearrangements. To understand their role in evolution, accurate IS taxonomy is essential. The IS 4 family is composed of ~70 elements and, like some other families, displays extremely elevated levels of internal divergence impeding its classification. The increasing availability of complete genome sequences provides a valuable source for the discovery of additional IS 4 elements. In this study, this genomic database was used to update the structural and functional definition of the IS 4 family.     

Tobacco transformants expressing the <Emphasis Type="Italic">Medicago truncatula</Emphasis> ornithine aminotransferase cDNA

The Medicago truncatula ornithine aminotransferase cDNA was cloned under the potent constitutive 35S RNA promoter of the cauliflower mosaic virus and transferred into the genome of tobacco Nicotiana tabacum SR1 plants. Transformed tobacco plants grew better in salinity stress, but did not differ in proline content under normal or stress conditions from control plants. It was assumed that the role of ornithine aminotransferase in the molecular mechanisms of stress resistance is not associated with additional proline synthesis.     

Use of fused <Emphasis Type="Italic">gfp</Emphasis> and <Emphasis Type="Italic">gus</Emphasis> reporters for the recovery of transformed <Emphasis Type="Italic">Medicago truncatula</Emphasis> somatic embryos without selective pressure

We developed an alternative methodology for in vitro selection of transgenic Medicago truncatula cv. Jemalong plants using a bifunctional construct in which the coding sequences for the green fluorescent protein (GFP) and the β-glucuronidase protein (GUS) are fused. An Agrobacterium-mediated transformation protocol was used followed by regeneration via somatic embryogenesis in the dark, to avoid the synthesis and the consequent autofluorescence of chlorophyll. This method is a clear advantage over antibiotic and herbicide selection in which survival of non-transformed tissue is commonly reported, with the reassurance that all the somatic embryos selected as GFP positive are transformed. This was subsequently corroborated by the detection of GUS activity in leaves, stems and roots of the regenerated plants. Without antibiotic selection, and performing the embryo induction in the dark, it was possible to attest the advantage of using GFP as an in vivo detectable reporter for early embryo selection. The fusion with the GUS coding sequence provided additional evidence for the transformation of the previously selected embryos.