Evidence of positive Darwinian selection in Omp85, a highly conserved bacterial outer membrane protein essential for cell viability

The role of selection in shaping diversity of natural M. tuberculosis populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to -9.5 × 10(-4) (95% CI -1.1 × 10(-3) to -6.8 × 10(-4)); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620-1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Quick diagnosis of human brain meningitis using omp85 gene amplicon as a genetic marker

The usual diagnosis of life-threatening human brain bacterial meningitis are expensive, time consuming or non-confirmatory. A quick PCR based diagnosis of meningitis in cerebrospinal fluids (CSF) using specific primers of virulent Omp85 gene of Neisseria meningitidis can detect as low as 1.0 ng of genomic DNA (G-DNA) in 80 min for confirmation of bacterial meningitis caused by N. meningitidis infection. The 257 bp amplicon of Omp85 gene does not show homology with other suspected pathogens in CSF and can be used as a specific genetic marker for diagnosis of the disease.

Differentiation of Xylella fastidiosa strains via multilocus sequence analysis of environmentally mediated genes (MLSA-E)

Isolates of the plant pathogen Xylella fastidiosa are genetically very similar, but studies on their biological traits have indicated differences in virulence and infection symptomatology. Taxonomic analyses have identified several subspecies, and phylogenetic analyses of housekeeping genes have shown broad host-based genetic differences; however, results are still inconclusive for genetic differentiation of isolates within subspecies. This study employs multilocus sequence analysis of environmentally mediated genes (MLSA-E; genes influenced by environmental factors) to investigate X. fastidiosa relationships and differentiate isolates with low genetic variability. Potential environmentally mediated genes, including host colonization and survival genes related to infection establishment, were identified a priori. The ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions (dN/dS) was calculated to select genes that may be under increased positive selection compared to previously studied housekeeping genes. Nine genes were sequenced from 54 X. fastidiosa isolates infecting different host plants across the United States. Results of maximum likelihood (ML) and Bayesian phylogenetic (BP) analyses are in agreement with known X. fastidiosa subspecies clades but show novel within-subspecies differentiation, including geographic differentiation, and provide additional information regarding host-based isolate variation and specificity. dN/dS ratios of environmentally mediated genes, though <1 due to high sequence similarity, are significantly greater than housekeeping gene dN/dS ratios and correlate with increased sequence variability. MLSA-E can more precisely resolve relationships between closely related bacterial strains with low genetic variability, such as X. fastidiosa isolates. Discovering the genetic relationships between X. fastidiosa isolates will provide new insights into the epidemiology of populations of X. fastidiosa, allowing improved disease management in economically important crops.

Crystallographic analysis of the C-terminal domain of the Escherichia coli lipoprotein BamC

In Gram-negative bacteria, the BAM complex catalyzes the essential process of assembling outer membrane proteins. The BAM complex in Escherichia coli consists of five proteins: one β-barrel membrane protein, BamA, and four lipoproteins, BamB, BamC, BamD and BamE. Here, the crystal structure of the C-terminal domain of E. coli BamC (BamC(C): Ala224-Ser343) refined to 1.5 Å resolution in space group H3 is reported. BamC(C) consists of a six-stranded antiparallel β-sheet, three α-helices and one 3(10)-helix. Sequence and surface analysis reveals that most of the conserved residues within BamC(C) are localized to form a continuous negatively charged groove that is involved in a major crystalline lattice contact in which a helix from a neighbouring BamC(C) binds against this surface. This interaction is topologically and architecturally similar to those seen in the substrate-binding grooves of other proteins with BamC-like folds. Taken together, these results suggest that an identified surface on the C-terminal domain of BamC may serve as an important protein-binding surface for interaction with other BAM-complex components or substrates.

Molecular evolution of candidate sour taste receptor gene PKD1L3 in mammals

The PKD1L3 gene encodes an ion channel protein that can interact with the PKD2L1 protein to form a candidate sour taste receptor. In the present study, we have analyzed the evolutionary patterns of PKD1L3 genes from 10 mammalian species. The results showed that PKD1L3 genes have evolved under a dominant purifying selection force. However, for some branches and sites, PKD1L3 genes were detected to have been operated by positive selection. Moreover, some of these positive evolutionary sites are likely to participate in acid stimulus recognition. In rodents, PKD1L3 genes evolved more rapidly than other mammalian lineages. Combined with other functional research reports, our results suggest that rodents may not be the most appropriate model for functional research on the PKD1L3 gene.

Crystal structure of Escherichia coli BamB, a lipoprotein component of the β-barrel assembly machinery complex

In Gram-negative bacteria, the BAM (β-barrel assembly machinery) complex catalyzes the essential process of assembling outer membrane proteins. The BAM complex in Escherichia coli consists of five proteins: one β-barrel membrane protein, BamA, and four lipoproteins, BamB, BamC, BamD, and BamE. Despite their role in outer membrane protein biogenesis, there is currently a lack of functional and structural information on the lipoprotein components of the BAM complex. Here, we report the first crystal structure of BamB, the largest and most functionally characterized lipoprotein component of the BAM complex. The crystal structure shows that BamB has an eight-bladed β-propeller structure, with four β-strands making up each blade. Mapping onto the structure the residues previously shown to be important for BamA interaction reveals that these residues, despite being far apart in the amino acid sequence, are localized to form a continuous solvent-exposed surface on one side of the β-propeller. Found on the same side of the β-propeller is a cluster of residues conserved among BamB homologs. Interestingly, our structural comparison study suggests that other proteins with a BamB-like fold often participate in protein or ligand binding, and that the binding interface on these proteins is located on the surface that is topologically equivalent to where the conserved residues and the residues that are important for BamA interaction are found on BamB. Our structural and bioinformatic analyses, together with previous biochemical data, provide clues to where the BamA and possibly a substrate interaction interface may be located on BamB.

Additive and synergistic bactericidal activity of antibodies directed against minor outer membrane proteins of Neisseria meningitidis

Neisseria meningitidis serogroup B is a major cause of bacterial meningitis in younger populations. The available vaccines are based on outer membrane vesicles obtained from wild-type strains. In children less than 2 years old they confer protection only against strains expressing homologous PorA, a major, variable outer membrane protein (OMP). We genetically modified a strain in order to eliminate PorA and to overproduce one or several minor and conserved OMPs. Using a mouse model mimicking children's PorA-specific bactericidal activity, it was demonstrated that overproduction of more than one minor OMP is required to elicit antibodies able to induce complement-mediated killing of strains expressing heterologous PorA. It is concluded that a critical density of bactericidal antibodies needs to be reached at the surface of meningococci to induce complement-mediated killing. With minor OMPs, this threshold is reached when more than one antigen is targeted, and this allows cross-protection.

Detecting the form of selection in the outer membrane protein C of Enterobacter aerogenes strains and Salmonella species

The types of selective pressure operating on the outer membrane protein C (ompC) of Enterobacter aerogenes strains, the causative agent for nosocomial infections, and Salmonella sp., the hazardous pathogen are investigated using the maximum likelihood-based codon substitution models. Although the rate of amino acid replacement to the silent substitution (omega) across the entire codon sites of ompC of E. aerogenes (omega=0.3194) and Salmonella sp. (omega=0.2047) indicate that the gene is subjected to purifying selection (i.e. omega<1), approximately 3.7% of ompC codon sites in E. aerogenes (omega=21.52) are under the influence of positive Darwinian selection (i.e. omega>1). Such contrast in the intensity of selective pressures in both pathogens could be associated with the differential response to the adverse environmental changes. In E. aerogenes, majority of the positively selected sites are located in the hypervariable cell-surface-exposed domains whereas the trans-membrane domains are functionally highly constrained.

Positively selected codons in immune-exposed loops of the vaccine candidate OMP-P1 of Haemophilus influenzae

The high levels of variation in surface epitopes can be considered as an evolutionary hallmark of immune selection. New computational tools enable analysis of this variation by identifying codons that exhibit high rates of amino acid changes relative to the synonymous substitution rate. In the outer membrane protein P1 of Haemophilus influenzae, a vaccine candidate for nontypeable strains, we identified four codons with this attribute in domains that did not correspond to known or assumed B- and T-cell epitopes of OMP-P1. These codons flank hypervariable domains and do not appear to be false positives as judged from parsimony and maximum likelihood analyses. Some closely spaced positively selected codons have been previously considered part of a transmembrane domain, which would render this region unsuited for inclusion in a vaccine. Secondary structure analysis, three-dimensional structural database searches, and homology modeling using FadL of E. coli as a structural homologue, however, revealed that all positively selected codons are located in or near extracellular looping domains. The spacing and level of diversity of these positively selected and exposed codons in OMP-P1 suggest that vaccine targets based on these and conserved flanking residues may provide broad coverage in H. influenzae.

Selection on protein-coding genes of natural cyanobacterial populations

We examined the distribution of synonymous and non-synonymous changes in 12 protein-coding genes of natural populations of cyanobacteria to infer changes in gene functionality. By comparing mutation distributions within and across species using the McDonald-Kreitman test, we found data sets to contain evidence for purifying selection (hetR of Trichodesmium, nifH of Cylindrospermopsis raceborskii and rpoC1 of Anabaena lemmermannii) and positive selection (kaiC of Microcoleus chthonoplastes and rbcX of Anabaena and Aphanizomenon sp.). Other genes from the same set of clonal isolates (petB and rbcL in M. chthonoplastes and Anabaena/Aphanizomenon, respectively) did not harbour evidence for either form of selection. The results of branch models of codon evolution agreed fully with the results of the McDonald-Kreitman test in terms of significance and absolute value of the dN/dS estimates. The high frequency of gene-specific mutation patterns and their association with branches that separate closely related cyanobacterial genera suggest that evolutionary tests are suited to uncover gene-specific selective differentiation in cyanobacterial genomes. At the same time, given the lack of information about the history of cyanobacteria, analysis of larger numbers of protein-coding genes of clonal cyanobacterial isolates will produce more detailed pictures of the effects of natural selection.

Functional and phylogenetic properties of the pore-forming beta-barrel transporters of the Omp85 family

beta-Barrel-shaped channels of the Omp85 family are involved in the translocation or assembly of proteins of bacterial, mitochondrial, and plastidic outer membranes. We have compared these proteins to understand the evolutionary development of the translocators. We have demonstrated that the proteins from proteobacteria and mitochondria have a pore diameter that is at least five times smaller than found for the Omp85 in cyanobacteria and plastids. This finding can explain why Omp85 from cyanobacteria (but not the homologous protein from proteobacteria) was remodeled to become the protein translocation pore after endosymbiosis. Further, the pore-forming region of the Omp85 proteins is restricted to the C terminus. Based on a phylogenetic analysis we have shown that the pore-forming domain displays a different evolutionary relationship than the N-terminal domain. In line with this, the affinity of the N-terminal domain to the C-terminal region of the Omp85 from plastids and cyanobacteria differs, even though the N-terminal domain is involved in gating of the pore in both groups. We have further shown that the N-terminal domain of nOmp85 takes part in homo-oligomerization. Thereby, the differences in the phylogeny of the two domains are explained by different functional constraints acting on the regions. The pore-forming domain, however, is further divided into two functional regions, where the distal C terminus itself forms a dimeric pore. Based on functional and phylogenetic analysis, we suggest an evolutionary scenario that explains the origin of the contemporary translocon.

Positive selection in the ComC-ComD system of Streptococcal Species

Competence-stimulating peptide (CSP) and ComD of the streptococcal species are a pheromone and its receptor, respectively, involved in the regulation of competence for natural genetic transformation. We show here that these molecules have undergone positive selection. This study is the first report of positive selection due to competition among bacterial populations.