Antibody profiling by proteome microarray reveals the immunogenicity of the attenuated smallpox vaccine modified vaccinia virus ankara is comparable to that of Dryvax

Modified vaccinia virus Ankara (MVA) is a highly attenuated vaccinia virus that is under consideration as an alternative to the conventional smallpox vaccine Dryvax. MVA was attenuated by extensive passage of vaccinia virus Ankara in chicken embryo fibroblasts. Several immunomodulatory genes and genes that influence host range are deleted or mutated, and replication is aborted in the late stage of infection in most nonavian cells. The effect of these mutations on immunogenicity is not well understood. Since the structural genes appear to be intact in MVA, it is hypothesized that critical targets for antibody neutralization have been retained. To test this, we probed microarrays of the Western Reserve (WR) proteome with sera from humans and macaques after MVA and Dryvax vaccination. As most protein sequences of MVA are 97 to 99% identical to those of other vaccinia virus strains, extensive binding cross-reactivity is expected, except for those deleted or truncated. Despite different hosts and immunization regimens, the MVA and Dryvax antibody profiles were broadly similar, with antibodies against membrane and core proteins being the best conserved. The responses to nonstructural proteins were less well conserved, although these are not expected to influence virus neutralization. The broadest antibody response was obtained for hyperimmune rabbits with WR, which is pathogenic in rabbits. These data indicate that, despite the mutations and deletions in MVA, its overall immunogenicity is broadly comparable to that of Dryvax, particularly at the level of antibodies to membrane proteins. The work supports other information suggesting that MVA may be a useful alternative to Dryvax.     

Construction, phenotypic analysis, and immunogenicity of a UL5/UL29 double deletion mutant of herpes simplex virus 2

A number of studies have shown that replication-defective mutant strains of herpes simplex virus (HSV) can induce protective immunity in animal systems against wild-type HSV challenge. However, all of those studies used viruses with single mutations. Because multiple, stable mutations provide optimal levels of safety for live vaccines, we felt that additional mutations needed to be engineered into a candidate vaccine strain for HSV-2 and genital herpes. We therefore isolated an HSV-2 strain with deletion mutations in two viral DNA replication protein genes, UL5 and UL29. The resulting double deletion mutant virus strain, dl5-29, fails to form plaques or to give any detectable single cycle yields in normal monkey or human cells. Nevertheless, dl5-29 expresses nearly the same pattern of gene products as the wild-type virus or the single mutant viruses and induces antibody titers in mice that are equivalent to those induced by single deletion mutant viruses. Therefore, it is feasible to isolate a mutant HSV strain with two mutations in essential genes and with an increased level of safety but which is still highly immunogenic.     

High fidelity of murine hepatitis virus replication is decreased in nsp14 exoribonuclease mutants

Replication fidelity of RNA virus genomes is constrained by the opposing necessities of generating sufficient diversity for adaptation and maintaining genetic stability, but it is unclear how the largest viral RNA genomes have evolved and are maintained under these constraints. A coronavirus (CoV) nonstructural protein, nsp14, contains conserved active-site motifs of cellular exonucleases, including DNA proofreading enzymes, and the severe acute respiratory syndrome CoV (SARS-CoV) nsp14 has 3'-to-5' exoribonuclease (ExoN) activity in vitro. Here, we show that nsp14 ExoN remarkably increases replication fidelity of the CoV murine hepatitis virus (MHV). Replacement of conserved MHV ExoN active-site residues with alanines resulted in viable mutant viruses with growth and RNA synthesis defects that during passage accumulated 15-fold more mutations than wild-type virus without changes in growth fitness. The estimated mutation rate for ExoN mutants was similar to that reported for other RNA viruses, whereas that of wild-type MHV was less than the established rates for RNA viruses in general, suggesting that CoVs with intact ExoN replicate with unusually high fidelity. Our results indicate that nsp14 ExoN plays a critical role in prevention or repair of nucleotide incorporation errors during genome replication. The established mutants are unique tools to test the hypothesis that high replication fidelity is required for the evolution and stability of large RNA genomes.     

Genomic Sequence and Virulence of Clonal Isolates of Vaccinia Virus Tiantan,the Chinese Smallpox Vaccine Strain

Despite the worldwide eradication of smallpox in 1979, the potential bioterrorism threat from variola virus and the ongoing use of vaccinia virus (VACV) as a vector for vaccine development argue for continued research on VACV. In China, the VACV Tiantan strain (TT) was used in the smallpox eradication campaign. Its progeny strain is currently being used to develop a human immunodeficiency virus (HIV) vaccine. Here we sequenced the full genomes of five TT clones isolated by plaque purification from the TT (752-1) viral stock. Phylogenetic analysis with other commonly used VACV strains showed that TT (752-1) and its clones clustered and exhibited higher sequence diversity than that found in Dryvax clones. The ∼190 kbp genomes of TT appeared to encode 273 open reading frames (ORFs). ORFs located in the middle of the genome were more conserved than those located at the two termini, where many virulence and immunomodulation associated genes reside. Several patterns of nucleotide changes including point mutations, insertions and deletions were identified. The polymorphisms in seven virulence-associated proteins and six immunomodulation-related proteins were analyzed. We also investigated the neuro- and skin- virulence of TT clones in mice and rabbits, respectively. The TT clones exhibited significantly less virulence than the New York City Board of Health (NYCBH) strain, as evidenced by less extensive weight loss and morbidity in mice as well as produced smaller skin lesions and lower incidence of putrescence in rabbits. The complete genome sequences, ORF annotations, and phenotypic diversity yielded from this study aid our understanding of the Chinese historic TT strain and are useful for HIV vaccine projects employing TT as a vector.     

Defective Interference in the Killer System of Saccharomyces cerevisiae

The K₁ killer virus (or plasmid) of Saccharomyces cerevisiae is a noninfectious double-stranded RNA genome found intracellularly packaged in an icosahedral capsid. This genome codes for a protein toxin and for resistance to that toxin. Defective interfering virus mutants are deletion derivatives of the killer virus double-stranded RNA genome; such mutants are called suppressive. Unlike strains carrying the wild-type genome, strains with these deletion derivatives are neither toxin producers nor toxin resistant. If both the suppressive and the wildtype virus are introduced into the same cell, most progeny become toxin-sensitive nonkillers (J. M. Somers, Genetics 74:571-579, 1973). Diploids formed by the mating of a killer with a suppressive strain were grown in liquid culture, and RNA was extracted from samples taken up to 41 generations after the mating. The ratio of killer RNA to suppressive RNA decreased with increasing generations; by 41 generations the killer RNA was barely detectable. The copy numbers of the suppressive genome and its parental killer were virtually the same in isogenic strains, as were the growth rates of diploid strains containing either virus alone. Therefore, suppressiveness, not being due to segregation or overgrowth by faster growing segregants, is likely due to preferential replication or maintenance of the suppressive genome. Three suppressive viruses, all derivatives of the same killer virus (T. K. Sweeney et al., Genetics 84:27-42, 1976), did not coexist stably. The evidence strongly indicates that the largest genome of the three slowly suppressed both of the smaller genomes, showing that larger genomes can suppress smaller ones and that suppression can occur between two suppressives. Of 48 isolates of strains carrying the suppressive viruses, 5 had newly detectable RNA species, all larger than the original suppressive genomes. At least seven genes necessary for maintenance of the wild-type killer virus (MAK genes) were needed by a suppressive mutant. No effect of ski mutations (affecting regulation of killer virus double-stranded RNA replication) on suppressiveness was observed.     

Role of biased hypermutation in evolution of subacute sclerosing panencephalitis virus from progenitor acute measles virus.

We identified an acute measles virus (Nagahata strain) closely related to a defective virus (Biken strain) isolated from a patient with subacute sclerosing panencephalitis (SSPE). The proteins of Nagahata strain measles virus are antigenically and electrophoretically similar to the proteins of Edmonston strain measles virus. However, the nucleotide sequence of the Nagahata matrix (M) gene is significantly different from the M genes of all the acute measles virus strains studied to date. The Nagahata M gene is strikingly similar to the M gene of Biken strain SSPE virus isolated several years later in the same locale. Eighty percent of the nucleotide differences between the Nagahata and Biken M genes are uridine-to-cytosine transitions known as biased hypermutation, which has been postulated to be caused by a cellular RNA-modifying activity. These biased mutations account for all but one of the numerous missense genetic changes predicted to cause amino acid substitutions. As a result, the Biken virus M protein loses conformation-specific epitopes that are conserved in the M proteins of Nagahata and Edmonston strain acute measles viruses. These conformation-specific epitopes are also absent in the cryptic M proteins encoded by the hypermutated M genes of two other defective SSPE viruses (Niigata and Yamagata strains). Nagahata-like sequences are found in the M genes of at least five other SSPE viruses isolated from three continents. These data indicate that Biken strain SSPE virus is derived from a progenitor closely resembling Nagahata strain acute measles virus and that biased hypermutation is largely responsible for the structural defects in the Biken virus M protein.     

Evolution of human influenza A viruses in nature: sequential mutations in the genomes of new H1N1

The genetic variation of the new pandemic H1N1 influenza A viruses isolated in 1977 was analyzed by two-dimensional oligonucleotide fingerprinting and RNA sequencing. Differences were observed in the fingerprints of the RNAs of these viruses, and analysis of the changes suggested that sequential mutations occurred in their genomes. Based on these data, a scheme is presented which proposes divergent evolution of strains from a common ancestry. Furthermore, it was found that mutations were not restricted to the genes coding for the hemagglutinin and the neuraminidase, but were scattered throughout the genome, suggesting that selective antibody pressure is not solely responsible for the emergence of genetic variants. Our data also strengthen the hypothesis that the new H1N1 influenza virus strains are derived from strains circulating in 1950.     

Antigenic and genetic variation in influenza A (H1N1) virus isolates recovered from a persistently infected immunodeficient child.

Antigenic and genetic variations have been analyzed in eight consecutive isolates recovered from a child with severe combined immunodeficiency syndrome persistently infected with naturally acquired type A (H1N1) influenza virus over a 10-month period. Hemagglutination inhibition reactions and T1 oligonucleotide fingerprinting demonstrated that these viruses were related to strains causing outbreaks in the United States at that time (1983 to 1984) but that antigenic and genetic differences between consecutive isolates could be detected. This variation between isolates was examined further by sequencing the RNAs encoding the HA1 region of the hemagglutinin (HA) and the nucleoprotein (NP) in five of the consecutive isolates. Multiple point mutations were detected in both genes, and a deletion of one amino acid was detected in the HA. Depending on the isolates compared, 5.8 x 10(-3) to 17 x 10(-3) substitutions per nucleotide site per year were detected in the RNAs encoding the HA1, and 3.5 x 10(-3) to 24 x 10(-3) substitutions per nucleotide site per year were detected in the NP gene. Fifty-four percent of the base changes in the HA1 and 73% in the NP led to amino acid substitutions. A progressive accumulation of mutations over time was not observed, suggesting that the genetic diversity of these viruses may best be interpreted as the result of shifts in the population equilibrium (quasi-species) of replicating variant genomes.     

Genome sequence analyses of Pseudomonas savastanoi pv. glycinea and subtractive hybridization-based comparative genomics with nine pseudomonads

Bacterial blight, caused by Pseudomonas savastanoi pv. glycinea (Psg), is a common disease of soybean. In an effort to compare a current field isolate with one isolated in the early 1960s, the genomes of two Psg strains, race 4 and B076, were sequenced using 454 pyrosequencing. The genomes of both Psg strains share more than 4,900 highly conserved genes, indicating very low genetic diversity between Psg genomes. Though conserved, genome rearrangements and recombination events occur commonly within the two Psg genomes. When compared to each other, 437 and 163 specific genes were identified in B076 and race 4, respectively. Most specific genes are plasmid-borne, indicating that acquisition and maintenance of plasmids may represent a major mechanism to change the genetic composition of the genome and even acquire new virulence factors. Type three secretion gene clusters of Psg strains are near identical with that of P. savastanoi pv. phaseolicola (Pph) strain 1448A and they shared 20 common effector genes. Furthermore, the coronatine biosynthetic cluster is present on a large plasmid in strain B076, but not in race 4. In silico subtractive hybridization-based comparative genomic analyses with nine sequenced phytopathogenic pseudomonads identified dozens of specific islands (SIs), and revealed that the genomes of Psg strains are more similar to those belonging to the same genomospecies such as Pph 1448A than to other phytopathogenic pseudomonads. The number of highly conserved genes (core genome) among them decreased dramatically when more genomes were included in the subtraction, suggesting the diversification of pseudomonads, and further indicating the genome heterogeneity among pseudomonads. However, the number of specific genes did not change significantly, suggesting these genes are indeed specific in Psg genomes. These results reinforce the idea of a species complex of P. syringae and support the reclassification of P. syringae into different species.     

Evidence for the common origin of viral and cellular sequences involved in sarcomagenic transformation.

The src genes of six different strains of avian sarcoma virus (ASV) were compared with those of a series of newly isolated sarcoma viruses, termed "recovery avian sarcoma viruses" (rASV's). The rASV's were isolated recently from chicken and quail tumors induced by transformation-defective (td) deletion mutants of Schmidt-Ruppin Rous sarcoma virus. The RNase T1-resistant oligonucleotide maps were constructed for the RNA genomes of different strains of ASV and td mutants. The src-specific sequences, characterized by RNase T1-resistant oligonucleotides ranging from 9 to 19 nucleotides long, were defined as those mapping between approximately 600 and 2,800 nucleotides from the 3' polyadenylate end of individual sarcoma viral RNAs, and missing in the corresponding td viral RNAs. Our results revealed that 12 src-specific oligonucleotides were highly conserved among several strains of ASV, including the rASV's, whereas certain strains of ASV were found to contain one to three characteristic src-specific oligonucleotides. We previously presented evidence supporting the idea that most of the src-specific sequences present in rASV RNAs are derived from cellular genetic information. Our present data indicate that the src genes of rASV's are closely related to other known ASVs. We conclude that the src genes of different strains of ASV and the cellular sarc sequences are of common origin, although some divergence has occurred among different viral src genes and related cellular sequences.     

Genetic heterogeneity in human T-cell leukemia/lymphoma virus type II.

DNA from the peripheral blood mononuclear cells of 17 different individuals infected with human T-cell lymphoma/leukemia virus type II (HTLV-II) was successfully amplified by the polymerase chain reaction (PCR) with the primer pair SK110/SK111. This primer pair is conserved among the pol genes of all primate T-cell lymphoma viruses (PTLV) and flanks a 140-bp fragment of DNA which, when used in comparative analyses, reflects the relative degree of diversity among PTLV genomes. Cloning, sequencing, and phylogenetic comparisons of these amplified 140-bp pol fragments indicated that there are at least two distinct genetic substrains of HTLV-II in the Western Hemisphere. These data were confirmed for selected isolates by performing PCR, cloning, and sequencing with to 10 additional primer pair-probe sets specific for different regions throughout the PTLV genome. HTLV-II isolates from Seminole, Guaymi, and Tobas Indians belong in the new substrain of HTLV-II, while the prototype MoT isolate defines the original substrain. There was greater diversity among HTLV-II New World strains than among HTLV-I New World strains. In fact, the heterogeneity among HTLV-II strains from the Western Hemisphere was similar to that observed in HTLV-I and simian T-cell lymphoma/leukemia virus type I isolates from around the world, including Japan, Africa, and Papua New Guinea. Given these geographic and anthropological considerations and assuming similar mutation rates and selective forces among the PTLV, these data suggest either that HTLV-II has existed for a long time in the indigenous Amerindian population or that HTLV-II isolates introduced into the New World were more heterogeneous than the HTLV-I strains introduced into the New World.     

Analysis of genome plasticity in pathogenic and commensal Escherichia coli isolates by use of DNA arrays

Genomes of prokaryotes differ significantly in size and DNA composition. Escherichia coli is considered a model organism to analyze the processes involved in bacterial genome evolution, as the species comprises numerous pathogenic and commensal variants. Pathogenic and nonpathogenic E. coli strains differ in the presence and absence of additional DNA elements contributing to specific virulence traits and also in the presence and absence of additional genetic information. To analyze the genetic diversity of pathogenic and commensal E. coli isolates, a whole-genome approach was applied. Using DNA arrays, the presence of all translatable open reading frames (ORFs) of nonpathogenic E. coli K-12 strain MG1655 was investigated in 26 E. coli isolates, including various extraintestinal and intestinal pathogenic E. coli isolates, 3 pathogenicity island deletion mutants, and commensal and laboratory strains. Additionally, the presence of virulence-associated genes of E. coli was determined using a DNA "pathoarray" developed in our laboratory. The frequency and distributional pattern of genomic variations vary widely in different E. coli strains. Up to 10% of the E. coli K-12-specific ORFs were not detectable in the genomes of the different strains. DNA sequences described for extraintestinal or intestinal pathogenic E. coli are more frequently detectable in isolates of the same origin than in other pathotypes. Several genes coding for virulence or fitness factors are also present in commensal E. coli isolates. Based on these results, the conserved E. coli core genome is estimated to consist of at least 3,100 translatable ORFs. The absence of K-12-specific ORFs was detectable in all chromosomal regions. These data demonstrate the great genome heterogeneity and genetic diversity among E. coli strains and underline the fact that both the acquisition and deletion of DNA elements are important processes involved in the evolution of prokaryotes.     

Identification of conserved lentiviral sequences as landmarks of genomic flexibility

Considering that recombinations produce quasispecies in lentivirus spreading, we identified and localized highly conserved sequences that may play an important role in viral ontology. Comparison of entire genomes, including 237 human, simian and non-primate mammal lentiviruses and 103 negative control viruses, led to identify 28 Conserved Lentiviral Sequences (CLSs). They were located mainly in the structural genes forming hot spots particularly in the gag and pol genes and to a lesser extent in LTRs and regulatory genes. The CLS pattern was the same throughout the different HIV-1 subtypes, except for some HIV-1-O strains. Only CLS 3 and 4 were detected in both negative control HTLV-1 oncornaviruses and D-particle-forming simian viruses, which are not immunodeficiency inducers and display a genetic stability. CLSs divided the virus genomes into domains allowing us to distinguish sequence families leading to the notion of 'species self' besides that of 'lentiviral self'. Most of acutely localized CLSs in HIV-1s (82%) corresponded to wide recombination segments being currently reported.     

Related haloarchaeal pleomorphic viruses contain different genome types

Archaeal viruses have been the subject of recent interest due to the diversity discovered in their virion architectures. Recently, a new group of haloarchaeal pleomorphic viruses has been discovered. It is distinctive in terms of the virion morphology and different genome types (ssDNA/dsDNA) harboured by rather closely related representatives. To date there are seven isolated viruses belonging to this group. Most of these share a cluster of five conserved genes, two of which encode major structural proteins. Putative proviruses and proviral remnants containing homologues of the conserved gene cluster were also identified suggesting a long-standing relationship of these viruses with their hosts. Comparative genomic analysis revealed three different ways of the genome organization, which possibly reflect different replication strategies employed by these viruses. The dsDNA genomes of two of these viruses were shown to contain single-strand interruptions. Further studies on one of the genomes suggested that the interruptions are located along the genome in a sequence-specific manner and exhibit polarity in distribution.     

The genomes of sheeppox and goatpox viruses

Sheeppox virus (SPPV) and goatpox virus (GTPV), members of the Capripoxvirus genus of the Poxviridae, are etiologic agents of important diseases of sheep and goats in northern and central Africa, southwest and central Asia, and the Indian subcontinent. Here we report the genomic sequence and comparative analysis of five SPPV and GTPV isolates, including three pathogenic field isolates and two attenuated vaccine viruses. SPPV and GTPV genomes are approximately 150 kbp and are strikingly similar to each other, exhibiting 96% nucleotide identity over their entire length. Wild-type genomes share at least 147 putative genes, including conserved poxvirus replicative and structural genes and genes likely involved in virulence and host range. SPPV and GTPV genomes are very similar to that of lumpy skin disease virus (LSDV), sharing 97% nucleotide identity. All SPPV and GTPV genes are present in LSDV. Notably in both SPPV and GTPV genomes, nine LSDV genes with likely virulence and host range functions are disrupted, including a gene unique to LSDV (LSDV132) and genes similar to those coding for interleukin-1 receptor, myxoma virus M003.2 and M004.1 genes (two copies each), and vaccinia virus F11L, N2L, and K7L genes. The absence of these genes in SPPV and GTPV suggests a significant role for them in the bovine host range. SPPV and GTPV genomes contain specific nucleotide differences, suggesting they are phylogenetically distinct. Relatively few genomic changes in SPPV and GTPV vaccine viruses account for viral attenuation, because they contain 71 and 7 genomic changes compared to their respective field strains. Notable genetic changes include mutation or disruption of genes with predicted functions involving virulence and host range, including two ankyrin repeat proteins in SPPV and three kelch-like proteins in GTPV. These comparative genomic data indicate the close genetic relationship among capripoxviruses, and they suggest that SPPV and GTPV are distinct and likely derived from an LSDV-like ancestor.     

Complete genome analysis of 33 ecologically and biologically diverse Rift Valley fever virus strains reveals widespread virus movement and low genetic diversity due to recent common ancestry

Rift Valley fever (RVF) virus is a mosquito-borne RNA virus responsible for large explosive outbreaks of acute febrile disease in humans and livestock in Africa with significant mortality and economic impact. The successful high-throughput generation of the complete genome sequence was achieved for 33 diverse RVF virus strains collected from throughout Africa and Saudi Arabia from 1944 to 2000, including strains differing in pathogenicity in disease models. While several distinct virus genetic lineages were determined, which approximately correlate with geographic origin, multiple exceptions indicative of long-distance virus movement have been found. Virus strains isolated within an epidemic (e.g., Mauritania, 1987, or Egypt, 1977 to 1978) exhibit little diversity, while those in enzootic settings (e.g., 1970s Zimbabwe) can be highly diverse. In addition, the large Saudi Arabian RVF outbreak in 2000 appears to have involved virus introduction from East Africa, based on the close ancestral relationship of a 1998 East African virus. Virus genetic diversity was low (~5%) and primarily involved accumulation of mutations at an average of 2.9 × 10⁻⁴ substitutions/site/year, although some evidence of RNA segment reassortment was found. Bayesian analysis of current RVF virus genetic diversity places the most recent common ancestor of these viruses in the late 1800s, the colonial period in Africa, a time of dramatic changes in agricultural practices and introduction of nonindigenous livestock breeds. In addition to insights into the evolution and ecology of RVF virus, these genomic data also provide a foundation for the design of molecular detection assays and prototype vaccines useful in combating this important disease.     

Gene amplification and expression by RNA viruses and potential for further application to plant gene transfer

The great majority of plant viruses encapsidate messenger-sense ssRNA and have no natural DNA phase in their life cycle. Despite their RNA nature, essentially any desired change can be introduced into such genomes by using recombinant DNA techniques with suitably constructed, expressible viral cDNA clones. For some viruses such as brome mosaic virus, these methods have been used to define the sequences controlling RNA-directed genomic RNA replication and the expression of internal genes via subgenomic mRNAs. The results suggest a surprising degree of genetic flexibility, which appears to be reflected in the varied gene complements and genetic organizations of presumably related plant and animal RNA viruses sharing conserved replication genes. Foreign genes inserted in such RNA virus genomes can be amplified and expressed to a high level in transfected plant cells. In addition to the potential use of such viruses as episomal expression vectors, it should be possible to couple the viral pathways of RNA-dependent RNA synthesis to amplify and to further regulate the expression of genes transformed into plant chromosomes.     

Characterization of genome structure of amphotropic and ecotropic wild mouse retroviruses.

We studied the RNA genomes of several wild mouse type C retroviruses by using RNase T1-oligonucleotide fingerprinting. The amphotropic and ecotropic viruses of field strain 1504 produced very similar oligonucleotide fingerprints, but each also had several unique oligonucleotides. All of these unique oligonucleotides were located in the env gene region and were probably responsible for the host range differences between these viruses, as well as the lymphomagenic and paralytogenic properties of the viruses. We obtained similar results with the amphotropic and ecotropic viruses of another field strain (4070), which was isolated from a mouse from a different trapping area. The amphotropic viruses of several field strains (strains 1504, 4070, and 1313) were more closely related than the ecotropic viruses of different strains (strains 1504, 4070, and 4996). These findings suggested that the genetic sequences of the amphotropic viruses are more conserved than those of ecotropic viruses isolated from the same wild mice.     

A phylogenomic analysis of Escherichia coli / Shigella group: implications of genomic features associated with pathogenicity and ecological adaptation

ABSTRACT: BACKGROUND: The Escherichia coli species contains a variety of commensal and pathogenic strains, and its intraspecific diversity is extraordinarily high. With the availability of an increasing number of E. coli strain genomes, a more comprehensive concept of their evolutionary history and ecological adaptation can be developed using phylogenomic analyses. In this study, we constructed two types of whole-genome phylogenies based on 34 E. coli strains using collinear genomic segments. The first phylogeny was based on the concatenated collinear regions shared by all of the studied genomes, and the second phylogeny was based on the variable collinear regions that are absent from at least one genome. Intuitively, the first phylogeny is likely to reveal the lineal evolutionary history among these strains (i.e., an evolutionary phylogeny), whereas the latter phylogeny is likely to reflect the whole-genome similarities of extant strains (i.e., a similarity phylogeny). RESULTS: Within the evolutionary phylogeny, the strains were clustered in accordance with known phylogenetic groups and phenotypes. When comparing evolutionary and similarity phylogenies, a concept emerges that Shigella may have originated from at least three distinct ancestors and evolved into a single clade. By scrutinizing the properties that are shared amongst Shigella strains but missing in other E. coli genomes, we found that the common regions of the Shigella genomes were mainly influenced by mobile genetic elements, implying that they may have experienced convergent evolution via horizontal gene transfer. Based on an inspection of certain key branches of interest, we identified several collinear regions that may be associated with the pathogenicity of specific strains. Moreover, by examining the annotated genes within these regions, further detailed evidence associated with pathogenicity was revealed. CONCLUSIONS: Collinear regions are reliable genomic features used for phylogenomic analysis among closely related genomes while linking the genomic diversity with phenotypic differences in a meaningful way. The pathogenicity of a strain may be associated with both the arrival of virulence factors and the modification of genomes via mutations. Such phylogenomic studies that compare collinear regions of whole genomes will help to better understand the evolution and adaptation of closely related microbes and E. coli in particular.