Jean Marie Rouillard

Yersinia pestis and the Plague of Justinian 541–543 AD: a genomic analysis

BACKGROUND Yersinia pestis has caused at least three human plague pandemics. The second (Black Death, 14-17th centuries) and third (19-20th centuries) have been genetically characterised, but there is only a limited understanding of the first pandemic, the Plague of Justinian (6-8th centuries). To address this gap, we sequenced and analysed draft genomes of Y pestis obtained from two individuals who died in the first pandemic. METHODS Teeth were removed from two individuals (known as A120 and A76) from the early medieval Aschheim-Bajuwarenring cemetery (Aschheim, Bavaria, Germany). We isolated DNA from the teeth using a modified phenol-chloroform method. We screened DNA extracts for the presence of the Y pestis-specific pla gene on the pPCP1 plasmid using primers and standards from an established assay, enriched the DNA, and then sequenced it. We reconstructed draft genomes of the infectious Y pestis strains, compared them with a database of genomes from 131 Y pestis strains from the second and third pandemics, and constructed a maximum likelihood phylogenetic tree. FINDINGS Radiocarbon dating of both individuals (A120 to 533 AD [plus or minus 98 years]; A76 to 504 AD [plus or minus 61 years]) places them in the timeframe of the first pandemic. Our phylogeny contains a novel branch (100% bootstrap at all relevant nodes) leading to the two Justinian samples. This branch has no known contemporary representatives, and thus is either extinct or unsampled in wild rodent reservoirs. The Justinian branch is interleaved between two extant groups, 0.ANT1 and 0.ANT2, and is distant from strains associated with the second and third pandemics. INTERPRETATION We conclude that the Y pestis lineages that caused the Plague of Justinian and the Black Death 800 years later were independent emergences from rodents into human beings. These results show that rodent species worldwide represent important reservoirs for the repeated emergence of diverse lineages of Y pestis into human populations. FUNDING McMaster University, Northern Arizona University, Social Sciences and Humanities Research Council of Canada, Canada Research Chairs Program, US Department of Homeland Security, US National Institutes of Health, Australian National Health and Medical Research Council.

Evolution combined with genomic study elucidates genetic bases of isobutanol tolerance in Escherichia coli

BackgroundIsobutanol is a promising next-generation biofuel with demonstrated high yield microbial production, but the toxicity of this molecule reduces fermentation volumetric productivity and final titer. Organic solvent tolerance is a complex, multigenic phenotype that has been recalcitrant to rational engineering approaches. We apply experimental evolution followed by genome resequencing and a gene expression study to elucidate genetic bases of adaptation to exogenous isobutanol stress.ResultsThe adaptations acquired in our evolved lineages exhibit antagonistic pleiotropy between minimal and rich medium, and appear to be specific to the effects of longer chain alcohols. By examining genotypic adaptation in multiple independent lineages, we find evidence of parallel evolution in marC, hfq, mdh, acrAB, gatYZABCD, and rph genes. Many isobutanol tolerant lineages show reduced RpoS activity, perhaps related to mutations in hfq or acrAB. Consistent with the complex, multigenic nature of solvent tolerance, we observe adaptations in a diversity of cellular processes. Many adaptations appear to involve epistasis between different mutations, implying a rugged fitness landscape for isobutanol tolerance. We observe a trend of evolution targeting post-transcriptional regulation and high centrality nodes of biochemical networks. Collectively, the genotypic adaptations we observe suggest mechanisms of adaptation to isobutanol stress based on remodeling the cell envelope and surprisingly, stress response attenuation.ConclusionsWe have discovered a set of genotypic adaptations that confer increased tolerance to exogenous isobutanol stress. Our results are immediately useful to further efforts to engineer more isobutanol tolerant host strains of E. coli for isobutanol production. We suggest that rpoS and post-transcriptional regulators, such as hfq, RNA helicases, and sRNAs may be interesting mutagenesis targets for future global phenotype engineering.

Versatile design and synthesis platform for visualizing genomes with Oligopaint FISH probes

A host of observations demonstrating the relationship between nuclear architecture and processes such as gene expression have led to a number of new technologies for interrogating chromosome positioning. Whereas some of these technologies reconstruct intermolecular interactions, others have enhanced our ability to visualize chromosomes in situ. Here, we describe an oligonucleotide- and PCR-based strategy for fluorescence in situ hybridization (FISH) and a bioinformatic platform that enables this technology to be extended to any organism whose genome has been sequenced. The oligonucleotide probes are renewable, highly efficient, and able to robustly label chromosomes in cell culture, fixed tissues, and metaphase spreads. Our method gives researchers precise control over the sequences they target and allows for single and multicolor imaging of regions ranging from tens of kilobases to megabases with the same basic protocol. We anticipate this technology will lead to an enhanced ability to visualize interphase and metaphase chromosomes.

Gene2Oligo: oligonucleotide design for in vitro gene synthesis.

There is substantial interest in implementing a bioinformatics tool that allows the design of oligonucleotides to support the development of in vitro gene synthesis. Current protocols to make long synthetic DNA molecules rely on the in vitro assembly of a set of short oligonucleotides, either by ligase chain reaction (LCR) or by assembly PCR. Ideally, such oligonucleotides should represent both strands of the final DNA molecule. They should be adjacent on the same strand and overlap the complementary oligonucleotides from the second strand to ensure good hybridization during assembly. This implies that the thermodynamic properties of each oligonucleotide have to be consistent across the set. Furthermore, any given oligonucleotide has to be totally specific to its target to avoid the creation of incorrectly assembled sequences. We have developed Gene2Oligo (http://berry.engin.umich.edu/gene2oligo/), a web-based tool that divides a long input DNA sequence into a set of adjacent oligonucleotides representing both DNA strands. The length of the oligonucleotides is dynamically optimized to ensure both the specificity and the uniform melting temperatures necessary for in vitro gene synthesis. We have successfully designed and used a set of oligonucleotides to synthesize the Saccharomyces cerevisiae cytochrome b5 by using both LCR and assembly PCR.

Biphenyl and Benzoate Metabolism in a Genomic Context: Outlining Genome-Wide Metabolic Networks in Burkholderia xenovorans LB400

ABSTRACT We designed and successfully implemented the use of in situ-synthesized 45-mer oligonucleotide DNA microarrays (XeoChips) for genome-wide expression profiling of Burkholderia xenovorans LB400, which is among the best aerobic polychlorinated biphenyl degraders known so far. We conducted differential gene expression profiling during exponential growth on succinate, benzoate, and biphenyl as sole carbon sources and investigated the transcriptome of early-stationary-phase cells grown on biphenyl. Based on these experiments, we outlined metabolic pathways and summarized other cellular functions in the organism relevant for biphenyl and benzoate degradation. All genes previously identified as being directly involved in biphenyl degradation were up-regulated when cells were grown on biphenyl compared to expression in succinate-grown cells. For benzoate degradation, however, genes for an aerobic coenzyme A activation pathway were up-regulated in biphenyl-grown cells, while the pathway for benzoate degradation via hydroxylation was up-regulated in benzoate-grown cells. The early-stationary-phase biphenyl-grown cells showed similar expression of biphenyl pathway genes, but a surprising up-regulation of C1 metabolic pathway genes was observed. The microarray results were validated by quantitative reverse transcription PCR with a subset of genes of interest. The XeoChips showed a chip-to-chip variation of 13.9%, compared to the 21.6% variation for spotted oligonucleotide microarrays, which is less variation than that typically reported for PCR product microarrays.

Integrating cancer genomics and proteomics in the post-genome era

The dawn of the post‐genome era is leading to extraordinary opportunities in biomedicine. Our group has embarked on a major effort to integrate genomics, transcriptomics and proteomics for the profiling of tumor tissues, an approach we refer to as operomics. Our major goals are the molecular classification of tumors and the identification of markers for the early detection of cancer. Molecular analyses of tumors rely on microdissected tissues, which are simultaneously investigated for genomic, transcriptomic and proteomic changes. Genomic alterations in tumor cells being investigated include deletions, amplifications and methylation changes across the entire genome as well as point mutations in specific genes. Expression analysis at the RNA level is being undertaken using oligonucleotide and cDNA based microarrays. An important aspect of our approach is the large‐scale identification and quantitative analysis of tumor proteins in whole cell lysates as well as in protein compartments. Protein separation strategies include two‐dimensional polyacrylamide gel electrophoresis and liquid chromatography. Specific protein subsets, of interest include membrane proteins, secreted proteins and antigenic proteins as sources of biomarkers for early detection of cancer. Our current approach is illustrated with findings stemming from our studies of human gliomas.

Amplification and overexpression of the IGF2 regulator PLAG1 in hepatoblastoma

There is evidence that 8q amplification is associated with poor prognosis in hepatoblastoma. A previous comparative genomic hybridization analysis identified a critical region in chromosomal bands 8q11.2–q13. Using restriction landmark genomic scanning in combination with a virtual genome scan, we showed that this region is delineated by sequences within contig NT_008183 of chromosomal subbands 8q11.22–q11.23. A real‐time PCR–based genomic copy number assay of 20 hepatoblastomas revealed gain or amplification in this critical chromosomal region in eight tumors. The expression of four genes and expressed sequence tags (ESTs) within this newly defined region was assayed by real‐time reverse transcriptase polymerase chain reaction (RT‐PCR) in four tumors with and six tumors without gain or amplification. The PLAG1 oncogene was found to be highly expressed in all but one tumor compared to normal liver tissue. Furthermore, quantitative RT‐PCR revealed that the expression level of the developmentally regulated transcription factor PLAG1 was 3–12 times greater in hepatoblastoma tumors and cell lines compared to age‐matched normal liver and comparable to the expression in fetal liver tissue. PLAG1 has been shown be a transcriptional activator of IGF2 in other tumor types. Using luciferase reporter assays, we demonstrated that PLAG1 transactivates transcription from the embryonic IGF2 promoter P3, also in hepatoblastoma cell lines. Thus, our results provide evidence that PLAG1 overexpression may be responsible for the frequently observed up‐regulation of IGF2 in hepatoblastoma and therefore may be implicated in the molecular pathogenesis of this childhood neoplasia.

A database of protein expression in lung cancer

We have developed a comprehensive approach to identifying molecular changes in lung cancer that includes both genomic and proteomic analyses. The related effort has produced a large amount of data pertaining to gene expression at the RNA and protein levels. As a result, we have constructed a database that contains protein expression data on lung cancer as well as other relevant data including DNA microarray derived data. A large number of proteins that are expressed in different types of lung cancer have been identified and have been correlated with the expression measures for their corresponding genes at the RNA level. The database is intended to facilitate our effort at developing novel classification schemes for lung cancer and the identification of novel markers for early diagnosis.

Ancient whole genome enrichment using baits built from modern DNA.

We report metrics from complete genome capture of nuclear DNA from extinct mammoths using biotinylated RNAs transcribed from an Asian elephant DNA extract. Enrichment of the nuclear genome ranged from 1.06- to 18.65-fold, to an apparent maximum threshold of ∼80% on-target. This projects an order of magnitude less costly complete genome sequencing from long-dead organisms, even when a reference genome is unavailable for bait design.

On-chip non-equilibrium dissociation curves and dissociation rate constants as methods to assess specificity of oligonucleotide probes.

Nucleic acid hybridization serves as backbone for many high-throughput systems for detection, expression analysis, comparative genomics and re-sequencing. Specificity of hybridization between probes and intended targets is always critical. Approaches to ensure and evaluate specificity include use of mismatch probes, obtaining dissociation curves rather than single temperature hybridizations, and comparative hybridizations. In this study, we quantify effects of mismatch type and position on intensity of hybridization signals and provide a new approach based on dissociation rate constants to evaluate specificity of hybridized signals in complex target mixtures. Using an extensive set of 18mer oligonucleotide probes on an in situ synthesized biochip platform, we demonstrate that mismatches in the center of the probe are more discriminating than mismatches toward the extremities of the probe and mismatches toward the attached end are less discriminating than those toward the loose end. The observed destabilizing effect of a mismatch type agreed in general with predictions using the nearest neighbor model. Use of a new parameter, specific dissociation temperature (Td-w, temperature of maximum specific dissociation rate constant), obtained from probe–target duplex dissociation profiles considerably improved the evaluation of specificity. These results have broad implications for hybridization data obtained from complex mixtures of nucleic acids.