Pablo Vinuesa

GET_HOMOLOGUES, a Versatile Software Package for Scalable and Robust Microbial Pangenome Analysis

ABSTRACT GET_HOMOLOGUES is an open-source software package that builds on popular orthology-calling approaches making highly customizable and detailed pangenome analyses of microorganisms accessible to nonbioinformaticians. It can cluster homologous gene families using the bidirectional best-hit, COGtriangles, or OrthoMCL clustering algorithms. Clustering stringency can be adjusted by scanning the domain composition of proteins using the HMMER3 package, by imposing desired pairwise alignment coverage cutoffs, or by selecting only syntenic genes. The resulting homologous gene families can be made even more robust by computing consensus clusters from those generated by any combination of the clustering algorithms and filtering criteria. Auxiliary scripts make the construction, interrogation, and graphical display of core genome and pangenome sets easy to perform. Exponential and binomial mixture models can be fitted to the data to estimate theoretical core genome and pangenome sizes, and high-quality graphics can be generated. Furthermore, pangenome trees can be easily computed and basic comparative genomics performed to identify lineage-specific genes or gene family expansions. The software is designed to take advantage of modern multiprocessor personal computers as well as computer clusters to parallelize time-consuming tasks. To demonstrate some of these capabilities, we survey a set of 50 Streptococcus genomes annotated in the Orthologous Matrix (OMA) browser as a benchmark case. The package can be downloaded at http://www.eead.csic.es/compbio/soft/gethoms.php and http://maya.ccg.unam.mx/soft/gethoms.php.

Multilocus sequence analysis for assessment of the biogeography and evolutionary genetics of four Bradyrhizobium species that nodulate soybeans on the asiatic continent.

ABSTRACT A highly supported maximum-likelihood species phylogeny for the genus Bradyrhizobium was inferred from a supermatrix obtained from the concatenation of partial atpD, recA, glnII, and rpoB sequences corresponding to 33 reference strains and 76 bradyrhizobia isolated from the nodules of Glycine max (soybean) trap plants inoculated with soil samples from Myanmar, India, Nepal, and Vietnam. The power of the multigene approach using multiple strains per species was evaluated in terms of overall tree resolution and phylogenetic congruence, representing a practical and portable option for bacterial molecular systematics. Potential pitfalls of the approach are highlighted. Seventy-five of the isolates could be classified as B. japonicum type Ia (USDA110/USDA122-like), B. liaoningense, B. yuanmingense, or B. elkanii, whereas one represented a novel Bradyrhizobium lineage. Most Nepalese B. japonicum Ia isolates belong to a highly epidemic clone closely related to strain USDA110. Significant phylogenetic evidence against the monophyly of the of B. japonicum I and Ia lineages was found. Analysis of their DNA polymorphisms revealed high population distances, significant genetic differentiation, and contrasting population genetic structures, suggesting that the strains in the Ia lineage are misclassified as B. japonicum. The DNA polymorphism patterns of all species conformed to the expectations of the neutral mutation and population equilibrium models and, excluding the B. japonicum Ia lineage, were consistent with intermediate recombination levels. All species displayed epidemic clones and had broad geographic and environmental distribution ranges, as revealed by mapping climate types and geographic origins of the isolates on the species tree.

Specific Detection of Bradyrhizobium and Rhizobium Strains Colonizing Rice (Oryza sativa) Roots by 16S-23S Ribosomal DNA Intergenic Spacer-Targeted PCR

ABSTRACT In addition to forming symbiotic nodules on legumes, rhizobial strains are members of soil or rhizosphere communities or occur as endophytes, e.g., in rice. Two rhizobial strains which have been isolated from root nodules of the aquatic legumes Aeschynomene fluminensis (IRBG271) and Sesbania aculeata(IRBG74) were previously found to promote rice growth. In addition to analyzing their phylogenetic positions, we assessed the suitability of the 16S-23S ribosomal DNA (rDNA) intergenic spacer (IGS) sequences for the differentiation of closely related rhizobial taxa and for the development of PCR protocols allowing the specific detection of strains in the environment. 16S rDNA sequence analysis (sequence identity, 99%) and phylogenetic analysis of IGS sequences showed that strain IRBG271 was related to but distinct from Bradyrhizobium elkanii. Rhizobium sp. (Sesbania) strain IRBG74 was located in theRhizobium-Agrobacterium cluster as a novel lineage according to phylogenetic 16S rDNA analysis (96.8 to 98.9% sequence identity with Agrobacterium tumefaciens;emended name, Rhizobium radiobacter). Strain IRBG74 harbored four copies of rRNA operons whose IGS sequences varied only slightly (2 to 9 nucleotides). The IGS sequence analyses allowed intraspecies differentiation, especially in the genusBradyrhizobium, as illustrated here for strains ofBradyrhizobium japonicum, B. elkanii, Bradyrhizobium liaoningense, and Bradyrhizobium sp. (Chamaecytisus) strain BTA-1. It also clearly differentiated fast-growing rhizobial species and strains, albeit with lower statistical significance. Moreover, the high sequence variability allowed the development of highly specific IGS-targeted nested-PCR assays. Strains IRBG74 and IRBG271 were specifically detected in complex DNA mixtures of numerous related bacteria and in the DNA of roots of gnotobiotically cultured or even of soil-grown rice plants after inoculation. Thus, IGS sequence analysis is an attractive technique for both microbial ecology and systematics.

Evolutionary genetics and biogeographic structure of Rhizobium gallicum sensu lato , a widely distributed bacterial symbiont of diverse legumes

We used phylogenetic and population genetics approaches to evaluate the importance of the evolutionary forces on shaping the genetic structure of Rhizobium gallicum and related species. We analysed 54 strains from several populations distributed in the Northern Hemisphere, using nucleotide sequences of three ‘core’ chromosomal genes (rrs, glnII and atpD) and two ‘auxiliary’ symbiotic genes (nifH and nodB) to elucidate the biogeographic history of the species and symbiotic ecotypes (biovarieties) within species. The analyses revealed that strains classified as Rhizobium mongolense and Rhizobium yanglingense belong to the chromosomal evolutionary lineage of R. gallicum and harbour symbiotic genes corresponding to a new biovar; we propose their reclassification as R. gallicum bv. orientale. The comparison of the chromosomal and symbiotic genes revealed evidence of lateral transfer of symbiotic information within and across species. Genetic differentiation analyses based on the chromosomal protein‐coding genes revealed a biogeographic pattern with three main populations, whereas the 16S rDNA sequences did not resolve that biogeographic pattern. Both the phylogenetic and population genetic analyses showed evidence of recombination at the rrs locus. We discuss our results in the light of the contrasting views of bacterial species expressed by microbial taxonomist and evolutionary biologists.

Rhizobium etli and Rhizobium gallicum Nodulate Common Bean (Phaseolus vulgaris) in a Traditionally Managed Milpa Plot in Mexico: Population Genetics and Biogeographic Implications

ABSTRACT The stability of the genetic structure of rhizobial populations nodulating Phaseolus vulgaris cultivated in a traditionally managed milpa plot in Mexico was studied over three consecutive years. The set of molecular markers analyzed (including partial rrs, glnII, nifH, and nodB sequences), along with host range experiments, placed the isolates examined in Rhizobium etli bv. phaseoli and Rhizobium gallicum bv. gallicum. Cluster analysis of multilocus enzyme electrophoresis and plasmid profile data separated the two species and identified numerically dominant clones within each of them. Population genetic analyses showed that there was high genetic differentiation between the two species and that there was low intrapopulation differentiation of the species over the 3 years. The results of linkage disequilibrium analyses are consistent with an epidemic genetic structure for both species, with frequent genetic exchange taking place within conspecific populations but not between the R. etli and R. gallicum populations. A subsample of isolates was selected and used for 16S ribosomal DNA PCR-restriction fragment length polymorphism analysis, nifH copy number determination, and host range experiments. Plasmid profiles and nifH hybridization patterns also revealed the occurrence of lateral plasmid transfer among distinct multilocus genotypes within species but not between species. Both species were recovered from nodules of the same plants, indicating that mechanisms other than host, spatial, or temporal isolation may account for the genetic barrier between the species. The biogeographic implications of finding an R. gallicum bv. gallicum population nodulating common bean in America are discussed.

Nodulation of Sesbania species by Rhizobium (Agrobacterium) strain IRBG74 and other rhizobia.

Summary Concatenated sequence analysis with 16S rRNA, rpoB and fusA genes identified a bacterial strain (IRBG74) isolated from root nodules of the aquatic legume Sesbania cannabina as a close relative of the plant pathogen Rhizobium radiobacter (syn. Agrobacterium tumefaciens). However, DNA:DNA hybridization with R. radiobacter, R. rubi, R. vitis and R. huautlense gave only 44%, 5%, 8% and 8% similarity respectively, suggesting that IRBG74 is potentially a new species. Additionally, it contained no vir genes and lacked tumour‐forming ability, but harboured a sym‐plasmid containing nifH and nodA genes similar to those in other Sesbania symbionts. Indeed, IRBG74 effectively nodulated S. cannabina and seven other Sesbania spp. that nodulate with Ensifer (Sinorhizobium)/Rhizobium strains with similar nodA genes to IRBG74, but not species that nodulate with Azorhizobium or Mesorhizobium. Light and electron microscopy revealed that IRBG74 infected Sesbania spp. via lateral root junctions under flooded conditions, but via root hairs under non‐flooded conditions. Thus, IRBG74 is the first confirmed legume‐nodulating symbiont from the Rhizobium (Agrobacterium) clade. Cross‐inoculation studies with various Sesbania symbionts showed that S. cannabina could form fully effective symbioses with strains in the genera Rhizobium and Ensifer, only ineffective ones with Azorhizobium strains, and either partially effective (Mesorhizobium huakii) or ineffective (Mesorhizobium plurifarium) symbioses with Mesorhizobium. These data are discussed in terms of the molecular phylogeny of Sesbania and its symbionts.

The lipid lysyl-phosphatidylglycerol is present in membranes of Rhizobium tropici CIAT899 and confers increased resistance to polymyxin B under acidic growth conditions

Lysyl-phosphatidylglycerol (LPG) is a well-known membrane lipid in several gram-positive bacteria but is almost unheard of in gram-negative bacteria. In Staphylococcus aureus, the gene product of mprF is responsible for LPG formation. Low pH-inducible genes, termed IpiA, have been identified in the gram-negative alpha-proteobacteria Rhizobium tropici and Sinorhizobium medicae in screens for acid-sensitive mutants and they encode homologs of MprF. An analysis of the sequenced bacterial genomes reveals that genes coding for homologs of MprF from S. aureus are present in several classes of organisms throughout the bacterial kingdom. In this study, we show that the expression of lpiA from R. tropici in the heterologous hosts Escherichia coli and Sinorhizobium meliloti causes formation of LPG. A wild-type strain of R. tropici forms LPG (about 1% of the total lipids) when the cells are grown in minimal medium at pH 4.5 but not when grown in minimal medium at neutral pH or in complex tryptone yeast (TY) medium at either pH. LPG biosynthesis does not occur when lpiA is deleted and is restored upon complementation of lpiA-deficient mutants with a functional copy of the lpiA gene. When grown in the low-pH medium, lpiA-deficient rhizobial mutants are over four times more susceptible to the cationic peptide polymyxin B than the wild type.

A Complete Set of Flagellar Genes Acquired by Horizontal Transfer Coexists with the Endogenous Flagellar System in Rhodobacter sphaeroides

Bacteria swim in liquid environments by means of a complex rotating structure known as the flagellum. Approximately 40 proteins are required for the assembly and functionality of this structure. Rhodobacter sphaeroides has two flagellar systems. One of these systems has been shown to be functional and is required for the synthesis of the well-characterized single subpolar flagellum, while the other was found only after the genome sequence of this bacterium was completed. In this work we found that the second flagellar system of R. sphaeroides can be expressed and produces a functional flagellum. In many bacteria with two flagellar systems, one is required for swimming, while the other allows movement in denser environments by producing a large number of flagella over the entire cell surface. In contrast, the second flagellar system of R. sphaeroides produces polar flagella that are required for swimming. Expression of the second set of flagellar genes seems to be positively regulated under anaerobic growth conditions. Phylogenic analysis suggests that the flagellar system that was initially characterized was in fact acquired by horizontal transfer from a gamma-proteobacterium, while the second flagellar system contains the native genes. Interestingly, other alpha-proteobacteria closely related to R. sphaeroides have also acquired a set of flagellar genes similar to the set found in R. sphaeroides, suggesting that a common ancestor received this gene cluster.

A comparative structural bioinformatics analysis of the insulin receptor family ectodomain based on phylogenetic information

The insulin receptor (IR), the insulin-like growth factor 1 receptor (IGF1R) and the insulin receptor-related receptor (IRR) are covalently-linked homodimers made up of several structural domains. The molecular mechanism of ligand binding to the ectodomain of these receptors and the resulting activation of their tyrosine kinase domain is still not well understood. We have carried out an amino acid residue conservation analysis in order to reconstruct the phylogeny of the IR Family. We have confirmed the location of ligand binding site 1 of the IGF1R and IR. Importantly, we have also predicted the likely location of the insulin binding site 2 on the surface of the fibronectin type III domains of the IR. An evolutionary conserved surface on the second leucine-rich domain that may interact with the ligand could not be detected. We suggest a possible mechanical trigger of the activation of the IR that involves a slight ‘twist’ rotation of the last two fibronectin type III domains in order to face the likely location of insulin. Finally, a strong selective pressure was found amongst the IRR orthologous sequences, suggesting that this orphan receptor has a yet unknown physiological role which may be conserved from amphibians to mammals.

A ClC Chloride Channel Homolog and Ornithine-Containing Membrane Lipids of Rhizobium tropici CIAT899 Are Involved in Symbiotic Efficiency and Acid Tolerance

Rhizobium tropici CIAT899 is highly tolerant to several environmental stresses and is a good competitor for nodule occupancy of common bean plants in acid soils. Random transposon mutagenesis was performed to identify novel genes of this strain involved in symbiosis and stress tolerance. Here, we present a genetic analysis of the locus disrupted by the Tn5 insertion in mutant 899-PV9, which lead to the discovery of sycA, a homolog of the ClC family of chloride channels and Cl-/H+ exchange transporters. A nonpolar deletion in this gene caused serious deficiencies in nodule development, nodulation competitiveness, and N2 fixation on Phaseolus vulgaris plants, probably due to its reduced ability to invade plant cells and to form stable symbiosomes, as judged by electron transmission microscopy. A second gene (olsC), found downstream of sycA, is homologous to aspartyl/asparaginyl beta-hydroxylases and modifies two species of ornithine-containing lipids in vivo, presumably by hydroxylation at a still-unknown position. A mutant carrying a nonpolar deletion in olsC is symbiotically defective, whereas overexpressed OlsC in the complemented strain provokes an acid-sensitive phenotype. This is the first report of a ClC homolog being essential for the establishment of a fully developed N2-fixing root nodule symbiosis and of a putative beta-hydroxylase that modifies ornithine-containing membrane lipids of R. tropici CIAT899, which, in turn, are contributing to symbiotic performance and acid tolerance.