Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons.

Phosphate homeostasis was studied in a monocotyledonous model plant through the characterization of the PHO1 gene family in rice (Oryza sativa). Bioinformatics and phylogenetic analysis showed that the rice genome has three PHO1 homologs, which cluster with the Arabidopsis (Arabidopsis thaliana) AtPHO1 and AtPHO1;H1, the only two genes known to be involved in root-to-shoot transfer of phosphate. In contrast to the Arabidopsis PHO1 gene family, all three rice PHO1 genes have a cis-natural antisense transcript located at the 5 ' end of the genes. Strand-specific quantitative reverse transcription-PCR analyses revealed distinct patterns of expression for sense and antisense transcripts for all three genes, both at the level of tissue expression and in response to nutrient stress. The most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. However, while the OsPHO1;2 sense transcript was relatively stable under various nutrient deficiencies, the antisense transcript was highly induced by inorganic phosphate (Pi) deficiency. Characterization of Ospho1;1 and Ospho1;2 insertion mutants revealed that only Ospho1;2 mutants had defects in Pi homeostasis, namely strong reduction in Pi transfer from root to shoot, which was accompanied by low-shoot and high-root Pi. Our data identify OsPHO1;2 as playing a key role in the transfer of Pi from roots to shoots in rice, and indicate that this gene could be regulated by its cis-natural antisense transcripts. Furthermore, phylogenetic analysis of PHO1 homologs in monocotyledons and dicotyledons revealed the emergence of a distinct clade of PHO1 genes in dicotyledons, which include members having roles other than long-distance Pi transport.

Molecular and biochemical characterisation of Trypanosoma cruzi phosphofructokinase

The characterisation of the gene encoding Trypanosoma cruzi CL Brener phosphofructokinase (PFK) and the biochemical properties of the expressed enzyme are reported here. In contradiction with previous reports, the PFK genes of CL Brener and YBM strain T. cruzi were found to be similar to their Leishmania mexicana and Trypanosoma brucei homologs in terms of both kinetic properties and size, with open reading frames encoding polypeptides with a deduced molecular mass of 53,483. The predicted amino acid sequence contains the C-terminal glycosome-targeting tripeptide SKL; this localisation was confirmed by immunofluorescence assays. In sequence comparisons with the genes of other eukaryotes, it was found that, despite being an adenosine triphosphate-dependent enzyme, T. cruzi PFK shows significant sequence similarity with inorganic pyrophosphate-dependent PFKs.

Specific duplication and dorsoventrally asymmetric expression patterns of Cycloidea-like genes in zygomorphic species of Ranunculaceae.

Floral bilateral symmetry (zygomorphy) has evolved several times independently in angiosperms from radially symmetrical (actinomorphic) ancestral states. Homologs of the Antirrhinum majus Cycloidea gene (Cyc) have been shown to control floral symmetry in diverse groups in core eudicots. In the basal eudicot family Ranunculaceae, there is a single evolutionary transition from actinomorphy to zygomorphy in the stem lineage of the tribe Delphinieae. We characterized Cyc homologs in 18 genera of Ranunculaceae, including the four genera of Delphinieae, in a sampling that represents the floral morphological diversity of this tribe, and reconstructed the evolutionary history of this gene family in Ranunculaceae. Within each of the two RanaCyL (Ranunculaceae Cycloidea-like) lineages previously identified, an additional duplication possibly predating the emergence of the Delphinieae was found, resulting in up to four gene copies in zygomorphic species. Expression analyses indicate that the RanaCyL paralogs are expressed early in floral buds and that the duration of their expression varies between species and paralog class. At most one RanaCyL paralog was expressed during the late stages of floral development in the actinomorphic species studied whereas all paralogs from the zygomorphic species were expressed, composing a species-specific identity code for perianth organs. The contrasted asymmetric patterns of expression observed in the two zygomorphic species is discussed in relation to their distinct perianth architecture.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Nematode selenoproteome: the use of the selenocysteine insertion system to decode one codon in an animal genome?

Selenocysteine (Sec) is co-translationally inserted into selenoproteins in response to codon UGA with the help of the selenocysteine insertion sequence (SECIS) element. The number of selenoproteins in animals varies, with humans having 25 and mice having 24 selenoproteins. To date, however, only one selenoprotein, thioredoxin reductase, has been detected in Caenorhabditis elegans, and this enzyme contains only one Sec. Here, we characterize the selenoproteomes of C.elegans and Caenorhabditis briggsae with three independent algorithms, one searching for pairs of homologous nematode SECIS elements, another searching for Cys- or Sec-containing homologs of potential nematode selenoprotein genes and the third identifying Sec-containing homologs of annotated nematode proteins. These methods suggest that thioredoxin reductase is the only Sec-containing protein in the C.elegans and C.briggsae genomes. In contrast, we identified additional selenoproteins in other nematodes. Assuming that Sec insertion mechanisms are conserved between nematodes and other eukaryotes, the data suggest that nematode selenoproteomes were reduced during evolution, and that in an extreme reduction case Sec insertion systems probably decode only a single UGA codon in C.elegans and C.briggsae genomes. In addition, all detected genes had a rare form of SECIS element containing a guanosine in place of a conserved adenosine present in most other SECIS structures, suggesting that in organisms with small selenoproteomes SECIS elements may change rapidly.

Structure and expression profile of the Arabidopsis PHO1 gene family indicates a broad role in inorganic phosphate homeostasis.

PHO1 has been recently identified as a protein involved in the loading of inorganic phosphate into the xylem of roots in Arabidopsis. The genome of Arabidopsis contains 11 members of the PHO1 gene family. The cDNAs of all PHO1 homologs have been cloned and sequenced. All proteins have the same topology and harbor a SPX tripartite domain in the N-terminal hydrophilic portion and an EXS domain in the C-terminal hydrophobic portion. The SPX and EXS domains have been identified in yeast (Saccharomyces cerevisiae) proteins involved in either phosphate transport or sensing or in sorting proteins to endomembranes. The Arabidopsis genome contains additional proteins of unknown function containing either a SPX or an EXS domain. Phylogenetic analysis indicated that the PHO1 family is subdivided into at least three clusters. Reverse transcription-PCR revealed a broad pattern of expression in leaves, roots, stems, and flowers for most genes, although two genes are expressed exclusively in flowers. Analysis of the activity of the promoter of all PHO1 homologs using promoter-beta-glucuronidase fusions revealed a predominant expression in the vascular tissues of roots, leaves, stems, or flowers. beta-Glucuronidase expression is also detected for several promoters in nonvascular tissue, including hydathodes, trichomes, root tip, root cortical/epidermal cells, and pollen grains. The expression pattern of PHO1 homologs indicates a likely role of the PHO1 proteins not only in the transfer of phosphate to the vascular cylinder of various tissues but also in the acquisition of phosphate into cells, such as pollen or root epidermal/cortical cells.

A coiled coil trigger site is essential for rapid binding of synaptobrevin to the SNARE acceptor complex.

Exocytosis from synaptic vesicles is driven by stepwise formation of a tight alpha-helical complex between the fusing membranes. The complex is composed of the three SNAREs: synaptobrevin 2, SNAP-25, and syntaxin 1a. An important step in complex formation is fast binding of vesicular synaptobrevin to the preformed syntaxin 1.SNAP-25 dimer. Exactly how this step relates to neurotransmitter release is not well understood. Here, we combined different approaches to gain insights into this reaction. Using computational methods, we identified a stretch in synaptobrevin 2 that may function as a coiled coil "trigger site." This site is also present in many synaptobrevin homologs functioning in other trafficking steps. Point mutations in this stretch inhibited binding to the syntaxin 1.SNAP-25 dimer and slowed fusion of liposomes. Moreover, the point mutations severely inhibited secretion from chromaffin cells. Altogether, this demonstrates that the trigger site in synaptobrevin is crucial for productive SNARE zippering.

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

ParA2, a Vibrio cholerae chromosome partitioning protein, forms left-handed helical filaments on DNA.

Most bacterial chromosomes contain homologs of plasmid partitioning (par) loci. These loci encode ATPases called ParA that are thought to contribute to the mechanical force required for chromosome and plasmid segregation. In Vibrio cholerae, the chromosome II (chrII) par locus is essential for chrII segregation. Here, we found that purified ParA2 had ATPase activities comparable to other ParA homologs, but, unlike many other ParA homologs, did not form high molecular weight complexes in the presence of ATP alone. Instead, formation of high molecular weight ParA2 polymers required DNA. Electron microscopy and three-dimensional reconstruction revealed that ParA2 formed bipolar helical filaments on double-stranded DNA in a sequence-independent manner. These filaments had a distinct change in pitch when ParA2 was polymerized in the presence of ATP versus in the absence of a nucleotide cofactor. Fitting a crystal structure of a ParA protein into our filament reconstruction showed how a dimer of ParA2 binds the DNA. The filaments formed with ATP are left-handed, but surprisingly these filaments exert no topological changes on the right-handed B-DNA to which they are bound. The stoichiometry of binding is one dimer for every eight base pairs, and this determines the geometry of the ParA2 filaments with 4.4 dimers per 120 A pitch left-handed turn. Our findings will be critical for understanding how ParA proteins function in plasmid and chromosome segregation.

Gut symbionts from distinct hosts exhibit genotoxic activity via divergent colibactin biosynthetic pathways.

Secondary metabolites produced by nonribosomal peptide synthetase (NRPS) or polyketide synthase (PKS) pathways are chemical mediators of microbial interactions in diverse environments. However, little is known about their distribution, evolution, and functional roles in bacterial symbionts associated with animals. A prominent example is "colibactin", a largely unknown family of secondary metabolites produced by Escherichia coli via a hybrid NRPS-PKS biosynthetic pathway, inflicting DNA damage upon eukaryotic cells and contributing to colorectal cancer and tumor formation in the mammalian gut. Thus far, homologs of this pathway have only been found in closely related Enterobacteriaceae, while a divergent variant of this gene cluster was recently discovered in a marine alphaproteobacterial Pseudovibrio strain. Herein, we sequenced the genome of Frischella perrara PEB0191, a bacterial gut symbiont of honey bees, and identified a homologous colibactin biosynthetic pathway related to those found in Enterobacteriaceae. We show that the colibactin genomic island (GI) has conserved gene synteny and biosynthetic module architecture across F. perrara, Enterobacteriaceae and the Pseudovibrio strain. Comparative metabolomics analyses of F. perrara and E. coli further reveal that these two bacteria produce related colibactin pathway-dependent metabolites. Finally, we demonstrate that F. perrara, like E. coli, causes DNA damage in eukaryotic cells in vitro in a colibactin pathway-dependent manner. Together, these results support that divergent variants of the colibactin biosynthetic pathway are widely distributed among bacterial symbionts, producing related secondary metabolites and likely endowing its producer with functional capabilities important for diverse symbiotic associations.

Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats.

BACKGROUND: Along the chromosome of the obligate intracellular bacteria Protochlamydia amoebophila UWE25, we recently described a genomic island Pam100G. It contains a tra unit likely involved in conjugative DNA transfer and lgrE, a 5.6-kb gene similar to five others of P. amoebophila: lgrA to lgrD, lgrF. We describe here the structure, regulation and evolution of these proteins termed LGRs since encoded by "Large G+C-Rich" genes. RESULTS: No homologs to the whole protein sequence of LGRs were found in other organisms. Phylogenetic analyses suggest that serial duplications producing the six LGRs occurred relatively recently and nucleotide usage analyses show that lgrB, lgrE and lgrF were relocated on the chromosome. The C-terminal part of LGRs is homologous to Leucine-Rich Repeats domains (LRRs). Defined by a cumulative alignment score, the 5 to 18 concatenated octacosapeptidic (28-meric) LRRs of LGRs present all a predicted alpha-helix conformation. Their closest homologs are the 28-residue RI-like LRRs of mammalian NODs and the 24-meres of some Ralstonia and Legionella proteins. Interestingly, lgrE, which is present on Pam100G like the tra operon, exhibits Pfam domains related to DNA metabolism. CONCLUSION: Comparison of the LRRs, enable us to propose a parsimonious evolutionary scenario of these domains driven by adjacent concatenations of LRRs. Our model established on bacterial LRRs can be challenged in eucaryotic proteins carrying less conserved LRRs, such as NOD proteins and Toll-like receptors.

The meiosis-specific recombinase hDmc1 forms ring structures and interacts with hRad51.

Eukaryotic cells encode two homologs of Escherichia coli RecA protein, Rad51 and Dmc1, which are required for meiotic recombination. Rad51, like E.coli RecA, forms helical nucleoprotein filaments that promote joint molecule and heteroduplex DNA formation. Electron microscopy reveals that the human meiosis-specific recombinase Dmc1 forms ring structures that bind single-stranded (ss) and double-stranded (ds) DNA. The protein binds preferentially to ssDNA tails and gaps in duplex DNA. hDmc1-ssDNA complexes exhibit an irregular, often compacted structure, and promote strand-transfer reactions with homologous duplex DNA. hDmc1 binds duplex DNA with reduced affinity to form nucleoprotein complexes. In contrast to helical RecA/Rad51 filaments, however, Dmc1 filaments are composed of a linear array of stacked protein rings. Consistent with the requirement for two recombinases in meiotic recombination, hDmc1 interacts directly with hRad51.

Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs.

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the "ortholog conjecture", or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins.

Fission yeast rad51 and dmc1, two efficient DNA recombinases forming helical nucleoprotein filaments.

Homologous recombination is important for the repair of double-strand breaks during meiosis. Eukaryotic cells require two homologs of Escherichia coli RecA protein, Rad51 and Dmc1, for meiotic recombination. To date, it is not clear, at the biochemical level, why two homologs of RecA are necessary during meiosis. To gain insight into this, we purified Schizosaccharomyces pombe Rad51 and Dmc1 to homogeneity. Purified Rad51 and Dmc1 form homo-oligomers, bind single-stranded DNA preferentially, and exhibit DNA-stimulated ATPase activity. Both Rad51 and Dmc1 promote the renaturation of complementary single-stranded DNA. Importantly, Rad51 and Dmc1 proteins catalyze ATP-dependent strand exchange reactions with homologous duplex DNA. Electron microscopy reveals that both S. pombe Rad51 and Dmc1 form nucleoprotein filaments. Rad51 formed helical nucleoprotein filaments on single-stranded DNA, whereas Dmc1 was found in two forms, as helical filaments and also as stacked rings. These results demonstrate that Rad51 and Dmc1 are both efficient recombinases in lower eukaryotes and reveal closer functional and structural similarities between the meiotic recombinase Dmc1 and Rad51. The DNA strand exchange activity of both Rad51 and Dmc1 is most likely critical for proper meiotic DNA double-strand break repair in lower eukaryotes.

Characterization of the PHO1 gene family and the responses to phosphate deficiency of Physcomitrella patens.

PHO1 was previously identified in Arabidopsis (Arabidopsis thaliana) as a protein involved in loading inorganic phosphate (Pi) into the xylem of roots and its expression was associated with the vascular cylinder. Seven genes homologous to AtPHO1 (PpPHO1;1-PpPHO1;7) have been identified in the moss Physcomitrella patens. The corresponding proteins harbor an SPX tripartite domain in the N-terminal hydrophilic portion and an EXS domain in the conserved C-terminal hydrophobic portion, both common features of the plant PHO1 family. Northern-blot analysis showed distinct expression patterns for the PpPHO1 genes, both at the tissue level and in response to phosphate deficiency. Transgenic P. patens expressing the beta-glucuronidase reporter gene under three different PpPHO1 promoters revealed distinct expression profiles in various tissues. Expression of PpPHO1;1 and PpPHO1;7 was specifically induced by Pi starvation. P. patens homologs to the Arabidopsis PHT1, DGD2, SQD1, and APS1 genes also responded to Pi deficiency by increased mRNA levels. Morphological changes associated with Pi deficiency included elongation of caulonemata with inhibition of the formation of side branches, resulting in colonies with greater diameter, but reduced mass compared to Pi-sufficient plants. Under Pi-deficient conditions, P. patens also increased the synthesis of ribonucleases and of an acid phosphatase, and increased the ratio of sulfolipids over phospholipids. These results indicate that P. patens and higher plants share some common strategies to adapt to Pi deficiency, although morphological changes are distinct, and that the PHO1 proteins are well conserved in bryophyte despite the lack of a developed vascular system.