Genome Sequence of Elizabethkingia anophelis Strain EaAs1, Isolated from the Asian Malaria Mosquito Anopheles stephensi.

We sequenced the genome of a strain of the Gram-negative bacterial species Elizabethkingia anophelis, which is an important component of the Anopheles mosquito microbiome. This genome sequence will add to the list of resources used to examine host-microbe interactions in mosquitoes.

Genome Sequences of Staphylococcus hominis Strains ShAs1, ShAs2, and ShAs3, Isolated from the Asian Malaria Mosquito Anopheles stephensi.

Staphylococcus hominis is a culturable component of the bacterial microbiome of Anopheles stephensi. Here, we present the annotated draft genome sequences of three S. hominis isolates from A. stephensi. These genomic resources will facilitate experiments to further our understanding of the role of bacteria in mosquito biology.

Genome Sequence of Stenotrophomonas maltophilia Strain SmAs1, Isolated From the Asian Malaria Mosquito Anopheles stephensi.

An isolate of Stenotrophomonas maltophilia was cultured from the Asian malaria vector Anopheles stephensi. Here, we present the annotated draft genome sequence of this S. maltophilia strain. This genomic resource will facilitate further characterization of bacteria associated with mosquitoes.

Decoupling of divergent gene regulation by sequence-specific DNA binding factors.

Divergent gene pairs (DGPs) are abundant in eukaryotic genomes. Since two genes in a DGP potentially share the same regulatory sequence, one might expect that they should be co-regulated. However, an inspection of yeast DGPs containing cell-cycle or stress response genes revealed that most DGPs are differentially-regulated. The mechanism underlying DGP differential regulation is not understood. Here, we showed that co- versus differential regulation cannot be explained by genetic features including promoter length, binding site orientation, TATA elements, nucleosome distribution, or presence of non-coding RNAs. Using time-lapse fluorescence microscopy, we carried out an in-depth study of a differentially regulated DGP, PFK26-MOB1. We found that their differential regulation is mainly achieved through two DNA-binding factors, Tbf1 and Mcm1. Similar to 'enhancer-blocking insulators' in higher eukaryotes, these factors shield the proximal promoter from the action of more distant transcription regulators. We confirmed the blockage function of Tbf1 using synthetic promoters. We further presented evidence that the blockage mechanism is widely used among genome-wide DGPs. Besides elucidating the DGP regulatory mechanism, our work revealed a novel class of insulators in yeast.

Complete Genome Sequence of SS52, a Strain of Escherichia coli O157:H7 Recovered from Supershedder Cattle.

Shiga toxin-producing Escherichia coli O157:H7 causes foodborne infections, and cattle are the primary reservoir. Some animals, known as supershedders, excrete orders of magnitude more E. coli O157:H7 in the feces than normal. Here, we report the complete genome sequence of the SS52 supershedder strain of E. coli O157:H7.

Comparative analysis of super-shedder strains of Escherichia coli O157:H7 reveals distinctive genomic features and a strongly aggregative adherent phenotype on bovine rectoanal junction squamous epithelial cells.

Shiga toxin-producing Escherichia coli O157:H7 (O157) are significant foodborne pathogens and pose a serious threat to public health worldwide. The major reservoirs of O157 are asymptomatic cattle which harbor the organism in the terminal recto-anal junction (RAJ). Some colonized animals, referred to as "super-shedders" (SS), are known to shed O157 in exceptionally large numbers (>104 CFU/g of feces). Recent studies suggest that SS cattle play a major role in the prevalence and transmission of O157, but little is known about the molecular mechanisms associated with super-shedding. Whole genome sequence analysis of an SS O157 strain (SS17) revealed a genome of 5,523,849 bp chromosome with 5,430 open reading frames and two plasmids, pO157 and pSS17, of 94,645 bp and 37,446 bp, respectively. Comparative analyses showed that SS17 is clustered with spinach-associated O157 outbreak strains, and belongs to the lineage I/II, clade 8, D group, and genotype 1, a subgroup of O157 with predicted hyper-virulence. A large number of non-synonymous SNPs and other polymorphisms were identified in SS17 as compared with other O157 strains (EC4115, EDL933, Sakai, TW14359), including in key adherence- and virulence-related loci. Phenotypic analyses revealed a distinctive and strongly adherent aggregative phenotype of SS17 on bovine RAJ stratified squamous epithelial (RSE) cells that was conserved amongst other SS isolates. Molecular genetic and functional analyses of defined mutants of SS17 suggested that the strongly adherent aggregative phenotype amongst SS isolates is LEE-independent, and likely results from a novel mechanism. Taken together, our study provides a rational framework for investigating the molecular mechanisms associated with SS, and strong evidence that SS O157 isolates have distinctive features and use a LEE-independent mechanism for hyper-adherence to bovine rectal epithelial cells.

Stochastic expression and epigenetic memory at the yeast HO promoter.

Eukaryotic gene regulation usually involves sequence-specific transcription factors and sequence-nonspecific cofactors. A large effort has been made to understand how these factors affect the average gene expression level among a population. However, little is known about how they regulate gene expression in individual cells. In this work, we address this question by mutating multiple factors in the regulatory pathway of the yeast HO promoter (HOpr) and probing the corresponding promoter activity in single cells using time-lapse fluorescence microscopy. We show that the HOpr fires in an "on/off" fashion in WT cells as well as in different genetic backgrounds. Many chromatin-related cofactors that affect the average level of HO expression do not actually affect the firing amplitude of the HOpr; instead, they affect the firing frequency among individual cell cycles. With certain mutations, the bimodal expression exhibits short-term epigenetic memory across the mitotic boundary. This memory is propagated in "cis" and reflects enhanced activator binding after a previous "on" cycle. We present evidence that the memory results from slow turnover of the histone acetylation marks.

Whole-genome sequencing reveals a link between ß-lactam resistance and synthetases of the alarmone (p)ppGpp in Staphylococcus aureus.

The overwhelming majority of methicillin-resistant Staphylococcus aureus (MRSA) clinical isolates exhibit a peculiar heterogeneous resistance to ß-lactam antibiotics: in cultures of such strains, the majority of cells display only a low level of methicillin resistance--often close to the MIC breakpoint of susceptible strains. Yet, in the same cultures, subpopulations of bacteria exhibiting very high levels of resistance are also present with variable frequencies, which are characteristic of the particular MRSA lineage. The mechanism of heterogeneous resistance is not understood. We describe here an experimental system for exploring the mechanism of heterogeneous resistance. Copies of the resistance gene mecA cloned into a temperature-sensitive plasmid were introduced into the fully sequenced methicillin-susceptible clinical isolate S. aureus strain 476. Transductants of strain 476 expressed methicillin resistance in a heterogeneous fashion: the great majority of cells showed only low MIC (0.75 µg/ml) for the antibiotic, but a minority population of highly resistant bacteria (MIC >300 µg/ml) was also present with a frequency of ∼10(-4). The genetic backgrounds of the majority and minority cells were compared by whole-genome sequencing: the only differences detectable were two point mutations in relA of the highly resistant minority population of bacteria. The relA gene codes for the synthesis of (p)ppGpp, an effector of the stringent stress response. Titration of (p)ppGpp showed increased amounts of this effector in the highly resistant cells. Involvement of (p)ppGpp synthesis genes may explain some of the perplexing aspects of ß-lactam resistance in MRSA, since many environmental and genetic changes can modulate cellular levels of (p)ppGpp.

Transcriptome-wide sequencing reveals numerous APOBEC1 mRNA-editing targets in transcript 3' UTRs.

Apolipoprotein B-editing enzyme, catalytic polypeptide-1 (APOBEC1) is a cytidine deaminase initially identified by its activity in converting a specific cytidine (C) to uridine (U) in apolipoprotein B (apoB) mRNA transcripts in the small intestine. Editing results in the translation of a truncated apoB isoform with distinct functions in lipid transport. To address the possibility that APOBEC1 edits additional mRNAs, we developed a transcriptome-wide comparative RNA sequencing (RNA-Seq) screen. We identified and validated 32 previously undescribed mRNA targets of APOBEC1 editing, all of which are located in AU-rich segments of transcript 3' untranslated regions (3' UTRs). Further analysis established several characteristic sequence features of editing targets, which were predictive for the identification of additional APOBEC1 substrates. The transcriptomics approach to RNA editing presented here dramatically expands the list of APOBEC1 mRNA editing targets and reveals a novel cellular mechanism for the modification of transcript 3' UTRs.

Tracking the in vivo evolution of multidrug resistance in Staphylococcus aureus by whole-genome sequencing.

The spread of multidrug-resistant Staphylococcus aureus (MRSA) strains in the clinical environment has begun to pose serious limits to treatment options. Yet virtually nothing is known about how resistance traits are acquired in vivo. Here, we apply the power of whole-genome sequencing to identify steps in the evolution of multidrug resistance in isogenic S. aureus isolates recovered periodically from the bloodstream of a patient undergoing chemotherapy with vancomycin and other antibiotics. After extensive therapy, the bacterium developed resistance, and treatment failed. Sequencing the first vancomycin susceptible isolate and the last vancomycin nonsusceptible isolate identified genome wide only 35 point mutations in 31 loci. These mutations appeared in a sequential order in isolates that were recovered at intermittent times during chemotherapy in parallel with increasing levels of resistance. The vancomycin nonsusceptible isolates also showed a 100-fold decrease in susceptibility to daptomycin, although this antibiotic was not used in the therapy. One of the mutated loci associated with decreasing vancomycin susceptibility (the vraR operon) was found to also carry mutations in six additional vancomycin nonsusceptible S. aureus isolates belonging to different genetic backgrounds and recovered from different geographic sites. As costs drop, whole-genome sequencing will become a useful tool in elucidating complex pathways of in vivo evolution in bacterial pathogens.

Genome wide identification of regulatory motifs in Bacillus subtilis.

BACKGROUND: To explain the vastly different phenotypes exhibited by the same organism under different conditions, it is essential that we understand how the organism's genes are coordinately regulated. While there are many excellent tools for predicting sequences encoding proteins or RNA genes, few algorithms exist to predict regulatory sequences on a genome wide scale with no prior information. RESULTS: To identify motifs involved in the control of transcription, an algorithm was developed that searches upstream of operons for improbably frequent dimers. The algorithm was applied to the B. subtilis genome, which is predicted to encode for approximately 200 DNA binding proteins. The dimers found to be over-represented could be clustered into 317 distinct groups, each thought to represent a class of motifs uniquely recognized by some transcription factor. For each cluster of dimers, a representative weight matrix was derived and scored over the regions upstream of the operons to predict the sites recognized by the cluster's factor, and a putative regulon of the operons immediately downstream of the sites was inferred. The distribution in number of operons per predicted regulon is comparable to that for well characterized transcription factors. The most highly over-represented dimers matched sigmaA, the T-box, and sigmaW sites. We have evidence to suggest that at least 52 of our clusters of dimers represent actual regulatory motifs, based on the groups' weight matrix matches to experimentally characterized sites, the functional similarity of the component operons of the groups' regulons, and the positional biases of the weight matrix matches. All predictions are assigned a significance value, and thresholds are set to avoid false positives. Where possible, we examine our false negatives, drawing examples from known regulatory motifs and regulons inferred from RNA expression data. CONCLUSIONS: We have demonstrated that in the case of B. subtilis our algorithm allows for the genome wide identification of regulatory sites. As well as recovering known sites, we predict new sites of yet uncharacterized factors. Results can be viewed at http://www.physics.rockefeller.edu/~mwangi/.