G. Brian Golding

Antibiotic resistance is ancient

The discovery of antibiotics more than 70 years ago initiated a period of drug innovation and implementation in human and animal health and agriculture. These discoveries were tempered in all cases by the emergence of resistant microbes. This history has been interpreted to mean that antibiotic resistance in pathogenic bacteria is a modern phenomenon; this view is reinforced by the fact that collections of microbes that predate the antibiotic era are highly susceptible to antibiotics. Here we report targeted metagenomic analyses of rigorously authenticated ancient DNA from 30,000-year-old Beringian permafrost sediments and the identification of a highly diverse collection of genes encoding resistance to β-lactam, tetracycline and glycopeptide antibiotics. Structure and function studies on the complete vancomycin resistance element VanA confirmed its similarity to modern variants. These results show conclusively that antibiotic resistance is a natural phenomenon that predates the modern selective pressure of clinical antibiotic use.

The closest BLAST hit is often not the nearest neighbor.

It is well known that basing phylogenetic reconstructions on uncorrected genetic distances can lead to errors in their reconstruction. Nevertheless, it is often common practice to report simply the most similar BLAST (Altschul et al. 1997) hit in genomic reports that discuss many genes (Ruepp et al. 2000; Freiberg et al. 1997). This is because BLAST hits can provide a rapid, efficient, and concise analysis of many genes at once. These hits are often interpreted to imply that the gene is most closely related to the gene or protein in the databases that returned the closest BLAST hit. Though these two may coincide, for many genes, particularly genes with few homologs, they may not be the same. There are a number of circumstances that can account for such limitations in accuracy (Eisen 2000). We stress here that genes appearing to be the most similar based on BLAST hits are often not each others closest relative phylogenetically. The extent to which this occurs depends on the availability of close relatives present in the databases. As an example we have chosen the analysis of the genomes of a crenarcheaota species Aeropyrum pernix, an organism with few close relatives fully sequenced, and Escherichia coli, an organism whose closest relative, Salmonella typhimurium, is completely sequenced.

A draft genome of Yersinia pestis from victims of the Black Death

Technological advances in DNA recovery and sequencing have drastically expanded the scope of genetic analyses of ancient specimens to the extent that full genomic investigations are now feasible and are quickly becoming standard. This trend has important implications for infectious disease research because genomic data from ancient microbes may help to elucidate mechanisms of pathogen evolution and adaptation for emerging and re-emerging infections. Here we report a reconstructed ancient genome of Yersinia pestis at 30-fold average coverage from Black Death victims securely dated to episodes of pestilence-associated mortality in London, England, 1348–1350. Genetic architecture and phylogenetic analysis indicate that the ancient organism is ancestral to most extant strains and sits very close to the ancestral node of all Y. pestis commonly associated with human infection. Temporal estimates suggest that the Black Death of 1347–1351 was the main historical event responsible for the introduction and widespread dissemination of the ancestor to all currently circulating Y. pestis strains pathogenic to humans, and further indicates that contemporary Y. pestis epidemics have their origins in the medieval era. Comparisons against modern genomes reveal no unique derived positions in the medieval organism, indicating that the perceived increased virulence of the disease during the Black Death may not have been due to bacterial phenotype. These findings support the notion that factors other than microbial genetics, such as environment, vector dynamics and host susceptibility, should be at the forefront of epidemiological discussions regarding emerging Y. pestis infections.

The origin of the eukaryotic cell

Molecular sequence data are beginning to provide important insights into the evolutionary origin of eukaryotic cells. Global phylogenies of numerous protein sequences indicate that the eukaryotic cell nucleus is a chimera, which has received major contributions from both a Gram-negative eubacterium and an archaebacterium. Recent studies also indicate that the formation of the nuclear envelope and the endoplasmic reticulum was accompanied by duplication of genes for the molecular chaperone proteins (e.g. hsp70, hsp90), which facilitate protein transport across membranes. Based on these observations, it is suggested that the ancestral eukaryotic cell arose by a unique endosymbiotic event involving engulfment of an eocyte archaebacterium by a Gram-negative eubacterial host.

Evolution of HSP70 Gene and Its Implications Regarding Relationships Between Archaebacteria, Eubacteria, and Eukaryotes

The 70-kDa heat-shock protein (HSP70) constitutes the most conserved protein present in all organisms that is known to date. Based on global alignment of HSP70 sequences from organisms representing all three domains, numerous sequence signatures that are specific for prokaryotic and eukaryotic homologs have been identified. HSP70s from the two archaebacterial species examined (viz., Halobacterium marismortui and Methanosarcina mazei) have been found to contain all eubacterial but no eukaryotic signature sequences. Based on several novel features of the HSP70 family of proteins (viz., presence of tandem repeats of a 9-amino-acid [a.a.] polypeptide sequence and structural similarity between the first and second quadrants of HSP70, homology of the N-terminal half of HSP70 to the bacterial MreB protein, presence of a conserved insert of 23–27 a.a. in all HSP70s except those from archaebacteria and gram-positive eubacteria) a model for the evolution of HSP70 gene from an early stage is proposed. The HSP70 homologs from archaebacteria and gram-positive bacteria lacking the insert in the N-terminal quadrants are indicated to be the ancestral form of the protein. Detailed phylogenetic analyses of HSP70 sequence data (viz., by bootstrap analyses, maximum parsimony, and maximum likelihood methods) provide evidence that archaebacteria are not monophyletic and show a close evolutionary linkage with the gram-positive eubacteria. These results do not support the traditional archaebacterial tree, where a close relationship between archaebacterial and eukaryotic homologs is observed. To explain the phylogenies based on HSP70 and other gene sequences, a model for the origin of eukaryotic cells involving fusion between archaebacteria and gram-negative eubacteria is proposed.

Transcriptional Profiling Implicates Novel Interactions between Abiotic Stress and Hormonal Responses in Thellungiella, a Close Relative of Arabidopsis

Thellungiella, an Arabidopsis (Arabidopsis thaliana)-related halophyte, is an emerging model species for studies designed to elucidate molecular mechanisms of abiotic stress tolerance. Using a cDNA microarray containing 3,628 unique sequences derived from previously described libraries of stress-induced cDNAs of the Yukon ecotype of Thellungiella salsuginea, we obtained transcript profiles of its response to cold, salinity, simulated drought, and rewatering after simulated drought. A total of 154 transcripts were differentially regulated under the conditions studied. Only six of these genes responded to all three stresses of drought, cold, and salinity, indicating a divergence among the end responses triggered by each of these stresses. Unlike in Arabidopsis, there were relatively few transcript changes in response to high salinity in this halophyte. Furthermore, the gene products represented among drought-responsive transcripts in Thellungiella associate a down-regulation of defense-related transcripts with exposure to water deficits. This antagonistic interaction between drought and biotic stress response may demonstrate Thellungiellas ability to respond precisely to environmental stresses, thereby conserving energy and resources and maximizing its survival potential. Intriguingly, changes of transcript abundance in response to cold implicate the involvement of jasmonic acid. While transcripts associated with photosynthetic processes were repressed by cold, physiological responses in plants developed at low temperature suggest a novel mechanism for photosynthetic acclimation. Taken together, our results provide useful starting points for more in-depth analyses of Thellungiellas extreme stress tolerance.

Pervasive Cryptic Epistasis in Molecular Evolution

The functional effects of most amino acid replacements accumulated during molecular evolution are unknown, because most are not observed naturally and the possible combinations are too numerous. We created 168 single mutations in wild-type Escherichia coli isopropymalate dehydrogenase (IMDH) that match the differences found in wild-type Pseudomonas aeruginosa IMDH. 104 mutant enzymes performed similarly to E. coli wild-type IMDH, one was functionally enhanced, and 63 were functionally compromised. The transition from E. coli IMDH, or an ancestral form, to the functional wild-type P. aeruginosa IMDH requires extensive epistasis to ameliorate the combined effects of the deleterious mutations. This result stands in marked contrast with a basic assumption of molecular phylogenetics, that sites in sequences evolve independently of each other. Residues that affect function are scattered haphazardly throughout the IMDH structure. We screened for compensatory mutations at three sites, all of which lie near the active site and all of which are among the least active mutants. No compensatory mutations were found at two sites indicating that a single site may engage in compound epistatic interactions. One complete and three partial compensatory mutations of the third site are remote and lie in a different domain. This demonstrates that epistatic interactions can occur between distant (>20Å) sites. Phylogenetic analysis shows that incompatible mutations were fixed in different lineages.

Evolution of allosteric control in glycogen phosphorylase

In relation to the primary sequence and three-dimensional structure of rabbit muscle glycogen phosphorylase, we have carried out a comparative sequence analysis of phosphorylases from human, rat, Dictyostelium, yeast, potato and Escherichia coli. Based on sequence similarity, a large region of the protein is shared by these enzymes extending from alpha-helix-1 to the last alpha-helix-33. Conserved residues are equally distributed between the N and C-terminal domains and occur primarily in buried residues. Phylogenetic analysis indicates that the two isozymes within either E. coli, potato or Dictyostelium are more closely related to each other than they are to other phosphorylases. Yeast phosphorylase is most closely related to the Dictyostelium isozymes. Mammalian muscle and brain isozymes are more closely related to each other than to the liver isozyme and the muscle isozyme is evolving at the slowest rate. All phosphorylases exhibit high conservation of active site and pyridoxal phosphate binding residues. Most phosphorylases also exhibit high conservation of sugar binding residues in the glycogen storage site. Phosphorylation and AMP binding site residues are poorly conserved in non-mammalian phosphorylases. In contrast, glucose-6-P binding residues are highly conserved in four of the seven non-mammalian enzymes. Analysis of interacting pairs of dimer contact residues indicates that they can be grouped into three relatively independent networks. One network contains phosphorylation and AMP binding residues and is poorly conserved in non-mammalian enzymes. A second network contains glucose-6-P binding residues and is highly conserved in enzymes containing a conserved glucose-6-P binding site. A third, conserved network contains residues within the tower helix and gate loop. A model for the evolution of allostery in phosphorylase is proposed, suggesting that glucose-6-P inhibition was an early control mechanism. The later creation of primarily distinct ligand binding sites for AMP/phosphorylation control may have allowed the establishment of a separate dimer contact network for propagating conformational changes leading to activation rather than inhibition of enzyme activity.

A Step Toward Barcoding Life: A Model-Based, Decision-Theoretic Method to Assign Genes to Preexisting Species Groups

A major part of the barcoding of life problem is assigning newly sequenced or sampled individuals to existing groups that are preidentified externally (by a taxonomist, for example). This problem involves evaluating the statistical evidence towards associating a sequence from a new individual with one group or another. The main concern of our current research is to perform this task in a fast and accurate manner. To accomplish this we have developed a model-based, decision-theoretic framework based on the coalescent theory. Under this framework, we utilized both distance and the posterior probability of a group, given the sequences from members of this group and the sequence from a newly sampled individual to assign this new individual. We believe that this approach makes efficient use of the available information in the data. Our preliminary results indicated that this approach is more accurate than using a simple measure of distance for assignment.

The role of laterally transferred genes in adaptive evolution

BackgroundBacterial genomes develop new mechanisms to tide them over the imposing conditions they encounter during the course of their evolution. Acquisition of new genes by lateral gene transfer may be one of the dominant ways of adaptation in bacterial genome evolution. Lateral gene transfer provides the bacterial genome with a new set of genes that help it to explore and adapt to new ecological niches.MethodsA maximum likelihood analysis was done on the five sequenced corynebacterial genomes to model the rates of gene insertions/deletions at various depths of the phylogeny.ResultsThe study shows that most of the laterally acquired genes are transient and the inferred rates of gene movement are higher on the external branches of the phylogeny and decrease as the phylogenetic depth increases. The newly acquired genes are under relaxed selection and evolve faster than their older counterparts. Analysis of some of the functionally characterised LGTs in each species has indicated that they may have a possible adaptive role.ConclusionThe five Corynebacterial genomes sequenced to date have evolved by acquiring between 8 – 14% of their genomes by LGT and some of these genes may have a role in adaptation.