Combined In Silico and In Vitro Study to Reveal the Structural Insights and Nucleotide-Binding Ability of the Transcriptional Regulator PehR from the Phytopathogen Ralstonia solanacearum.

The transcriptional regulator PehR regulates the synthesis of the extracellular plant cell wall-degrading enzyme polygalacturonase, which is essential in the bacterial wilt of plants caused by one of the most devastating plant phytopathogens, Ralstonia solanacearum. The bacterium has a wide global distribution infecting many different plant species, resulting in massive agricultural and economic losses. Because the PehR molecular structure has not yet been determined and the structural consequences of PehR on ligand binding have not been thoroughly investigated, we have used an in silico approach combined with in vitro experiments for the first time to characterize the PehR regulator from a local isolate (Tezpur, Assam, India) of the phytopathogenic bacterium R. solanacearum F1C1. In this study, an in silico approach was employed to model the 3D structure of the PehR regulator, followed by the binding analysis of different ligands against this regulatory protein. Molecular docking studies suggest that ATP has the highest binding affinity for the PehR regulator. By using molecular dynamics (MD) simulation analysis, involving root-mean-square deviation, root-mean-square fluctuations, hydrogen bonding, radius of gyration, solvent-accessible surface area, and principal component analysis, it was possible to confirm the sudden conformational changes of the PehR regulator caused by the presence of ATP. We used an in vitro approach to further validate the formation of the PehR-ATP complex. In this approach, recombinant DNA technology was used to clone, express, and purify the gene encoding the PehR regulator from R. solanacearum F1C1. Purified PehR was used in ATP-binding experiments using fluorescence spectroscopy and Fourier transform infrared spectroscopy, the outcomes of which showed a potent binding to ATP. The putative PehR-ATP-binding analysis revealed the importance of the amino acids Lys190, Glu191, Arg192, Arg375, and Asp378 for the ATP-binding process, but further study is required to confirm this. It will be simpler to comprehend the catalytic mechanisms of a crucial PehR regulator process in R. solanacearum with the aid of the ATP-binding process hints provided by these structural biology applications.

Microliter spotting and micro-colony observation: A rapid and simple approach for counting bacterial colony forming units.

For enumerating viable bacteria, traditional dilution plating to count colony forming units (CFUs) has always been the preferred method in microbiology owing to its simplicity, albeit being laborious and time-consuming. Similar CFU counts can be obtained by quantifying growing micro-colonies in conjunction with the benefits of a microscope. Here, we employed a simple method of five to ten microliter spotting of a diluted bacterial culture multiple times on a single Petri dish followed by determining CFU by counting micro-colonies using a phase-contrast microscope. In this method, the CFU of an Escherichia coli culture can be estimated within a four-hour period after spotting. Further, within a ten-hour period after spotting, CFU in a culture of Ralstonia solanacearum, a bacterium with a generation time of around 2 h, can be estimated. The CFU number determined by micro-colonies observed for 106-fold dilutions or lower is similar to that obtained by the dilution plating method for 107-fold dilutions or lower. Micro-colony numbers observed in the early hours of growth (2 h in case of E. coli and 8 h in case of R. solanacearum) were found to remain consistent at later hours (4 h in case of E. coli and 10 h in case of R. solanacearum), where the visibility of the colonies was better due to a noticeable increase in the size of the colonies. This suggested that micro-colonies observed in the early hours indeed represent the bacterial number in the culture. Practical applications to this counting method were employed in studying the rifampicin-resistant mutation rate as well as performing a fluctuation test in E. coli. The spotting method described here to enumerate bacterial CFU results in reduction of labour, time and resources.

Machine learning approach identifies prominent codons from different degenerate groups influencing gene expression in bacteria.

Unequal usage of synonymous codons is known as codon usage bias (CUB), which is generally different between the high-expression genes (HEG) and low-expression genes (LEG) in organisms is not yet adequately reported across different bacteria. In this study, a machine learning-based approach was implemented initially to find out codons that are significantly different between the HEG and LEG in Escherichia coli. It identified Cys codons such as UGU and UGC, Lys codons such as AAA and AAG that were least influenced by gene expression. Codons such as UCU (Ser), CUG (Leu), GGG (Gly), CGG (Arg) etc. were identified to be influenced maximum by the gene expression. The study was extended to analyze codon usage in 683 other bacterial species. Cys (UGU/UGC) and Ser (AGU/AGC) codons were identified being the least different between the two groups of genes across these bacterial species. Codons such as CGA, CUG, GGG, GCC, ACC, AUA, and AUC were identified to be influenced by the gene expression across majority of these species. This study supports the role of CUB on gene expression across bacteria and demonstrates a commonality among bacteria regarding behavior of certain codons with regard to gene expression.

Incorporation of transition to transversion ratio and nonsense mutations, improves the estimation of the number of synonymous and non-synonymous sites in codons.

A common approach to estimate the strength and direction of selection acting on protein coding sequences is to calculate the dN/dS ratio. The method to calculate dN/dS has been widely used by many researchers and many critical reviews have been made on its application after the proposition by Nei and Gojobori in 1986. However, the method is still evolving considering the non-uniform substitution rates and pretermination codons. In our study of SNPs in 586 genes across 156 Escherichia coli strains, synonymous polymorphism in 2-fold degenerate codons were higher in comparison to that in 4-fold degenerate codons, which could be attributed to the difference between transition (Ti) and transversion (Tv) substitution rates where the average rate of a transition is four times more than that of a transversion in general. We considered both the Ti/Tv ratio, and nonsense mutation in pretermination codons, to improve estimates of synonymous (S) and non-synonymous (NS) sites. The accuracy of estimating dN/dS has been improved by considering the Ti/Tv ratio and nonsense substitutions in pretermination codons. We showed that applying the modified approach based on Ti/Tv ratio and pretermination codons results in higher values of dN/dS in 29 common genes of equal reading-frames between E. coli and Salmonella enterica. This study emphasizes the robustness of amino acid composition with varying codon degeneracy, as well as the pretermination codons when calculating dN/dS values.

Glycome Profiling and Bioprospecting Potential of the Himalayan Buddhist Handmade Paper of Tawang Region of Arunachal Pradesh.

The paper and pulp industry (PPI) is one of the largest industries that contribute to the growing economy of the world. While wood remains the primary raw material of the PPIs, the demand for paper has also grown alongside the expanding global population, leading to deforestation and ecological imbalance. Wood-based paper production is associated with enormous utilization of water resources and the release of different wastes and untreated sludge that degrades the quality of the environment and makes it unsafe for living creatures. In line with this, the indigenous handmade paper making from the bark of Daphne papyracea, Wall. ex G. Don by the Monpa tribe of Arunachal Pradesh, India is considered as a potential alternative to non-wood fiber. This study discusses the species distribution modeling of D. papyracea, community-based production of the paper, and glycome profiling of the paper by plant cell wall glycan-directed monoclonal antibodies. The algorithms used for ecological and geographical modeling indicated the maximum predictive distribution of the plant toward the western parts of Arunachal Pradesh. It was also found that the suitable distribution of D. papyracea was largely affected by the precipitation and temperature variables. Plant cell walls are primarily made up of cellulose, hemicellulose, lignin, pectin, and glycoproteins. Non-cellulosic cell wall glycans contribute significantly to various physical properties such as density, crystallinity, and tensile strength of plant cell walls. Therefore, a detailed analysis of non-cellulosic cell wall glycan through glycome profiling and glycosyl residue composition analysis is important for the polymeric composition and commercial processing of D. papyracea paper. ELISA-based glycome profiling results demonstrated that major classes of cell wall glycans such as xylan, arabinogalactans, and rhamnogalacturonan-I were present on D. papyracea paper. The presence of these polymers in the Himalayan Buddhist handmade paper of Arunachal Pradesh is correlated with its high tensile strength. The results of this study imply that non-cellulosic cell wall glycans are required for the production of high-quality paper. To summarize, immediate action is required to strengthen the centuries-old practice of handmade paper, which can be achieved through education, workshops, technical know-how, and effective marketing aid to entrepreneurs.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria.

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

Immunoinformatics mapping of potential epitopes in SARS-CoV-2 structural proteins.

All approved coronavirus disease 2019 (COVID-19) vaccines in current use are safe, effective, and reduce the risk of severe illness. Although data on the immunological presentation of patients with COVID-19 is limited, increasing experimental evidence supports the significant contribution of B and T cells towards the resolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Despite the availability of several COVID-19 vaccines with high efficacy, more effective vaccines are still needed to protect against the new variants of SARS-CoV-2. Employing a comprehensive immunoinformatic prediction algorithm and leveraging the genetic closeness with SARS-CoV, we have predicted potential immune epitopes in the structural proteins of SARS-CoV-2. The S and N proteins of SARS-CoV-2 and SARS-CoVs are main targets of antibody detection and have motivated us to design four multi-epitope vaccines which were based on our predicted B- and T-cell epitopes of SARS-CoV-2 structural proteins. The cardinal epitopes selected for the vaccine constructs are predicted to possess antigenic, non-allergenic, and cytokine-inducing properties. Additionally, some of the predicted epitopes have been experimentally validated in published papers. Furthermore, we used the C-ImmSim server to predict effective immune responses induced by the epitope-based vaccines. Taken together, the immune epitopes predicted in this study provide a platform for future experimental validations which may facilitate the development of effective vaccine candidates and epitope-based serological diagnostic assays.

Structural insight into locked nucleic acid based novel antisense modifications: A DFT calculations at monomer and MD simulations at oligomer level.

In the present study, five novel LNA based antisense modifications have been proposed. A conformational search was carried out using TANGO, followed by geometry optimization using MOPAC. Based on their electronic energies the most stable conformation for each modification was identified. Further, DFT based full geometry optimization on the most stable conformations at the gas phase B3LYP/6-31G(d,p) using a Gaussian03 and single point energy calculations on the optimized structures at the solvent phase B3LYP/6-311G(d,p) level of theory were done to derive their quantum chemical descriptors using the Gaussian09. A comparison of global reactivity descriptors confirmed that the LNA based modifications were the most reactive. Base-pair stability was recorded by observing the binding energies and base-pairing conformations of modified GC base pairs at the B3LYP/6-311G(d,p) level of theory. Molecular dynamics simulations have been performed at the oligomer duplex level by incorporating individual modifications on 20-mer RNA-RNA duplexes using AMBER16. Free energy calculations of duplex structures suggested that incorporation of A2 modification into the RNA-RNA duplex increased the duplex binding affinity similar to LNA. Whereas, the A3 modification showed less binding compared to LNA but improved binding compared to MOE. This computational approach using quantum chemical methods may be very useful to propose better modifications than the existing ones before performing the experiments in the area of antisense technology.

The Entner-Doudoroff and Nonoxidative Pentose Phosphate Pathways Bypass Glycolysis and the Oxidative Pentose Phosphate Pathway in Ralstonia solanacearum.

In Ralstonia solanacearum, a devastating phytopathogen whose metabolism is poorly understood, we observed that the Entner-Doudoroff (ED) pathway and nonoxidative pentose phosphate pathway (non-OxPPP) bypass glycolysis and OxPPP under glucose oxidation. Evidence derived from 13C stable isotope feeding and genome annotation-based comparative metabolic network analysis supported the observations. Comparative metabolic network analysis derived from the currently available 53 annotated R. solanacearum strains, including a recently reported strain (F1C1), representing the four phylotypes, confirmed the lack of key genes coding for phosphofructokinase (pfk-1) and phosphogluconate dehydrogenase (gnd) enzymes that are relevant for glycolysis and OxPPP, respectively. R. solanacearum F1C1 cells fed with [13C]glucose (99% [1-13C]glucose or 99% [1,2-13C]glucose or 40% [13C6]glucose) followed by gas chromatography-mass spectrometry (GC-MS)-based labeling analysis of fragments from amino acids, glycerol, and ribose provided clear evidence that rather than glycolysis and the OxPPP, the ED pathway and non-OxPPP are the main routes sustaining metabolism in R. solanacearum The 13C incorporation in the mass ions of alanine (m/z 260 and m/z 232), valine (m/z 288 and m/z 260), glycine (m/z 218), serine (m/z 390 and m/z 362), histidine (m/z 440 and m/z 412), tyrosine (m/z 466 and m/z 438), phenylalanine (m/z 336 and m/z 308), glycerol (m/z 377), and ribose (m/z 160) mapped the pathways supporting the observations. The outcomes help better define the central carbon metabolic network of R. solanacearum that can be integrated with 13C metabolic flux analysis as well as flux balance analysis studies for defining the metabolic phenotypes.IMPORTANCE Understanding the metabolic versatility of Ralstonia solanacearum is important, as it regulates the trade-off between virulence and metabolism (1, 2) in a wide range of plant hosts. Due to a lack of clear evidence until this work, several published research papers reported on the potential roles of glycolysis and the oxidative pentose phosphate pathway (OxPPP) in R. solanacearum (3, 4). This work provided evidence from 13C stable isotope feeding and genome annotation-based comparative metabolic network analysis that the Entner-Doudoroff pathway and non-OxPPP bypass glycolysis and OxPPP during the oxidation of glucose, a component of the host xylem pool that serves as a potential carbon source (5). The outcomes help better define the central carbon metabolic network of R. solanacearum that can be integrated with 13C metabolic flux analysis as well as flux balance analysis studies for defining the metabolic phenotypes. The study highlights the need to critically examine phytopathogens whose metabolism is poorly understood.

Essential role of the ESX-3 associated eccD3 locus in maintaining the cell wall integrity of Mycobacterium smegmatis.

Mycobacterial pathogens have evolved a unique secretory apparatus called the Type VII secretion system (T7SS) which comprises of five gene clusters designated as ESX1, ESX2, ESX3, ESX4, and ESX5. Of these the ESX3 T7SS plays an important role in the regulatory uptake of iron from the environment, thereby enabling the bacteria to establish successful infection in the host. However, ESX3 secretion system is conserved among all the mycobacterial species including the fast-growing nonpathogenic species M. smegmatis. Although the function of ESX3 T7SS is known to be absolutely critical for establishing infection by M. tuberculosis, its conserved nature in all the pathogenic and nonpathogenic mycobacterial species intrigues to explore the additional functional roles in Mycobacterium species through which potent targets for drugs can be identified and developed. In the present study, we investigated the possible role of EccD3, a transmembrane protein of the ESX3 T7SS in M. smegmatis by deleting the entire eccD3 gene by efficient allelic exchange method. The preliminary investigations through the creation of knockout mutant of the eccD3 gene indicate that this secretory apparatus has an important role in maintaining the cell wall integrity which was evident from the abnormal colony morphology, lack of biofilm formation and difference in cell wall permeability.

An Innovative Approach to Study Ralstonia solanacearum Pathogenicity in 6 to 7 Days Old Tomato Seedlings by Root Dip Inoculation.

Ralstonia solanacearum (F1C1) is a Gram-negative plant pathogenic bacterium that causes lethal wilt disease in a wide range of plant species. This pathogen is very well known for its unpredictable behavior during infection and wilting its host. Because of its mysterious infection behavior, virulence and pathogenicity standardization are still a big challenge in the case of R. solanacearum. Here, we report an innovative pathogenicity assay of R. solanacearum (F1C1) in the early stage of tomato seedlings by root dip inoculation. In this assay, we employed 6-7days old tomato seedlings for infection grown under nutrients free and gnotobiotic condition. After that, pathogenicity assay was performed by maintaining the inoculated seedlings in 1.5 or 2 ml sterile microfuge tubes. During infection, wilting symptom starts appearing from ~48 h post inoculation and the pathogenicity assay gets completed within seven days of post inoculation. This method is rapid, consistent as well as less resource dependent in terms of labor, space and cost to screen large numbers of plants. Hence, this newly developed assay is an easy and useful approach to study pathogen virulence functions and its interaction with the host plant during wilting and disease progression at the seedling stage.

Codon degeneracy and amino acid abundance influence the measures of codon usage bias: improved Nc (N̂c ) and ENCprime (N̂'c ) measures.

Effective number of codons (N^c) and its variant N^'c (effective number of codons prime) are the two widely used methods for measuring unequal usage of synonymous codons in coding sequences, known as the codon usage bias (CUB). The mathematical formula used in calculating N^c and N^'c values is giving inappropriate measures of CUB in case of low abundance of amino acids. In addition, the magnitude of error also varies according to codon degeneracy. In this study, a modified formula for N^c and N^'c has been developed to measure the CUB more accurately. Online implementations of the modified formula are available in the web portal at http://agnigarh.tezu.ernet.in/~ssankar/cub.php.

Disruption of comA homolog in Ralstonia solanacearum does not impair its twitching motility.

Ralstonia solanacearum is an important phyto-pathogenic bacterium. The bacterium exhibits type IV pili meditated twitching motility that has been implicated in the process of natural transformation in it. A comA gene homolog, alike in several other naturally competent bacteria, has been already reported in this bacterium. However, there are no report of direct link between comA and twitching motility during the natural transformation process in this pathogen. In order to figure out any connection between comA and twitching motility, we created an insertion mutation in comA gene homolog of R. solanacearum F1C1 strain. As anticipated, the insertion mutant (CBRS01 strain) was inefficient for natural transformation. CBRS01 strain was found to be proficient for twitching motility alike the wild-type F1C1. This is interesting since recent findings of Salzer et al. (2016;Environ Microbiol;18:65-74) showed deficiency of twitching motility due to comEC gene (comA homolog) mutation in another naturally competent Gram-negative bacterium Thermus thermophilus. Additionally, we also found CBRS01 strain to be proficient for extracellular cellulase activity and virulence on tomato seedlings. Our findings in this work indicate that an R. solanacearum strain inefficient in undergoing natural transformation can, however, be proficient in exhibiting twitching motility.

Comparative analysis of codon usage bias in Crenarchaea and Euryarchaea genome reveals differential preference of synonymous codons to encode highly expressed ribosomal and RNA polymerase proteins.

The present study was undertaken to investigate the pattern of optimal codon usage in Archaea. Comparative analysis was executed to understand the pattern of codon usage bias between the high expression genes (HEG) and the whole genomes in two Archaeal phyla, Crenarchaea and Euryarchaea. The G+C% of the HEG was found to be less in comparison to the genome G+C% in Crenarchaea, whereas reverse was the case in Euryarchaea. The preponderance of U/A ending codons that code for HEG in Crenarchaea was in sharp contrast to the C/G ended ones in Euryarchaea. The analysis revealed prevalence of Uending codons even within theWWY(nucleotide ambiguity code) families in Crenarchaea vis-à-vis Euryarchaea, bacteria and Eukarya. No plausible interpretation of the observed disparity could be made either in the context of tRNA gene composition or genome G+C%. The results in this study attested that the preferential biasness for codons in HEG of Crenarchaea might be different from Euryarchaea. The main highlights are (i) varied CUB in the HEG and in the whole genomes in Euryarchaea and Crenarchaea. (ii) Crenarchaea was found to have some unusual optimal codons (OCs) compared to other organisms. (iii) G+C% (and GC3) of the HEG were different from the genome G+C% in the two phyla. (iv) Genome G+C% and tRNA gene number failed to explain CUB in Crenarchaea. (v) Translational selection is possibly responsible for A+T rich OCs in Crenarchaea.

Discrepancy among the synonymous codons with respect to their selection as optimal codon in bacteria.

The different triplets encoding the same amino acid, termed as synonymous codons, are not equally abundant in a genome. Factors such as G + C% and tRNA are known to influence their abundance in a genome. However, the order of the nucleotide in each codon per se might also be another factor impacting on its abundance values. Of the synonymous codons for specific amino acids, some are preferentially used in the high expression genes that are referred to as the 'optimal codons' (OCs). In this study, we compared OCs of the 18 amino acids in 221 species of bacteria. It is observed that there is amino acid specific influence for the selection of OCs. There is also influence of phylogeny in the choice of OCs for some amino acids such as Glu, Gln, Lys and Leu. The phenomenon of codon bias is also supported by the comparative studies of the abundance values of the synonymous codons with same G + C. It is likely that the order of the nucleotides in the triplet codon is also perhaps involved in the phenomenon of codon usage bias in organisms.

Selection on GGU and CGU codons in the high expression genes in bacteria.

The fourfold degenerate site (FDS) in coding sequences is important for studying the effect of any selection pressure on codon usage bias (CUB) because nucleotide substitution per se is not under any such pressure at the site due to the unaltered amino acid sequence in a protein. We estimated the frequency variation of nucleotides at the FDS across the eight family boxes (FBs) defined as Um(g), the unevenness measure of a gene g. The study was made in 545 species of bacteria. In many bacteria, the Um(g) correlated strongly with Nc'-a measure of the CUB. Analysis of the strongly correlated bacteria revealed that the U-ending codons (GGU, CGU) were preferred to the G-ending codons (GGG, CGG) in Gly and Arg FBs even in the genomes with G+C % higher than 65.0. Further evidence suggested that these codons can be used as a good indicator of selection pressure on CUB in genomes with higher G+C %.

Constraint on di-nucleotides by codon usage bias in bacterial genomes.

It has been reported earlier that the relative di-nucleotide frequency (RDF) in different parts of a genome is similar while the frequency is variable among different genomes. So RDF is termed as genome signature in bacteria. It is not known if the constancy in RDF is governed by genome wide mutational bias or by selection. Here we did comparative analysis of RDF between the inter-genic and the coding sequences in seventeen bacterial genomes, whose gene expression data was available. The constraint on di-nucleotides was found to be higher in the coding sequences than that in the inter-genic regions and the constraint at the 2nd codon position was more than that in the 3rd position within a genome. Further analysis revealed that the constraint on di-nucleotides at the 2nd codon position is greater in the high expression genes (HEG) than that in the whole genomes as well as in the low expression genes (LEG). We analyzed RDF at the 2nd and the 3rd codon positions in simulated coding sequences that were computationally generated by keeping the codon usage bias (CUB) according to genome G+C composition and the sequence of amino acids unaltered. In the simulated coding sequences, the constraint observed was significantly low and no significant difference was observed between the HEG and the LEG in terms of di-nucleotide constraint. This indicated that the greater constraint on di-nucleotides in the HEG was due to the stronger selection on CUB in these genes in comparison to the LEG within a genome. Further, we did comparative analyses of the RDF in the HEG rpoB and rpoC of 199 bacteria, which revealed a common pattern of constraints on di-nucleotides at the 2nd codon position across these bacteria. To validate the role of CUB on di-nucleotide constraint, we analyzed RDF at the 2nd and the 3rd codon positions in simulated rpoB/rpoC sequences. The analysis revealed that selection on CUB is an important attribute for the constraint on di-nucleotides at these positions in bacterial genomes. We believe that this study has come with major findings of the role of CUB on di-nucleotide constraint in bacterial genomes.

gltB/D mutants of Xanthomonas oryzae pv. oryzae are virulence deficient.

Xanthomonas oryzae pv. oryzae (Xoo) causes bacterial blight, a serious disease of rice. Upon clip inoculation of rice leaves, Xoo causes typical V-shaped lesions whose leading edge moves through the mid-veinal region. We have isolated a virulence deficient mutant of Xoo, referred to as BXO808 that causes limited lesions which primarily extend through the side-veinal regions of rice leaves. Functional complementation studies identified a clone, pSR19, from a cosmid genomic library that restored wild-type virulence and lesion phenotype to BXO808. Transposon mutagenesis of the pSR19 clone, marker exchange experiments, and targeted mutagenesis, revealed that the BXO808 phenotype is due to mutation in the gltB/D genes of Xoo, which encode glutamate synthase subunits α and ß, respectively. The gltB/D mutants that were generated in this study also exhibited virulence deficiency, an altered lesion phenotype and growth deficiency on minimal medium with low levels of ammonium as a sole nitrogen source. This is the first report that mutations in the gltB/D genes of Xoo cause virulence deficiency.

Transfer RNA gene numbers may not be completely responsible for the codon usage bias in asparagine, isoleucine, phenylalanine, and tyrosine in the high expression genes in bacteria.

It is generally believed that the effect of translational selection on codon usage bias is related to the number of transfer RNA genes in bacteria, which is more with respect to the high expression genes than the whole genome. Keeping this in the background, we analyzed codon usage bias with respect to asparagine, isoleucine, phenylalanine, and tyrosine amino acids. Analysis was done in seventeen bacteria with the available gene expression data and information about the tRNA gene number. In most of the bacteria, it was observed that codon usage bias and tRNA gene number were not in agreement, which was unexpected. We extended the study further to 199 bacteria, limiting to the codon usage bias in the two highly expressed genes rpoB and rpoC which encode the RNA polymerase subunits ß and ß', respectively. In concordance with the result in the high expression genes, codon usage bias in rpoB and rpoC genes was also found to not be in agreement with tRNA gene number in many of these bacteria. Our study indicates that tRNA gene numbers may not be the sole determining factor for translational selection of codon usage bias in bacterial genomes.