SNP Discovery in European Anchovy (Engraulis encrasicolus, L) by High-Throughput Transcriptome and Genome Sequencing

12 p.

The Demand for Eastern Oysters, Crassostrea virginica, from the Gulf of Mexico in the Presence of Vibrio vulnificus

California, in response to health concerns, initiated a program on 1 March 1991 which required anyone selling eastern oysters, Crassostrea virginica, from the Gulf of Mexico area to notify potential consumers that there was a risk in consuming them raw. This mandatory warning, followed shortly thereafter by a similar warning in other states, including Louisiana and Florida, received extensive press cover-age throughout the country and particularly in the Gulf area. This paper examines the extent to which the demand for Gulf-area oysters has been reduced as a result of mandatory warning labels and negative publicity. In general, the results suggest that since 1991 the “summer” dockside price has been reduced by about 50% as a result of warning labels and associated negative publicity, while the “winter” dockside price has been reduced by about 30%.

Predicting the distribution of Vibrio vulnificus in Chesapeake Bay

Vibrio vulnificus is a gram-negative pathogenic bacterium endemic to coastal waters worldwide, and a leading cause of seafood related mortality. Because of human health concerns, understanding the ecology of the species and potentially predicting its distribution is of great importance. We evaluated and applied a previously published qPCR assay to water samples (n = 235) collected from the main-stem of the Chesapeake Bay (2007 – 2008) by Maryland and Virginia State water quality monitoring programs. Results confirmed strong relationships between the likelihood of Vibrio vulnificus presence and both temperature and salinity that were used to develop a logistic regression model. The habitat model demonstrated a high degree of concordance (93%), and robustness as subsequent bootstrapping (n=1000) did not change model output (P > 0.05). We forced this empirical habitat model with temperature and salinity predictions generated by a regional hydrodynamic modeling system to demonstrate its utility in future pathogen forecasting efforts in the Chesapeake Bay.

Elucidation of hepatitis C virus transmission and early diversification by single genome sequencing.

A precise molecular identification of transmitted hepatitis C virus (HCV) genomes could illuminate key aspects of transmission biology, immunopathogenesis and natural history. We used single genome sequencing of 2,922 half or quarter genomes from plasma viral RNA to identify transmitted/founder (T/F) viruses in 17 subjects with acute community-acquired HCV infection. Sequences from 13 of 17 acute subjects, but none of 14 chronic controls, exhibited one or more discrete low diversity viral lineages. Sequences within each lineage generally revealed a star-like phylogeny of mutations that coalesced to unambiguous T/F viral genomes. Numbers of transmitted viruses leading to productive clinical infection were estimated to range from 1 to 37 or more (median = 4). Four acutely infected subjects showed a distinctly different pattern of virus diversity that deviated from a star-like phylogeny. In these cases, empirical analysis and mathematical modeling suggested high multiplicity virus transmission from individuals who themselves were acutely infected or had experienced a virus population bottleneck due to antiviral drug therapy. These results provide new quantitative and qualitative insights into HCV transmission, revealing for the first time virus-host interactions that successful vaccines or treatment interventions will need to overcome. Our findings further suggest a novel experimental strategy for identifying full-length T/F genomes for proteome-wide analyses of HCV biology and adaptation to antiviral drug or immune pressures.

Comparative analysis of Klebsiella pneumoniae genomes identifies a phospholipase D family protein as a novel virulence factor

BACKGROUND: Klebsiella pneumoniae strains are pathogenic to animals and humans, in which they are both a frequent cause of nosocomial infections and a re-emerging cause of severe community-acquired infections. K. pneumoniae isolates of the capsular serotype K2 are among the most virulent. In order to identify novel putative virulence factors that may account for the severity of K2 infections, the genome sequence of the K2 reference strain Kp52.145 was determined and compared to two K1 and K2 strains of low virulence and to the reference strains MGH 78578 and NTUH-K2044.

RESULTS: In addition to diverse functions related to host colonization and virulence encoded in genomic regions common to the four strains, four genomic islands specific for Kp52.145 were identified. These regions encoded genes for the synthesis of colibactin toxin, a putative cytotoxin outer membrane protein, secretion systems, nucleases and eukaryotic-like proteins. In addition, an insertion within a type VI secretion system locus included sel1 domain containing proteins and a phospholipase D family protein (PLD1). The pld1 mutant was avirulent in a pneumonia model in mouse. The pld1 mRNA was expressed in vivo and the pld1 gene was associated with K. pneumoniae isolates from severe infections. Analysis of lipid composition of a defective E. coli strain complemented with pld1 suggests an involvement of PLD1 in cardiolipin metabolism.

CONCLUSIONS: Determination of the complete genome of the K2 reference strain identified several genomic islands comprising putative elements of pathogenicity. The role of PLD1 in pathogenesis was demonstrated for the first time and suggests that lipid metabolism is a novel virulence mechanism of K. pneumoniae.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled 'reservoir' host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers. © 2012 Biek et al.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

Prediction of Staphylococcus aureus antimicrobial resistance by whole-genome sequencing

Whole-genome sequencing (WGS) could potentially provide a single platform for extracting all the information required to predict an organism’s phenotype. However, its ability to provide accurate predictions has not yet been demonstrated in large independent studies of specific organisms. In this study, we aimed to develop a genotypic prediction method for antimicrobial susceptibilities. The whole genomes of 501 unrelated Staphylococcus aureus isolates were sequenced, and the assembled genomes were interrogated using BLASTn for a panel of known resistance determinants (chromosomal mutations and genes carried on plasmids). Results were compared with phenotypic susceptibility testing for 12 commonly used antimicrobial agents (penicillin, methicillin, erythromycin, clindamycin, tetracycline, ciprofloxacin, vancomycin, trimethoprim, gentamicin, fusidic acid, rifampin, and mupirocin) performed by the routine clinical laboratory. We investigated discrepancies by repeat susceptibility testing and manual inspection of the sequences and used this information to optimize the resistance determinant panel and BLASTn algorithm. We then tested performance of the optimized tool in an independent validation set of 491 unrelated isolates, with phenotypic results obtained in duplicate by automated broth dilution (BD Phoenix) and disc diffusion. In the validation set, the overall sensitivity and specificity of the genomic prediction method were 0.97 (95% confidence interval [95% CI], 0.95 to 0.98) and 0.99 (95% CI, 0.99 to 1), respectively, compared to standard susceptibility testing methods. The very major error rate was 0.5%, and the major error rate was 0.7%. WGS was as sensitive and specific as routine antimicrobial susceptibility testing methods. WGS is a promising alternative to culture methods for resistance prediction in S. aureus and ultimately other major bacterial pathogens.

Genome sequencing: a grin without a cat

In some types of unicellular algae, the chloroplasts have their own nucleus — a legacy of the time when the chloroplast was a free-living cell. The sequence of the genome in one such nucleus is now revealed.

Construction of a river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel and preliminary RH mapping of chromosomes 3 and 10

The buffalo (Bubalus bubalis) not only is a useful source of milk, it also provides meat and works as a natural source of labor and biogas. To establish a project for buffalo genome mapping a 5,000-rad whole genome radiation hybrid panel was constructed for river buffalo and used to build preliminary RH maps from two chromosomes (BBU 3 and BBU10). The preliminary maps contain 66 markers, including coding genes, cattle ESTs and microsatellite loci. The RH maps presented here are the starting point for mapping additional loci, in particular, genes and expressed sequence tags that will allow detailed comparative maps between buffalo, cattle and other species to be constructed. A large quantity of DNA has been prepared from the cell lines forming the RH panel reported here and will be made publicly available to the international community both for the study of chromosome evolution and for the improvement of traits important to the role of buffalo in animal agriculture.

Construction of a river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel and preliminary RH mapping of chromosomes 3 and 10

The buffalo (Bubalus bubalis) is a source of milk and meat, and also serves as a draft animal. In this study, a 5000-rad whole-genome radiation hybrid (RH) panel for river buffalo was constructed and used to build preliminary RH maps for BBU3 and BBU10 chromosomes. The preliminary maps contain 66 markers, including coding genes, cattle expressed sequence tags (ESTs) and microsatellite loci. The RH maps presented here are the starting point for mapping additional loci that will allow detailed comparative maps between buffalo, cattle and other species whose genomes may be mapped in the future. A large quantity of DNA has been prepared from the cell lines forming the river buffalo RH panel and will be made publicly available to the international community both for the study of chromosome evolution and for the improvement of traits important to the role of buffalo in animal agriculture.