Evolutionary, structural and functional relationships revealed by comparative analysis of syntenic genes in Rhizobiales

Determination of the metabolic pathways for degradation of naphthalene and pyrene in Amycolatopsis sp. Poz14

Polycyclic aromatic hydrocarbons (PAHs) constitute important soil contaminants derived from petroleum. Poz14 strain can degrade pyrene and naphthalene. Its genome presented 9333 genes, among them those required for PAHs degradation. By phylogenomic analysis, the strain might be assigned to Amycolatopsis nivea. The strain was grown in glucose, pyrene, and naphthalene to compare their proteomes; 180 proteins were detected in total, and 90 of them were exclusives for xenobiotic conditions. Functions enriched with the xenobiotics belonged to transcription, translation, modification of proteins and transport of inorganic ions. Enriched pathways were pentose phosphate, proteasome and RNA degradation; in contrast, in glucose were glycolysis/gluconeogenesis and glyoxylate cycle. Proteins proposed to participate in the upper PAHs degradation were multicomponent oxygenase complexes, Rieske oxygenases, and dioxygenases; in the lower pathways were ortho-cleavage of catechol, phenylacetate, phenylpropionate, benzoate, and anthranilate. The catechol dioxygenase activity was measured and found increased when the strain was grown in naphthalene. Amycolatopsis sp. Poz14 genome and proteome revealed the PAHs degradation pathways and functions helping to contend the effects of such process.

Optogenetics in Sinorhizobium meliloti Enables Spatial Control of Exopolysaccharide Production and Biofilm Structure

Microorganisms play a vital role in shaping the soil environment and enhancing plant growth by interacting with plant root systems. Because of the vast diversity of cell types involved, combined with dynamic and spatial heterogeneity, identifying the causal contribution of a defined factor, such as a microbial exopolysaccharide (EPS), remains elusive. Synthetic approaches that enable orthogonal control of microbial pathways are a promising means to dissect such complexity. Here we report the implementation of a synthetic, light-activated, transcriptional control platform using the blue-light responsive DNA binding protein EL222 in the nitrogen fixing soil bacterium Sinorhizobium meliloti. By fine-tuning the system, we successfully achieved optical control of an EPS production pathway without significant basal expression under noninducing (dark) conditions. Optical control of EPS recapitulated important behaviors such as a mucoid plate phenotype and formation of structured biofilms, enabling spatial control of biofilm structures in S. meliloti. The successful implementation of optically controlled gene expression in S. meliloti enables systematic investigation of how genotype and microenvironmental factors together shape phenotype in situ.

Synthesis of Rhizobial Exopolysaccharides and Their Importance for Symbiosis with Legume Plants

Rhizobia dwell and multiply in the soil and represent a unique group of bacteria able to enter into a symbiotic interaction with plants from the Fabaceae family and fix atmospheric nitrogen inside de novo created plant organs, called nodules. One of the key determinants of the successful interaction between these bacteria and plants are exopolysaccharides, which represent species-specific homo- and heteropolymers of different carbohydrate units frequently decorated by non-carbohydrate substituents. Exopolysaccharides are typically built from repeat units assembled by the Wzx/Wzy-dependent pathway, where individual subunits are synthesized in conjunction with the lipid anchor undecaprenylphosphate (und-PP), due to the activity of glycosyltransferases. Complete oligosaccharide repeat units are transferred to the periplasmic space by the activity of the Wzx flippase, and, while still being anchored in the membrane, they are joined by the polymerase Wzy. Here we have focused on the genetic control over the process of exopolysaccharides (EPS) biosynthesis in rhizobia, with emphasis put on the recent advancements in understanding the mode of action of the key proteins operating in the pathway. A role played by exopolysaccharide in Rhizobium-legume symbiosis, including recent data confirming the signaling function of EPS, is also discussed.

Genomic studies of nitrogen-fixing rhizobial strains from Phaseolus vulgaris seeds and nodules

Background

Rhizobia are soil bacteria that establish symbiotic relationships with legumes and fix nitrogen in root nodules. We recently reported that several nitrogen-fixing rhizobial strains, belonging to Rhizobium phaseoli, R. trifolii, R. grahamii and Sinorhizobium americanum, were able to colonize Phaseolus vulgaris (common bean) seeds. To gain further insight into the traits that support this ability, we analyzed the genomic sequences and proteomes of R. phaseoli (CCGM1) and S. americanum (CCGM7) strains from seeds and compared them with those of the closely related strains CIAT652 and CFNEI73, respectively, isolated only from nodules.

Results

In a fine structural study of the S. americanum genomes, the chromosomes, megaplasmids and symbiotic plasmids were highly conserved and syntenic, with the exception of the smaller plasmid, which appeared unrelated. The symbiotic tract of CCGM7 appeared more disperse, possibly due to the action of transposases. The chromosomes of seed strains had less transposases and strain-specific genes. The seed strains CCGM1 and CCGM7 shared about half of their genomes with their closest strains (3353 and 3472 orthologs respectively), but a large fraction of the rest also had homology with other rhizobia. They contained 315 and 204 strain-specific genes, respectively, particularly abundant in the functions of transcription, motility, energy generation and cofactor biosynthesis. The proteomes of seed and nodule strains were obtained and showed a particular profile for each of the strains. About 82 % of the proteins in the comparisons appeared similar. Forty of the most abundant proteins in each strain were identified; these proteins in seed strains were involved in stress responses and coenzyme and cofactor biosynthesis and in the nodule strains mainly in central processes. Only 3 % of the abundant proteins had hypothetical functions.

Conclusions

Functions that were enriched in the genomes and proteomes of seed strains possibly participate in the successful occupancy of the new niche. The genome of the strains had features possibly related to their presence in the seeds. This study helps to understand traits of rhizobia involved in seed adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3053-z) contains supplementary material, which is available to authorized users.

Species in Wolbachia? Proposal for the designation of 'Candidatus Wolbachia bourtzisii', 'Candidatus Wolbachia onchocercicola', 'Candidatus Wolbachia blaxteri', 'Candidatus Wolbachia brugii', 'Candidatus Wolbachia taylori', 'Candidatus Wolbachia collembolicola' and 'Candidatus Wolbachia multihospitum' for the different species within Wolbachia supergroups

Wolbachia are highly extended bacterial endosymbionts that infect arthropods and filarial nematodes and produce contrasting phenotypes on their hosts. Wolbachia taxonomy has been understudied. Currently, Wolbachia strains are classified into phylogenetic supergroups. Here we applied phylogenomic analyses to study Wolbachia evolutionary relationships and examined metrics derived from their genome sequences such as average nucleotide identity (ANI), in silico DNA-DNA hybridization (DDH), G+C content, and synteny to shed light on the taxonomy of these bacteria. Draft genome sequences of strains wDacA and wDacB obtained from the carmine cochineal insect Dactylopius coccus were included. Although all analyses indicated that each Wolbachia supergroup represents a distinct evolutionary lineage, we found that some of the analyzed supergroups showed enough internal heterogeneity to be considered as assemblages of more than one species. Thus, supergroups would represent supraspecific groupings. Consequently, Wolbachia pipientis nomen species would apply only to strains of supergroup B and we propose the designation of 'Candidatus Wolbachia bourtzisii', 'Candidatus Wolbachia onchocercicola', 'Candidatus Wolbachia blaxterii', 'Candidatus Wolbachia brugii', 'Candidatus Wolbachia taylorii', 'Candidatus Wolbachia collembolicola' and 'Candidatus Wolbachia multihospitis' for other supergroups.

Characterization of soil bacterial assemblies in Brazilian savanna-like vegetation reveals acidobacteria dominance

The Brazilian Cerrado is the second largest biome in Brazil and is considered a biodiversity hotspot. In this work, we compared the bacterial communities in Cerrado soil associated with four types of native vegetation (Cerrado Denso, Cerrado sensu stricto, Campo Sujo, and Mata de Galeria) by ribosomal RNA intergenic spacer analysis, terminal fragment restriction length polymorphism and pyrosequencing. The fingerprinting results were very similar. The bacterial communities of Cerrado Denso and Cerrado sensu stricto grouped together and were distinct from those in Campo Sujo and Mata de Galeria. Pyrosequencing generated approximately 40,000 16S rRNA gene sequences per sample and allowed the identification of 17 phyla in soil samples under Cerrado vegetation. Acidobacteria were dominant in all areas studied with a relative frequency of 40-47 %, followed closely by Proteobacteria accounting for 34-40 % of the sequences. Results from all molecular techniques used suggested that the bacterial communities of Cerrado sensu stricto and Cerrado Denso are very similar to each other, while Campo Sujo forms a separate group, and Mata de Galeria is the most distinct with higher species richness. This is the first extensive study of native Cerrado soil microbiota, an important but endangered biome.

Rhizobial extrachromosomal replicon variability, stability and expression in natural niches

In bacteria, niche adaptation may be determined by mobile extrachromosomal elements. A remarkable characteristic of Rhizobium and Ensifer (Sinorhizobium) but also of Agrobacterium species is that almost half of the genome is contained in several large extrachromosomal replicons (ERs). They encode a plethora of functions, some of them required for bacterial survival, niche adaptation, plasmid transfer or stability. In spite of this, plasmid loss is common in rhizobia upon subculturing. Rhizobial gene-expression studies in plant rhizospheres with novel results from transcriptomic analysis of Rhizobium phaseoli in maize and Phaseolus vulgaris roots highlight the role of ERs in natural niches and allowed the identification of common extrachromosomal genes expressed in association with plant rootlets and the replicons involved.

Conservation of gene order and content in the circular chromosomes of 'Candidatus Liberibacter asiaticus' and other Rhizobiales

'Ca. Liberibacter asiaticus,' an insect-vectored, obligate intracellular bacterium associated with citrus-greening disease, also called "HLB," is a member of the Rhizobiales along with nitrogen-fixing microsymbionts Sinorhizobium meliloti and Bradyrhizobium japonicum, plant pathogen Agrobacterium tumefaciens and facultative intracellular mammalian pathogen Bartonella henselae. Comparative analyses of their circular chromosomes identified 514 orthologous genes shared among all five species. Shared among all five species are 50 identical blocks of microsyntenous orthologous genes (MOGs), containing a total of 283 genes. While retaining highly conserved genomic blocks of microsynteny, divergent evolution, horizontal gene transfer and niche specialization have disrupted macrosynteny among the five circular chromosomes compared. Highly conserved microsyntenous gene clusters help define the Rhizobiales, an order previously defined by 16S RNA gene similarity and herein represented by the three families: Bartonellaceae, Bradyrhizobiaceae and Rhizobiaceae. Genes without orthologs in the other four species help define individual species. The circular chromosomes of each of the five Rhizobiales species examined had genes lacking orthologs in the other four species. For example, 63 proteins are encoded by genes of 'Ca. Liberibacter asiaticus' not shared with other members of the Rhizobiales. Of these 63 proteins, 17 have predicted functions related to DNA replication or RNA transcription, and some of these may have roles related to low genomic GC content. An additional 17 proteins have predicted functions relevant to cellular processes, particularly modifications of the cell surface. Seventeen unshared proteins have specific metabolic functions including a pathway to synthesize cholesterol encoded by a seven-gene operon. The remaining 12 proteins encoded by 'Ca. Liberibacter asiaticus' genes not shared with other Rhizobiales are of bacteriophage origin. 'Ca. Liberibacter asiaticus' shares 11 genes with only Sinorhizobium meliloti and 12 genes are shared with only Bartonella henselae.

The genetics of symbiotic nitrogen fixation: comparative genomics of 14 rhizobia strains by resolution of protein clusters

The symbiotic relationship between legumes and nitrogen fixing bacteria is critical for agriculture, as it may have profound impacts on lowering costs for farmers, on land sustainability, on soil quality, and on mitigation of greenhouse gas emissions. However, despite the importance of the symbioses to the global nitrogen cycling balance, very few rhizobial genomes have been sequenced so far, although there are some ongoing efforts in sequencing elite strains. In this study, the genomes of fourteen selected strains of the order Rhizobiales, all previously fully sequenced and annotated, were compared to assess differences between the strains and to investigate the feasibility of defining a core ‘symbiome’—the essential genes required by all rhizobia for nodulation and nitrogen fixation. Comparison of these whole genomes has revealed valuable information, such as several events of lateral gene transfer, particularly in the symbiotic plasmids and genomic islands that have contributed to a better understanding of the evolution of contrasting symbioses. Unique genes were also identified, as well as omissions of symbiotic genes that were expected to be found. Protein comparisons have also allowed the identification of a variety of similarities and differences in several groups of genes, including those involved in nodulation, nitrogen fixation, production of exopolysaccharides, Type I to Type VI secretion systems, among others, and identifying some key genes that could be related to host specificity and/or a better saprophytic ability. However, while several significant differences in the type and number of proteins were observed, the evidence presented suggests no simple core symbiome exists. A more abstract systems biology concept of nitrogen fixing symbiosis may be required. The results have also highlighted that comparative genomics represents a valuable tool for capturing specificities and generalities of each genome.

Sequence variability of Rhizobiales orthologs and relationship with physico-chemical characteristics of proteins

Background

Chromosomal orthologs can reveal the shared ancestral gene set and their evolutionary trends. Additionally, physico-chemical properties of encoded proteins could provide information about functional adaptation and ecological niche requirements.

Results

We analyzed 7080 genes (five groups of 1416 orthologs each) from Rhizobiales species (S. meliloti, R. etli, and M. loti, plant symbionts; A. tumefaciens, a plant pathogen; and B. melitensis, an animal pathogen). We evaluated their phylogenetic relationships and observed three main topologies. The first, with closer association of R. etli to A. tumefaciens; the second with R. etli closer to S. meliloti; and the third with A. tumefaciens and S. meliloti as the closest pair. This was not unusual, given the close relatedness of these three species. We calculated the synonymous (dS) and nonsynonymous (dN) substitution rates of these orthologs, and found that informational and metabolic functions showed relatively low dN rates; in contrast, genes from hypothetical functions and cellular processes showed high dN rates. An alternative measure of sequence variability, percentage of changes by species, was used to evaluate the most specific proportion of amino acid residues from alignments. When dN was compared with that measure a high correlation was obtained, revealing that much of evolutive information was extracted with the percentage of changes by species at the amino acid level. By analyzing the sequence variability of orthologs with a set of five properties (polarity, electrostatic charge, formation of secondary structures, molecular volume, and amino acid composition), we found that physico-chemical characteristics of proteins correlated with specific functional roles, and association of species did not follow their typical phylogeny, probably reflecting more adaptation to their life styles and niche preferences. In addition, orthologs with low dN rates had residues with more positive values of polarity, volume and electrostatic charge.

Conclusions

These findings revealed that even when orthologs perform the same function in each genomic background, their sequences reveal important evolutionary tendencies and differences related to adaptation.

This article was reviewed by: Dr. Purificación López-García, Prof. Jeffrey Townsend (nominated by Dr. J. Peter Gogarten), and Ms. Olga Kamneva.

Comparison of the 'Ca. Liberibacter asiaticus' genome adapted for an intracellular lifestyle with other members of the Rhizobiales

An intracellular plant pathogen 'Candidatus Liberibacter asiaticus,' a member of the Rhizobiales, is related to Sinorhizobium meliloti, Bradyrhizobium japonicum, nitrogen fixing endosymbionts, Agrobacterium tumefaciens, a plant pathogen, and Bartonella henselae, an intracellular mammalian pathogen. Whole chromosome comparisons identified at least 50 clusters of conserved orthologous genes found on the chromosomes of all five metabolically diverse species. The intracellular pathogens 'Ca. Liberibacter asiaticus' and Bartonella henselae have genomes drastically reduced in gene content and size as well as a relatively low content of guanine and cytosine. Codon and amino acid preferences that emphasize low guanosine and cytosine usage are globally employed in these genomes, including within regions of microsynteny and within signature sequences of orthologous proteins. The length of orthologous proteins is generally conserved, but not their isoelectric points, consistent with extensive amino acid substitutions to accommodate selection for low GC content. The 'Ca. Liberibacter asiaticus' genome apparently has all of the genes required for DNA replication present in Sinorhizobium meliloti except it has only two, rather than three RNaseH genes. The gene set required for DNA repair has only one rather than ten DNA ligases found in Sinorhizobium meliloti, and the DNA PolI of 'Ca. Liberibacter asiaticus' lacks domains needed for excision repair. Thus the ability of 'Ca. Liberibacter asiaticus' to repair mutations in its genome may be impaired. Both 'Ca. Liberibacter asiaticus and Bartonella henselae lack enzymes needed for the metabolism of purines and pyrimidines, which must therefore be obtained from the host. The 'Ca. Liberibacter asiaticus' genome also has a greatly reduced set of sigma factors used to control transcription, and lacks sigma factors 24, 28 and 38. The 'Ca. Liberibacter asiaticus' genome has all of the hallmarks of a reduced genome of a pathogen adapted to an intracellular lifestyle.

Intragenomic diversity of Rhizobium leguminosarum bv. trifolii clover nodule isolates

BACKGROUND

Soil bacteria from the genus Rhizobium are characterized by a complex genomic architecture comprising chromosome and large plasmids. Genes responsible for symbiotic interactions with legumes are usually located on one of the plasmids, named the symbiotic plasmid (pSym). The plasmids have a great impact not only on the metabolic potential of rhizobia but also underlie genome rearrangements and plasticity.

RESULTS

Here, we analyzed the distribution and sequence variability of markers located on chromosomes and extrachromosomal replicons of Rhizobium leguminosarum bv. trifolii strains originating from nodules of clover grown in the same site in cultivated soil. First, on the basis of sequence similarity of repA and repC replication genes to the respective counterparts of chromids reported in R. leguminosarum bv. viciae 3841 and R. etli CFN42, chromid-like replicons were distinguished from the pool of plasmids of the nodule isolates studied. Next, variability of the gene content was analyzed in the different genome compartments, i.e., the chromosome, chromid-like and 'other plasmids'. The stable and unstable chromosomal and plasmid genes were detected on the basis of hybridization data. Displacement of a few unstable genes between the chromosome, chromid-like and 'other plasmids', as well as loss of some markers was observed in the sampled strains. Analyses of chosen gene sequences allowed estimation of the degree of their adaptation to the three genome compartments as well as to the host.

CONCLUSIONS

Our results showed that differences in distribution and sequence divergence of plasmid and chromosomal genes can be detected even within a small group of clover nodule isolates recovered from clovers grown at the same site. Substantial divergence of genome organization could be detected especially taking into account the content of extrachromosomal DNA. Despite the high variability concerning the number and size of plasmids among the studied strains, conservation of the location as well as dynamic distribution of the individual genes (especially replication genes) of a particular genome compartment were demonstrated. The sequence divergence of particular genes may be affected by their location in the given genome compartment. The 'other plasmid' genes are less adapted to the host genome than the chromosome and chromid-like genes.

argC Orthologs from Rhizobiales show diverse profiles of transcriptional efficiency and functionality in Sinorhizobium meliloti

Several factors can influence ortholog replacement between closely related species. We evaluated the transcriptional expression and metabolic performance of ortholog substitution complementing a Sinorhizobium meliloti argC mutant with argC from Rhizobiales (Agrobacterium tumefaciens, Rhizobium etli, and Mesorhizobium loti). The argC gene is necessary for the synthesis of arginine, an amino acid that is central to protein and cellular metabolism. Strains were obtained carrying plasmids with argC orthologs expressed under the speB and argC (S. meliloti) and lac (Escherichia coli) promoters. Complementation analysis was assessed by growth, transcriptional activity, enzymatic activity, mRNA levels, specific detection of ArgC proteomic protein, and translational efficiency. The argC orthologs performed differently in each complementation, reflecting the diverse factors influencing gene expression and the ability of the ortholog product to function in a foreign metabolic background. Optimal complementation was directly related to sequence similarity with S. meliloti, and was inversely related to species signature, with M. loti argC showing the poorest performance, followed by R. etli and A. tumefaciens. Different copy numbers of genes and amounts of mRNA and protein were produced, even with genes transcribed from the same promoter, indicating that coding sequences play a role in the transcription and translation processes. These results provide relevant information for further genomic analyses and suggest that orthologous gene substitutions between closely related species are not completely functionally equivalent.

Genomic and evolutionary comparisons of diazotrophic and pathogenic bacteria of the order Rhizobiales

BACKGROUND

Species belonging to the Rhizobiales are intriguing and extensively researched for including both bacteria with the ability to fix nitrogen when in symbiosis with leguminous plants and pathogenic bacteria to animals and plants. Similarities between the strategies adopted by pathogenic and symbiotic Rhizobiales have been described, as well as high variability related to events of horizontal gene transfer. Although it is well known that chromosomal rearrangements, mutations and horizontal gene transfer influence the dynamics of bacterial genomes, in Rhizobiales, the scenario that determine pathogenic or symbiotic lifestyle are not clear and there are very few studies of comparative genomic between these classes of prokaryotic microorganisms trying to delineate the evolutionary characterization of symbiosis and pathogenesis.

RESULTS

Non-symbiotic nitrogen-fixing bacteria and bacteria involved in bioremediation closer to symbionts and pathogens in study may assist in the origin and ancestry genes and the gene flow occurring in Rhizobiales. The genomic comparisons of 19 species of Rhizobiales, including nitrogen-fixing, bioremediators and pathogens resulted in 33 common clusters to biological nitrogen fixation and pathogenesis, 15 clusters exclusive to all nitrogen-fixing bacteria and bacteria involved in bioremediation, 13 clusters found in only some nitrogen-fixing and bioremediation bacteria, 01 cluster exclusive to some symbionts, and 01 cluster found only in some pathogens analyzed. In BBH performed to all strains studied, 77 common genes were obtained, 17 of which were related to biological nitrogen fixation and pathogenesis. Phylogenetic reconstructions for Fix, Nif, Nod, Vir, and Trb showed possible horizontal gene transfer events, grouping species of different phenotypes.

CONCLUSIONS

The presence of symbiotic and virulence genes in both pathogens and symbionts does not seem to be the only determinant factor for lifestyle evolution in these microorganisms, although they may act in common stages of host infection. The phylogenetic analysis for many distinct operons involved in these processes emphasizes the relevance of horizontal gene transfer events in the symbiotic and pathogenic similarity.

Conserved symbiotic plasmid DNA sequences in the multireplicon pangenomic structure of Rhizobium etli

Strains of the same bacterial species often show considerable genomic variation. To examine the extent of such variation in Rhizobium etli, the complete genome sequence of R. etli CIAT652 and the partial genomic sequences of six additional R. etli strains having different geographical origins were determined. The sequences were compared with each other and with the previously reported genome sequence of R. etli CFN42. DNA sequences common to all strains constituted the greater part of these genomes and were localized in both the chromosome and large plasmids. About 700 to 1,000 kb of DNA that did not match sequences of the complete genomes of strains CIAT652 and CFN42 was unique to each R. etli strain. These sequences were distributed throughout the chromosome as individual genes or chromosomal islands and in plasmids, and they encoded accessory functions, such as transport of sugars and amino acids, or secondary metabolism; they also included mobile elements and hypothetical genes. Sequences corresponding to symbiotic plasmids showed high levels of nucleotide identity (about 98 to 99%), whereas chromosomal sequences and the sequences with matches to other plasmids showed lower levels of identity (on average, about 90 to 95%). We concluded that R. etli has a pangenomic structure with a core genome composed of both chromosomal and plasmid sequences, including a highly conserved symbiotic plasmid, despite the overall genomic divergence.

Metabolic capacity of Sinorhizobium (Ensifer) meliloti strains as determined by phenotype MicroArray analysis

Sinorhizobium meliloti is a soil bacterium that fixes atmospheric nitrogen in plant roots. The high genetic diversity of its natural populations has been the subject of extensive analysis. Recent genomic studies of several isolates revealed a high content of variable genes, suggesting a correspondingly large phenotypic differentiation among strains of S. meliloti. Here, using the Phenotype MicroArray (PM) system, hundreds of different growth conditions were tested in order to compare the metabolic capabilities of the laboratory reference strain Rm1021 with those of four natural S. meliloti isolates previously analyzed by comparative genomic hybridization (CGH). The results of PM analysis showed that most phenotypic differences involved carbon source utilization and tolerance to osmolytes and pH, while fewer differences were scored for nitrogen, phosphorus, and sulfur source utilization. Only the variability of the tested strain in tolerance to sodium nitrite and ammonium sulfate of pH 8 was hypothesized to be associated with the genetic polymorphisms detected by CGH analysis. Colony and cell morphologies and the ability to nodulate Medicago truncatula plants were also compared, revealing further phenotypic diversity. Overall, our results suggest that the study of functional (phenotypic) variability of S. meliloti populations is an important and complementary step in the investigation of genetic polymorphism of rhizobia and may help to elucidate rhizobial evolutionary dynamics, including adaptation to diverse environments.

Analysis of plasmid genes by phylogenetic profiling and visualization of homology relationships using Blast2Network

BACKGROUND

Phylogenetic methods are well-established bioinformatic tools for sequence analysis, allowing to describe the non-independencies of sequences because of their common ancestor. However, the evolutionary profiles of bacterial genes are often complicated by hidden paralogy and extensive and/or (multiple) horizontal gene transfer (HGT) events which make bifurcating trees often inappropriate. In this context, plasmid sequences are paradigms of network-like relationships characterizing the evolution of prokaryotes. Actually, they can be transferred among different organisms allowing the dissemination of novel functions, thus playing a pivotal role in prokaryotic evolution. However, the study of their evolutionary dynamics is complicated by the absence of universally shared genes, a prerequisite for phylogenetic analyses.

RESULTS

To overcome such limitations we developed a bioinformatic package, named Blast2Network (B2N), allowing the automatic phylogenetic profiling and the visualization of homology relationships in a large number of plasmid sequences. The software was applied to the study of 47 completely sequenced plasmids coming from Escherichia, Salmonella and Shigella spps.

CONCLUSION

The tools implemented by B2N allow to describe and visualize in a new way some of the evolutionary features of plasmid molecules of Enterobacteriaceae; in particular it helped to shed some light on the complex history of Escherichia, Salmonella and Shigella plasmids and to focus on possible roles of unannotated proteins.The proposed methodology is general enough to be used for comparative genomic analyses of bacteria.

Construction of signature-tagged mutant library in Mesorhizobium loti as a powerful tool for functional genomics

Rhizobia are nitrogen-fixing soil bacteria that establish endosymbiosis with some leguminous plants. The completion of several rhizobial genome sequences provides opportunities for genome-wide functional studies of the physiological roles of many rhizobial genes. In order to carry out genome-wide phenotypic screenings, we have constructed a large mutant library of the nitrogen-fixing symbiotic bacterium, Mesorhizobium loti, by transposon mutagenesis. Transposon insertion mutants were generated using the signature-tagged mutagenesis (STM) technique and a total of 29 330 independent mutants were obtained. Along with the collection of transposon mutants, we have determined the transposon insertion sites for 7892 clones, and confirmed insertions in 3680 non-redundant M. loti genes (50.5% of the total number of M. loti genes). Transposon insertions were randomly distributed throughout the M. loti genome without any bias toward G+C contents of insertion target sites and transposon plasmids used for the mutagenesis. We also show the utility of STM mutants by examining the specificity of signature tags and test screenings for growth- and nodulation-deficient mutants. This defined mutant library allows for genome-wide forward- and reverse-genetic functional studies of M. loti and will serve as an invaluable resource for researchers to further our understanding of rhizobial biology.

On the nature of fur evolution: a phylogenetic approach in Actinobacteria

Background

An understanding of the evolution of global transcription regulators is essential for comprehending the complex networks of cellular metabolism that have developed among related organisms. The fur gene encodes one of those regulators – the ferric uptake regulator Fur – widely distributed among bacteria and known to regulate different genes committed to varied metabolic pathways. On the other hand, members of the Actinobacteria comprise an ecologically diverse group of bacteria able to inhabit various natural environments, and for which relatively little is currently understood concerning transcriptional regulation.

Results

BLAST analyses revealed the presence of more than one fur homologue in most members of the Actinobacteria whose genomes have been fully sequenced. We propose a model to explain the evolutionary history of fur within this well-known bacterial phylum: the postulated scenario includes one duplication event from a primitive regulator, which probably had a broad range of co-factors and DNA-binding sites. This duplication predated the appearance of the last common ancestor of the Actinobacteria, while six other duplications occurred later within specific groups of organisms, particularly in two genera: Frankia and Streptomyces. The resulting paralogues maintained main biochemical properties, but became specialised for regulating specific functions, coordinating different metal ions and binding to unique DNA sequences. The presence of syntenic regions surrounding the different fur orthologues supports the proposed model, as do the evolutionary distances and topology of phylogenetic trees built using both Neighbor-Joining and Maximum-Likelihood methods.

Conclusion

The proposed fur evolutionary model, which includes one general duplication and two in-genus duplications followed by divergence and specialization, explains the presence and diversity of fur genes within the Actinobacteria. Although a few rare horizontal gene transfer events have been reported, the model is consistent with the view of gene duplication as a main force of microbial genomes evolution. The parallel study of Fur phylogeny across diverse organisms offers a solid base to guide functional studies and allows the comparison between response mechanisms in relation with the surrounding environment. The survey of regulators among related genomes provides a relevant tool for understanding the evolution of one of the first lines of cellular adaptability, control of DNA transcription.

Ferrichrome utilization in a mesorhizobial population: microevolution of a three-locus system

The ability to utilize the siderophore ferrichrome as an iron source was found to be a variable trait in a field population of mesorhizobia. To investigate the genetic basis of this variation, genes required for ferrichrome utilization (fhu genes) were characterized in Mesorhizobium strain R88B, an Fhu(+) member of the population. Functional fhu genes were present at three loci. Two genes of the ferrichrome ABC transporter, fhuBD, were identified at an fhu1 locus downstream of the symbiosis island that was integrated at the phe-tRNA gene. The fhuA gene encoding the ferrichrome outer membrane receptor was located in the fhu2 locus together with non-functional fhuDB genes, while the fhuC gene encoding the ATPase required for ferrichrome transport was part of the fhu3 locus that included genes required to form a functional TonB complex. None of the fhu genes were present in the sequenced Mesorhizobium loti strain MAFF303099. Comparisons with MAFF303099 suggested that the fhu2 and fhu3 loci evolved through small-scale (< 5 kb) acquisitions and deletions. Despite their independent origins, the three fhu loci were coordinately regulated in response to iron availability. Within the mesorhizobial population, the ability to utilize ferrichrome was most strongly correlated with the presence of the fhuA gene. We hypothesize that the ferrichrome transport system evolved through cycles of gene acquisition and deletion, with the positive selection pressure of an iron-poor or siderophore-rich environment being offset by the negative pressure of the outer membrane receptor being a target for phage.

Secrets of soil survival revealed by the genome sequence of Arthrobacter aurescens TC1

Arthrobacter sp. strains are among the most frequently isolated, indigenous, aerobic bacterial genera found in soils. Member of the genus are metabolically and ecologically diverse and have the ability to survive in environmentally harsh conditions for extended periods of time. The genome of Arthrobacter aurescens strain TC1, which was originally isolated from soil at an atrazine spill site, is composed of a single 4,597,686 basepair (bp) circular chromosome and two circular plasmids, pTC1 and pTC2, which are 408,237 bp and 300,725 bp, respectively. Over 66% of the 4,702 open reading frames (ORFs) present in the TC1 genome could be assigned a putative function, and 13.2% (623 genes) appear to be unique to this bacterium, suggesting niche specialization. The genome of TC1 is most similar to that of Tropheryma, Leifsonia, Streptomyces, and Corynebacterium glutamicum, and analyses suggest that A. aurescens TC1 has expanded its metabolic abilities by relying on the duplication of catabolic genes and by funneling metabolic intermediates generated by plasmid-borne genes to chromosomally encoded pathways. The data presented here suggest that Arthrobacter's environmental prevalence may be due to its ability to survive under stressful conditions induced by starvation, ionizing radiation, oxygen radicals, and toxic chemicals.

Soil systems contain the greatest diversity of microorganisms on earth, with 5,000–10,000 species of microorganism per gram of soil. Arthrobacter sp. strains have a primitive life cycle and are among the most frequently isolated, indigenous soil bacteria, found in common and deep subsurface soils, arctic ice, and environments contaminated with industrial chemicals and radioactive materials. To better understand how these bacteria survive in environmentally harsh conditions, the authors used a structural genomics approach to identify genes involved in soil survival of Arthrobacter aurescens strain TC1, a bacterium originally isolated for its ability to degrade the herbicide atrazine. They found that the genome of this bacterium comprises a single circular chromosome and two plasmids that encode for a large number proteins involved in stress responses due to starvation, desiccation, oxygen radicals, and toxic chemicals. A. aurescens' metabolic versatility is in part due to the presence of duplicated catabolic genes and its ability to funnel plasmid-derived intermediates into chromosomally encoded pathways. Arthrobacter's array of genes that allow for survival in stressful conditions and its ability to produce a temperature-tolerant “cyst”-like resting cell render this soil microorganism able to survive and prosper in a variety of environmental conditions.

Syntenic arrangements of the surface polysaccharide biosynthesis genes in Rhizobium leguminosarum

We applied a genomic approach in the identification of genes required for the biosynthesis of different polysaccharides in Rhizobium leguminosarum bv. trifolii TA1 (RtTA1). Pulsed-field gel electrophoresis analyses of undigested genomic DNA revealed that the RtTA1 genome is partitioned into a chromosome and four large plasmids. The combination of sequencing of RtTA1 library BAC clones and PCR amplification of polysaccharide genes from the RtTA1 genome led to the identification of five large regions and clusters, as well as many separate potential polysaccharide biosynthesis genes dispersed in the genome. We observed an apparent abundance of genes possibly linked to lipopolysaccharide biosynthesis. All RtTA1 polysaccharide biosynthesis regions showed a high degree of conserved synteny between R. leguminosarum bv. viciae and/or Rhizobium etli. A majority of the genes displaying a conserved order also showed high sequence identity levels.

Systematic analysis of head-to-head gene organization: evolutionary conservation and potential biological relevance

Several “head-to-head” (or “bidirectional”) gene pairs have been studied in individual experiments, but genome-wide analysis of this gene organization, especially in terms of transcriptional correlation and functional association, is still insufficient. We conducted a systematic investigation of head-to-head gene organization focusing on structural features, evolutionary conservation, expression correlation and functional association. Of the present 1,262, 1,071, and 491 head-to-head pairs identified in human, mouse, and rat genomes, respectively, pairs with 1– to 400–base pair distance between transcription start sites form the majority (62.36%, 64.15%, and 55.19% for human, mouse, and rat, respectively) of each dataset, and the largest group is always the one with a transcription start site distance of 101 to 200 base pairs. The phylogenetic analysis among Fugu, chicken, and human indicates a negative selection on the separation of head-to-head genes across vertebrate evolution, and thus the ancestral existence of this gene organization. The expression analysis shows that most of the human head-to-head genes are significantly correlated, and the correlation could be positive, negative, or alternative depending on the experimental conditions. Finally, head-to-head genes statistically tend to perform similar functions, and gene pairs associated with the significant cofunctions seem to have stronger expression correlations. The findings indicate that the head-to-head gene organization is ancient and conserved, which subjects functionally related genes to correlated transcriptional regulation and thus provides an exquisite mechanism of transcriptional regulation based on gene organization. These results have significantly expanded the knowledge about head-to-head gene organization. Supplementary materials for this study are available at http://www.scbit.org/h2h.

It was commonly assumed that higher eukaryotic genomes are loosely organized and genes are interspersed in the whole genome sequences. However, experiments have continuously identified eukaryotic head-to-head gene pairs with genes located closely next to each other, possibly sharing a same promoter; and preliminary genomic surveys have even proved head-to-head gene pair to be a common feature of human genome. The authors report a systematic investigation of head-to-head gene pairs in terms of the genomic structure, evolutionary conservation, expressional correlation, and functional association. The authors first identified some common structural and distributional patterns in three representative mammalian genomes: human, mouse, and rat. Then, through comparative analyses between human, chicken, and Fugu, they observed a conservation tendency of head-to-head gene pairs in vertebrates. Finally, interactive analyses of expressional and functional association yielded some interesting results, including the significant expression correlation of head-to-head genes, especially for the pairs with significant functional association. The main conclusion of this paper is that the head-to-head gene organization is ancient and conserved, subjecting functionally related genes to coregulated transcription. Lists of head-to-head gene pairs in human, mouse, rat, chicken, and Fugu are provided, while some individual pairs in need of further in-depth investigations are highlighted.