Yuriy Fofanov

Molecular complexity of successive bacterial epidemics deconvoluted by comparative pathogenomics

Understanding the fine-structure molecular architecture of bacterial epidemics has been a long-sought goal of infectious disease research. We used short-read-length DNA sequencing coupled with mass spectroscopy analysis of SNPs to study the molecular pathogenomics of three successive epidemics of invasive infections involving 344 serotype M3 group A Streptococcus in Ontario, Canada. Sequencing the genome of 95 strains from the three epidemics, coupled with analysis of 280 biallelic SNPs in all 344 strains, revealed an unexpectedly complex population structure composed of a dynamic mixture of distinct clonally related complexes. We discovered that each epidemic is dominated by micro- and macrobursts of multiple emergent clones, some with distinct strain genotype–patient phenotype relationships. On average, strains were differentiated from one another by only 49 SNPs and 11 insertion-deletion events (indels) in the core genome. Ten percent of SNPs are strain specific; that is, each strain has a unique genome sequence. We identified nonrandom temporal–spatial patterns of strain distribution within and between the epidemic peaks. The extensive full-genome data permitted us to identify genes with significantly increased rates of nonsynonymous (amino acid-altering) nucleotide polymorphisms, thereby providing clues about selective forces operative in the host. Comparative expression microarray analysis revealed that closely related strains differentiated by seemingly modest genetic changes can have significantly divergent transcriptomes. We conclude that enhanced understanding of bacterial epidemics requires a deep-sequencing, geographically centric, comparative pathogenomics strategy.

Where does Neisseria

BackgroundPathogenicity islands (PAIs) or genomic islands (GEIs) are considered to be the result of a recent horizontal transfer. Detecting PAIs/GEIs as well as their putative source can provide insight into the organism’s pathogenicity within its host. Previously we introduced a tool called S-plot which provides a visual representation of the variation in compositional properties across and between genomic sequences. Utilizing S-plot and new functionality developed here, we examined 18 publicly available Neisseria genomes, including strains of both pathogenic and non-pathogenic species, in order to identify regions of unusual compositional properties (RUCPs) using both a sliding window as well as a gene-by-gene approach.ResultsNumerous GEIs and PAIs were identified including virulence genes previously found within the pathogenic Neisseria species. While some genes were conserved amongst all species, only pathogenic species, or an individual species, a number of genes were detected that are unique to an individual strain. While the majority of such genes have an origin unknown, a number of putative sources including pathogenic and capsule-containing bacteria were determined, indicative of gene exchange between Neisseria spp. and other bacteria within their microhabitat. Furthermore, we uncovered evidence that both N. meningitidis and N. gonorrhoeae have separately acquired DNA from their human host. Data suggests that all three Neisseria species have received horizontally transferred elements post-speciation.ConclusionsUsing this approach, we were able to not only find previously identified regions of virulence but also new regions which may be contributing to the virulence of the species. This comparative analysis provides a means for tracing the evolutionary history of the acquisition of foreign DNA within this genus. Looking specifically at the RUCPs present within the 18 genomes considered, a stronger similarity between N. meningitidis and N. lactamica is observed, suggesting that N. meningitidis arose before N. gonorrhoeae.

Gene expression profile of human lymphoid CEM cells sensitive and resistant to glucocorticoid-evoked apoptosis

Three closely related clones of leukemic lymphoid CEM cells were compared for their gene expression responses to the glucocorticoid dexamethasone (Dex). All three contained receptors for Dex, but only two responded by undergoing apoptosis. After a time of exposure to Dex that ended late in the interval preceding onset of apoptosis, gene microarray analyses were carried out. The results indicate that the expression of a limited, distinctive set of genes was altered in the two apoptosis-prone clones, not in the resistant clone. That clone showed altered expression of different sets of genes, suggesting that a molecular switch converted patterns of gene expression between the two phenotypes: apoptosis-prone and apoptosis-resistant. The results are consistent with the hypothesis that altered expression of a distinctive network of genes after glucocorticoid administration ultimately triggers apoptosis of leukemic lymphoid cells. The altered genes identified provide new foci for study of their role in cell death.

Characterization of Escherichia coli MG1655 grown in a low-shear modeled microgravity environment

BackgroundExtra-cellular shear force is an important environmental parameter that is significant both medically and in the space environment. Escherichia coli cells grown in a low-shear modeled microgravity (LSMMG) environment produced in a high aspect rotating vessel (HARV) were subjected to transcriptional and physiological analysis.ResultsAerobic LSMMG cultures were grown in rich (LB) and minimal (MOPS + glucose) medium with a normal gravity vector HARV control. Reproducible changes in transcription were seen, but no specific LSMMG responsive genes were identified. Instead, absence of shear and a randomized gravity vector appears to cause local extra-cellular environmental changes, which elicit reproducible cellular responses. In minimal media, the majority of the significantly up- or down-regulated genes of known function were associated with the cell envelope. In rich medium, most LSMMG down-regulated genes were involved in translation. No observable changes in post-culture stress responses and antibiotic sensitivity were seen in cells immediately after exposure to LSMMG. Comparison with earlier studies of Salmonella enterica serovar Typhimurium conducted under similar growth conditions, revealed essentially no similarity in the genes that were significantly up- or down-regulated.ConclusionComparison of these results to previous studies suggests that different organisms may dramatically differ in their responses to medically significant low-shear and space environments. Depending on their specific response, some organisms, such as Salmonella, may become preadapted in a manner that predisposes them to increased virulence.

Gene networks in glucocorticoid-evoked apoptosis of leukemic cells.

To discover the genes responsible for the apoptosis evoked by glucocorticoids in leukemic lymphoid cells, we have begun gene array analysis on microchips. Three clones of CEM cells were compared: C7-14, C1-15 and C1-6. C7-14 and C1-15 are subclones from the original clones C7 (sensitive to apoptosis by glucocorticoids) and C1 (resistant). C1-6 is a spontaneous revertant to sensitivity from the C1 clone. Previously we presented data on the sets of genes whose expression is altered in these cell clones after 20 h exposure to dexamethasone (Dex). The two sensitive clones, which respond by undergoing apoptosis starting about 24h after Dex is added, both showed >2.5-fold induction of 39 genes and 2-fold reduction of expressed levels from 21 genes. C1-15, the resistant clone, showed alterations in a separate set of genes. In this paper, we present further analysis of the data on genes regulated in these cell clones after 20 h Dex and compare them with the genes regulated after 12h Dex. Some, but not all the genes found altered at 20 h are altered at 12h, consistent with our hypothesis that sequential gene regulation eventually provokes full apoptosis. We also compare the levels of basal gene expression in the three clones. At the basal level no single gene stands out, but small sets of genes differ >2-fold in basal expression between the two sensitive and the resistant clone. A number of the genes basally higher in the resistant clone are potentially anti-apoptotic. This is consistent with our hypothesis that the resistant cells have undergone a general shift in gene expression.

PIQA: pipeline for Illumina G1 genome analyzer data quality assessment

Summary: PIQA is a quality analysis pipeline designed to examine genomic reads produced by Next Generation Sequencing technology (Illumina G1 Genome Analyzer). A short statistical summary, as well as tile-by-tile and cycle-by-cycle graphical representation of clusters density, quality scores and nucleotide frequencies allow easy identification of various technical problems including defective tiles, mistakes in sample/library preparations and abnormalities in the frequencies of appearance of sequenced genomic reads. PIQA is written in the R statistical programming language and is compatible with bustard, fastq and scarf Illumina G1 Genome Analyzer data formats. Availability: The PIQA pipeline, installation instructions and examples are available at the supplementary web site (http://bioinfo.uh.edu/PIQA). Contact: yfofanov@bioinfo.uh.edu

Ehrlichia chaffeensis TRP120 binds a G+C-rich motif in host cell DNA and exhibits eukaryotic transcriptional activator function.

ABSTRACT Ehrlichia chaffeensis is an obligately intracellular bacterium that modulates host cell gene transcription in the mononuclear phagocyte, but the host gene targets and mechanisms involved in transcriptional modulation are not well-defined. In this study, we identified a novel tandem repeat DNA-binding domain in the E. chaffeensis 120-kDa tandem repeat protein (TRP120) that directly binds host cell DNA. TRP120 was observed by immunofluorescent microscopy in the nucleus of E. chaffeensis-infected host cells and was detected in nuclear extracts by Western immunoblotting with TRP120-specific antibody. The TRP120 binding sites and associated host cell target genes were identified using high-throughput deep sequencing (Illumina) of immunoprecipitated DNA (chromatin immunoprecipitation and high-throughput DNA sequencing). Multiple em motif elicitation (MEME) analysis of the most highly enriched TRP120-bound sequences revealed a G+C-rich DNA motif, and recombinant TRP120 specifically bound synthetic oligonucleotides containing the motif. TRP120 target gene binding sites were mapped most frequently to intersecting regions (intron/exon; 49%) but were also identified in upstream regulatory regions (25%) and downstream locations (26%). Genes targeted by TRP120 were most frequently associated with transcriptional regulation, signal transduction, and apoptosis. TRP120 targeted inflammatory chemokine genes, CCL2, CCL20, and CXCL11, which were strongly upregulated during E. chaffeensis infection and were also upregulated by direct transfection with recombinant TRP120. This study reveals that TRP120 is a novel DNA-binding protein that is involved in a host gene transcriptional regulation strategy.

Whole Genome Sequencing of Mutation Accumulation Lines Reveals a Low Mutation Rate in the Social Amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10−9, with a Poisson confidence interval of 4.1×10−9 − 9.5×10−9, per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10−11, with a Poisson confidence interval ranging from 7.4×10−13 to 1.6×10−10, is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.

Metagenomic Analysis of the Airborne Environment in Urban Spaces

The organisms in aerosol microenvironments, especially densely populated urban areas, are relevant to maintenance of public health and detection of potential epidemic or biothreat agents. To examine aerosolized microorganisms in this environment, we performed sequencing on the material from an urban aerosol surveillance program. Whole metagenome sequencing was applied to DNA extracted from air filters obtained during periods from each of the four seasons. The composition of bacteria, plants, fungi, invertebrates, and viruses demonstrated distinct temporal shifts. Bacillus thuringiensis serovar kurstaki was detected in samples known to be exposed to aerosolized spores, illustrating the potential utility of this approach for identification of intentionally introduced microbial agents. Together, these data demonstrate the temporally dependent metagenomic complexity of urban aerosols and the potential of genomic analytical techniques for biosurveillance and monitoring of threats to public health.

Human-blind probes and primers for dengue virus identification: Exhaustive analysis of subsequences present in the human and 83 dengue genome sequences

Reliable detection and identification of pathogens in complex biological samples, in the presence of contaminating DNA from a variety of sources, is an important and challenging diagnostic problem for the development of field tests. The problem is compounded by the difficulty of finding a single, unique genomic sequence that is present simultaneously in all genomes of a species of closely related pathogens and absent in the genomes of the host or the organisms that contribute to the sample background. Here we describe ‘host‐blind probe design’– a novel strategy of designing probes based on highly frequent genomic signatures found in the pathogen genomes of interest but absent from the host genome. Upon hybridization, an array of such informative probes will produce a unique pattern that is a genetic fingerprint for each pathogen strain. This multiprobe approach was applied to 83 dengue virus genome sequences, available in public databases, to design and perform in silico microarray experiments. The resulting patterns allow one to unequivocally distinguish the four major serotypes, and within each serotype to identify the most similar strain among those that have been completely sequenced. In an environment where dengue is indigenous, this would allow investigators to determine if a particular isolate belongs to an ongoing outbreak or is a previously circulating version. Using our probe set, the probability that misdiagnosis at the serotype level would occur is ≈ 1 : 10150.