Examination of the Enterotoxigenic Escherichia coli Population Structure during Human Infection.

UNLABELLED: Enterotoxigenic E. coli (ETEC) can cause severe diarrhea and death in children in developing countries; however, bacterial diversity in natural infection is uncharacterized. In this study, we explored the natural population variation of ETEC from individuals with cholera-like diarrhea. Genomic sequencing and comparative analysis of multiple ETEC isolates from twelve cases of severe diarrhea demonstrated clonal populations in the majority of subjects (10/12). In contrast, a minority of individuals (2/12) yielded phylogenomically divergent ETEC isolates. Detailed examination revealed that isolates also differed in virulence factor content. These genomic data suggest that severe, cholera-like ETEC infections are largely caused by a clonal population of organisms within individual patients. Additionally, the isolation of similar clones from geographically and temporally dispersed cases with similar clinical presentations suggests that some isolates are particularly suited for virulence. The identification of multiple genomically diverse isolates with variable virulence factor profiles from a single subject highlights the dynamic nature of ETEC, as well as a potential weakness in the examination of cultures obtained from a single colony in clinical settings. These findings have implications for vaccine design and provide a framework for the study of population variation in other human pathogens. IMPORTANCE: Enterotoxigenic Escherichia coli (ETEC) has been identified as one of the major causes of diarrheal diseases in children as well as travelers. It has been previously appreciated that this pathogenic variant of E. coli is diverse, both at the genomic level, as defined with multilocus sequence typing, and with regard to the presence or absence of virulence factors within clonal groups. Using whole-genome sequencing and comparative analysis, we identified and characterized diverse enterotoxigenic E. coli isolates from individual patients. In 17% of patients, we identified multiple distinct ETEC isolates, each with unique genomic features and in some cases diverse virulence factor profiles. These studies ascertained that any one person may be colonized by multiple pathogenic ETEC isolates, which may impact how we think about the development of vaccines and therapeutics against these organisms.

Standardized metadata for human pathogen/vector genomic sequences.

High throughput sequencing has accelerated the determination of genome sequences for thousands of human infectious disease pathogens and dozens of their vectors. The scale and scope of these data are enabling genotype-phenotype association studies to identify genetic determinants of pathogen virulence and drug/insecticide resistance, and phylogenetic studies to track the origin and spread of disease outbreaks. To maximize the utility of genomic sequences for these purposes, it is essential that metadata about the pathogen/vector isolate characteristics be collected and made available in organized, clear, and consistent formats. Here we report the development of the GSCID/BRC Project and Sample Application Standard, developed by representatives of the Genome Sequencing Centers for Infectious Diseases (GSCIDs), the Bioinformatics Resource Centers (BRCs) for Infectious Diseases, and the U.S. National Institute of Allergy and Infectious Diseases (NIAID), part of the National Institutes of Health (NIH), informed by interactions with numerous collaborating scientists. It includes mapping to terms from other data standards initiatives, including the Genomic Standards Consortium's minimal information (MIxS) and NCBI's BioSample/BioProjects checklists and the Ontology for Biomedical Investigations (OBI). The standard includes data fields about characteristics of the organism or environmental source of the specimen, spatial-temporal information about the specimen isolation event, phenotypic characteristics of the pathogen/vector isolated, and project leadership and support. By modeling metadata fields into an ontology-based semantic framework and reusing existing ontologies and minimum information checklists, the application standard can be extended to support additional project-specific data fields and integrated with other data represented with comparable standards. The use of this metadata standard by all ongoing and future GSCID sequencing projects will provide a consistent representation of these data in the BRC resources and other repositories that leverage these data, allowing investigators to identify relevant genomic sequences and perform comparative genomics analyses that are both statistically meaningful and biologically relevant.

High-level relatedness among Mycobacterium abscessus subsp. massiliense strains from widely separated outbreaks.

Three recently sequenced strains isolated from patients during an outbreak of Mycobacterium abscessus subsp. massiliense infections at a cystic fibrosis center in the United States were compared with 6 strains from an outbreak at a cystic fibrosis center in the United Kingdom and worldwide strains. Strains from the 2 cystic fibrosis outbreaks showed high-level relatedness with each other and major-level relatedness with strains that caused soft tissue infections during an epidemic in Brazil. We identified unique single-nucleotide polymorphisms in cystic fibrosis and soft tissue outbreak strains, separate single-nucleotide polymorphisms only in cystic fibrosis outbreak strains, and unique genomic traits for each subset of isolates. Our findings highlight the necessity of identifying M. abscessus to the subspecies level and screening all cystic fibrosis isolates for relatedness to these outbreak strains. We propose 2 diagnostic strategies that use partial sequencing of rpoB and secA1 genes and a multilocus sequence typing protocol.

Molecular signature of high yield (growth) influenza a virus reassortants prepared as candidate vaccine seeds.

BACKGROUND: Human influenza virus isolates generally grow poorly in embryonated chicken eggs. Hence, gene reassortment of influenza A wild type (wt) viruses is performed with a highly egg adapted donor virus, A/Puerto Rico/8/1934 (PR8), to provide the high yield reassortant (HYR) viral 'seeds' for vaccine production. HYR must contain the hemagglutinin (HA) and neuraminidase (NA) genes of wt virus and one to six 'internal' genes from PR8. Most studies of influenza wt and HYRs have focused on the HA gene. The main objective of this study is the identification of the molecular signature in all eight gene segments of influenza A HYR candidate vaccine seeds associated with high growth in ovo. METHODOLOGY: The genomes of 14 wt parental viruses, 23 HYRs (5 H1N1; 2, 1976 H1N1-SOIV; 2, 2009 H1N1pdm; 2 H2N2 and 12 H3N2) and PR8 were sequenced using the high-throughput sequencing pipeline with big dye terminator chemistry. RESULTS: Silent and coding mutations were found in all internal genes derived from PR8 with the exception of the M gene. The M gene derived from PR8 was invariant in all 23 HYRs underlining the critical role of PR8 M in high yield phenotype. None of the wt virus derived internal genes had any silent change(s) except the PB1 gene in X-157. The highest number of recurrent silent and coding mutations was found in NS. With respect to the surface antigens, the majority of HYRs had coding mutations in HA; only 2 HYRs had coding mutations in NA. SIGNIFICANCE: In the era of application of reverse genetics to alter influenza A virus genomes, the mutations identified in the HYR gene segments associated with high growth in ovo may be of great practical benefit to modify PR8 and/or wt virus gene sequences for improved growth of vaccine 'seed' viruses.

Genomic insights into the emerging human pathogen Mycobacterium massiliense.

Mycobacterium massiliense (Mycobacterium abscessus group) is an emerging pathogen causing pulmonary disease and skin and soft tissue infections. We report the genome sequence of the type strain CCUG 48898.

Genomic diversity of 2010 Haitian cholera outbreak strains.

The millions of deaths from cholera during the past 200 y, coupled with the morbidity and mortality of cholera in Haiti since October 2010, are grim reminders that Vibrio cholerae, the etiologic agent of cholera, remains a scourge. We report the isolation of both V. cholerae O1 and non-O1/O139 early in the Haiti cholera epidemic from samples collected from victims in 18 towns across eight Arrondissements of Haiti. The results showed two distinct populations of V. cholerae coexisted in Haiti early in the epidemic. As non-O1/O139 V. cholerae was the sole pathogen isolated from 21% of the clinical specimens, its role in this epidemic, either alone or in concert with V. cholerae O1, cannot be dismissed. A genomic approach was used to examine similarities and differences among the Haitian V. cholerae O1 and V. cholerae non-O1/O139 strains. A total of 47 V. cholerae O1 and 29 V. cholerae non-O1/O139 isolates from patients and the environment were sequenced. Comparative genome analyses of the 76 genomes and eight reference strains of V. cholerae isolated in concurrent epidemics outside Haiti and 27 V. cholerae genomes available in the public database demonstrated substantial diversity of V. cholerae and ongoing flux within its genome.

PB2 residue 158 is a pathogenic determinant of pandemic H1N1 and H5 influenza a viruses in mice.

Influenza A viruses are human and animal pathogens that cause morbidity and mortality, which range from mild to severe. The 2009 H1N1 pandemic was caused by the emergence of a reassortant H1N1 subtype (H1N1pdm) influenza A virus containing gene segments that originally circulated in human, avian, and swine virus reservoirs. The molecular determinants of replication and pathogenesis of H1N1pdm viruses in humans and other mammals are poorly understood. Therefore, we set out to elucidate viral determinants critical to the pathogenesis of this novel reassortant using a mouse model. We found that a glutamate-to-glycine substitution at residue 158 of the PB2 gene (PB2-E158G) increased the morbidity and mortality of the parental H1N1pdm virus. Results from mini-genome replication assays in human cells and virus titration in mouse tissues demonstrated that PB2-E158G is a pathogenic determinant, because it significantly increases viral replication rates. The virus load in PB2-E158G-infected mouse lungs was 1,300-fold higher than that of the wild-type virus. Our data also show that PB2-E158G had a much stronger influence on the RNA replication and pathogenesis of H1N1pdm viruses than PB2-E627K, which is a known pathogenic determinant. Remarkably, PB2-E158G substitutions also altered the pathotypes of two avian H5 viruses in mice, indicating that this residue impacts genetically divergent influenza A viruses and suggesting that this region of PB2 could be a new antiviral target. Collectively, the data presented in this study demonstrate that PB2-E158G is a novel pathogenic determinant of influenza A viruses in the mouse model. We speculate that PB2-E158G may be important in the adaptation of avian PB2 genes to other mammals, and BLAST sequence analysis identified a naturally occurring human H1N1pdm isolate that has this substitution. Therefore, future surveillance efforts should include scrutiny of this region of PB2 because of its potential impact on pathogenesis.

Unseasonal transmission of H3N2 influenza A virus during the swine-origin H1N1 pandemic.

The initial wave of swine-origin influenza A virus (pandemic H1N1/09) in the United States during the spring and summer of 2009 also resulted in an increased vigilance and sampling of seasonal influenza viruses (H1N1 and H3N2), even though they are normally characterized by very low incidence outside of the winter months. To explore the nature of virus evolution during this influenza "off-season," we conducted a phylogenetic analysis of H1N1 and H3N2 sequences sampled during April to June 2009 in New York State. Our analysis revealed that multiple lineages of both viruses were introduced and cocirculated during this time, as is typical of influenza virus during the winter. Strikingly, however, we also found strong evidence for the presence of a large transmission chain of H3N2 viruses centered on the south-east of New York State and which continued until at least 1 June 2009. These results suggest that the unseasonal transmission of influenza A viruses may be more widespread than is usually supposed.

The early diversification of influenza A/H1N1pdm.

Background Since its initial detection in April 2009, the A/H1N1pdm influenza virus has spread rapidly in humans, with over 5,700 human deaths. However, little is known about the evolutionary dynamics of H1N1pdm and its geographic and temporal diversification.Methods Phylogenetic analysis was conducted upon the concatenated coding regions of whole-genome sequences from 290 H1N1pdm isolates sampled globally between April 1 - July 9, 2009, including relatively large samples from the US states of Wisconsin and New York. Results At least 7 phylogenetically distinct viral clades have disseminated globally and co-circulated in localities that experienced multiple introductions of H1N1pdm. The epidemics in New York and Wisconsin were dominated by two different clades, both phylogenetically distinct from the viruses first identified in California and Mexico, suggesting an important role for founder effects in determining local viral population structures. Conclusions Determining the global diversity of H1N1pdm is central to understanding the evolution and spatial spread of the current pandemic, and to predict its future impact on human populations. Our results indicate that H1N1pdm has already diversified into distinct viral lineages with defined spatial patterns.

Evolutionary dynamics of human rotaviruses: balancing reassortment with preferred genome constellations.

Group A human rotaviruses (RVs) are a major cause of severe gastroenteritis in infants and young children. Yet, aside from the genes encoding serotype antigens (VP7; G-type and VP4; P-type), little is known about the genetic make-up of emerging and endemic human RV strains. To gain insight into the diversity and evolution of RVs circulating at a single location over a period of time, we sequenced the eleven-segmented, double-stranded RNA genomes of fifty-one G3P[8] strains collected from 1974 to 1991 at Children's Hospital National Medical Center, Washington, D. C. During this period, G1P[8] strains typically dominated, comprising on average 56% of RV infections each year in hospitalized children. A notable exception was in the 1976 and 1991 winter seasons when the incidence of G1P[8] infections decreased dramatically, a trend that correlated with a significant increase in G3P[8] infections. Our sequence analysis indicates that the 1976 season was characterized by the presence of several genetically distinct, co-circulating clades of G3P[8] viruses, which contained minor but significant differences in their encoded proteins. These 1976 lineages did not readily exchange gene segments with each other, but instead remained stable over the course of the season. In contrast, the 1991 season contained a single major clade, whose genome constellation was similar to one of the 1976 clades. The 1991 clade may have gained a fitness advantage after reassorting with as of yet unidentified RV strain(s). This study reveals for the first time that genetically distinct RV clades of the same G/P-type can co-circulate and cause disease. The findings from this study also suggest that, although gene segment exchange occurs, most reassortant strains are replaced over time by lineages with preferred genome constellations. Elucidation of the selective pressures that favor maintenance of RVs with certain sets of genes may be necessary to anticipate future vaccine needs.

Insights into the Musa genome: syntenic relationships to rice and between Musa species.

BACKGROUND: Musa species (Zingiberaceae, Zingiberales) including bananas and plantains are collectively the fourth most important crop in developing countries. Knowledge concerning Musa genome structure and the origin of distinct cultivars has greatly increased over the last few years. Until now, however, no large-scale analyses of Musa genomic sequence have been conducted. This study compares genomic sequence in two Musa species with orthologous regions in the rice genome. RESULTS: We produced 1.4 Mb of Musa sequence from 13 BAC clones, annotated and analyzed them along with 4 previously sequenced BACs. The 443 predicted genes revealed that Zingiberales genes share GC content and distribution characteristics with eudicot and Poaceae genomes. Comparison with rice revealed microsynteny regions that have persisted since the divergence of the Commelinid orders Poales and Zingiberales at least 117 Mya. The previously hypothesized large-scale duplication event in the common ancestor of major cereal lineages within the Poaceae was verified. The divergence time distributions for Musa-Zingiber (Zingiberaceae, Zingiberales) orthologs and paralogs provide strong evidence for a large-scale duplication event in the Musa lineage after its divergence from the Zingiberaceae approximately 61 Mya. Comparisons of genomic regions from M. acuminata and M. balbisiana revealed highly conserved genome structure, and indicated that these genomes diverged circa 4.6 Mya. CONCLUSION: These results point to the utility of comparative analyses between distantly-related monocot species such as rice and Musa for improving our understanding of monocot genome evolution. Sequencing the genome of M. acuminata would provide a strong foundation for comparative genomics in the monocots. In addition a genome sequence would aid genomic and genetic analyses of cultivated Musa polyploid genotypes in research aimed at localizing and cloning genes controlling important agronomic traits for breeding purposes.

Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy.

We sequenced 2.2 Mb representing triplicated genome segments of Brassica oleracea, which are each paralogous with one another and homologous with a segmentally duplicated region of the Arabidopsis thaliana genome. Sequence annotation identified 177 conserved collinear genes in the B. oleracea genome segments. Analysis of synonymous base substitution rates indicated that the triplicated Brassica genome segments diverged from a common ancestor soon after divergence of the Arabidopsis and Brassica lineages. This conclusion was corroborated by phylogenetic analysis of protein families. Using A. thaliana as an outgroup, 35% of the genes inferred to be present when genome triplication occurred in the Brassica lineage have been lost, most likely via a deletion mechanism, in an interspersed pattern. Genes encoding proteins involved in signal transduction or transcription were not found to be significantly more extensively retained than those encoding proteins classified with other functions, but putative proteins predicted in the A. thaliana genome were underrepresented in B. oleracea. We identified one example of gene loss from the Arabidopsis lineage. We found evidence for the frequent insertion of gene fragments of nuclear genomic origin and identified four apparently intact genes in noncollinear positions in the B. oleracea and A. thaliana genomes.