Co-carriage of Plasmid NDM and Chromosomal KPC in Klebsiella pneumoniae ST255 Human Wound Isolate in Brazil.

Multidrug-resistant K. pneumoniae is one of the main causes of hospital-acquired infections worldwide and frequently carries antimicrobial resistance genes in moving elements. In this study, we described a K. pneumoniae clinical isolate carrying simultaneous chromosomal blaKPC, and plasmid-mediated blaNDM and blaOXA-9. The isolate is multidrug-resistant and belongs to ST 225. While blaKPC were identified in the chromosome, the blaNDM was mediated by IncFII(K) plasmid and the blaOXA-9, in a IncFIB(K) plasmid. The blaKPC context was composed by Tn4401 transposon and two insertion sequences ISKpn6 and ISKpn7. The co-production of diverse ß-lactamases brings an alert about a new adaptive profile of K. pneumoniae strains and their dissemination in the hospital-acquired infectious.

Insights into agar and secondary metabolite pathways from the genome of the red alga Gracilaria domingensis (Rhodophyta, Gracilariales).

Gracilariales is a clade of florideophycean red macroalgae known for being the main source of agar. We present a de novo genome assembly and annotation of Gracilaria domingensis, an agarophyte alga with flattened thallus widely distributed along Central and South American Atlantic intertidal zones. In addition to structural analysis, an organizational comparison was done with other Rhodophyta genomes. The nuclear genome has 78 Mbp, with 11,437 predicted coding genes, 4,075 of which did not have hits in sequence databases. We also predicted 1,567 noncoding RNAs, distributed in 14 classes. The plastid and mitochondrion genome structures were also obtained. Genes related to agar synthesis were identified. Genes for type II galactose sulfurylases could not be found. Genes related to ascorbate synthesis were found. These results suggest an intricate connection of cell wall polysaccharide synthesis and the redox systems through the use of L-galactose in Rhodophyta. The genome of G. domingensis should be valuable to phycological and aquacultural research, as it is the first tropical and Western Atlantic red macroalgal genome to be sequenced.

A transcriptome-based signature of pathological angiogenesis predicts breast cancer patient survival.

The specific genes and molecules that drive physiological angiogenesis differ from those involved in pathological angiogenesis, suggesting distinct mechanisms for these seemingly related processes. Unveiling genes and pathways preferentially associated with pathologic angiogenesis is key to understanding its mechanisms, thereby facilitating development of novel approaches to managing angiogenesis-dependent diseases. To better understand these different processes, we elucidated the transcriptome of the mouse retina in the well-accepted oxygen-induced retinopathy (OIR) model of pathological angiogenesis. We identified 153 genes changed between normal and OIR retinas, which represent a molecular signature relevant to other angiogenesis-dependent processes such as cancer. These genes robustly predict the survival of breast cancer patients, which was validated in an independent 1,000-patient test cohort (40% difference in 15-year survival; p = 2.56 x 10-21). These results suggest that the OIR model reveals key genes involved in pathological angiogenesis, and these may find important applications in stratifying tumors for treatment intensification or for angiogenesis-targeted therapies.

A ligand motif enables differential vascular targeting of endothelial junctions between brain and retina.

Endothelial heterogeneity has important implications in health and disease. Molecular markers selectively expressed in the vasculature of different organs and tissues are currently being explored in targeted therapies with promising results in preclinical and clinical studies. Noteworthy is the role that combinatorial approaches such as phage display have had in identifying such markers by using phage as nanoparticles and surrogates for billions of different peptides, screening noninvasively the vascular lumen for binding sites. Here, we show that a new peptide motif that emerged from such combinatorial screening of the vasculature binds selectively to blood vessels in the brain in vivo but not to vessels in other organs. Peptides containing a conserved motif in which amino acids Phenylalanine-Arginine-Tryptophan (FRW) predominate could be visualized by transmission electron microscopy bound to the junctions between endothelial cells in all areas of the brain, including the optic nerve, but not in other barrier-containing tissues, such as intestines and testis. Remarkably, peptides containing the motif do not bind to vessels in the retina, implying an important molecular difference between these two vascular barriers. Furthermore, the peptide allows for in vivo imaging, demonstrating that new tools for studying and imaging the brain are likely to emerge from this motif.

Origin and diversification of Xanthomonas citri subsp. citri pathotypes revealed by inclusive phylogenomic, dating, and biogeographic analyses.

BACKGROUND: Xanthomonas citri subsp. citri pathotypes cause bacterial citrus canker, being responsible for severe agricultural losses worldwide. The A pathotype has a broad host spectrum, while A* and Aw are more restricted both in hosts and in geography. Two previous phylogenomic studies led to contrasting well-supported clades for sequenced genomes of these pathotypes. No extensive biogeographical or divergence dating analytic approaches have been so far applied to available genomes. RESULTS: Based on a larger sampling of genomes than in previous studies (including six new genomes sequenced by our group, adding to a total of 95 genomes), phylogenomic analyses resulted in different resolutions, though overall indicating that A + AW is the most likely true clade. Our results suggest the high degree of recombination at some branches and the fast diversification of lineages are probable causes for this phylogenetic blurring effect. One of the genomes analyzed, X. campestris pv. durantae, was shown to be an A* strain; this strain has been reported to infect a plant of the family Verbenaceae, though there are no reports of any X. citri subsp. citri pathotypes infecting any plant outside the Citrus genus. Host reconstruction indicated the pathotype ancestor likely had plant hosts in the family Fabaceae, implying an ancient jump to the current Rutaceae hosts. Extensive dating analyses indicated that the origin of X. citri subsp. citri occurred more recently than the main phylogenetic splits of Citrus plants, suggesting dispersion rather than host-directed vicariance as the main driver of geographic expansion. An analysis of 120 pathogenic-related genes revealed pathotype-associated patterns of presence/absence. CONCLUSIONS: Our results provide novel insights into the evolutionary history of X. citri subsp. citri as well as a sound phylogenetic foundation for future evolutionary and genomic studies of its pathotypes.

Molecular Adaptations of Bacterial Mercuric Reductase to the Hypersaline Kebrit Deep in the Red Sea.

The hypersaline Kebrit Deep brine pool in the Red Sea is characterized by high levels of toxic heavy metals. Here, we describe two structurally related mercuric reductases (MerAs) from this site which were expressed in Escherichia coli Sequence similarities suggest that both genes are derived from proteobacteria, most likely the Betaproteobacteria or Gammaproteobacteria We show that one of the enzymes (K35NH) is strongly inhibited by NaCl, while the other (K09H) is activated in a NaCl-dependent manner. We infer from this difference that the two forms might support the detoxification of mercury in bacterial microorganisms that employ the compatible solutes and salt-in strategies, respectively. Three-dimensional structure modeling shows that all amino acid substitutions unique to each type are located outside the domain responsible for formation of the active MerA homodimer, and the vast majority of these are found on the surface of the molecule. Moreover, K09H exhibits the predominance of acidic over hydrophobic side chains that is typical of halophilic salt-dependent proteins. These findings enhance our understanding of how selection pressures imposed by two environmental stressors have endowed MerA enzymes with catalytic properties that can potentially function in microorganisms that utilize distinct mechanisms for osmotic balance in hypersaline environments.IMPORTANCE Analysis of two structurally homologous but catalytically distinct mercuric reductases from the Kebrit Deep brine in the Red Sea sheds light on the adaptations that enable microorganisms to cope simultaneously with extreme salinity and toxic mercury compounds. One is strongly inhibited by high NaCl concentrations, while the other exhibits NaCl-dependent activation. Their different activity profiles imply that they may derive from bacterial microorganisms that utilize compatible solutes and salt-in strategies, respectively, to maintain osmotic balance. Three-dimensional modeling reveals that regions not involved in formation of the active homodimer are conserved between the two. However, in the NaCl-dependent form, distinct amino acid substitutions are found in areas that are critical for stability in high salt. The work provides insights into how two environmental stressors have shaped the structure of orthologous enzymes through selection and adaptation, enabling them to retain their catalytic function in what may be very different cellular contexts.

IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses.

Viruses represent the most abundant life forms on the planet. Recent experimental and computational improvements have led to a dramatic increase in the number of viral genome sequences identified primarily from metagenomic samples. As a result of the expanding catalog of metagenomic viral sequences, there exists a need for a comprehensive computational platform integrating all these sequences with associated metadata and analytical tools. Here we present IMG/VR (https://img.jgi.doe.gov/vr/), the largest publicly available database of 3908 isolate reference DNA viruses with 264 413 computationally identified viral contigs from >6000 ecologically diverse metagenomic samples. Approximately half of the viral contigs are grouped into genetically distinct quasi-species clusters. Microbial hosts are predicted for 20 000 viral sequences, revealing nine microbial phyla previously unreported to be infected by viruses. Viral sequences can be queried using a variety of associated metadata, including habitat type and geographic location of the samples, or taxonomic classification according to hallmark viral genes. IMG/VR has a user-friendly interface that allows users to interrogate all integrated data and interact by comparing with external sequences, thus serving as an essential resource in the viral genomics community.

Proteomics-based identification of differentially abundant proteins reveals adaptation mechanisms of Xanthomonas citri subsp. citri during Citrus sinensis infection.

BACKGROUND: Xanthomonas citri subsp. citri (Xac) is the causal agent of citrus canker. A proteomic analysis under in planta infectious and non-infectious conditions was conducted in order to increase our knowledge about the adaptive process of Xac during infection. RESULTS: For that, a 2D-based proteomic analysis of Xac at 1, 3 and 5 days after inoculation, in comparison to Xac growth in NB media was carried out and followed by MALDI-TOF-TOF identification of 124 unique differentially abundant proteins. Among them, 79 correspond to up-regulated proteins in at least one of the three stages of infection. Our results indicate an important role of proteins related to biofilm synthesis, lipopolysaccharides biosynthesis, and iron uptake and metabolism as possible modulators of plant innate immunity, and revealed an intricate network of proteins involved in reactive oxygen species adaptation during Plants` Oxidative Burst response. We also identified proteins previously unknown to be involved in Xac-Citrus interaction, including the hypothetical protein XAC3981. A mutant strain for this gene has proved to be non-pathogenic in respect to classical symptoms of citrus canker induced in compatible plants. CONCLUSIONS: This is the first time that a protein repertoire is shown to be active and working in an integrated manner during the infection process in a compatible host, pointing to an elaborate mechanism for adaptation of Xac once inside the plant.

A novel mercuric reductase from the unique deep brine environment of Atlantis II in the Red Sea.

A unique combination of physicochemical conditions prevails in the lower convective layer (LCL) of the brine pool at Atlantis II (ATII) Deep in the Red Sea. With a maximum depth of over 2000 m, the pool is characterized by acidic pH (5.3), high temperature (68 °C), salinity (26%), low light levels, anoxia, and high concentrations of heavy metals. We have established a metagenomic dataset derived from the microbial community in the LCL, and here we describe a gene for a novel mercuric reductase, a key component of the bacterial detoxification system for mercuric and organomercurial species. The metagenome-derived gene and an ortholog from an uncultured soil bacterium were synthesized and expressed in Escherichia coli. The properties of their products show that, in contrast to the soil enzyme, the ATII-LCL mercuric reductase is functional in high salt, stable at high temperatures, resistant to high concentrations of Hg(2+), and efficiently detoxifies Hg(2+) in vivo. Interestingly, despite the marked functional differences between the orthologs, their amino acid sequences differ by less than 10%. Site-directed mutagenesis and kinetic analysis of the mutant enzymes, in conjunction with three-dimensional modeling, have identified distinct structural features that contribute to extreme halophilicity, thermostability, and high detoxification capacity, suggesting that these were acquired independently during the evolution of this enzyme. Thus, our work provides fundamental structural insights into a novel protein that has undergone multiple biochemical and biophysical adaptations to promote the survival of microorganisms that reside in the extremely demanding environment of the ATII-LCL.

Sequence verification of synthetic DNA by assembly of sequencing reads.

Gene synthesis attempts to assemble user-defined DNA sequences with base-level precision. Verifying the sequences of construction intermediates and the final product of a gene synthesis project is a critical part of the workflow, yet one that has received the least attention. Sequence validation is equally important for other kinds of curated clone collections. Ensuring that the physical sequence of a clone matches its published sequence is a common quality control step performed at least once over the course of a research project. GenoREAD is a web-based application that breaks the sequence verification process into two steps: the assembly of sequencing reads and the alignment of the resulting contig with a reference sequence. GenoREAD can determine if a clone matches its reference sequence. Its sophisticated reporting features help identify and troubleshoot problems that arise during the sequence verification process. GenoREAD has been experimentally validated on thousands of gene-sized constructs from an ORFeome project, and on longer sequences including whole plasmids and synthetic chromosomes. Comparing GenoREAD results with those from manual analysis of the sequencing data demonstrates that GenoREAD tends to be conservative in its diagnostic. GenoREAD is available at www.genoread.org.

The complete genome sequence of 'Candidatus Liberibacter americanus', associated with Citrus huanglongbing.

Liberibacter spp. form a Rhizobiaceae clade of phloem-limited pathogens of limited host range. Two obligately parasitic species have been sequenced: 'Candidatus Liberibacter asiaticus', which causes citrus huanglongbing (HLB) worldwide, and 'Ca. L. solanacearum', which causes potato "zebra chip" disease. A third (proposed) species, Liberibacter crescens, was isolated from mountain papaya, grown in axenic culture, and sequenced. In an effort to identify common host determinants, the complete genomic DNA sequence of a second HLB species, 'Ca. L. americanus' strain 'São Paulo' was determined. The circular genome of 1,195,201 bp had an average 31.12% GC content and 983 predicted protein encoding genes, 800 (81.4%) of which had a predicted function. There were 658 genes common to all sequenced Liberibacter spp. and only 8 genes common to 'Ca. L. americanus' and 'Ca. L. asiaticus' but not found in 'Ca. L. solanacearum'. Surprisingly, most of the lipopolysaccharide biosynthetic genes were missing from the 'Ca. L. americanus' genome, as well as OmpA and a key regulator of flagellin, all indicating a 'Ca. L. americanus' strategy of avoiding production of major pathogen-associated molecular patterns present in 'Ca. L. asiaticus' and 'Ca. L. solanacearum'. As with 'Ca. L. asiaticus', one of two 'Ca. L. americanus' prophages replicated as an excision plasmid and carried potential lysogenic conversion genes that appeared fragmentary or degenerated in 'Ca. L. solanacearum'.

Vascular ligand-receptor mapping by direct combinatorial selection in cancer patients.

Molecules differentially expressed in blood vessels among organs or between damaged and normal tissues, are attractive therapy targets; however, their identification within the human vasculature is challenging. Here we screened a peptide library in cancer patients to uncover ligand-receptors common or specific to certain vascular beds. Surveying ~2.35 x 10(6) motifs recovered from biopsies yielded a nonrandom distribution, indicating that systemic tissue targeting is feasible. High-throughput analysis by similarity search, protein arrays, and affinity chromatography revealed four native ligand-receptors, three of which were previously unrecognized. Two are shared among multiple tissues (integrin α4/annexin A4 and cathepsin B/apolipoprotein E3) and the other two have a restricted and specific distribution in normal tissue (prohibitin/annexin A2 in white adipose tissue) or cancer (RAGE/leukocyte proteinase-3 in bone metastases). These findings provide vascular molecular markers for biotechnology and medical applications.

Analysis of twelve genomes of the bacterium Kerstersia gyiorum from brown-throated sloths (Bradypus variegatus), the first from a non-human host.

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates.

Global proteome of the saprophytic strain Leptospira biflexa and comparative analysis with pathogenic strain Leptospira interrogans uncover new pathogenesis mechanisms.

Leptospira is a genus of bacteria that includes free-living saprophytic species found in water or soil, and pathogenic species, which are the etiologic agents of leptospirosis. Besides all the efforts, there are only a few proteins described as virulence factors in the pathogenic strain L. interrogans. This work aims to perform L. biflexa serovar Patoc1 strain Paris global proteome and to compare with the proteome database of pathogenic L. interrogans serovar Copenhageni strain Fiocruz L1-130. We identified a total of 2327 expressed proteins of L. biflexa by mass spectrometry. Using the Get Homologues software with the global proteome of L. biflexa and L. interrogans, we found orthologous proteins classified into conserved, low conserved, and specific proteins. Comparative bioinformatic analyses were performed to understand the biological functions of the proteins, subcellular localization, the presence of signal peptide, structural domains, and motifs using public softwares. These results lead to the selection of 182 low conserved within the saprophyte, and 176 specific proteins of L. interrogans. It is anticipated that these findings will indicate further studies to uncover virulence factors in the pathogenic strain. This work presents for the first time the global proteome of saprophytic strain L. biflexa serovar Patoc, strain Patoc1. SIGNIFICANCE: The comparative analysis established an array of specific proteins in pathogenic strain that will narrow down the identification of immune protective proteins that will help fight leptospirosis.

Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range.

BACKGROUND: Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain A(w)12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. RESULTS: To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. CONCLUSIONS: Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

The plant pathogen Pseudomonas syringae pv. tomato is genetically monomorphic and under strong selection to evade tomato immunity.

Recently, genome sequencing of many isolates of genetically monomorphic bacterial human pathogens has given new insights into pathogen microevolution and phylogeography. Here, we report a genome-based micro-evolutionary study of a bacterial plant pathogen, Pseudomonas syringae pv. tomato. Only 267 mutations were identified between five sequenced isolates in 3,543,009 nt of analyzed genome sequence, which suggests a recent evolutionary origin of this pathogen. Further analysis with genome-derived markers of 89 world-wide isolates showed that several genotypes exist in North America and in Europe indicating frequent pathogen movement between these world regions. Genome-derived markers and molecular analyses of key pathogen loci important for virulence and motility both suggest ongoing adaptation to the tomato host. A mutational hotspot was found in the type III-secreted effector gene hopM1. These mutations abolish the cell death triggering activity of the full-length protein indicating strong selection for loss of function of this effector, which was previously considered a virulence factor. Two non-synonymous mutations in the flagellin-encoding gene fliC allowed identifying a new microbe associated molecular pattern (MAMP) in a region distinct from the known MAMP flg22. Interestingly, the ancestral allele of this MAMP induces a stronger tomato immune response than the derived alleles. The ancestral allele has largely disappeared from today's Pto populations suggesting that flagellin-triggered immunity limits pathogen fitness even in highly virulent pathogens. An additional non-synonymous mutation was identified in flg22 in South American isolates. Therefore, MAMPs are more variable than expected differing even between otherwise almost identical isolates of the same pathogen strain.

Reconciliation of sequence data and updated annotation of the genome of Agrobacterium tumefaciens C58, and distribution of a linear chromosome in the genus Agrobacterium.

Two groups independently sequenced the Agrobacterium tumefaciens C58 genome in 2001. We report here consolidation of these sequences, updated annotation, and additional analysis of the evolutionary history of the linear chromosome, which is apparently limited to the biovar I group of Agrobacterium.

Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni.

Schistosoma mansoni is the primary causative agent of schistosomiasis, which affects 200 million individuals in 74 countries. We generated 163,000 expressed-sequence tags (ESTs) from normalized cDNA libraries from six selected developmental stages of the parasite, resulting in 31,000 assembled sequences and 92% sampling of an estimated 14,000 gene complement. By analyzing automated Gene Ontology assignments, we provide a detailed view of important S. mansoni biological systems, including characterization of metazoa-specific and eukarya-conserved genes. Phylogenetic analysis suggests an early divergence from other metazoa. The data set provides insights into the molecular mechanisms of tissue organization, development, signaling, sexual dimorphism, host interactions and immune evasion and identifies novel proteins to be investigated as vaccine candidates and potential drug targets.

SIS: a program to generate draft genome sequence scaffolds for prokaryotes.

BACKGROUND: Decreasing costs of DNA sequencing have made prokaryotic draft genome sequences increasingly common. A contig scaffold is an ordering of contigs in the correct orientation. A scaffold can help genome comparisons and guide gap closure efforts. One popular technique for obtaining contig scaffolds is to map contigs onto a reference genome. However, rearrangements that may exist between the query and reference genomes may result in incorrect scaffolds, if these rearrangements are not taken into account. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. RESULTS: We present a linear-time algorithm that can generate a set of contig scaffolds for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, even though in this general case there is no guarantee that all scaffolds in the scaffold set will be correct. We compare the performance of sis, the program that implements the algorithm, to seven other scaffold-generating programs. The results of our tests show that sis has overall better performance. CONCLUSIONS: sis is a new easy-to-use tool to generate contig scaffolds, available both as stand-alone and as a web server. The good performance of sis in our tests adds evidence that large-scale inversions are widespread in prokaryotic genomes.

Distribution of nitrogen fixation and nitrogenase-like sequences amongst microbial genomes.

BACKGROUND: The metabolic capacity for nitrogen fixation is known to be present in several prokaryotic species scattered across taxonomic groups. Experimental detection of nitrogen fixation in microbes requires species-specific conditions, making it difficult to obtain a comprehensive census of this trait. The recent and rapid increase in the availability of microbial genome sequences affords novel opportunities to re-examine the occurrence and distribution of nitrogen fixation genes. The current practice for computational prediction of nitrogen fixation is to use the presence of the nifH and/or nifD genes. RESULTS: Based on a careful comparison of the repertoire of nitrogen fixation genes in known diazotroph species we propose a new criterion for computational prediction of nitrogen fixation: the presence of a minimum set of six genes coding for structural and biosynthetic components, namely NifHDK and NifENB. Using this criterion, we conducted a comprehensive search in fully sequenced genomes and identified 149 diazotrophic species, including 82 known diazotrophs and 67 species not known to fix nitrogen. The taxonomic distribution of nitrogen fixation in Archaea was limited to the Euryarchaeota phylum; within the Bacteria domain we predict that nitrogen fixation occurs in 13 different phyla. Of these, seven phyla had not hitherto been known to contain species capable of nitrogen fixation. Our analyses also identified protein sequences that are similar to nitrogenase in organisms that do not meet the minimum-gene-set criteria. The existence of nitrogenase-like proteins lacking conserved co-factor ligands in both diazotrophs and non-diazotrophs suggests their potential for performing other, as yet unidentified, metabolic functions. CONCLUSIONS: Our predictions expand the known phylogenetic diversity of nitrogen fixation, and suggest that this trait may be much more common in nature than it is currently thought. The diverse phylogenetic distribution of nitrogenase-like proteins indicates potential new roles for anciently duplicated and divergent members of this group of enzymes.