CiteULike: Xavier's library [81 articles]

Analyzing patterns of microbial evolution using the mauve genome alignment system.

AE Darling — 2008-04-08T08:45:04-00:00

Methods in molecular biology (Clifton, N.J.), Vol. 396 (2007), pp. 135-152.

During the course of evolution, genomes can undergo large-scale mutation events such as rearrangement and lateral transfer. Such mutations can result in significant variations in gene order and gene content among otherwise closely related organisms. The Mauve genome alignment system can successfully identify such rearrangement and lateral transfer events in comparisons of multiple microbial genomes even under high levels of recombination. This chapter outlines the main features of Mauve and provides examples that describe how to use Mauve to conduct a rigorous multiple genome comparison and study evolutionary patterns.

Genome sequencing shows that European isolates of Francisella tularensis subspecies tularensis are almost identical to US laboratory strain Schu S4.

RR Chaudhuri — 2008-03-07T16:30:09-00:00

PLoS ONE, Vol. 2, No. 4. (2007)

BACKGROUND: Francisella tularensis causes tularaemia, a life-threatening zoonosis, and has potential as a biowarfare agent. F. tularensis subsp. tularensis, which causes the most severe form of tularaemia, is usually confined to North America. However, a handful of isolates from this subspecies was obtained in the 1980s from ticks and mites from Slovakia and Austria. Our aim was to uncover the origins of these enigmatic European isolates. METHODOLOGY/PRINCIPAL FINDINGS: We determined the complete genome sequence of FSC198, a European isolate of F. tularensis subsp. tularensis, by whole-genome shotgun sequencing and compared it to that of the North American laboratory strain Schu S4. Apparent differences between the two genomes were resolved by re-sequencing discrepant loci in both strains. We found that the genome of FSC198 is almost identical to that of Schu S4, with only eight SNPs and three VNTR differences between the two sequences. Sequencing of these loci in two other European isolates of F. tularensis subsp. tularensis confirmed that all three European isolates are also closely related to, but distinct from Schu S4. CONCLUSIONS/SIGNIFICANCE: The data presented here suggest that the Schu S4 laboratory strain is the most likely source of the European isolates of F. tularensis subsp. tularensis and indicate that anthropogenic activities, such as movement of strains or animal vectors, account for the presence of these isolates in Europe. Given the highly pathogenic nature of this subspecies, the possibility that it has become established wild in the heartland of Europe carries significant public health implications.

Complete Genomic Characterization of a Pathogenic A.II Strain of Francisella tularensis Subspecies tularensis.

SM Beckstrom-Sternberg — 2008-03-07T16:29:56-00:00

PLoS ONE, Vol. 2, No. 9. (2007)

Francisella tularensis is the causative agent of tularemia, which is a highly lethal disease from nature and potentially from a biological weapon. This species contains four recognized subspecies including the North American endemic F. tularensis subsp. tularensis (type A), whose genetic diversity is correlated with its geographic distribution including a major population subdivision referred to as A.I and A.II. The biological significance of the A.I - A.II genetic differentiation is unknown, though there are suggestive ecological and epidemiological correlations. In order to understand the differentiation at the genomic level, we have determined the complete sequence of an A.II strain (WY96-3418) and compared it to the genome of Schu S4 from the A.I population. We find that this A.II genome is 1,898,476 bp in size with 1,820 genes, 1,303 of which code for proteins. While extensive genomic variation exists between "WY96" and Schu S4, there is only one whole gene difference. This one gene difference is a hypothetical protein of unknown function. In contrast, there are numerous SNPs (3,367), small indels (1,015), IS element differences (7) and large chromosomal rearrangements (31), including both inversions and translocations. The rearrangement borders are frequently associated with IS elements, which would facilitate intragenomic recombination events. The pathogenicity island duplicated regions (DR1 and DR2) are essentially identical in WY96 but vary relative to Schu S4 at 60 nucleotide positions. Other potential virulence-associated genes (231) varied at 559 nucleotide positions, including 357 non-synonymous changes. Molecular clock estimates for the divergence time between A.I and A.II genomes for different chromosomal regions ranged from 866 to 2131 years before present. This paper is the first complete genomic characterization of a member of the A.II clade of Francisella tularensis subsp. tularensis.

Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human-pathogenic strains

Laurence Rohmer — 2007-06-06T00:28:02-00:00

Genome Biology, Vol. 8 (05 June 2007), R102.

The complete genome sequence of Francisella tularensis, the causative agent of tularemia.

P Larsson — 2005-02-11T16:10:31-00:00

Nat Genet, Vol. 37, No. 2. (February 2005), pp. 153-159.

Francisella tularensis is one of the most infectious human pathogens known. In the past, both the former Soviet Union and the US had programs to develop weapons containing the bacterium. We report the complete genome sequence of a highly virulent isolate of F. tularensis (1,892,819 bp). The sequence uncovers previously uncharacterized genes encoding type IV pili, a surface polysaccharide and iron-acquisition systems. Several virulence-associated genes were located in a putative pathogenicity island, which was duplicated in the genome. More than 10% of the putative coding sequences contained insertion-deletion or substitution mutations and seemed to be deteriorating. The genome is rich in IS elements, including IS630 Tc-1 mariner family transposons, which are not expected in a prokaryote. We used a computational method for predicting metabolic pathways and found an unexpectedly high proportion of disrupted pathways, explaining the fastidious nutritional requirements of the bacterium. The loss of biosynthetic pathways indicates that F. tularensis is an obligate host-dependent bacterium in its natural life cycle. Our results have implications for our understanding of how highly virulent human pathogens evolve and will expedite strategies to combat them.

Chromosome rearrangement and diversification of Francisella tularensis revealed by the type B (OSU18) genome sequence.

JF Petrosino — 2008-03-07T16:29:28-00:00

J Bacteriol, Vol. 188, No. 19. (October 2006), pp. 6977-6985.

The gamma-proteobacterium Francisella tularensis is one of the most infectious human pathogens, and the highly virulent organism F. tularensis subsp. tularensis (type A) and less virulent organism F. tularensis subsp. holarctica (type B) are most commonly associated with significant disease in humans and animals. Here we report the complete genome sequence and annotation for a low-passage type B strain (OSU18) isolated from a dead beaver found near Red Rock, Okla., in 1978. A comparison of the F. tularensis subsp. holarctica sequence with that of F. tularensis subsp. tularensis strain Schu4 (P. Larsson et al., Nat. Genet. 37:153-159, 2005) highlighted genetic differences that may underlie different pathogenicity phenotypes and the evolutionary relationship between type A and type B strains. Despite extensive DNA sequence identity, the most significant difference between type A and type B isolates is the striking amount of genomic rearrangement that exists between the strains. All but two rearrangements can be attributed to homologous recombination occurring between two prominent insertion elements, ISFtu1 and ISFtu2. Numerous pseudogenes have been found in the genomes and are likely contributors to the difference in virulence between the strains. In contrast, no rearrangements have been observed between the OSU18 genome and the genome of the type B live vaccine strain (LVS), and only 448 polymorphisms have been found within non-transposase-coding sequences whose homologs are intact in OSU18. Nonconservative differences between the two strains likely include the LVS attenuating mutation(s).

Dissemination in Japanese hospitals of strains of Staphylococcus aureus heterogeneously resistant to vancomycin.

K Hiramatsu — 2008-03-06T11:11:16-00:00

Lancet, Vol. 350, No. 9092. (6 December 1997), pp. 1670-1673.

BACKGROUND: Since the discovery of the vancomycin-resistant Staphylococcus aureus (VRSA) strain Mu50 (minimum inhibitory concentration [MIC] 8 mg/L), there has been concern about the potential spread of such strains throughout Japanese hospitals. Two important questions need to be answered: (1) what is the prevalence of VRSA, and (2) by what mechanism does vancomycin resistance occur. METHODS: The vancomycin susceptibilities of three methicillin-resistant S aureus (MRSA) strains (Mu50, Mu3, and H1) and the methicillin-susceptible S aureus type strain FDA209P were compared by MIC determinations and population analysis. Mu3 (MIC 3 mg/L) was isolated from the sputum of a patient with pneumonia after surgery who had failed vancomycin therapy. H1 (MIC 2 mg/L), which is a representative vancomycin-susceptible MRSA strain, was isolated from a patient with pneumonia who responded favourably to vancomycin therapy. Subclones of Mu3 with increased resistance against vancomycin were selected with serial concentrations of vancomycin and their MICs were determined. The prevalence of VRSA and Mu3-like strains in Japanese hospitals was estimated by population analysis from 1149 clinical MRSA isolates obtained from 203 hospitals throughout Japan. The genetic traits of the Mu3 and Mu50 strains were compared with clonotypes of MRSA from around the world. FINDINGS: Mu3 and Mu50 had an identical pulsed-field gel electrophoresis banding pattern. When grown in a drug-free medium, Mu3 produced subpopulation of cells with varying degrees of vancomycin resistance, thus demonstrating natural heterogeneity, or variability, in susceptibility to vancomycin. In the presence of vancomycin, Mu3 produced subclones with resistance roughly proportional to the concentrations of vancomycin used. Selection of Mu3 with 8 mg/L or more of vancomycin gave rise to subclones with vancomycin resistance equal to that of Mu50 (MIC 8 mg/L) at a frequency of 1/1,000,000. During screening of Japanese MRSA strains, no strain of VRSA additional to Mu50 was found. The prevalence of MRSA isolates heterogeneously resistant to vancomycin was 20% in Juntendo University Hospital, 9.3% in the other seven university hospitals, and 1.3% in non-university hospitals or clinics. INTERPRETATION: Heterogeneously resistant VRSA is a preliminary stage that allows development into VRSA upon exposure to vancomycin. Heterogeneously resistant VRSA was found in hospitals throughout Japan. This finding could explain, at least partly, the frequent therapeutic failure of MRSA infection with vancomycin in Japan.

Escherichia coli K-12: a cooperatively developed annotation snapshot2005

Monica Riley — 2006-01-25T20:51:14-00:00

Nucleic Acids Research, Vol. 34, No. 1. (2006), pp. 1-9.

Genome dynamics and diversity of Shigella species, the etiologic agents of bacillary dysentery

Fan Yang — 2005-11-23T12:22:32-00:00

Nucleic Acids Research, Vol. 33, No. 19. (2005), pp. 6445-6458.

Genome sequence of Staphylococcus aureus strain Newman and comparative analysis of staphylococcal genomes: polymorphism and evolution of two major pathogenicity islands.

T Baba — 2008-03-04T15:55:28-00:00

J Bacteriol, Vol. 190, No. 1. (January 2008), pp. 300-310.

Strains of Staphylococcus aureus, an important human pathogen, display up to 20% variability in their genome sequence, and most sequence information is available for human clinical isolates that have not been subjected to genetic analysis of virulence attributes. S. aureus strain Newman, which was also isolated from a human infection, displays robust virulence properties in animal models of disease and has already been extensively analyzed for its molecular traits of staphylococcal pathogenesis. We report here the complete genome sequence of S. aureus Newman, which carries four integrated prophages, as well as two large pathogenicity islands. In agreement with the view that S. aureus Newman prophages contribute important properties to pathogenesis, fewer virulence factors are found outside of the prophages than for the highly virulent strain MW2. The absence of drug resistance genes reflects the general antibiotic-susceptible phenotype of S. aureus Newman. Phylogenetic analyses reveal clonal relationships between the staphylococcal strains Newman, COL, NCTC8325, and USA300 and a greater evolutionary distance to strains MRSA252, MW2, MSSA476, N315, Mu50, JH1, JH9, and RF122. However, polymorphism analysis of two large pathogenicity islands distributed among these strains shows that the two islands were acquired independently from the evolutionary pathway of the chromosomal backbones of staphylococcal genomes. Prophages and pathogenicity islands play central roles in S. aureus virulence and evolution.

Genome and virulence determinants of high virulence community-acquired MRSA.

T Baba — 2008-03-04T15:55:21-00:00

Lancet, Vol. 359, No. 9320. (25 May 2002), pp. 1819-1827.

BACKGROUND: A new type of meticillin-resistant Staphylococcus aureus (MRSA), designated community-acquired MRSA, is becoming increasingly noticeable in the community, some strains of which cause fatal infections in otherwise healthy individuals. By contrast with hospital-acquired MRSA, community-acquired MRSA is more susceptible to non b-lactam antibiotics. We investigated the high virulence potential of certain strains of this bacterium. METHODS: We ascertained the whole genome sequence of MW2, a strain of community-acquired MRSA, by shotgun cloning and sequencing. MW2 caused fatal septicaemia and septic arthritis in a 16-month-old girl in North Dakota, USA, in 1998. The genome of this strain was compared with those of hospital-acquired MRSA strains, including N315 and Mu50. FINDINGS: Meticillin resistance gene (mecA) in MW2 was carried by a novel allelic form (type IVa) of staphylococcal cassette chromosome mec (SCCmec), by contrast with type II in N315 and Mu50. Type IVa SCCmec did not carry any of the multiple antibiotic resistance genes reported in type II SCCmec. By contrast, 19 additional virulence genes were recorded in the MW2 genome. All but two of these virulence genes were noted in four of the seven genomic islands of MW2. INTERPRETATION: MW2 carried a range of virulence and resistance genes that was distinct from those displayed on the chromosomes of extant S aureus strains. Most genes were carried by specific allelic forms of genomic islands in the MW2 chromosome. The combination of allelic forms of genomic islands is the genetic basis that determines the pathogenicity of medically important phenotypes of S aureus, including those of community-acquired MRSA strains.

Whole genome sequencing of meticillin-resistant Staphylococcus aureus.

M Kuroda — 2008-03-04T15:55:14-00:00

Lancet, Vol. 357, No. 9264. (21 April 2001), pp. 1225-1240.

BACKGROUND: Staphylococcus aureus is one of the major causes of community-acquired and hospital-acquired infections. It produces numerous toxins including superantigens that cause unique disease entities such as toxic-shock syndrome and staphylococcal scarlet fever, and has acquired resistance to practically all antibiotics. Whole genome analysis is a necessary step towards future development of countermeasures against this organism. METHODS: Whole genome sequences of two related S aureus strains (N315 and Mu50) were determined by shot-gun random sequencing. N315 is a meticillin-resistant S aureus (MRSA) strain isolated in 1982, and Mu50 is an MRSA strain with vancomycin resistance isolated in 1997. The open reading frames were identified by use of GAMBLER and GLIMMER programs, and annotation of each was done with a BLAST homology search, motif analysis, and protein localisation prediction. FINDINGS: The Staphylococcus genome was composed of a complex mixture of genes, many of which seem to have been acquired by lateral gene transfer. Most of the antibiotic resistance genes were carried either by plasmids or by mobile genetic elements including a unique resistance island. Three classes of new pathogenicity islands were identified in the genome: a toxic-shock-syndrome toxin island family, exotoxin islands, and enterotoxin islands. In the latter two pathogenicity islands, clusters of exotoxin and enterotoxin genes were found closely linked with other gene clusters encoding putative pathogenic factors. The analysis also identified 70 candidates for new virulence factors. INTERPRETATION: The remarkable ability of S aureus to acquire useful genes from various organisms was revealed through the observation of genome complexity and evidence of lateral gene transfer. Repeated duplication of genes encoding superantigens explains why S aureus is capable of infecting humans of diverse genetic backgrounds, eliciting severe immune reactions. Investigation of many newly identified gene products, including the 70 putative virulence factors, will greatly improve our understanding of the biology of staphylococci and the processes of infectious diseases caused by S aureus.

Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus

Sarah Highlander — 2007-11-07T21:43:29-00:00

BMC Microbiology, Vol. 7 (06 November 2007), 99.

Complete genome sequence of USA300, an epidemic clone of community-acquired meticillin-resistant Staphylococcus aureus.

BA Diep — 2008-03-04T15:55:01-00:00

Lancet, Vol. 367, No. 9512. (4 March 2006), pp. 731-739.

BACKGROUND: USA300, a clone of meticillin-resistant Staphylococcus aureus, is a major source of community-acquired infections in the USA, Canada, and Europe. Our aim was to sequence its genome and compare it with those of other strains of S aureus to try to identify genes responsible for its distinctive epidemiological and virulence properties. METHODS: We ascertained the genome sequence of FPR3757, a multidrug resistant USA300 strain, by random shotgun sequencing, then compared it with the sequences of ten other staphylococcal strains. FINDINGS: Compared with closely related S aureus, we noted that almost all of the unique genes in USA300 clustered in novel allotypes of mobile genetic elements. Some of the unique genes are involved in pathogenesis, including Panton-Valentine leucocidin and molecular variants of enterotoxin Q and K. The most striking feature of the USA300 genome is the horizontal acquisition of a novel mobile genetic element that encodes an arginine deiminase pathway and an oligopeptide permease system that could contribute to growth and survival of USA300. We did not detect this element, termed arginine catabolic mobile element (ACME), in other S aureus strains. We noted a high prevalence of ACME in S epidermidis, suggesting not only that ACME transfers into USA300 from S epidermidis, but also that this element confers a selective advantage to this ubiquitous commensal of the human skin. INTERPRETATION: USA300 has acquired mobile genetic elements that encode resistance and virulence determinants that could enhance fitness and pathogenicity.

Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance.

MT Holden — 2008-03-04T15:54:52-00:00

Proc Natl Acad Sci U S A, Vol. 101, No. 26. (29 June 2004), pp. 9786-9791.

Staphylococcus aureus is an important nosocomial and community-acquired pathogen. Its genetic plasticity has facilitated the evolution of many virulent and drug-resistant strains, presenting a major and constantly changing clinical challenge. We sequenced the approximately 2.8-Mbp genomes of two disease-causing S. aureus strains isolated from distinct clinical settings: a recent hospital-acquired representative of the epidemic methicillin-resistant S. aureus EMRSA-16 clone (MRSA252), a clinically important and globally prevalent lineage; and a representative of an invasive community-acquired methicillin-susceptible S. aureus clone (MSSA476). A comparative-genomics approach was used to explore the mechanisms of evolution of clinically important S. aureus genomes and to identify regions affecting virulence and drug resistance. The genome sequences of MRSA252 and MSSA476 have a well conserved core region but differ markedly in their accessory genetic elements. MRSA252 is the most genetically diverse S. aureus strain sequenced to date: approximately 6% of the genome is novel compared with other published genomes, and it contains several unique genetic elements. MSSA476 is methicillin-susceptible, but it contains a novel Staphylococcal chromosomal cassette (SCC) mec-like element (designated SCC(476)), which is integrated at the same site on the chromosome as SCCmec elements in MRSA strains but encodes a putative fusidic acid resistance protein. The crucial role that accessory elements play in the rapid evolution of S. aureus is clearly illustrated by comparing the MSSA476 genome with that of an extremely closely related MRSA community-acquired strain; the differential distribution of large mobile elements carrying virulence and drug-resistance determinants may be responsible for the clinically important phenotypic differences in these strains.

Insights on evolution of virulence and resistance from the complete genome analysis of an early methicillin-resistant Staphylococcus aureus strain and a biofilm-producing methicillin-resistant Staphylococcus epidermidis strain.

SR Gill — 2008-03-04T15:53:56-00:00

J Bacteriol, Vol. 187, No. 7. (April 2005), pp. 2426-2438.

Staphylococcus aureus is an opportunistic pathogen and the major causative agent of numerous hospital- and community-acquired infections. Staphylococcus epidermidis has emerged as a causative agent of infections often associated with implanted medical devices. We have sequenced the approximately 2.8-Mb genome of S. aureus COL, an early methicillin-resistant isolate, and the approximately 2.6-Mb genome of S. epidermidis RP62a, a methicillin-resistant biofilm isolate. Comparative analysis of these and other staphylococcal genomes was used to explore the evolution of virulence and resistance between these two species. The S. aureus and S. epidermidis genomes are syntenic throughout their lengths and share a core set of 1,681 open reading frames. Genome islands in nonsyntenic regions are the primary source of variations in pathogenicity and resistance. Gene transfer between staphylococci and low-GC-content gram-positive bacteria appears to have shaped their virulence and resistance profiles. Integrated plasmids in S. epidermidis carry genes encoding resistance to cadmium and species-specific LPXTG surface proteins. A novel genome island encodes multiple phenol-soluble modulins, a potential S. epidermidis virulence factor. S. epidermidis contains the cap operon, encoding the polyglutamate capsule, a major virulence factor in Bacillus anthracis. Additional phenotypic differences are likely the result of single nucleotide polymorphisms, which are most numerous in cell envelope proteins. Overall differences in pathogenicity can be attributed to genome islands in S. aureus which encode enterotoxins, exotoxins, leukocidins, and leukotoxins not found in S. epidermidis.

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

MM Mwangi — 2007-07-03T11:46:27-00:00

Proc Natl Acad Sci U S A, Vol. 104, No. 22. (29 May 2007), pp. 9451-9456.

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.

Complete genome sequence of Shigella flexneri5b and comparison with Shigella flexneri2a

Huan Nie — 2006-07-07T04:05:00-00:00

BMC Genomics, Vol. 7 (06 July 2006), 173.

Genome sequence of Shigella flexneri 2a: insights into pathogenicity through comparison with genomes of Escherichia coli K12 and O157.

Q Jin — 2008-03-04T15:43:19-00:00

Nucleic Acids Res, Vol. 30, No. 20. (15 October 2002), pp. 4432-4441.

We have sequenced the genome of Shigella flexneri serotype 2a, the most prevalent species and serotype that causes bacillary dysentery or shigellosis in man. The whole genome is composed of a 4 607 203 bp chromosome and a 221 618 bp virulence plasmid, designated pCP301. While the plasmid shows minor divergence from that sequenced in serotype 5a, striking characteristics of the chromosome have been revealed. The S.flexneri chromosome has, astonishingly, 314 IS elements, more than 7-fold over those possessed by its close relatives, the non-pathogenic K12 strain and enterohemorrhagic O157:H7 strain of Escherichia coli. There are 13 translocations and inversions compared with the E.coli sequences, all involve a segment larger than 5 kb, and most are associated with deletions or acquired DNA sequences, of which several are likely to be bacteriophage-transmitted pathogenicity islands. Furthermore, S.flexneri, resembling another human-restricted enteric pathogen, Salmonella typhi, also has hundreds of pseudogenes compared with the E.coli strains. All of these could be subjected to investigations towards novel preventative and treatment strategies against shigellosis.

Complete genome sequence and comparative genomics of Shigella flexneri serotype 2a strain 2457T.

J Wei — 2008-03-04T15:43:11-00:00

Infect Immun, Vol. 71, No. 5. (May 2003), pp. 2775-2786.

We determined the complete genome sequence of Shigella flexneri serotype 2a strain 2457T (4,599,354 bp). Shigella species cause >1 million deaths per year from dysentery and diarrhea and have a lifestyle that is markedly different from those of closely related bacteria, including Escherichia coli. The genome exhibits the backbone and island mosaic structure of E. coli pathogens, albeit with much less horizontally transferred DNA and lacking 357 genes present in E. coli. The strain is distinctive in its large complement of insertion sequences, with several genomic rearrangements mediated by insertion sequences, 12 cryptic prophages, 372 pseudogenes, and 195 S. flexneri-specific genes. The 2457T genome was also compared with that of a recently sequenced S. flexneri 2a strain, 301. Our data are consistent with Shigella being phylogenetically indistinguishable from E. coli. The S. flexneri-specific regions contain many genes that could encode proteins with roles in virulence. Analysis of these will reveal the genetic basis for aspects of this pathogenic organism's distinctive lifestyle that have yet to be explained.

Genome dynamics and diversity of Shigella species, the etiologic agents of bacillary dysentery.

F Yang — 2007-09-22T14:17:17-00:00

Nucleic Acids Res, Vol. 33, No. 19. (2005), pp. 6445-6458.

The Shigella bacteria cause bacillary dysentery, which remains a significant threat to public health. The genus status and species classification appear no longer valid, as compelling evidence indicates that Shigella, as well as enteroinvasive Escherichia coli, are derived from multiple origins of E.coli and form a single pathovar. Nevertheless, Shigella dysenteriae serotype 1 causes deadly epidemics but Shigella boydii is restricted to the Indian subcontinent, while Shigella flexneri and Shigella sonnei are prevalent in developing and developed countries respectively. To begin to explain these distinctive epidemiological and pathological features at the genome level, we have carried out comparative genomics on four representative strains. Each of the Shigella genomes includes a virulence plasmid that encodes conserved primary virulence determinants. The Shigella chromosomes share most of their genes with that of E.coli K12 strain MG1655, but each has over 200 pseudogenes, 300 approximately 700 copies of insertion sequence (IS) elements, and numerous deletions, insertions, translocations and inversions. There is extensive diversity of putative virulence genes, mostly acquired via bacteriophage-mediated lateral gene transfer. Hence, via convergent evolution involving gain and loss of functions, through bacteriophage-mediated gene acquisition, IS-mediated DNA rearrangements and formation of pseudogenes, the Shigella spp. became highly specific human pathogens with variable epidemiological and pathological features.

Extensive mosaic structure revealed by the complete genome sequence of uropathogenic Escherichia coli.

RA Welch — 2005-10-27T20:40:51-00:00

Proc Natl Acad Sci U S A, Vol. 99, No. 26. (24 December 2002), pp. 17020-17024.

We present the complete genome sequence of uropathogenic Escherichia coli, strain CFT073. A three-way genome comparison of the CFT073, enterohemorrhagic E. coli EDL933, and laboratory strain MG1655 reveals that, amazingly, only 39.2% of their combined (nonredundant) set of proteins actually are common to all three strains. The pathogen genomes are as different from each other as each pathogen is from the benign strain. The difference in disease potential between O157:H7 and CFT073 is reflected in the absence of genes for type III secretion system or phage- and plasmid-encoded toxins found in some classes of diarrheagenic E. coli. The CFT073 genome is particularly rich in genes that encode potential fimbrial adhesins, autotransporters, iron-sequestration systems, and phase-switch recombinases. Striking differences exist between the large pathogenicity islands of CFT073 and two other well-studied uropathogenic E. coli strains, J96 and 536. Comparisons indicate that extraintestinal pathogenic E. coli arose independently from multiple clonal lineages. The different E. coli pathotypes have maintained a remarkable synteny of common, vertically evolved genes, whereas many islands interrupting this common backbone have been acquired by different horizontal transfer events in each strain.

How to become a uropathogen: comparative genomic analysis of extraintestinal pathogenic Escherichia coli strains.

E Brzuszkiewicz — 2008-03-03T18:21:45-00:00

Proc Natl Acad Sci U S A, Vol. 103, No. 34. (22 August 2006), pp. 12879-12884.

Uropathogenic Escherichia coli (UPEC) strain 536 (O6:K15:H31) is one of the model organisms of extraintestinal pathogenic E. coli (ExPEC). To analyze this strain's genetic basis of urovirulence, we sequenced the entire genome and compared the data with the genome sequence of UPEC strain CFT073 (O6:K2:H1) and to the available genomes of nonpathogenic E. coli strain MG1655 (K-12) and enterohemorrhagic E. coli. The genome of strain 536 is approximately 292 kb smaller than that of strain CFT073. Genomic differences between both UPEC are mainly restricted to large pathogenicity islands, parts of which are unique to strain 536 or CFT073. Genome comparison underlines that repeated insertions and deletions in certain parts of the genome contribute to genome evolution. Furthermore, 427 and 432 genes are only present in strain 536 or in both UPEC, respectively. The majority of the latter genes is encoded within smaller horizontally acquired DNA regions scattered all over the genome. Several of these genes are involved in increasing the pathogens' fitness and adaptability. Analysis of virulence-associated traits expressed in the two UPEC O6 strains, together with genome comparison, demonstrate the marked genetic and phenotypic variability among UPEC. The ability to accumulate and express a variety of virulence-associated genes distinguishes ExPEC from many commensals and forms the basis for the individual virulence potential of ExPEC. Accordingly, instead of a common virulence mechanism, different ways exist among ExPEC to cause disease.

The genome sequence of avian pathogenic Escherichia coli strain O1:K1:H7 shares strong similarities with human extraintestinal pathogenic E. coli genomes.

TJ Johnson — 2008-03-03T18:21:30-00:00

J Bacteriol, Vol. 189, No. 8. (April 2007), pp. 3228-3236.

Escherichia coli strains that cause disease outside the intestine are known as extraintestinal pathogenic E. coli (ExPEC) and include human uropathogenic E. coli (UPEC) and avian pathogenic E. coli (APEC). Regardless of host of origin, ExPEC strains share many traits. It has been suggested that these commonalities may enable APEC to cause disease in humans. Here, we begin to test the hypothesis that certain APEC strains possess potential to cause human urinary tract infection through virulence genotyping of 1,000 APEC and UPEC strains, generation of the first complete genomic sequence of an APEC (APEC O1:K1:H7) strain, and comparison of this genome to all available human ExPEC genomic sequences. The genomes of APEC O1 and three human UPEC strains were found to be remarkably similar, with only 4.5% of APEC O1's genome not found in other sequenced ExPEC genomes. Also, use of multilocus sequence typing showed that some of the sequenced human ExPEC strains were more like APEC O1 than other human ExPEC strains. This work provides evidence that at least some human and avian ExPEC strains are highly similar to one another, and it supports the possibility that a food-borne link between some APEC and UPEC strains exists. Future studies are necessary to assess the ability of APEC to overcome the hurdles necessary for such a food-borne transmission, and epidemiological studies are required to confirm that such a phenomenon actually occurs.

Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12.

T Hayashi — 2006-07-18T18:54:29-00:00

DNA Res, Vol. 8, No. 1. (28 February 2001), pp. 11-22.

Escherichia coli O157:H7 is a major food-borne infectious pathogen that causes diarrhea, hemorrhagic colitis, and hemolytic uremic syndrome. Here we report the complete chromosome sequence of an O157:H7 strain isolated from the Sakai outbreak, and the results of genomic comparison with a benign laboratory strain, K-12 MG1655. The chromosome is 5.5 Mb in size, 859 Kb larger than that of K-12. We identified a 4.1-Mb sequence highly conserved between the two strains, which may represent the fundamental backbone of the E. coli chromosome. The remaining 1.4-Mb sequence comprises of O157:H7-specific sequences, most of which are horizontally transferred foreign DNAs. The predominant roles of bacteriophages in the emergence of O157:H7 is evident by the presence of 24 prophages and prophage-like elements that occupy more than half of the O157:H7-specific sequences. The O157:H7 chromosome encodes 1632 proteins and 20 tRNAs that are not present in K-12. Among these, at least 131 proteins are assumed to have virulence-related functions. Genome-wide codon usage analysis suggested that the O157:H7-specific tRNAs are involved in the efficient expression of the strain-specific genes. A complete set of the genes specific to O157:H7 presented here sheds new insight into the pathogenicity and the physiology of O157:H7, and will open a way to fully understand the molecular mechanisms underlying the O157:H7 infection.

Genome sequence of enterohaemorrhagic Escherichia coli O157:H7

Nicole Perna — 2006-05-31T15:04:05-00:00

Nature, Vol. 409, No. 6819. (25 January 2001), pp. 529-533.

Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: A comparative genomics approach.

Swaine L Chen — 2006-04-11T19:33:43-00:00

Proc Natl Acad Sci U S A (3 April 2006)

Escherichia coli is a model laboratory bacterium, a species that is widely distributed in the environment, as well as a mutualist and pathogen in its human hosts. As such, E. coli represents an attractive organism to study how environment impacts microbial genome structure and function. Uropathogenic E. coli (UPEC) must adapt to life in several microbial communities in the human body, and has a complex life cycle in the bladder when it causes acute or recurrent urinary tract infection (UTI). Several studies designed to identify virulence factors have focused on genes that are uniquely represented in UPEC strains, whereas the role of genes that are common to all E. coli has received much less attention. Here we describe the complete 5,065,741-bp genome sequence of a UPEC strain recovered from a patient with an acute bladder infection and compare it with six other finished E. coli genome sequences. We searched 3,470 ortholog sets for genes that are under positive selection only in UPEC strains. Our maximum likelihood-based analysis yielded 29 genes involved in various aspects of cell surface structure, DNA metabolism, nutrient acquisition, and UTI. These results were validated by resequencing a subset of the 29 genes in a panel of 50 urinary, periurethral, and rectal E. coli isolates from patients with UTI. These studies outline a computational approach that may be broadly applicable for studying strain-specific adaptation and pathogenesis in other bacteria.

The complete genome sequence of Escherichia coli K-12.

FR Blattner — 2005-07-21T08:30:29-00:00

Science, Vol. 277, No. 5331. (5 September 1997), pp. 1453-1474.

The 4,639,221-base pair sequence of Escherichia coli K-12 is presented. Of 4288 protein-coding genes annotated, 38 percent have no attributed function. Comparison with five other sequenced microbes reveals ubiquitous as well as narrowly distributed gene families; many families of similar genes within E. coli are also evident. The largest family of paralogous proteins contains 80 ABC transporters. The genome as a whole is strikingly organized with respect to the local direction of replication; guanines, oligonucleotides possibly related to replication and recombination, and most genes are so oriented. The genome also contains insertion sequence (IS) elements, phage remnants, and many other patches of unusual composition indicating genome plasticity through horizontal transfer.

Genome sequence of a serotype M3 strain of group A Streptococcus: phage-encoded toxins, the high-virulence phenotype, and clone emergence.

SB Beres — 2008-03-03T18:04:37-00:00

Proc Natl Acad Sci U S A, Vol. 99, No. 15. (23 July 2002), pp. 10078-10083.

Genome sequences are available for many bacterial strains, but there has been little progress in using these data to understand the molecular basis of pathogen emergence and differences in strain virulence. Serotype M3 strains of group A Streptococcus (GAS) are a common cause of severe invasive infections with unusually high rates of morbidity and mortality. To gain insight into the molecular basis of this high-virulence phenotype, we sequenced the genome of strain MGAS315, an organism isolated from a patient with streptococcal toxic shock syndrome. The genome is composed of 1,900,521 bp, and it shares approximately 1.7 Mb of related genetic material with genomes of serotype M1 and M18 strains. Phage-like elements account for the great majority of variation in gene content relative to the sequenced M1 and M18 strains. Recombination produces chimeric phages and strains with previously uncharacterized arrays of virulence factor genes. Strain MGAS315 has phage genes that encode proteins likely to contribute to pathogenesis, such as streptococcal pyrogenic exotoxin A (SpeA) and SpeK, streptococcal superantigen (SSA), and a previously uncharacterized phospholipase A(2) (designated Sla). Infected humans had anti-SpeK, -SSA, and -Sla antibodies, indicating that these GAS proteins are made in vivo. SpeK and SSA were pyrogenic and toxic for rabbits. Serotype M3 strains with the phage-encoded speK and sla genes increased dramatically in frequency late in the 20th century, commensurate with the rise in invasive disease caused by M3 organisms. Taken together, the results show that phage-mediated recombination has played a critical role in the emergence of a new, unusually virulent clone of serotype M3 GAS.

Molecular genetic anatomy of inter- and intraserotype variation in the human bacterial pathogen group A Streptococcus

Stephen Beres — 2006-12-18T17:25:03-00:00

PNAS, Vol. 103, No. 18. (2 May 2006), pp. 7059-7064.

In recent years we have studied the relationship between strain genotypes and patient phenotypes in group A Streptococcus (GAS), a model human bacterial pathogen that causes extensive morbidity and mortality worldwide. We have concentrated our efforts on serotype M3 organisms because these strains are common causes of pharyngeal and invasive infections, produce unusually severe invasive infections, and can exhibit epidemic behavior. Our studies have been hindered by the lack of genome-scale phylogenies of multiple GAS strains and whole-genome sequences of multiple serotype M3 strains recovered from individuals with defined clinical phenotypes. To remove some of these impediments, we sequenced to closure the genome of four additional GAS strains and conducted comparative genomic resequencing of 12 contemporary serotype M3 strains representing distinct genotypes and phenotypes. Serotype M3 strains are a single phylogenetic lineage. Strains from asymptomatic throat carriers were significantly less virulent for mice than sterile-site isolates and evolved to a less virulent phenotype by multiple genetic pathways. Strain persistence or extinction between epidemics was strongly associated with presence or absence, respectively, of the prophage encoding streptococcal pyrogenic exotoxin A. A serotype M3 clone significantly underrepresented among necrotizing fasciitis cases has a unique frameshift mutation that truncates MtsR, a transcriptional regulator controlling expression of genes encoding iron-acquisition proteins. Expression microarray analysis of this clone confirmed significant alteration in expression of genes encoding iron metabolism proteins. Our analysis provided unprecedented detail about the molecular anatomy of bacterial strain genotype-patient phenotype relationships. 10.1073/pnas.0510279103

Complete genome of acute rheumatic fever-associated serotype M5 Streptococcus pyogenes strain manfredo.

MT Holden — 2008-03-03T18:04:11-00:00

J Bacteriol, Vol. 189, No. 4. (February 2007), pp. 1473-1477.

Comparisons of the 1.84-Mb genome of serotype M5 Streptococcus pyogenes strain Manfredo with previously sequenced genomes emphasized the role of prophages in diversification of S. pyogenes and the close relationship between strain Manfredo and MGAS8232, another acute rheumatic fever-associated strain.

Progress toward characterization of the group A Streptococcus metagenome: complete genome sequence of a macrolide-resistant serotype M6 strain.

DJ Banks — 2008-03-03T18:03:57-00:00

J Infect Dis, Vol. 190, No. 4. (15 August 2004), pp. 727-738.

We describe the genome sequence of a macrolide-resistant strain (MGAS10394) of serotype M6 group A Streptococcus (GAS). The genome is 1,900,156 bp in length, and 8 prophage-like elements or remnants compose 12.4% of the chromosome. A 8.3-kb prophage remnant encodes the SpeA4 variant of streptococcal pyrogenic exotoxin A. The genome of strain MGAS10394 contains a chimeric genetic element composed of prophage genes and a transposon encoding the mefA gene conferring macrolide resistance. This chimeric element also has a gene encoding a novel surface-exposed protein (designated "R6 protein"), with an LPKTG cell-anchor motif located at the carboxyterminus. Surface expression of this protein was confirmed by flow cytometry. Humans with GAS pharyngitis caused by serotype M6 strains had antibody against the R6 protein present in convalescent, but not acute, serum samples. Our studies add to the theme that GAS prophage-encoded extracellular proteins contribute to host-pathogen interactions in a strain-specific fashion.

Genome sequence of an M3 strain of Streptococcus pyogenes reveals a large-scale genomic rearrangement in invasive strains and new insights into phage evolution.

I Nakagawa — 2008-03-03T18:03:41-00:00

Genome Res, Vol. 13, No. 6A. (June 2003), pp. 1042-1055.

Group Astreptococcus (GAS) is a gram-positive bacterial pathogen that causes various suppurative infections and nonsuppurative sequelae. Since the late 1980s, streptococcal toxic-shock like syndrome (STSS) and severe invasive GAS infections have been reported globally. Here we sequenced the genome of serotype M3 strain SSI-1, isolated from an STSS patient in Japan, and compared it with those of other GAS strains. The SSI-1 genome is composed of 1,884,275 bp, and 1.7 Mb of the sequence is highly conserved relative to strain SF370 (serotype M1) and MGAS8232 (serotype M18), and almost completely conserved relative to strain MGAS315 (serotype M3). However, a large genomic rearrangement has been shown to occur across the replication axis between the homologous rrn-comX1 regions and between two prophage-coding regions across the replication axis. Atotal of 1 Mb of chromosomal DNA is inverted across the replication axis. Interestingly, the recombinations between the prophage regions are within the phage genes, and the genes encoding superantigens and mitogenic factors are interchanged between two prophages. This genomic rearrangement occurs in 65% of clinical isolates (64/94) collected after 1990, whereas it is found in only 25% of clinical isolates (7/28) collected before 1985. These observations indicate that streptococcal phages represent important plasticity regions in the GAS chromosome where recombination between homologous phage genes can occur and result not only in new phage derivatives, but also in large chromosomal rearrangements.

Genome sequence of a serotype M28 strain of group a streptococcus: potential new insights into puerperal sepsis and bacterial disease specificity.

NM Green — 2008-03-03T18:03:28-00:00

J Infect Dis, Vol. 192, No. 5. (1 September 2005), pp. 760-770.

Puerperal sepsis, a major cause of death of young women in Europe in the 1800s, was due predominantly to the gram-positive pathogen group A Streptococcus. Studies conducted during past decades have shown that serotype M28 strains are the major group A Streptococcus organisms responsible for many of these infections. To begin to increase our understanding of their enrichment in puerperal sepsis, we sequenced the genome of a genetically representative strain. This strain has genes encoding a novel array of prophage virulence factors, cell-surface proteins, and other molecules likely to contribute to host-pathogen interactions. Importantly, genes for 7 inferred extracellular proteins are encoded by a 37.4-kb foreign DNA element that is shared with group B Streptococcus and is present in all serotype M28 strains. Proteins encoded by the 37.4-kb element were expressed extracellularly and in human infections. Acquisition of foreign genes has helped create a disease-specialist clone of this pathogen.

Evolutionary origin and emergence of a highly successful clone of serotype M1 group a Streptococcus involved multiple horizontal gene transfer events.

P Sumby — 2008-03-03T18:03:14-00:00

J Infect Dis, Vol. 192, No. 5. (1 September 2005), pp. 771-782.

To better understand the molecular events involved in the origin of new pathogenic bacteria, we studied the evolution of a highly virulent clone of serotype M1 group A Streptococcus (GAS). Genomic, DNA-DNA microarray, and single-nucleotide polymorphism analyses indicated that this clone evolved through a series of horizontal gene transfer events that involved (1) the acquisition of prophages encoding streptococcal pyrogenic exotoxin A and extracellular DNases and (2) the reciprocal recombination of a 36-kb chromosomal region encoding the extracellular toxins NAD+-glycohydrolase (NADase) and streptolysin O (SLO). These gene transfer events were associated with significantly increased production of SLO and NADase. Virtual identity in the 36-kb region present in contemporary serotype M1 and M12 isolates suggests that a serotype M12 strain served as the donor of this region. Multiple horizontal gene transfer events were a crucial factor in the evolutionary origin and emergence of a very abundant contemporary clone of serotype M1 GAS.

Genome sequence and comparative microarray analysis of serotype M18 group A Streptococcus strains associated with acute rheumatic fever outbreaks.

JC Smoot — 2008-03-03T18:02:26-00:00

Proc Natl Acad Sci U S A, Vol. 99, No. 7. (2 April 2002), pp. 4668-4673.

Acute rheumatic fever (ARF), a sequelae of group A Streptococcus (GAS) infection, is the most common cause of preventable childhood heart disease worldwide. The molecular basis of ARF and the subsequent rheumatic heart disease are poorly understood. Serotype M18 GAS strains have been associated for decades with ARF outbreaks in the U.S. As a first step toward gaining new insight into ARF pathogenesis, we sequenced the genome of strain MGAS8232, a serotype M18 organism isolated from a patient with ARF. The genome is a circular chromosome of 1,895,017 bp, and it shares 1.7 Mb of closely related genetic material with strain SF370 (a sequenced serotype M1 strain). Strain MGAS8232 has 178 ORFs absent in SF370. Phages, phage-like elements, and insertion sequences are the major sources of variation between the genomes. The genomes of strain MGAS8232 and SF370 encode many of the same proven or putative virulence factors. Importantly, strain MGAS8232 has genes encoding many additional secreted proteins involved in human-GAS interactions, including streptococcal pyrogenic exotoxin A (scarlet fever toxin) and two uncharacterized pyrogenic exotoxin homologues, all phage-associated. DNA microarray analysis of 36 serotype M18 strains from diverse localities showed that most regions of variation were phages or phage-like elements. Two epidemics of ARF occurring 12 years apart in Salt Lake City, UT, were caused by serotype M18 strains that were genetically identical, or nearly so. Our analysis provides a critical foundation for accelerated research into ARF pathogenesis and a molecular framework to study the plasticity of GAS genomes.

Complete genome sequence of an M1 strain of Streptococcus pyogenes.

JJ Ferretti — 2008-03-03T18:01:59-00:00

Proc Natl Acad Sci U S A, Vol. 98, No. 8. (10 April 2001), pp. 4658-4663.

The 1,852,442-bp sequence of an M1 strain of Streptococcus pyogenes, a Gram-positive pathogen, has been determined and contains 1,752 predicted protein-encoding genes. Approximately one-third of these genes have no identifiable function, with the remainder falling into previously characterized categories of known microbial function. Consistent with the observation that S. pyogenes is responsible for a wider variety of human disease than any other bacterial species, more than 40 putative virulence-associated genes have been identified. Additional genes have been identified that encode proteins likely associated with microbial "molecular mimicry" of host characteristics and involved in rheumatic fever or acute glomerulonephritis. The complete or partial sequence of four different bacteriophage genomes is also present, with each containing genes for one or more previously undiscovered superantigen-like proteins. These prophage-associated genes encode at least six potential virulence factors, emphasizing the importance of bacteriophages in horizontal gene transfer and a possible mechanism for generating new strains with increased pathogenic potential.

An introduction to MCMC for machine learning

C Andrieu — 2006-08-24T15:04:55-00:00

(2003)

This purpose of this introductory paper is threefold. First, it introduces the Monte Carlo method with emphasis on probabilistic machine learning. Second, it reviews the main building blocks of modern Markov chain Monte Carlo simulation, thereby providing and introduction to the remaining papers of this special issue. Lastly, it discusses new interesting research horizons.

Monte Carlo Strategies in Scientific Computing

Jun Liu — 2006-07-31T12:04:14-00:00

(17 October 2002)

A large number of scientists and engineers employ Monte Carlo simulation and related global optimization techniques (such as simulated annealing) as an essential tool in their work. For such scientists, there is a need to keep up to date with several recent advances in Monte Carlo methodologies such as cluster methods, data- augmentation, simulated tempering and other auxiliary variable methods. There is also a trend in moving towards a population-based approach. All these advances in one way or another were motivated by the need to sample from very complex distribution for which traditional methods would tend to be trapped in local energy minima. It is our aim to provide a self-contained and up to date treatment of the Monte Carlo method to this audience. The Monte Carlo method is a computer-based statistical sampling approach for solving numerical problems concerned with a complex system. The methodology was initially developed in the field of statistical physics during the early days of electronic computing (1945-55) and has now been adopted by researchers in almost all scientific fields. The fundamental idea for constructing Markov chain based Monte Carlo algorithms was introduced in the 1950s. This idea was later extended to handle more and more complex physical systems. In the 1980s, statisticians and computer scientists developed Monter Carlo-based algorithms for a wide variety of integration and optimization tasks. In the 1990s, the method began to play an important role in computational biology. Over the past fifty years, reasearchers in diverse scientific fields have studied the Monte Carlo method and contributed to its development. Today, a large number of scientisits and engineers employ Monte Carlo techniques as an essential tool in their work. For such scientists, there is a need to keep up-to-date with recent advances in Monte Carlo methodologies.

Markov Chain Monte Carlo in Practice

WR Gilks — 2006-01-19T12:28:51-00:00

(01 December 1995)

General state-space Markov chain theory has evolved to make it both more accessible and more powerful. Markov Chain Monte Carlo in Practice introduces MCMC methods and their applications while also providing some theoretical background. Considering the broad audience, the editors emphasize practice rather than theory and keep the technical content to a minimum. They offer step-by-step instructions for using the methods presented and show the importance of MCMC in real applications with examples ranging from the simple to the more complex in fields such as archaeology, astronomy, biostatistics, genetics, epidemiology, and image analysis.

Monte Carlo sampling methods using Markov chains and their applications

WK Hastings — 2006-12-27T04:52:05-00:00

Biometrika, Vol. 57, No. 1. (1 April 1970), pp. 97-109.

A generalization of the sampling method introduced by Metropolis et al. (1953) is presented along with an exposition of the relevant theory, techniques of application and methods and difficulties of assessing the error in Monte Carlo estimates. Examples of the methods, including the generation of random orthogonal matrices and potential applications of the methods to numerical problems arising in statistics, are discussed. 10.1093/biomet/57.1.97

Markov Chain Monte Carlo Method and Its Application

Stephen Brooks — 2006-09-05T17:41:17-00:00

The Markov chain Monte Carlo (MCMC) method, as a computer-intensive statistical tool, has enjoyed an enormous upsurge in interest over the last few years. This paper provides a simple, comprehensive and tutorial review of some of the most common areas of research in this field. We begin by discussing how MCMC algorithms can be constructed from standard building-blocks to produce Markov chains with the desired stationary distribution. We also motivate and discuss more complex ideas that have been proposed in the literature, such as continuous time and dimension jumping methods. We discuss some implementational issues associated with MCMC methods. We take a look at the arguments for and against multiple replications, consider how long chains should be run for and how to determine suitable starting points. We also take a look at graphical models and how graphical approaches can be used to simplify MCMC implementation. Finally, we present a couple of examples, which we use as case-studies to highlight some of the points made earlier in the text. In particular, we use a simple changepoint model to illustrate how to tackle a typical Bayesian modelling problem via the MCMC method, before using mixture model problems to provide illustrations of good sampler output and of the implementation of a reversible jump MCMC algorithm.

Bayesian Computation and Stochastic Systems

Julian Besag — 2006-10-24T18:03:57-00:00

Statistical Science, Vol. 10, No. 1. (1995), pp. 3-41.

Markov chain Monte Carlo (MCMC) methods have been used extensively in statistical physics over the last 40 years, in spatial statistics for the past 20 and in Bayesian image analysis over the last decade. In the last five years, MCMC has been introduced into significance testing, general Bayesian inference and maximum likelihood estimation. This paper presents basic methodology of MCMC, emphasizing the Bayesian paradigm, conditional probability and the intimate relationship with Markov random fields in spatial statistics. Hastings algorithms are discussed, including Gibbs, Metropolis and some other variations. Pairwise difference priors are described and are used subsequently in three Bayesian applications, in each of which there is a pronounced spatial or temporal aspect to the modeling. The examples involve logistic regression in the presence of unobserved covariates and ordinal factors; the analysis of agricultural field experiments, with adjustment for fertility gradients; and processing of low-resolution medical images obtained by a gamma camera. Additional methodological issues arise in each of these applications and in the Appendices. The paper lays particular emphasis on the calculation of posterior probabilities and concurs with others in its view that MCMC facilitates a fundamental breakthrough in applied Bayesian modeling.