Direct Amplification, Sequencing and Profiling of Chlamydia trachomatis Strains in Single and Mixed Infection Clinical Samples

Sequencing bacterial genomes from DNA isolated directly from clinical samples offers the promise of rapid and precise acquisition of informative genetic information. In the case of Chlamydia trachomatis, direct sequencing is particularly desirable because it obviates the requirement for culture in mammalian cells, saving time, cost and the possibility of missing low abundance strains. In this proof of concept study, we developed methodology that would allow genome-scale direct sequencing, using a multiplexed microdroplet PCR enrichment technology to amplify a 100 kb region of the C. trachomatis genome with 500 1.1–1.3 kb overlapping amplicons (5-fold amplicon redundancy). We integrated comparative genomic data into a pipeline to preferentially select conserved sites for amplicon design. The 100 kb target region could be amplified from clinical samples, including remnants from diagnostics tests, originating from the cervix, urethra and urine, For rapid analysis of these data, we developed a framework for whole-genome based genotyping called binstrain. We used binstrain to estimate the proportion of SNPs originating from 14 C. trachomatis reference serotype genomes in each sample. Direct DNA sequencing methods such as the one described here may have an important role in understanding the biology of C. trachomatis mixed infections and the natural genetic variation of the species within clinically relevant ecological niches.

[1]  N. Nagelkerke,et al.  The epidemiology of Chlamydia trachomatis within a sexually transmitted diseases core group. , 1996, The Journal of infectious diseases.

[2]  K. Millman,et al.  Molecular and mutation trends analyses of omp1 alleles for serovar E of Chlamydia trachomatis. Implications for the immunopathogenesis of disease. , 1997, The Journal of clinical investigation.

[3]  R. W. Davis,et al.  Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis. , 1998, Science.

[4]  S. Salzberg,et al.  Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. , 2000, Nucleic acids research.

[5]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[6]  D. Dean,et al.  Evidence for long-term cervical persistence of Chlamydia trachomatis by omp1 genotyping. , 2000, The Journal of infectious diseases.

[7]  F. Blattner,et al.  Mauve: multiple alignment of conserved genomic sequence with rearrangements. , 2004, Genome research.

[8]  M. P. Cummings.  PHYLIP (Phylogeny Inference Package) , 2004 .

[9]  S. Franceschi,et al.  The natural course of Chlamydia trachomatis infection in asymptomatic Colombian women: a 5-year follow-up study. , 2005, The Journal of infectious diseases.

[10]  W. Bruno,et al.  Evolution of Chlamydia trachomatis diversity occurs by widespread interstrain recombination involving hotspots. , 2006, Genome research.

[11]  Ziheng Yang.  PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[12]  Jack A. M. Leunissen,et al.  Turning CFCs into salt. , 1996, Nucleic Acids Res..

[13]  I. Goodhead,et al.  Chlamydia trachomatis: genome sequence analysis of lymphogranuloma venereum isolates. , 2007, Genome research.

[14]  Joseph A. Cook.  Eliminating blinding trachoma. , 2008, The New England journal of medicine.

[15]  Bill W Colston,et al.  High-throughput quantitative polymerase chain reaction in picoliter droplets. , 2008, Analytical chemistry.

[16]  R. Demars,et al.  Interstrain Gene Transfer in Chlamydia trachomatis In Vitro: Mechanism and Significance , 2007, Journal of bacteriology.

[17]  D. Dean,et al.  Multiple Chlamydiaceae Species in Trachoma: Implications for Disease Pathogenesis and Control , 2008, PLoS medicine.

[18]  W. Bruno,et al.  Predicting Phenotype and Emerging Strains among Chlamydia trachomatis Infections , 2009, Emerging infectious diseases.

[19]  P. Keim,et al.  Humans and evolutionary and ecological forces shaped the phylogeography of recently emerged diseases , 2009, Nature Reviews Microbiology.

[20]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[21]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[22]  D. Rockey,et al.  Genome Sequencing of Recent Clinical Chlamydia trachomatis Strains Identifies Loci Associated with Tissue Tropism and Regions of Apparent Recombination , 2010, Infection and Immunity.

[23]  T. Read,et al.  Hypervirulent Chlamydia trachomatis Clinical Strain Is a Recombinant between Lymphogranuloma Venereum (L2) and D Lineages , 2011, mBio.

[24]  R. Skilton,et al.  Development of a Transformation System for Chlamydia trachomatis: Restoration of Glycogen Biosynthesis by Acquisition of a Plasmid Shuttle Vector , 2011, PLoS pathogens.

[25]  Matthias Meyer,et al.  A draft genome of Yersinia pestis from victims of the Black Death , 2011, Nature.

[26]  Kevin Galinsky,et al.  Hybrid selection for sequencing pathogen genomes from clinical samples , 2011, Genome Biology.

[27]  N. Verlander,et al.  Genital chlamydial infection: association between clinical features, organism genotype and load. , 2011, Journal of medical microbiology.

[28]  T. Read,et al.  Interplay of recombination and selection in the genomes of Chlamydia trachomatis , 2011, Biology Direct.

[29]  B. Moumen,et al.  Targeted genome enrichment for efficient purification of endosymbiont DNA from host DNA , 2012, Symbiosis.

[30]  T. Read,et al.  Population genomics of Chlamydia trachomatis: insights on drift, selection, recombination, and population structure. , 2012, Molecular biology and evolution.

[31]  T. Read,et al.  Genome-wide recombination in Chlamydia trachomatis , 2012, Nature Genetics.

[32]  Katherine H. Huang,et al.  A framework for human microbiome research , 2012, Nature.

[33]  S. Boyer,et al.  Spider: An R package for the analysis of species identity and evolution, with particular reference to DNA barcoding , 2012, Molecular ecology resources.

[34]  Justin Zobel,et al.  Short read sequence typing (SRST): multi-locus sequence types from short reads , 2012, BMC Genomics.

[35]  Julian Parkhill,et al.  Whole genome analysis of diverse Chlamydia trachomatis strains identifies phylogenetic relationships masked by current clinical typing , 2012, Nature Genetics.

[36]  Leping Li,et al.  ART: a next-generation sequencing read simulator , 2012, Bioinform..

[37]  R. Valdivia,et al.  Virulence determinants in the obligate intracellular pathogen Chlamydia trachomatis revealed by forward genetic approaches , 2012, Proceedings of the National Academy of Sciences.

[38]  Julian Parkhill,et al.  Whole-genome sequences of Chlamydia trachomatis directly from clinical samples without culture , 2013, Genome research.

[39]  T. Putman,et al.  Culture-independent sequence analysis of Chlamydia trachomatis in urogenital specimens identifies regions of recombination and in-patient sequence mutations. , 2013, Microbiology.