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.

Escherichia coli strains that allow antibiotic-free plasmid selection and maintenance by repressor titration

We report the construction of two novel Escherichia coli strains (DH1lacdapD and DH1lacP2dapD) that facilitate the antibiotic-free selection and stable maintenance of recombinant plasmids in complex media. They contain the essential chromosomal gene, dapD, under the control of the lac operator/promoter. Unless supplemented with IPTG (which induces expression of dapD) or DAP, these cells lyse. However, when the strains are transformed with a multicopy plasmid containing the lac operator, the operator competitively titrates the LacI repressor and allows expression of dapD from the lac promoter. Thus transformants can be isolated and propagated simply by their ability to grow on any medium by repressor titration selection. No antibiotic resistance genes or other protein expressing sequences are required on the plasmid, and antibiotics are not necessary for plasmid selection, making these strains a valuable tool for therapeutic DNA and recombinant protein production. We describe the construction of these strains and demonstrate plasmid selection and maintenance by repressor titration, using the new pORT plasmid vectors designed to facilitate recombinant DNA exploitation.

Structural variation of the pseudoautosomal region between and within inbred mouse strains.

The pseudoautosomal region (PAR) is a segment of shared homology between the sex chromosomes. Here we report additional probes for this region of the mouse genome. Genetic and fluorescence in situ hybridization analyses indicate that one probe, PAR-4, hybridizes to the pseudoautosomal telomere and a minor locus at the telomere of chromosome 9 and that a PCR assay based on the PAR-4 sequence amplifies only the pseudoautosomal locus (DXYHgu1). The region detected by PAR-4 is structurally unstable; it shows polymorphism both between mouse strains and between animals of the same inbred strain, which implies an unusually high mutation rate. Variation occurs in the region adjacent to a (TTAGGG)n array. Two pseudoautosomal probes can also hybridize to the distal telomeres of chromosomes 9 and 13, and all three telomeres contain DXYMov15. The similarity between these telomeres may reflect ancestral telomere-telomere exchange.

Receptor usage and differential downregulation of CD46 by measles virus wild-type and vaccine strains.

Recently, two cell surface molecules, CD46 and moesin, have been found to be functionally associated with measles virus (MV) infectivity of cells. We investigated the receptor usage of MV wild-type, subacute sclerosing panencephalitis, and vaccine strains and their effect on the down-regulation of CD46 after infection. We found that the infection of human cell lines with all 19 MV strains tested was inhibitable with antibodies against CD46. In contrast, not all strains of MV led to the downregulation of CD46 following infection. The group of CD46 non-downregulating strains comprised four lymphotropic wild-type isolates designated AB, DF, DL, and WTF. Since the downregulation of CD46 is caused by interaction with newly synthesized MV hemagglutinin (MV-H), we tested the capability of recombinant MV-H proteins to downregulate CD46. Recombinant MV-H proteins of MV strains Edmonston, Halle, and CM led to the down-regulation of CD46, whereas those of DL and WTF did not. This observed differential downregulation by different MV strains has profound consequences, since lack of CD46 on the cell surface leads to susceptibility of cells to complement lysis. These results suggest that lymphotropic wild-type strains of MV which do not downregulate CD46 may have an advantage for replication in vivo. The relatively weak immune response against attenuated vaccine strains of MV compared with wild-type strains might be related to this phenomenon.

Physical interaction between components of DNA mismatch repair and nucleotide excision repair

Nucleotide excision repair (NER) and DNA mismatch repair are required for some common processes although the biochemical basis for this requirement is unknown. Saccharomyces cerevisiae RAD14 was identified in a two-hybrid screen using MSH2 as “bait,” and pairwise interactions between MSH2 and RAD1, RAD2, RAD3, RAD10, RAD14, and RAD25 subsequently were demonstrated by two-hybrid analysis. MSH2 coimmunoprecipitated specifically with epitope-tagged versions of RAD2, RAD10, RAD14, and RAD25. MSH2 and RAD10 were found to interact in msh3 msh6 and mlh1 pms1 double mutants, suggesting a direct interaction with MSH2. Mutations in MSH2 increased the UV sensitivity of NER-deficient yeast strains, and msh2 mutations were epistatic to the mutator phenotype observed in NER-deficient strains. These data suggest that MSH2 and possibly other components of DNA mismatch repair exist in a complex with NER proteins, providing a biochemical and genetical basis for these proteins to function in common processes.

The p20 and Ded1 proteins have antagonistic roles in eIF4E-dependent translation in Saccharomyces cerevisiae

The translation initiation factor eIF4E mediates the binding of the small ribosomal subunit to the cap structure at the 5′ end of the mRNA. In Saccharomyces cerevisiae, the cap-binding protein eIF4E is mainly associated with eIF4G, forming the cap-binding complex eIF4F. Other proteins are detected upon purification of the complex on cap-affinity columns. Among them is p20, a protein of unknown function encoded by the CAF20 gene. Here, we show a negative regulatory role for the p20 protein in translation initiation. Deletion of CAF20 partially suppresses mutations in translation initiation factors. Overexpression of the p20 protein results in a synthetic enhancement of translation mutation phenotypes. Similar effects are observed for mutations in the DED1 gene, which we have isolated as a multicopy suppressor of a temperature-sensitive eIF4E mutation. The DED1 gene encodes a putative RNA helicase of the DEAD-box family. The analyses of its suppressor activity, of polysome profiles of ded1 mutant strains, and of synthetic lethal interactions with different translation mutants indicate that the Ded1 protein has a role in translation initiation in S. cerevisiae.

The mating-type and pathogenicity locus of the fungus Ustilago hordei spans a 500-kb region

The fungal pathogen Ustilago hordei causes the covered smut disease of barley and oats. Mating and pathogenicity in this fungus are controlled by the MAT locus, which contains two distinct gene complexes, a and b. In this study, we tagged the a and b regions with the recognition sequence for the restriction enzyme I-SceI and determined that the distance between the complexes is 500 kb in a MAT-1 strain and 430 kb in a MAT-2 strain. Characterization of the organization of the known genes within the a and b gene complexes provided evidence for nonhomology and sequence inversion between MAT-1 and MAT-2. Antibiotic-resistance markers also were used to tag the a gene complex in MAT-1 strains (phleomycin) and the b gene complex in MAT-2 strains (hygromycin). Crosses were performed with these strains and progeny resistant to both antibiotics were recovered at a very low frequency, suggesting that recombination is suppressed within the MAT region. Overall, the chromosome homologues carrying the MAT locus of U. hordei share features with primitive sex chromosomes, with the added twist that the MAT locus also controls pathogenicity.

The prion model for [URE3] of yeast: Spontaneous generation and requirements for propagation

The genetic properties of the non-Mendelian element, [URE3], suggest that it is a prion (infectious protein) form of Ure2p, a mediator of nitrogen regulation in Saccharomyces cerevisiae. Into a ure2Δ strain (necessarily lacking [URE3]), we introduced a plasmid overproducing Ure2p. This induced the frequent “spontaneous generation” of [URE3], with properties identical to the original [URE3]. Altering the translational frame only in the prion-inducing domain of URE2 shows that it is Ure2 protein (and not URE2 RNA) that induces appearance of [URE3]. The proteinase K-resistance of Ure2p is unique to [URE3] strains and is not seen in nitrogen regulation of normal strains. The prion-inducing domain of Ure2p (residues 1–65) can propagate [URE3] in the absence of the C-terminal part of the molecule. In contrast, the C-terminal part of Ure2p cannot be converted to the prion (inactive) form without the prion-inducing domain covalently attached. These experiments support the prion model for [URE3] and extend our understanding of its propagation.

Effects of Genome Position and the DNA Damage Checkpoint on the Structure and Frequency of sod2 Gene Amplification in Fission Yeast

The Schizosaccharomyces pombe sod2 gene, located near the telomere on the long arm of chromosome I, encodes a Na+ (or Li+)/H+ antiporter. Amplification of sod2 has previously been shown to confer resistance to LiCl. We analyzed 20 independent LiCl-resistant strains and found that the only observed mechanism of resistance is amplification of sod2. The amplicons are linear, extrachromosomal elements either 225 or 180 kb long, containing both sod2 and telomere sequences. To determine whether proximity to a telomere is necessary for sod2 amplification, a strain was constructed in which the gene was moved to the middle of the same chromosomal arm. Selection of LiCl-resistant strains in this genetic background also yielded amplifications of sod2, but in this case the amplified DNA was exclusively chromosomal. Thus, proximity to a telomere is not a prerequisite for gene amplification in S. pombe but does affect the mechanism. Relative to wild-type cells, mutants with defects in the DNA damage aspect of the rad checkpoint control pathway had an increased frequency of sod2 amplification, whereas mutants defective in the S-phase completion checkpoint did not. Two models for generating the amplified DNA are presented.

Regional and strain-specific gene expression mapping in the adult mouse brain

To determine the genetic causes and molecular mechanisms responsible for neurobehavioral differences in mice, we used highly parallel gene expression profiling to detect genes that are differentially expressed between the 129SvEv and C57BL/6 mouse strains at baseline and in response to seizure. In addition, we identified genes that are differentially expressed in specific brain regions. We found that approximately 1% of expressed genes are differentially expressed between strains in at least one region of the brain and that the gene expression response to seizure is significantly different between the two inbred strains. The results lead to the identification of differences in gene expression that may account for distinct phenotypes in inbred strains and the unique functions of specific brain regions.

Live attenuated HIV vaccines: Predicting the tradeoff between efficacy and safety

The utility of live attenuated vaccines for controlling HIV epidemics is being debated. Live attenuated HIV vaccines (LAHVs) could be extremely effective in protecting against infection with wild-type strains, but may not be completely safe as the attenuated strain could cause AIDS in some vaccinated individuals. We present a theoretical framework for evaluating the consequences of the tradeoff between vaccine efficacy (in terms of preventing new infections with wild-type strains) and safety (in terms of vaccine-induced AIDS deaths). We use our framework to predict, for Zimbabwe and Thailand, the epidemiological impact of 1,000 different (specified by efficacy and safety characteristics) LAHVs. We predict that paradoxically: (i) in Zimbabwe (where transmission is high) LAHVs would significantly decrease the AIDS death rate, but (ii) in Thailand (where transmission is low) exactly the same vaccines (in terms of efficacy and safety characteristics) would increase the AIDS death rate. Our results imply that a threshold transmission rate exists that determines whether any given LAHV has a beneficial or a detrimental impact. We also determine the vaccine perversity point, which is defined in terms of the fraction of vaccinated individuals who progress to AIDS as a result of the vaccine strain. Vaccination with any LAHV that causes more than 5% of vaccinated individuals to progress to AIDS in 25 years would, even 50 years later, lead to perversity (i.e., increase the annual AIDS death rate) in Thailand; these same vaccines would lead to decreases in the annual AIDS death rate in Zimbabwe.

Nitrogen management and the future of food: Lessons from the management of energy and carbon

The food system dominates anthropogenic disruption of the nitrogen cycle by generating excess fixed nitrogen. Excess fixed nitrogen, in various guises, augments the greenhouse effect, diminishes stratospheric ozone, promotes smog, contaminates drinking water, acidifies rain, eutrophies bays and estuaries, and stresses ecosystems. Yet, to date, regulatory efforts to limit these disruptions largely ignore the food system. There are many parallels between food and energy. Food is to nitrogen as energy is to carbon. Nitrogen fertilizer is analogous to fossil fuel. Organic agriculture and agricultural biotechnology play roles analogous to renewable energy and nuclear power in political discourse. Nutrition research resembles energy end-use analysis. Meat is the electricity of food. As the agriculture and food system evolves to contain its impacts on the nitrogen cycle, several lessons can be extracted from energy and carbon: (i) set the goal of ecosystem stabilization; (ii) search the entire production and consumption system (grain, livestock, food distribution, and diet) for opportunities to improve efficiency; (iii) implement cap-and-trade systems for fixed nitrogen; (iv) expand research at the intersection of agriculture and ecology, and (v) focus on the food choices of the prosperous. There are important nitrogen-carbon links. The global increase in fixed nitrogen may be fertilizing the Earth, transferring significant amounts of carbon from the atmosphere to the biosphere, and mitigating global warming. A modern biofuels industry someday may produce biofuels from crop residues or dedicated energy crops, reducing the rate of fossil fuel use, while losses of nitrogen and other nutrients are minimized.

Genetic rearrangements in the rfb regions of Vibrio cholerae O1 and O139.

The recent emergence of a pathogenic new non-O1 serotype (O139) of Vibrio cholerae has led to numerous studies in an attempt to identify the origins of this new strain. Our studies indicate that O139 strains have clear differences in the surface polysaccharides when compared with O1 strains: the lipopolysaccharide can be described as semi-rough. Southern hybridization with the O1 rfb region demonstrates that O139 strains no longer contain any of the rfb genes required for the synthesis of the O1 O-antigen or its modification and also lack at least 6 kb of additional contiguous DNA. However, O139 strains have retained rfaD and have a single open reading frame closely related to three small open reading frames of the O1 rfb region. This region is closely related to the H-repeat of Escherichia coli and to the transposases of a number of insertion sequence elements and has all the features of an insertion sequence element that has been designated VcIS1. Transposon insertion mutants defective in O139 O-antigen (and capsule) biosynthesis map to the same fragment as VcIS1. Preliminary sequence data of complementing clones indicate that this DNA encodes a galactosyl-transferase and other enzymes for the utilization of galactose in polysaccharide biosynthesis. We propose a mechanism by which both the Ogawa serotype of O1 strains and the O139 serotype strains may have evolved.

Two-domain reconstitution of a functional protein histidine kinase

In prokaryotes, in the absence of protein serine/threonine/tyrosine kinases, protein histidine kinases play a major role in signal transduction involved in cellular adaptation to various environmental changes and stresses. Histidine kinases phosphorylate their cognate response regulators at a specific aspartic acid residue with ATP in response to particular environmental signals. In this His-Asp phosphorelay signal transduction system, it is still unknown how the histidine kinase exerts its enzymatic function. Here we demonstrate that the cytoplasmic kinase domain of EnvZ, a transmembrane osmosensor of Escherichia coli can be further divided into two distinct functional subdomains: subdomain A [EnvZ(C)⋅(223–289); 67 residues] and subdomain B [EnvZ(C)⋅(290–450); 161 residues]. Subdomain A, with a high helical content, contains the autophosphorylation site, H–243, and forms a stable dimer having the recognition site for OmpR, the cognate response regulator of EnvZ. Subdomain B, an α/β-protein, exists as a monomer. When mixed, the two subdomains reconstitute the kinase function to phosphorylate subdomain A at His-243 in the presence of ATP. Subsequently, the phosphorylated subdomain A is able to transfer its phosphate group to OmpR. The two-domain structure of this histidine kinase provides an insight into the structural arrangement of the enzyme and its transphosphorylation mechanism.

Mutagenicity in Escherichia coli of the major DNA adduct derived from the endogenous mutagen malondialdehyde

The spectrum of mutations induced by the naturally occurring DNA adduct pyrimido[1,2-α]purin-10(3H)-one (M1G) was determined by site-specific approaches using M13 vectors replicated in Escherichia coli. M1G was placed at position 6256 in the (−)-strand of M13MB102 by ligating the oligodeoxynucleotide 5′-GGT(M1G)TCCG-3′ into a gapped-duplex derivative of the vector. Unmodified and M1G-modified genomes containing either a cytosine or thymine at position 6256 of the (+)-strand were transformed into repair-proficient and repair-deficient E. coli strains, and base pair substitutions were quantitated by hybridization analysis. Modified genomes containing a cytosine opposite M1G resulted in roughly equal numbers of M1G→A and M1G→T mutations with few M1G→C mutations. The total mutation frequency was ≈1%, which represents a 500-fold increase in mutations compared with unmodified M13MB102. Transformation of modified genomes containing a thymine opposite M1G allowed an estimate to be made of the ability of M1G to block replication. The (−)-strand was replicated >80% of the time in the unadducted genome but only 20% of the time when M1G was present. Correction of the mutation frequency for the strand bias of replication indicated that the actual frequency of mutations induced by M1G was 18%. Experiments using E. coli with different genetic backgrounds indicated that the SOS response enhances the mutagenicity of M1G and that M1G is a substrate for repair by the nucleotide excision repair complex. These studies indicate that M1G, which is present endogenously in DNA of healthy human beings, is a strong block to replication and an efficient premutagenic lesion.