Whole-genome phylogenies of the family Bacillaceae and expansion of the sigma factor gene family in the Bacillus cereus species-group.

BACKGROUND: The Bacillus cereus sensu lato group consists of six species (B. anthracis, B. cereus, B. mycoides, B. pseudomycoides, B. thuringiensis, and B. weihenstephanensis). While classical microbial taxonomy proposed these organisms as distinct species, newer molecular phylogenies and comparative genome sequencing suggests that these organisms should be classified as a single species (thus, we will refer to these organisms collectively as the Bc species-group). How do we account for the underlying similarity of these phenotypically diverse microbes? It has been established for some time that the most rapidly evolving and evolutionarily flexible portions of the bacterial genome are regulatory sequences and transcriptional networks. Other studies have suggested that the sigma factor gene family of these organisms has diverged and expanded significantly relative to their ancestors; sigma factors are those portions of the bacterial transcriptional apparatus that control RNA polymerase recognition for promoter selection. Thus, examining sigma factor divergence in these organisms would concurrently examine both regulatory sequences and transcriptional networks important for divergence. We began this examination by comparison to the sigma factor gene set of B. subtilis. RESULTS: Phylogenetic analysis of the Bc species-group utilizing 157 single-copy genes of the family Bacillaceae suggests that several taxonomic revisions of the genus Bacillus should be considered. Within the Bc species-group there is little indication that the currently recognized species form related sub-groupings, suggesting that they are members of the same species. The sigma factor gene family encoded by the Bc species-group appears to be the result of a dynamic gene-duplication and gene-loss process that in previous analyses underestimated the true heterogeneity of the sigma factor content in the Bc species-group. CONCLUSIONS: Expansion of the sigma factor gene family appears to have preferentially occurred within the extracytoplasmic function (ECF) sigma factor genes, while the primary alternative (PA) sigma factor genes are, in general, highly conserved with those found in B. subtilis. Divergence of the sigma-controlled transcriptional regulons among various members of the Bc species-group likely has a major role in explaining the diversity of phenotypic characteristics seen in members of the Bc species-group.

Accelerated evolution of the electron transport chain in anthropoid primates.

Mitochondria are both the power plant of the cell and a central integrator of signals that govern the lifespan, replication and death of the cell. Perhaps as a consequence, genes that encode components of the mitochondrial electron transport chain (ETC) are generally conserved. Therefore, it is surprising that many of these genes in anthropoid primates (New World monkeys, Old World monkeys and apes, including humans) have been major targets of darwinian positive selection. Sequence comparisons have provided evidence that marked increases of non-synonymous substitution rates occurred in anthropoid ETC genes that encode subunits of Complex III and IV, and the electron carrier molecule cytochrome c (CYC). Two important questions are: (i) how has evolution altered ETC function? and; (ii) how might functional changes in the ETC be linked to evolution of an expanded neocortical brain?

A third isoform of cytochrome c oxidase subunit VIII is present in mammals.

The terminal enzyme of the mitochondrial respiratory chain, cytochrome c oxidase (COX), contains three mitochondrial and ten nuclear encoded subunits in mammals. Three of the nuclear subunits (VIa, VIIa, and VIII) have muscle and non-muscle-specific isoforms, subunit IV contains a lung-specific isoform, and subunit VIb contains a testes-specific isoform. For subunit VIII, the smallest nuclear encoded COX polypeptide, we have now found a third gene (COX 8-3), which has been identified in human, lemur, rat, and mouse, suggesting that it is present in a broad range of Eutherian mammals. Sequence similarity and gene structure support the homology of COX8-3 to the other subunit VIII isoforms, indicating that all three are the product of gene duplications. COX VIII-3 protein is mitochondrially-targeted, as shown by a fluorescent COX VIII3/DsRed fusion protein. Both the mitochondrial targeting and its sequence conservation suggest that COXVIII-3 functions as part of the COX holoenzyme and could have a tissue-specific role, as is the case for the other two isoforms. Questions remain about where COX8-3 is predominantly expressed. However, detection of full-length cDNAs, lower levels of sequence divergence at the first and second codon positions compared to the third, and a conserved gene structure indicate that COX VIII-3 is an expressed gene whose origin dates to at least 91 million years ago.

Amino acid replacement is rapid in primates for the mature polypeptides of COX subunits, but not for their targeting presequences.

We examined inferred amino acid replacements for 16 genes that encode the proteins of the cytochrome c oxidase (COX) holoenzyme in eight vertebrate species. Phylogeny-based analysis revealed that the human lineage (primates) has had an unusually large, statistically significant, number of amino acid replacements in the mature protein coding region of these genes. This finding is similar to earlier observations of an accelerated non-synonymous substitution rate for some lineages of primates for COX1, COX2, COX4, and COX7AH. In contrast, the mitochondrial targeting presequences of these same proteins have not undergone a concomitant rate change. This more comprehensive analysis suggests that COX5A, COX6B, COX6C, COX7C, and COX8L have also undergone an acceleration in amino acid replacement rates in anthropoid primates. Some of these rate accelerations (e.g. in COX5A and COX7C) are so pronounced that non-human mammalian sequences are more similar to sequences from Xenopus or zebrafish than they are to human. Since the functions of the targeting and mature proteins of these polypeptides are different, the mature portions of these genes are likely to have undergone a functionally significant change that is adaptive in nature.

Sequence variation in the herpes simplex virus U(S)1 ocular virulence determinant.

PURPOSE: The herpes simplex virus type 1 (HSV-1) U(S)1 gene encodes host-range and ocular virulence determinants. Mutations in U(S)1 affecting virulence are known in strain OD4, but the genomic variation across several strains is not known. The goal was to determine the degree of sequence variation in the gene from several ocular HSV isolates. METHODS: The U(S)1 gene from six ocular HSV-1 isolates, as well as strains KOS and F, were sequenced, and bioinformatics analyses were applied to the data. RESULTS: Strains 17, F, CJ394, and CJ311 had identical amino acid sequences. With the other strains, most of the variability was concentrated in the amino-terminal third of the protein. MEME analysis identified a 63-residue core sequence (motif 1) present in all α-herpesvirus U(S)1 homologs that were located in a region identified as structured. Ten amino acids were absolutely conserved in all the α-herpesvirus U(S)1 homologs and were all located in the central core. Consensus-binding motifs for cyclin-dependent kinases and pocket proteins were also identified. CONCLUSIONS: These results suggest that significant sequence variation exists in the U(S)1 gene, that the α22 protein contains a conserved central core region with structurally variable regions at the amino- and carboxyl termini, that 10 amino acids are conserved in α-herpes U(S)1 homologs, and that additional host proteins may interact with the HSV-1 U(S)1 and U(S)1.5 proteins. This information will be valuable in designing further studies on structure-function relationships and on the role these play in host-range determination and keratitis.

Rapid electrostatic evolution at the binding site for cytochrome c on cytochrome c oxidase in anthropoid primates.

Cytochrome c (CYC) oxidase (COX), a multisubunit enzyme that functions in mitochondrial aerobic energy production, catalyzes the transfer of electrons from CYC to oxygen and participates in creating the electrochemical gradient used for ATP synthesis. Modeling three-dimensional structural data on COX and CYC reveals that 57 of the >1,500 COX residues can be implicated in binding CYC. Because of the functional importance of the transfer of electrons to oxygen, it might be expected that natural selection would drastically constrain amino acid replacement rates of CYC and COX. Instead, in anthropoid primates, although not in other mammals, CYC and COX show markedly accelerated amino acid replacement rates, with the COX acceleration being much greater at the positions that bind CYC than at those that do not. Specifically, in the anthropoid lineage descending from the last common ancestor of haplorhines (tarsiers and anthropoids) to that of anthropoids (New World monkeys and catarrhines) and that of catarrhines (Old World monkeys and apes, including humans), a minimum of 27 of the 57 COX amino acid residues that bind CYC were replaced, most frequently from electrostatically charged to noncharged residues. Of the COX charge-bearing residues involved in binding CYC, half (11 of 22) have been replaced with uncharged residues. CYC residues that interact with COX residues also frequently changed, but only two of the CYC changes altered charge. We suggest that reducing the electrostatic interaction between COX and CYC was part of the adaptive evolution underlying the emergence of anthropoid primates.

Rapid nonsynonymous evolution of the iron-sulfur protein in anthropoid primates.

Cytochrome c (CYC) and 9 of the 13 subunits of cytochrome c oxidase (complex IV; COX) were previously shown to have accelerated rates of nonsynonymous substitution in anthropoid primates. Cytochrome b, the mtDNA encoded subunit of ubiquinol-cytochrome c reductase (complex III), also showed an accelerated nonsynonymous substitution rate in anthropoid primates but rate information about the nuclear encoded subunits of complex III has been lacking. We now report that phylogenetic and relative rates analysis of a nuclear encoded catalytically active subunit of complex III, the iron-sulfur protein (ISP), shows an accelerated rate of amino acid replacement similar to cytochrome b. Because both ISP and subunit 9, whose function is not directly related to electron transport, are produced by cleavage into two subunits of the initial translation product of a single gene, it is probable that these two subunits of complex III have essentially identical underlying rates of mutation. Nevertheless, we find that the catalytically active ISP has an accelerated rate of amino acid replacement in anthropoid primates whereas the catalytically inactive subunit 9 does not.

Retention of a duplicate gene through changes in subcellular targeting: an electron transport protein homologue localizes to the golgi.

Cytochrome c oxidase (COX), the terminal enzyme complex of the electron transport chain, contains 13 subunits, 3 encoded by mitochondrial DNA and 10 by nuclear. Several of the nuclear subunits, including subunit VIIa, are known to have two tissue- and development-specific isoforms in mammals. A recently identified third member of the gene family, COX7AR, encodes a protein previously thought to function in mitochondria. However, observation of fluorescent pCOX7AR C-terminal fusion proteins in HeLa cells showed that pCOX7AR is localized to the Golgi apparatus. Sequence analyses indicate that the duplication of COX7AR occurred prior to the origin of the Euteleostomi (bony vertebrates) and that pCOX7AR is more highly conserved than the two other isoforms. These results indicate that, after gene duplication and modification of the mitochondrial targeting signal, pCOX7AR was evolutionarily altered to a new and apparently important function in the Golgi. These results also suggest that predictions of function from homology can be misleading and show that specialization and modification of subcellular localization are similar to cis-element subfunctionalization. In cis-element subfunctionalization, complementary null mutations occur to the cis-elements of the descendents of a gene duplication, causing both descendent genes to be obligate. In the process described in this paper, which could be termed subcellular subfunctionalization, complementary null mutations can occur to the subcellular localization signals of the descendants of a gene duplication, causing both descendent genes to be similarly obligate. Noncomplementary null mutations could also uncover an alternate localization, which is the more likely case for pCOX7AR.

Adaptive evolution of cytochrome c oxidase subunit VIII in anthropoid primates.

Cytochrome c oxidase (COX) is a 13-subunit protein complex that catalyzes the last step in mitochondrial electron transfer in mammals. Of the 10 subunits encoded by nuclear DNA (three are mtDNA products), some are expressed as tissue- and/or development-specific isoforms. For COX subunit VIII, previous work showed that expression of the contractile muscle-specific isoform gene, COX8H, is absent in humans and Old World monkeys, and the other isoform gene, COX8L, is expressed ubiquitously. Here, we show that COX8H is transcribed in most primate clades, but its expression is absent in catarrhines, that is, in Old World monkeys and hominids (apes, including humans), having become a pseudogene in the stem of the catarrhines. The ubiquitously expressed isoform, COX8L, underwent nonsynonymous rate acceleration and elevation in the ratio of nonsynonymous/synonymous changes in the stem of anthropoid primates (New World monkeys and catarrhines), possibly setting the stage for loss of the heart-type (H) isoform. The most rapidly evolving region of VIII-L is one that interacts with COX I, suggesting that the changes are functionally coadaptive. Because accelerated rates of nonsynonymous substitutions in anthropoids such as observed for COX8L are also shown by genes for at least 13 other electron transport chain components, these encoded amino acid replacements may be viewed as part of a series of coadaptive changes that optimized the anthropoid biochemical machinery for aerobic energy metabolism. We argue that these changes were linked to the evolution of an expanded neocortex in anthropoid primates.