Michael H. Holmes

We present a draft sequence of the genome of Aedes aegypti, the primary vector for yellow fever and dengue fever, which at approximately 1376 million base pairs is about 5 times the size of the genome of the malaria vector Anopheles gambiae. Nearly 50% of the Ae. aegypti genome consists of transposable elements. These contribute to a factor of approximately 4 to 6 increase in average gene length and in sizes of intergenic regions relative to An. gambiae and Drosophila melanogaster. Nonetheless, chromosomal synteny is generally maintained among all three insects, although conservation of orthologous gene order is higher (by a factor of approximately 2) between the mosquito species than between either of them and the fruit fly. An increase in genes encoding odorant binding, cytochrome P450, and cuticle domains relative to An. gambiae suggests that members of these protein families underpin some of the biological differences between the two mosquito species.

Abstract Motivation: Most sequence comparison methods assume that the data being compared are trustworthy, but this is not the case with raw DNA sequences obtained from automatic sequencing machines. Nevertheless, sequence comparisons need to be done on them in order to remove vector splice sites and contaminants. This step is necessary before other genomic data processing stages can be carried out, such as fragment assembly or EST clustering. A specialized tool is therefore needed to solve this apparent dilemma. Results: We have designed and implemented a program that specifically addresses the problem. This program, called Lucy, has been in use since 1998 at The Institute for Genomic Research (TIGR). During this period, many rounds of experience-driven modifications were made to Lucy to improve its accuracy and its ability to deal with extremely difficult input cases. We believe we have finally obtained a useful program which strikes a delicate balance among the many issues involved in the raw sequence cleaning problem, and we wish to share it with the research community. Availability: Lucy is available directly from TIGR (http://www.tigr.org/softlab). Academic users can download Lucy after accepting a free academic use license. Business users may need to pay a license fee to use Lucy for commercial purposes. Contact: Questions regarding the quality assessment module of Lucy should be directed to Michael Holmes (mholmes@tigr.org). Questions regarding other aspects of Lucy should be directed to Hui-Hsien Chou (hhchou@iastate.edu).

We describe here a clinical daptomycin treatment failure in a patient with recurrent methicillin-resistant Staphylococcus aureus (MRSA) bacteremia in whom daptomycin was administered after a failed empirical treatment course with vancomycin and piperacillin-tazobactam. We had the opportunity to compare the genome sequences of an isogenic pair of daptomycin-susceptible and -resistant MRSA isolates obtained before and after initiation of daptomycin therapy, respectively. The genotype of both isolates was USA800, ST5, SCCmec type IV, agr type II. There was no increase in cell wall thickness in the daptomycin-resistant strain despite having decreased susceptibility to both vancomycin and daptomycin. By comparing the genome sequences by pyrosequencing, we identified a polymorphism (S337L) in the tenth transmembrane segment of the multiple peptide resistance factor, MprF, encoding lysyl phosphatidylglycerol transferase. This enzyme has been shown previously to promote repulsion of daptomycin at the cell surface by addition of positively charged lysine to phosphatidylglycerol. Also, the hlb open reading frame (ORF) encoding the β-toxin was interrupted by a prophage in the daptomycin-susceptible strain; this phage was missing in the daptomycin-resistant isolate and the hlb ORF was restored. Loss of the phage in the resistant isolate also resulted in loss of the virulence factor genes clpP, scn, and sak. This is the first study to use pyrosequencing to compare the genomes of a daptomycin-susceptible/resistant MRSA isolate pair obtained during failed daptomycin therapy in humans.

BACKGROUND: A low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates. RESULTS: We have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2). CONCLUSION: Whole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.