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.

A common evolutionary origin for mitochondria and hydrogenosomes.

Trichomonads are among the earliest eukaryotes to diverge from the main line of eukaryotic descent. Keeping with their ancient nature, these facultative anaerobic protists lack two "hallmark" organelles found in most eukaryotes: mitochondria and peroxisomes. Trichomonads do, however, contain an unusual organelle involved in carbohydrate metabolism called the hydrogenosome. Like mitochondria, hydrogenosomes are double-membrane bounded organelles that produce ATP using pyruvate as the primary substrate. Hydrogenosomes are, however, markedly different from mitochondria as they lack DNA, cytochromes and the citric acid cycle. Instead, they contain enzymes typically found in anaerobic bacteria and are capable of producing molecular hydrogen. We show here that hydrogenosomes contain heat shock proteins, Hsp70, Hsp60, and Hsp10, with signature sequences that are conserved only in mitochondrial and alpha-Gram-negative purple bacterial Hsps. Biochemical analysis of hydrogenosomal Hsp60 shows that the mature protein isolated from the organelle lacks a short, N-terminal sequence, similar to that observed for most nuclear-encoded mitochondrial matrix proteins. Moreover, phylogenetic analyses of hydrogenosomal Hsp70, Hsp60, and Hsp10 show that these proteins branch within a monophyletic group composed exclusively of mitochondrial homologues. These data establish that mitochondria and hydrogenosomes have a common eubacterial ancestor and imply that the earliest-branching eukaryotes contained the endosymbiont that gave rise to mitochondria in higher eukaryotes.

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.

A common origin for woody Sonchus and five related genera in the Macaronesian islands: molecular evidence for extensive radiation.

Woody Sonchus and five related genera (Babcockia, Taeckholmia, Sventenia, Lactucosonchus, and Prenanthes) of the Macaronesian islands have been regarded as an outstanding example of adaptive radiation in angiosperms. Internal transcribed spacer region of the nuclear rDNA (ITS) sequences were used to demonstrate that, despite the extensive morphological and ecological diversity of the plants, the entire alliance in insular Macaronesia has a common origin. The sequence data place Lactucosonchus as sister group to the remainder of the alliance and also indicate that four related genera are in turn sister groups to subg. Dendrosonchus and Taeckholmia. This implies that the woody members of Sonchus were derived from an ancestor similar to allied genera now present on the Canary Islands. It is also evident that the alliance probably occurred in the Canary Islands during the late Miocene or early Pliocene. A rapid radiation of major lineages in the alliance is consistent with an unresolved polytomy near the base and low ITS sequence divergence. Increase of woodiness is concordant with other insular endemics and refutes the relictural nature of woody Sonchus in the Macaronesian islands.

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.

Studies on the biosynthesis of taxol: the taxane carbon skeleton is not of mevalonoid origin.

A cell culture of Taxus chinensis was established to produce the diterpene 2alpha,5alpha,10beta,14beta-tetra-acetoxy4 ++ +(20),11-taxadiene (taxuyunnanine C) in 2.6% (dry weight) yield. The incorporation of [U-13C6]glucose, [1-13C]glucose, and [1,2-13C2]acetate into this diterpene was analyzed by NMR spectroscopy. Label from [1,2-13C2]acetate was diverted to the four acetyl groups of taxuyunnanine C, but not to the taxane ring system. Label from [1-13C]glucose and [U-13C6]glucose was efficiently incorporated into both the taxane ring system and the acetyl groups. The four isoprenoid moieties of the diterpene showed identical labeling patterns. The analysis of long-range 13C13C couplings in taxuyunnanine C obtained from an experiment with [U-13C6]glucose documents the involvement of an intramolecular rearrangement in the biosynthesis of the isoprenoid precursor. The labeling patterns are inconsistent with the mevalonate pathway. The taxoid data share important features with the alternative pathway of isoprenoid biosynthesis operating in certain eubacteria Rohmer, M., Knani, M., Simonin, P., Sutter, B. & Sahm, H. (1993) Biochem. J. 295, 517-524].

Monophyletic origin and unique dispersal patterns of domestic fowls.

With the aim of elucidating in greater detail the genealogical origin of the present domestic fowls of the world, we have determined mtDNA sequences of the D-loop regions for a total of 21 birds, of which 12 samples belong to red junglefowl (Gallus gallus) comprising three subspecies (six Gallus gallus gallus, three Gallus gallus spadiceus, and three Gallus gallus bankiva) and nine represent diverse domestic breeds (Gallus gallus domesticus). We also sequenced four green junglefowl (Gallus varius), two Lafayette's junglefowl (Gallus lafayettei), and one grey junglefowl (Gallus sonneratii). We then constructed a phylogenetic tree for these birds by the use of nucleotide sequences, choosing the Japanese quail (Coturnix coturnix japonica) as an outgroup. We found that a continental population of G. g. gallus was the real matriarchic origin of all the domestic poultries examined in this study. It is also of particular interest that there were no discernible differences among G. gallus subspecies; G. g. bankiva was a notable exception. This was because G. g. spadiceus and a continental population of G. g. gallus formed a single cluster in the phylogenetic tree. G. g. bankiva, on the other hand, was a distinct entity, thus deserving its subspecies status. It implies that a continental population of G. g. gallus sufficed as the monophyletic ancestor of all domestic breeds. We also discussed a possible significance of the initial dispersal pattern of the present domestic fowls, using the phylogenetic tree.

Ancient single origin for Malagasy primates.

We report new evidence that bears decisively on a long-standing controversy in primate systematics. DNA sequence data for the complete cytochrome b gene, combined with an expanded morphological data set, confirm the results of a previous study and again indicate that all extant Malagasy lemurs originated from a single common ancestor. These results, as well as those from other genetic studies, call for a revision of primate classifications in which the dwarf and mouse lemurs are placed within the Afro-Asian lorisiforms. The phylogenetic results, in agreement with paleocontinental data, indicate an African origin for the common ancestor of lemurs and lorises (the Strepsirrhini). The molecular data further suggest the surprising conclusion that lemurs began evolving independently by the early Eocene at the latest. This indicates that the Malagasy primate lineage is more ancient than generally thought and places the split between the two strepsirrhine lineages well before the appearance of known Eocene fossil primates. We conclude that primate origins were marked by rapid speciation and diversification sometime before the late Paleocene.

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.

Antisense globin RNA in mouse erythroid tissues: structure, origin, and possible function.

The aim of the experiments described in this paper was to test for the presence of antisense globin RNA in mouse erythroid tissues and, if found, to characterize these molecules. The present study made use of a multistep procedure in which a molecular tag is attached to cellular RNA by ligation with a defined ribooligonucleotide. The act of ligation preserves the termini of RNA molecules, which become the junctions between cellular RNAs and the ligated ribooligonucleotide. It also unambiguously preserves the identity of cellular RNA as a sense or antisense molecule through all subsequent manipulations. Using this approach, we identified and characterized antisense beta-globin RNA in erythroid spleen cells and reticulocytes from anemic mice. We show in this paper that the antisense globin RNA is fully complementary to spliced globin mRNA, indicative of the template/transcript relationship. It terminates at the 5' end with a uridylate stretch, reflecting the presence of poly(A) at the 3' end of the sense globin mRNA. With respect to the structure of their 3' termini, antisense globin RNA can be divided into three categories: full-size molecules corresponding precisely to globin mRNA, truncated molecules lacking predominantly 14 3'-terminal nucleotides, and extended antisense RNA containing 17 additional 3'-terminal nucleotides. The full-size antisense globin RNA contains two 14-nt-long complementary sequences within its 3'-terminal segment corresponding to the 5'-untranslated region of globin mRNA. This, together with the nature of the predominant truncation, suggests a mechanism by which antisense RNA might give rise to new sense-strand globin mRNA.

In vivo protein-DNA interactions at human DNA replication origin.

Protein-DNA interactions were studied in vivo at the region containing a human DNA replication origin, located at the 3' end of the lamin B2 gene and partially overlapping the promoter of another gene, located downstream. DNase I treatment of nuclei isolated from both exponentially growing and nonproliferating HL-60 cells showed that this region has an altered, highly accessible, chromatin structure. High-resolution analysis of protein-DNA interactions in a 600-bp area encompassing the origin was carried out by the in vivo footprinting technique based on the ligation-mediated polymerase chain reaction. In growing HL-60 cells, footprints at sequences homologous to binding sites for known transcription factors (members of the basic-helix-loop-helix family, nuclear respiratory factor 1, transcription factor Sp1, and upstream binding factor) were detected in the region corresponding to the promoter of the downstream gene. Upon conversion of cells to a nonproliferative state, a reduction in the intensity of these footprints was observed that paralleled the diminished transcriptional activity of the genomic area. In addition to these protections, in close correspondence to the replication initiation site, a prominent footprint was detected that extended over 70 nucleotides on one strand only. This footprint was absent from nonproliferating HL-60 cells, indicating that this specific protein-DNA interaction might be involved in the process of origin activation.

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.