Evolutionary genomics of Staphylococcus aureus: Insights into the origin of methicillin-resistant strains and the toxic shock syndrome epidemic

An emerging theme in medical microbiology is that extensive variation exists in gene content among strains of many pathogenic bacterial species. However, this topic has not been investigated on a genome scale with strains recovered from patients with well-defined clinical conditions. Staphylococcus aureus is a major human pathogen and also causes economically important infections in cows and sheep. A DNA microarray representing >90% of the S. aureus genome was used to characterize genomic diversity, evolutionary relationships, and virulence gene distribution among 36 strains of divergent clonal lineages, including methicillin-resistant strains and organisms causing toxic shock syndrome. Genetic variation in S. aureus is very extensive, with ≈22% of the genome comprised of dispensable genetic material. Eighteen large regions of difference were identified, and 10 of these regions have genes that encode putative virulence factors or proteins mediating antibiotic resistance. We find that lateral gene transfer has played a fundamental role in the evolution of S. aureus. The mec gene has been horizontally transferred into distinct S. aureus chromosomal backgrounds at least five times, demonstrating that methicillin-resistant strains have evolved multiple independent times, rather than from a single ancestral strain. This finding resolves a long-standing controversy in S. aureus research. The epidemic of toxic shock syndrome that occurred in the 1970s was caused by a change in the host environment, rather than rapid geographic dissemination of a new hypervirulent strain. DNA microarray analysis of large samples of clinically characterized strains provides broad insights into evolution, pathogenesis, and disease emergence.

The Ku-like protein from Saccharomyces cerevisiae is required in vitro for the assembly of a stable multiprotein complex at a eukaryotic origin of replication.

We have previously shown that three distinct DNA-binding activities, in crude form, are necessary for the ATP-dependent assembly of a specific and stable multiprotein complex at a yeast origin of replication. Here we show the purification of one of these DNA binding activities, referred to as origin binding factor 2 (OBF2). The purified protein is a heterodimer composed of two polypeptides with molecular mass values of 65 and 80 kDa as determined by SDS/PAGE. Purified OBF2 not only binds DNA but also supports the formation of a protein complex at essential sequences within the ARS121 origin of replication. Interestingly, OBF2 binds tightly and nonspecifically to both duplex DNA and single-stranded DNA. The interaction with duplex DNA occurs at the termini. N-terminal sequencing of the 65-kDa subunit has revealed that this polypeptide is identical to the previously identified HDF1 peptide, a yeast homolog of the small subunit of the mammalian Ku autoantigen. Although the potential involvement of Ku in DNA metabolic events has been proposed, this is the first requirement for a Ku-like protein in the assembly of a protein complex at essential sequences within a eukaryotic origin of replication.

Class 3 Hox genes in insects and the origin of zen.

We have cloned, from a beetle and a locust, genes that are homologous to the class 3 Hox genes of vertebrates. Outside the homeobox they share sequence motifs with the Drosophila zerknüllt (zen) and z2 genes, and like zen, are expressed only in extraembryonic membranes. We conclude that the zen genes of Drosophila derive from a Hox class 3 sequence that formed part of the common ancestral Hox cluster, but that in insects this (Hox) gene has lost its role in patterning the anterio-posterior axis of the embryo, and acquired a new function. In the lineage leading to Drosophila, the zen genes have diverged particularly rapidly.

Exon shuffling and the origin of the mitochondrial targeting function in plant cytochrome c1 precursor.

Since most of the examples of "exon shuffling" are between vertebrate genes, the view is often expressed that exon shuffling is limited to the evolutionarily recent lineage of vertebrates. Although exon shuffling in plants has been inferred from the analysis of intron phases of plant genes [Long, M., Rosenberg, C. & Gilbert, W. (1995) Proc. Natl. Acad. Sci. USA 92, 12495-12499] and from the comparison of two functionally unknown sunflower genes [Domon, C. & Steinmetz, A. (1994) Mol. Gen. Genet. 244, 312-317], clear cases of exon shuffling in plant genes remain to be uncovered. Here, we report an example of exon shuffling in two important nucleus-encoded plant genes: cytosolic glyceraldehyde-3-phosphate dehydrogenase (cytosolic GAPDH or GapC) and cytochrome c1 precursor. The intron-exon structures of the shuffled region indicate that the shuffling event took place at the DNA sequence level. In this case, we can establish a donor-recipient relationship for the exon shuffling. Three amino terminal exons of GapC have been donated to cytochrome c1, where, in a new protein environment, they serve as a source of the mitochondrial targeting function. This finding throws light upon an old important but unsolved question in gene evolution: the origin of presequences or transit peptides that generally exist in nucleus-encoded organelle genes.

The root of the universal tree and the origin of eukaryotes based on elongation factor phylogeny.

The genes for the protein synthesis elongation factors Tu (EF-Tu) and G (EF-G) are the products of an ancient gene duplication, which appears to predate the divergence of all extant organismal lineages. Thus, it should be possible to root a universal phylogeny based on either protein using the second protein as an outgroup. This approach was originally taken independently with two separate gene duplication pairs, (i) the regulatory and catalytic subunits of the proton ATPases and (ii) the protein synthesis elongation factors EF-Tu and EF-G. Questions about the orthology of the ATPase genes have obscured the former results, and the elongation factor data have been criticized for inadequate taxonomic representation and alignment errors. We have expanded the latter analysis using a broad representation of taxa from all three domains of life. All phylogenetic methods used strongly place the root of the universal tree between two highly distinct groups, the archaeons/eukaryotes and the eubacteria. We also find that a combined data set of EF-Tu and EF-G sequences favors placement of the eukaryotes within the Archaea, as the sister group to the Crenarchaeota. This relationship is supported by bootstrap values of 60-89% with various distance and maximum likelihood methods, while unweighted parsimony gives 58% support for archaeal monophyly.

Dynamic friction and the origin of the complexity of earthquake sources.

We study a simple antiplane fault of finite length embedded in a homogeneous isotropic elastic solid to understand the origin of seismic source heterogeneity in the presence of nonlinear rate- and state-dependent friction. All the mechanical properties of the medium and friction are assumed homogeneous. Friction includes a characteristic length that is longer than the grid size so that our models have a well-defined continuum limit. Starting from a heterogeneous initial stress distribution, we apply a slowly increasing uniform stress load far from the fault and we simulate the seismicity for a few 1000 events. The style of seismicity produced by this model is determined by a control parameter associated with the degree of rate dependence of friction. For classical friction models with rate-independent friction, no complexity appears and seismicity is perfectly periodic. For weakly rate-dependent friction, large ruptures are still periodic, but small seismicity becomes increasingly nonstationary. When friction is highly rate-dependent, seismicity becomes nonperiodic and ruptures of all sizes occur inside the fault. Highly rate-dependent friction destabilizes the healing process producing premature healing of slip and partial stress drop. Partial stress drop produces large variations in the state of stress that in turn produce earthquakes of different sizes. Similar results have been found by other authors using the Burridge and Knopoff model. We conjecture that all models in which static stress drop is only a fraction of the dynamic stress drop produce stress heterogeneity.

Skeletal matrices, muci, and the origin of invertebrate calcification.

The sudden appearance of calcified skeletons among many different invertebrate taxa at the Precambrian-Cambrian transition may have required minor reorganization of preexisting secretory functions. In particular, features of the skeletal organic matrix responsible for regulating crystal growth by inhibition may be derived from mucous epithelial excretions. The latter would have prevented spontaneous calcium carbonate overcrusting of soft tissues exposed to the highly supersaturated Late Proterozoic ocean [Knoll, A. H., Fairchild, I. J. & Swett, K. (1993) Palaios 8, 512-525], a putative function for which we propose the term "anticalcification." We tested this hypothesis by comparing the serological properties of skeletal water-soluble matrices and mucous excretions of three invertebrates--the scleractinian coral Galaxea fascicularis and the bivalve molluscs Mytilus edulis and Mercenaria mercenaria. Crossreactivities recorded between muci and skeletal water-soluble matrices suggest that these different secretory products have a high degree of homology. Furthermore, freshly extracted muci of Mytilus were found to inhibit calcium carbonate precipitation in solution.

On the origin of the order Artiodactyla.

The first known members of the order Artiodactyla appeared suddenly throughout the Holarctic region at the beginning of the Eocene. They are characterized by distinctive cursorial skeletal specializations. Owing to their abrupt appearance and the lack of transitional forms, the origin of the order is problematic. Descent from a "condylarth," specifically the arctocyonid Chriacus, has been suggested based on dental resemblances, but until now postcranial anatomy seemed to preclude close relationship between Arctocyonidae and Artiodactyla. A middle Paleocene specimen of a small arctocyonid (?Chriacus) reported here is much more similar to the oldest artiodactyl, Diacodexis, in the derived condition of the hindlimb, reviving the possibility that Artiodactyla evolved from an arctocyonid.

Archaeal-eubacterial mergers in the origin of Eukarya: phylogenetic classification of life.

A symbiosis-based phylogeny leads to a consistent, useful classification system for all life. "Kingdoms" and "Domains" are replaced by biological names for the most inclusive taxa: Prokarya (bacteria) and Eukarya (symbiosis-derived nucleated organisms). The earliest Eukarya, anaerobic mastigotes, hypothetically originated from permanent whole-cell fusion between members of Archaea (e.g., Thermoplasma-like organisms) and of Eubacteria (e.g., Spirochaeta-like organisms). Molecular biology, life-history, and fossil record evidence support the reunification of bacteria as Prokarya while subdividing Eukarya into uniquely defined subtaxa: Protoctista, Animalia, Fungi, and Plantae.

A single-stranded DNA binding protein binds the origin of replication of the duplex kinetoplast DNA.

Replication of the kinetoplast DNA (kDNA) minicircle of trypanosomatids initiates at a conserved 12-nt sequence, 5'-GGGGTTGGTGTA-3', termed the universal minicircle sequence (UMS). A sequence-specific single-stranded DNA-binding protein from Crithidia fasciculata binds the heavy strand of the 12-mer UMS. Whereas this UMS-binding protein (UMSBP) does not bind a duplex UMS dodecamer, it binds the double-stranded kDNA minicircle as well as a duplex minicircle fragment containing the origin-associated UMS. Binding of the minicircle origin region by the single-stranded DNA binding protein suggested the local unwinding of the DNA double helix at this site. Modification of thymine residues at this site by KMnO4 revealed that the UMS resides within an unwound or otherwise sharply distorted DNA at the minicircle origin region. Computer analysis predicts the sequence-directed curving of the minicircle origin region. Electrophoresis of a minicircle fragment containing the origin region in polyacrylamide gels revealed a significantly lower electrophoretic mobility than expected from its length. The fragment anomalous electrophoretic mobility is displayed only in its native conformation and is dependent on temperature and gel porosity, indicating the local curving of the DNA double helix. We suggest that binding of UMSBP at the minicircle origin of replication is possible through local unwinding of the DNA double helix at the UMS site. It is hypothesized here that this local melting is initiated through the untwisting of unstacked dinucleotide sequences at the bent origin site.

Genetic absolute dating based on microsatellites and the origin of modern humans.

We introduce a new genetic distance for microsatellite loci, incorporating features of the stepwise mutation model, and test its performance on microsatellite polymorphisms in humans, chimpanzees, and gorillas. We find that it performs well in determining the relations among the primates, but less well than other distance measures (not based on the stepwise mutation model) in determining the relations among closely related human populations. However, the deepest split in the human phylogeny seems to be accurately reconstructed by the new distance and separates African and non-African populations. The new distance is independent of population size and therefore allows direct estimation of divergence times if the mutation rate is known. Based on 30 microsatellite polymorphisms and a recently reported average mutation rate of 5.6 x 10(-4) at 15 dinucleotide microsatellites, we estimate that the deepest split in the human phylogeny occurred about 156,000 years ago. Unlike most previous estimates, ours requires no external calibration of the rate of molecular evolution. We can use such calibrations, however, to test our estimate.

Considerations on the folding topology and evolutionary origin of cadherin domains.

Cell-cell adhesion in zonula adherens and desmosomal junctions is mediated by cadherins, and recent crystal structures of the first domain from murine N-cadherin provide a plausible molecular basis for this adhesive action. A structure-based sequence analysis of this adhesive domain indicates that its fold is common to all extracellular cadherin domains. The cadherin folding topology is also shown to be similar to immunoglobulin-like domains and to other Greek-key beta-sandwich structures, as diverse as domains from plant cytochromes, bacterial cellulases, and eukaryotic transcription factors. Sequence similarities between cadherins and these other molecules are very low, however, and intron patterns are also different. On balance, independent origins for a favorable folding topology seem more likely than evolutionary divergence from an ancestor common to cadherins and immunoglobulins.

Evolutionary origin of Plasmodium and other Apicomplexa based on rRNA genes.

We have explored the evolutionary history of the Apicomplexa and two related protistan phyla, Dinozoa and Ciliophora, by comparing the nucleotide sequences of small subunit ribosomal RNA genes. We conclude that the Plasmodium lineage, to which the malarial parasites belong, diverged from other apicomplexan lineages (piroplasmids and coccidians) several hundred million years ago, perhaps even before the Cambrian. The Plasmodium radiation, which gave rise to several species parasitic to humans, occurred approximately 129 million years ago; Plasmodium parasitism of humans has independently arisen several times. The origin of apicomplexans (Plasmodium), dinoflagellates, and ciliates may be > 1 billion years old, perhaps older than the three multicellular kingdoms of animals, plants, and fungi. Digenetic parasitism independently evolved several times in the Apicomplexa.

In vitro cleavage and joining at the viral origin of replication by the replication initiator protein of tomato yellow leaf curl virus.

Replication of the single-stranded DNA genome of geminiviruses occurs via a double-stranded intermediate that is subsequently used as a template for rolling-circle replication of the viral strand. Only one of the proteins encoded by the virus, here referred to as replication initiator protein (Rep protein), is indispensable for replication. We show that the Rep protein of tomato yellow leaf curl virus initiates viral-strand DNA synthesis by introducing a nick in the plus strand within the nonanucleotide 1TAATATT decreases 8AC, identical among all geminiviruses. After cleavage, the Rep protein remains bound to the 5' end of the cleaved strand. In addition, we show that the Rep protein has a joining activity, suggesting that it acts as a terminase, thus resolving the nascent viral single strand into genome-sized units.