Two-component signal transduction systems and regulation of virulence factors in <Emphasis Type="Italic">Xanthomonas</Emphasis>: a perspective

Two-component signal transduction systems (TCSTSs), consisting of a histidine kinase and a response regulator, play a critical role in regulating virulence gene expression in Gram-negative phytopathogenic bacteria Xanthomonas spp.. To date, 12 TCSTS genes have been identified, accounting for approximately 10% of the TCSTS genes in each genome that have been experimentally identified to be related to pathogenesis. These TCSTSs modulate the expression of a number of virulence factors through diverse molecular mechanisms such as interacting with DNA, protein-binding and involvement in second messenger metabolism, which generates a high level of regulatory versatility. Here we summarize the current knowledge in this field and discuss the emerging themes and remaining questions that are important in deciphering the signaling network of TCSTSs in Xanthomonas.     

A genomic analysis of two-component signal transduction in Streptococcus pneumoniae

A genomics-based approach was used to identify the entire gene complement of putative two-component signal transduction systems (TCSTSs) in Streptococcus pneumoniae. A total of 14 open reading frames (ORFs) were identified as putative response regulators, 13 of which were adjacent to genes encoding probable histidine kinases. Both the histidine kinase and response regulator proteins were categorized into subfamilies on the basis of phylogeny. Through a systematic programme of mutagenesis, the importance of each novel TCSTS was determined with respect to viability and pathogenicity. One TCSTS was identified that was essential for the growth of S. pneumoniaeThis locus was highly homologous to the yycFG gene pair encoding the essential response regulator/histidine kinase proteins identified in Bacillus subtilis and Staphylococcus aureus. Separate deletions of eight other loci led in each case to a dramatic attenuation of growth in a mouse respiratory tract infection model, suggesting that these signal transduction systems are important for the in vivo adaptation and pathogenesis of S. pneumoniae. The identification of conserved TCSTSs important for both pathogenicity and viability in a Gram-positive pathogen highlights the potential of two-component signal transduction as a multicomponent target for antibacterial drug discovery.     

Genome-scale mutagenesis and phenotypic characterization of two-component signal transduction systems in Xanthomonas campestris pv. campestris ATCC 33913

The gram-negative bacterium Xanthomonas campestris pv. campestris is the causal agent of black rot disease of cruciferous plants. Its genome encodes a large repertoire of two-component signal transduction systems (TCSTSs), which consist of histidine kinases and response regulators (RR) to monitor and respond to environmental stimuli. To investigate the biological functions of these TCSTS genes, we aimed to inactivate all 54 RR genes in X. campestris pv. campestris ATCC 33913, and successfully generated 51 viable mutants using the insertion inactivation method. Plant inoculation identified two novel response regulator genes (XCC1958 and XCC3107) that are involved in virulence of this strain. Genetic complementation demonstrated that XCC3107, designated as vgrR (virulence and growth regulator), also affects bacterial growth and activity of extracellular proteases. In addition, we assessed the survival of these mutants under various stresses, including osmotic stress, high sodium concentration, heat shock, and sodium dodecyl sulfate exposure, and identified a number of genes that may be involved in the general stress response of X. campestris pv. campestris. Mutagenesis and phenotypic characterization of RR genes in this study will facilitate future studies on signaling networks in this important phytopathogenic bacterium.     

Two-Component Signal Transduction Systems of Desulfovibrio vulgaris: Structural and Phylogenetic Analysis and Deduction of Putative Cognate Pairs

A large number of two-component signal transduction systems (TCSTS) including 59 putative sensory histidine kinases (HK) and 55 response regulators (RR) were identified from the Desulfovibrio vulgaris genome. In this study, the structural and phylogenetic analyses of all putative TCSTSs in D. vulgaris were performed. The results showed that D. vulgaris contained 21 hybrid-type HKs, implying that multiple-step phosphorelay may be a common signal transduction mechanism in D. vulgaris. Despite the low sequence similarity that restricted the resolution of the phylogenetic analyses, most TCSTS components of D. vulgaris were found clustered into several subfamilies previously recognized in Escherichia coli and Bacillus subtilis. An attempt was made in this investigation to identify the possible cognate HK-RR pairs not linked on the chromosome in D. vulgaris based on similar expression patterns in response to various environmental and genetic changes. Expression levels of all HK and RR genes were measured using whole-genome microarrays. Five groups of HK-RR genes not linked on the chromosome were identified as possible cognate pairs in D. vulgaris. The results provided a preliminary list of possible cognate HK-RR pairs and constitute a basis for further exploration of interaction and physiological function of TCSTSs in D. vulgaris. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Gail Simmons]     

Relationship between genome similarity and DNA-DNA hybridization among closely related bacteria

DNA-DNA hybridization has been established as an important technology in bacterial species taxonomy and phylogenetic analysis. In this study, we analyzed how the efficiency with which the genomic DNA from one species hybridizes to the genomic DNA of another species (DNA-DNA hybridization) in microarray analysis relates to the similarity between two genomes. We found that the predicted DNA-DNA hybridization based on genome sequence similarity correlated well with the experimentally determined microarray hybridization. Between closely related strains, significant numbers of highly divergent genes (<55% identity) and/or the accumulation of mismatches between conserved genes lowered the DNA-DNA hybridization signal, and this reduced the hybridization signals to below 70% for even bacterial strains with over 97% 16S rRNA gene identity. In addition, our results also suggest that a DNA-DNA hybridization signal intensity of over 40% indicates that two genomes at least shared 30% conserved genes (>60% gene identity). This study may expand our knowledge of DNA-DNA hybridization based on genomic sequence similarity comparison and further provide insights for bacterial phylogeny analyses.     

Non-gamma-proteobacteria gene islands contribute to the Xanthomonas genome

Horizontal gene transfer, a process through which genomes acquire sequences from distantly related organisms, is believed to be a major source of genetic diversity in bacteria. A central question concerning the impact of gene transfer on bacterial genome evolution is the proportion of horizontally transferred sequences within genomes. Through BLAST search, we found that the genomes of two phytopathogens, Xanthomonas campestris pv. campestris and Xanthomonas axonopodis pv. citri, have close to 40% of the genes with the highest similarity to genes from phylogenetically distant organisms (non-gamma-proteobacteria). Most of these genes are found to be contiguous in the genome, forming genome islands, which may have been transferred from other organisms. Overall, the total number of genes within genome islands corresponds to almost one quarter of the entire xanthomonad genomes. Interestingly, many of the genes in these islands are functionally related to plant pathogenesis and virulence. Thus, these results suggest that horizontally transferred genes are clustered in the genome, and may facilitate fitness in new environments, as in the case of plant-bacteria interaction.     

Building interacting partner predictors using co-varying residue pairs between histidine kinase and response regulator pairs of 48 bacterial two-component systems

The two‐component system (TCS) is a signal transduction system that involves a histidine kinase (HK) and a response regulator (RR). Although up to hundreds of TCSs may operate in parallel in a bacterial cell, the high‐fidelity of a TCS signaling is well maintained, minimizing irrelevant crosstalk between TCSs. When a HK gene and a RR gene in a given TCS system exist in neighboring positions, it is almost certain that their protein products (i.e., HK and RR) are interacting partners. However, large bacterial genomes often have multiple HK genes and/or cognate RR genes that are not neighboring positions. In many partially assembled genomes, some HK genes and RR genes belong to different contigs. In these cases, it is not clear which HK(s) and RR(s) interact. By combining information‐theoretic and graph‐theoretic approaches, we developed a computational method identifying co‐evolving residue pairs between HKs and cognate RRs and predicting the interacting HK:RR pairs for each TCS. In addition, we built a TCSppWWW webserver ( http://compath.org/platcom/tcs ) that takes query sequences of pairing candidates and predicts their HK:RR pairing using precomputed models. The current release of TCSppWWW provides predictors for 48 TCSs using over 20,000 protein sequences from about 900 bacterial genomes. Three different types of predictors using Random Forest, RBF Network, and Naïve Bayes are provided. Once a set of HK and RR candidate sequences are submitted, TCSppWWW aligns query sequences to the precomputed multiple sequence alignment of HK:RR pairs, extracts co‐evolving column positions, then returns prediction results with prediction margin and additional information. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

New features of Asian Crassostrea oyster mitochondrial genomes: A novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.     

Chromosome assignment of six dog genes by FISH, and correlation with dog-human Zoo-FISH data

Cross-species chromosome painting analyses have recently demonstrated the presence of regions of conserved synteny between the human and domestic dog genomes, aiding the search for candidate genes for inherited traits. Concerted efforts to subchromosomally assign substantial numbers of dog gene sequences are now needed in order to refine these comparative data, both in terms of marker density and resolution. We have developed novel PCR markers representing three dog genes (ALB, FOS, HNRPA2B1) for which no sequence or mapping data were previously available, to our knowledge. These, in addition to three gene markers previously described (ALDOA, RPE65, VCAM1), were used to isolate and chromosomally assign corresponding large insert genomic clones by fluorescence in situ hybridization (FISH). Chromosome assignments for these six dog genes are discussed in terms of those of the human orthologues, and correlated with existing comparative mapping information, identifying one apparent exception to existing Zoo-FISH data, and aiding refinement of the boundaries of conserved chromosome segments in both genomes.     

Global analysis of the eukaryotic pathways and networks regulated by Salmonella typhimurium in mouse intestinal infection in vivo

Background

Citrus canker is a disease that has severe economic impact on the citrus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri, whose genome sequence was made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C.

Results

We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on type III secreted effector repertoires. In addition to pthA, already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI, are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content differences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, type IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP, which codes for a natriuretic protein.

Conclusion

We have uncovered numerous commonalities and differences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. The gained knowledge will be instrumental for improving citrus canker control.     

Gene identification in novel eukaryotic genomes by self-training algorithm

Finding new protein-coding genes is one of the most important goals of eukaryotic genome sequencing projects. However, genomic organization of novel eukaryotic genomes is diverse and ab initio gene finding tools tuned up for previously studied species are rarely suitable for efficacious gene hunting in DNA sequences of a new genome. Gene identification methods based on cDNA and expressed sequence tag (EST) mapping to genomic DNA or those using alignments to closely related genomes rely either on existence of abundant cDNA and EST data and/or availability on reference genomes. Conventional statistical ab initio methods require large training sets of validated genes for estimating gene model parameters. In practice, neither one of these types of data may be available in sufficient amount until rather late stages of the novel genome sequencing. Nevertheless, we have shown that gene finding in eukaryotic genomes could be carried out in parallel with statistical models estimation directly from yet anonymous genomic DNA. The suggested method of parallelization of gene prediction with the model parameters estimation follows the path of the iterative Viterbi training. Rounds of genomic sequence labeling into coding and non-coding regions are followed by the rounds of model parameters estimation. Several dynamically changing restrictions on the possible range of model parameters are added to filter out fluctuations in the initial steps of the algorithm that could redirect the iteration process away from the biologically relevant point in parameter space. Tests on well-studied eukaryotic genomes have shown that the new method performs comparably or better than conventional methods where the supervised model training precedes the gene prediction step. Several novel genomes have been analyzed and biologically interesting findings are discussed. Thus, a self-training algorithm that had been assumed feasible only for prokaryotic genomes has now been developed for ab initio eukaryotic gene identification.     

The genome of Xanthomonas campestris pv. campestris B100 and its use for the reconstruction of metabolic pathways involved in xanthan biosynthesis

The complete genome sequence of the Xanthomonas campestris pv. campestris strain B100 was established. It consisted of a chromosome of 5,079,003bp, with 4471 protein-coding genes and 62 RNA genes. Comparative genomics showed that the genes required for the synthesis of xanthan and xanthan precursors were highly conserved among three sequenced X. campestris pv. campestris genomes, but differed noticeably when compared to the remaining four Xanthomonas genomes available. For the xanthan biosynthesis genes gumB and gumK earlier translational starts were proposed, while gumI and gumL turned out to be unique with no homologues beyond the Xanthomonas genomes sequenced. From the genomic data the biosynthesis pathways for the production of the exopolysaccharide xanthan could be elucidated. The first step of this process is the uptake of sugars serving as carbon and energy sources wherefore genes for 15 carbohydrate import systems could be identified. Metabolic pathways playing a role for xanthan biosynthesis could be deduced from the annotated genome. These reconstructed pathways concerned the storage and metabolization of the imported sugars. The recognized sugar utilization pathways included the Entner-Doudoroff and the pentose phosphate pathway as well as the Embden-Meyerhof pathway (glycolysis). The reconstruction indicated that the nucleotide sugar precursors for xanthan can be converted from intermediates of the pentose phosphate pathway, some of which are also intermediates of glycolysis or the Entner-Doudoroff pathway. Xanthan biosynthesis requires in particular the nucleotide sugars UDP-glucose, UDP-glucuronate, and GDP-mannose, from which xanthan repeat units are built under the control of the gum genes. The updated genome annotation data allowed reconsidering and refining the mechanistic model for xanthan biosynthesis.     

Intraphylum diversity and complex evolution of cyanobacterial aminoacyl-tRNA synthetases

A comparative genomic analysis of 35 cyanobacterial strains has revealed that the gene complement of aminoacyl-tRNA synthetases (AARSs) and routes for aminoacyl-tRNA synthesis may differ among the species of this phylum. Several genes encoding AARS paralogues were identified in some genomes. In-depth phylogenetic analysis was done for each of these proteins to gain insight into their evolutionary history. GluRS, HisRS, ArgRS, ThrRS, CysRS, and Glu-Q-RS showed evidence of a complex evolutionary course as indicated by a number of inconsistencies with our reference tree for cyanobacterial phylogeny. In addition to sequence data, support for evolutionary hypotheses involving horizontal gene transfer or gene duplication events was obtained from other observations including biased sequence conservation, the presence of indels (insertions or deletions), or vestigial traces of ancestral redundant genes. We present evidences for a novel protein domain with two putative transmembrane helices recruited independently by distinct AARS in particular cyanobacteria.     

An HMM-Based Comparative Genomic Framework for Detecting Introgression in Eukaryotes

One outcome of interspecific hybridization and subsequent effects of evolutionary forces is introgression, which is the integration of genetic material from one species into the genome of an individual in another species. The evolution of several groups of eukaryotic species has involved hybridization, and cases of adaptation through introgression have been already established. In this work, we report on PhyloNet-HMM—a new comparative genomic framework for detecting introgression in genomes. PhyloNet-HMM combines phylogenetic networks with hidden Markov models (HMMs) to simultaneously capture the (potentially reticulate) evolutionary history of the genomes and dependencies within genomes. A novel aspect of our work is that it also accounts for incomplete lineage sorting and dependence across loci. Application of our model to variation data from chromosome 7 in the mouse (Mus musculus domesticus) genome detected a recently reported adaptive introgression event involving the rodent poison resistance gene Vkorc1, in addition to other newly detected introgressed genomic regions. Based on our analysis, it is estimated that about 9% of all sites within chromosome 7 are of introgressive origin (these cover about 13 Mbp of chromosome 7, and over 300 genes). Further, our model detected no introgression in a negative control data set. We also found that our model accurately detected introgression and other evolutionary processes from synthetic data sets simulated under the coalescent model with recombination, isolation, and migration. Our work provides a powerful framework for systematic analysis of introgression while simultaneously accounting for dependence across sites, point mutations, recombination, and ancestral polymorphism.     

A novel group of type I polyketide synthases (PKS) in animals and the complex phylogenomics of PKSs

Type I polyketide synthases (PKSs), and related fatty acid synthases (FASs), represent a large group of proteins encoded by a diverse gene family that occurs in eubacteria and eukaryotes (mainly in fungi). Collectively, enzymes encoded by this gene family produce a wide array of polyketide compounds that encompass a broad spectrum of biological activity including antibiotic, antitumor, antifungal, immunosuppressive, and predator defense functional roles. We employed a phylogenomics approach to estimate relationships among members of this gene family from eubacterial and eukaryotic genomes. Our results suggest that some animal genomes (sea urchins, birds, and fish) possess a previously unidentified group of pks genes, in addition to possessing fas genes used in fatty acid metabolism. These pks genes in the chicken, fish, and sea urchin genomes do not appear to be closely related to any other animal or fungal genes, and instead are closely related to pks genes from the slime mold Dictyostelium and eubacteria. Continued accumulation of genome sequence data from diverse animal lineages is required to clarify whether the presence of these (non-fas) pks genes in animal genomes owes their origins to horizontal gene transfer (from eubacterial or Dictostelium genomes) or to more conventional patterns of vertical inheritance coupled with massive gene loss in several animal lineages. Additionally, results of our broad-scale phylogenetic analyses bolster the support for previous hypotheses of horizontal gene transfer of pks genes from bacterial to fungal and protozoan lineages.