A Retrospective Study on Genetic Heterogeneity within Treponema Strains: Subpopulations Are Genetically Distinct in a Limited Number of Positions

Background

Pathogenic uncultivable treponemes comprise human and animal pathogens including agents of syphilis, yaws, bejel, pinta, and venereal spirochetosis in rabbits and hares. A set of 10 treponemal genome sequences including those of 4 Treponema pallidum ssp. pallidum (TPA) strains (Nichols, DAL-1, Mexico A, SS14), 4 T. p. ssp. pertenue (TPE) strains (CDC-2, Gauthier, Samoa D, Fribourg-Blanc), 1 T. p. ssp. endemicum (TEN) strain (Bosnia A) and one strain (Cuniculi A) of Treponema paraluisleporidarum ecovar Cuniculus (TPLC) were examined with respect to the presence of nucleotide intrastrain heterogeneous sites.

Methodology/Principal Findings

The number of identified intrastrain heterogeneous sites in individual genomes ranged between 0 and 7. Altogether, 23 intrastrain heterogeneous sites (in 17 genes) were found in 5 out of 10 investigated treponemal genomes including TPA strains Nichols (n = 5), DAL-1 (n = 4), and SS14 (n = 7), TPE strain Samoa D (n = 1), and TEN strain Bosnia A (n = 5). Although only one heterogeneous site was identified among 4 tested TPE strains, 16 such sites were identified among 4 TPA strains. Heterogeneous sites were mostly strain-specific and were identified in four tpr genes (tprC, GI, I, K), in genes involved in bacterial motility and chemotaxis (fliI, cheC-fliY), in genes involved in cell structure (murC), translation (prfA), general and DNA metabolism (putative SAM dependent methyltransferase, topA), and in seven hypothetical genes.

Conclusions/Significance

Heterogeneous sites likely represent both the selection of adaptive changes during infection of the host as well as an ongoing diversifying evolutionary process.     

Whole Genome Sequence of Treponema pallidum ssp. pallidum,Strain Mexico A,Suggests Recombination between Yaws and Syphilis Strains

Background

Treponema pallidum ssp. pallidum (TPA), the causative agent of syphilis, and Treponema pallidum ssp. pertenue (TPE), the causative agent of yaws, are closely related spirochetes causing diseases with distinct clinical manifestations. The TPA Mexico A strain was isolated in 1953 from male, with primary syphilis, living in Mexico. Attempts to cultivate TPA Mexico A strain under in vitro conditions have revealed lower growth potential compared to other tested TPA strains.

Methodology/Principal Findings

The complete genome sequence of the TPA Mexico A strain was determined using the Illumina sequencing technique. The genome sequence assembly was verified using the whole genome fingerprinting technique and the final sequence was annotated. The genome size of the Mexico A strain was determined to be 1,140,038 bp with 1,035 predicted ORFs. The Mexico A genome sequence was compared to the whole genome sequences of three TPA (Nichols, SS14 and Chicago) and three TPE (CDC-2, Samoa D and Gauthier) strains. No large rearrangements in the Mexico A genome were found and the identified nucleotide changes occurred most frequently in genes encoding putative virulence factors. Nevertheless, the genome of the Mexico A strain, revealed two genes (TPAMA_0326 (tp92) and TPAMA_0488 (mcp2-1)) which combine TPA- and TPE- specific nucleotide sequences. Both genes were found to be under positive selection within TPA strains and also between TPA and TPE strains.

Conclusions/Significance

The observed mosaic character of the TPAMA_0326 and TPAMA_0488 loci is likely a result of inter-strain recombination between TPA and TPE strains during simultaneous infection of a single host suggesting horizontal gene transfer between treponemal subspecies.     

Whole genome sequences of three Treponema pallidum ssp. pertenue strains: yaws and syphilis treponemes differ in less than 0.2% of the genome sequence

Background

The yaws treponemes, Treponema pallidum ssp. pertenue (TPE) strains, are closely related to syphilis causing strains of Treponema pallidum ssp. pallidum (TPA). Both yaws and syphilis are distinguished on the basis of epidemiological characteristics, clinical symptoms, and several genetic signatures of the corresponding causative agents.

Methodology/Principal Findings

To precisely define genetic differences between TPA and TPE, high-quality whole genome sequences of three TPE strains (Samoa D, CDC-2, Gauthier) were determined using next-generation sequencing techniques. TPE genome sequences were compared to four genomes of TPA strains (Nichols, DAL-1, SS14, Chicago). The genome structure was identical in all three TPE strains with similar length ranging between 1,139,330 bp and 1,139,744 bp. No major genome rearrangements were found when compared to the four TPA genomes. The whole genome nucleotide divergence (d_A) between TPA and TPE subspecies was 4.7 and 4.8 times higher than the observed nucleotide diversity (π) among TPA and TPE strains, respectively, corresponding to 99.8% identity between TPA and TPE genomes. A set of 97 (9.9%) TPE genes encoded proteins containing two or more amino acid replacements or other major sequence changes. The TPE divergent genes were mostly from the group encoding potential virulence factors and genes encoding proteins with unknown function.

Conclusions/Significance

Hypothetical genes, with genetic differences, consistently found between TPE and TPA strains are candidates for syphilitic treponemes virulence factors. Seventeen TPE genes were predicted under positive selection, and eleven of them coded either for predicted exported proteins or membrane proteins suggesting their possible association with the cell surface. Sequence changes between TPE and TPA strains and changes specific to individual strains represent suitable targets for subspecies- and strain-specific molecular diagnostics.     

Evolutionary Processes in the Emergence and Recent Spread of the Syphilis Agent,Treponema pallidum

The incidence of syphilis has risen worldwide in the last decade in spite of being an easily treated infection. The causative agent of this sexually transmitted disease is the bacterium Treponema pallidum subspecies pallidum (TPA), very closely related to subsp. pertenue (TPE) and endemicum (TEN), responsible for the human treponematoses yaws and bejel, respectively. Although much focus has been placed on the question of the spatial and temporary origins of TPA, the processes driving the evolution and epidemiological spread of TPA since its divergence from TPE and TEN are not well understood. Here, we investigate the effects of recombination and selection as forces of genetic diversity and differentiation acting during the evolution of T. pallidum subspecies. Using a custom-tailored procedure, named phylogenetic incongruence method, with 75 complete genome sequences, we found strong evidence for recombination among the T. pallidum subspecies, involving 12 genes and 21 events. In most cases, only one recombination event per gene was detected and all but one event corresponded to intersubspecies transfers, from TPE/TEN to TPA. We found a clear signal of natural selection acting on the recombinant genes, which is more intense in their recombinant regions. The phylogenetic location of the recombination events detected and the functional role of the genes with signals of positive selection suggest that these evolutionary processes had a key role in the evolution and recent expansion of the syphilis bacteria and significant implications for the selection of vaccine candidates and the design of a broadly protective syphilis vaccine.     

Whole Genome Sequence of the Treponema Fribourg-Blanc: Unspecified Simian Isolate Is Highly Similar to the Yaws Subspecies

Background

Unclassified simian strain Treponema Fribourg-Blanc was isolated in 1966 from baboons (Papio cynocephalus) in West Africa. This strain was morphologically indistinguishable from T. pallidum ssp. pallidum or ssp. pertenue strains, and it was shown to cause human infections.

Methodology/Principal Findings

To precisely define genetic differences between Treponema Fribourg-Blanc (unclassified simian isolate, FB) and T. pallidum ssp. pertenue strains (TPE), a high quality sequence of the whole Fribourg-Blanc genome was determined with 454-pyrosequencing and Illumina sequencing platforms. Combined average coverage of both methods was greater than 500×. Restriction target sites (n = 1,773), identified in silico, of selected restriction enzymes within the Fribourg-Blanc genome were verified experimentally and no discrepancies were found. When compared to the other three sequenced TPE genomes (Samoa D, CDC-2, Gauthier), no major genome rearrangements were found. The Fribourg-Blanc genome clustered with other TPE strains (especially with the TPE CDC-2 strain), while T. pallidum ssp. pallidum strains clustered separately as well as the genome of T. paraluiscuniculi strain Cuniculi A. Within coding regions, 6 deletions, 5 insertions and 117 substitutions differentiated Fribourg-Blanc from other TPE genomes.

Conclusions/Significance

The Fribourg-Blanc genome showed similar genetic characteristics as other TPE strains. Therefore, we propose to rename the unclassified simian isolate to Treponema pallidum ssp. pertenue strain Fribourg-Blanc. Since the Fribourg-Blanc strain was shown to cause experimental infection in human hosts, non-human primates could serve as possible reservoirs of TPE strains. This could considerably complicate recent efforts to eradicate yaws. Genetic differences specific for Fribourg-Blanc could then contribute for identification of cases of animal-derived yaws infections.     

On the origin of the treponematoses: a phylogenetic approach

Background

Since the first recorded epidemic of syphilis in 1495, controversy has surrounded the origins of the bacterium Treponema pallidum subsp. pallidum and its relationship to the pathogens responsible for the other treponemal diseases: yaws, endemic syphilis, and pinta. Some researchers have argued that the syphilis-causing bacterium, or its progenitor, was brought from the New World to Europe by Christopher Columbus and his men, while others maintain that the treponematoses, including syphilis, have a much longer history on the European continent.

Methodology/Principal Findings

We applied phylogenetics to this problem, using data from 21 genetic regions examined in 26 geographically disparate strains of pathogenic Treponema. Of all the strains examined, the venereal syphilis-causing strains originated most recently and were more closely related to yaws-causing strains from South America than to other non-venereal strains. Old World yaws-causing strains occupied a basal position on the tree, indicating that they arose first in human history, and a simian strain of T. pallidum was found to be indistinguishable from them.

Conclusions/Significance

Our results lend support to the Columbian theory of syphilis''s origin while suggesting that the non-sexually transmitted subspecies arose earlier in the Old World. This study represents the first attempt to address the problem of the origin of syphilis using molecular genetics, as well as the first source of information regarding the genetic make-up of non-venereal strains from the Western hemisphere.     

Whole genome sequences of Treponema pallidum subsp. endemicum isolated from Cuban patients: The non-clonal character of isolates suggests a persistent human infection rather than a single outbreak

Bejel (endemic syphilis) is a neglected non-venereal disease caused by Treponema pallidum subsp. endemicum (TEN). Although it is mostly present in hot, dry climates, a few cases have been found outside of these areas. The aim of this work was the sequencing and analysis of TEN isolates obtained from “syphilis patients” in Cuba, which is not considered an endemic area for bejel. Genomes were obtained by pool segment genome sequencing or direct sequencing methods, and the bioinformatics analysis was performed according to an established pipeline. We obtained four genomes with 100%, 81.7%, 52.6%, and 21.1% breadth of coverage, respectively. The sequenced genomes revealed a non-clonal character, with nucleotide variability ranging between 0.2–10.3 nucleotide substitutions per 100 kbp among the TEN isolates. Nucleotide changes affected 27 genes, and the analysis of the completely sequenced genome also showed a recombination event between tprC and tprI, in TP0488 as well as in the intergenic region between TP0127–TP0129. Despite limitations in the quality of samples affecting breadth of sequencing coverage, the determined non-clonal character of the isolates suggests a persistent infection in the Cuban population rather than a single outbreak caused by imported case.     

Syphilis Screening among 27,150 Pregnant Women in South Chinese Rural Areas Using Point-of-Care Tests

Objectives

To determine the prevalence and correlates of syphilis among pregnant women in rural areas of South China.

Methods

Point-of-care syphilis testing was provided at 71 health facilities in less developed, rural areas of Guangdong Province. Positive samples were confirmed at a local referral center by toluidine red unheated serum tests (TRUST) and Treponema pallidum particle agglutination (TPPA) tests.

Results

Altogether 27,150 pregnant women in rural Guangdong were screened for syphilis. 106 (0.39%) syphilis cases were diagnosed, of which 78 (73.6%) received treatment for syphilis. Multivariate analysis revealed that older pregnant women (31–35 years old, aOR 2.7, 95% CI 0.99–7.32; older than 35 years old, aOR 5.9, 95% CI 2.13–16.34) and those with a history of adverse pregnant outcomes (aOR 3.64, 95% CI 2.30–5.76) were more likely to be infected with syphilis.

Conclusions

A high prevalence of syphilis exists among pregnant women living in rural areas of South China. Enhanced integration of syphilis screening with other routine women''s health services (OB GYN, family planning) may be useful for controlling China''s syphilis epidemic.     

Evaluation of Genome Sequencing Quality in Selected Plant Species Using Expressed Sequence Tags

Background

With the completion of genome sequencing projects for more than 30 plant species, large volumes of genome sequences have been produced and stored in online databases. Advancements in sequencing technologies have reduced the cost and time of whole genome sequencing enabling more and more plants to be subjected to genome sequencing. Despite this, genome sequence qualities of multiple plants have not been evaluated.

Methodology/Principal Finding

Integrity and accuracy were calculated to evaluate the genome sequence quality of 32 plants. The integrity of a genome sequence is presented by the ratio of chromosome size and genome size (or between scaffold size and genome size), which ranged from 55.31% to nearly 100%. The accuracy of genome sequence was presented by the ratio between matched EST and selected ESTs where 52.93% ∼ 98.28% and 89.02% ∼ 98.85% of the randomly selected clean ESTs could be mapped to chromosome and scaffold sequences, respectively. According to the integrity, accuracy and other analysis of each plant species, thirteen plant species were divided into four levels. Arabidopsis thaliana, Oryza sativa and Zea mays had the highest quality, followed by Brachypodium distachyon, Populus trichocarpa, Vitis vinifera and Glycine max, Sorghum bicolor, Solanum lycopersicum and Fragaria vesca, and Lotus japonicus, Medicago truncatula and Malus × domestica in that order. Assembling the scaffold sequences into chromosome sequences should be the primary task for the remaining nineteen species. Low GC content and repeat DNA influences genome sequence assembly.

Conclusion

The quality of plant genome sequences was found to be lower than envisaged and thus the rapid development of genome sequencing projects as well as research on bioinformatics tools and the algorithms of genome sequence assembly should provide increased processing and correction of genome sequences that have already been published.     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

Genome Sequence of the Versatile Fish Pathogen Edwardsiella tarda Provides Insights into its Adaptation to Broad Host Ranges and Intracellular Niches

Background

Edwardsiella tarda is the etiologic agent of edwardsiellosis, a devastating fish disease prevailing in worldwide aquaculture industries. Here we describe the complete genome of E. tarda, EIB202, a highly virulent and multi-drug resistant isolate in China.

Methodology/Principal Findings

E. tarda EIB202 possesses a single chromosome of 3,760,463 base pairs containing 3,486 predicted protein coding sequences, 8 ribosomal rRNA operons, and 95 tRNA genes, and a 43,703 bp conjugative plasmid harboring multi-drug resistant determinants and encoding type IV A secretion system components. We identified a full spectrum of genetic properties related to its genome plasticity such as repeated sequences, insertion sequences, phage-like proteins, integrases, recombinases and genomic islands. In addition, analysis also indicated that a substantial proportion of the E. tarda genome might be devoted to the growth and survival under diverse conditions including intracellular niches, with a large number of aerobic or anaerobic respiration-associated proteins, signal transduction proteins as well as proteins involved in various stress adaptations. A pool of genes for secretion systems, pili formation, nonfimbrial adhesions, invasions and hemagglutinins, chondroitinases, hemolysins, iron scavenging systems as well as the incomplete flagellar biogenesis might feature its surface structures and pathogenesis in a fish body.

Conclusion/Significance

Genomic analysis of the bacterium offered insights into the phylogeny, metabolism, drug-resistance, stress adaptation, and virulence characteristics of this versatile pathogen, which constitutes an important first step in understanding the pathogenesis of E. tarda to facilitate construction of a practical effective vaccine used for combating fish edwardsiellosis.     

Molecular typing of Treponema pallidum: a systematic review and meta-analysis

Background

Syphilis is resurgent in many regions of the world. Molecular typing is a robust tool for investigating strain diversity and epidemiology. This study aimed to review original research on molecular typing of Treponema pallidum (T. pallidum) with three objectives: (1) to determine specimen types most suitable for molecular typing; (2) to determine T. pallidum subtype distribution across geographic areas; and (3) to summarize available information on subtypes associated with neurosyphilis and macrolide resistance.

Methodology/Principal Findings

Two researchers independently searched five databases from 1998 through 2010, assessed for eligibility and study quality, and extracted data. Search terms included “Treponema pallidum,” or “syphilis,” combined with the subject headings “molecular,” “subtyping,” “typing,” “genotype,” and “epidemiology.” Sixteen eligible studies were included. Publication bias was not statistically significant by the Begg rank correlation test. Medians, inter-quartile ranges, and 95% confidence intervals were determined for DNA extraction and full typing efficiency. A random-effects model was used to perform subgroup analyses to reduce obvious between-study heterogeneity. Primary and secondary lesions and ear lobe blood specimens had an average higher yield of T. pallidum DNA (83.0% vs. 28.2%, χ² = 247.6, p<0.001) and an average higher efficiency of full molecular typing (80.9% vs. 43.1%, χ² = 102.3, p<0.001) compared to plasma, whole blood, and cerebrospinal fluid. A pooled analysis of subtype distribution based on country location showed that 14d was the most common subtype, and subtype distribution varied across geographic areas. Subtype data associated with macrolide resistance and neurosyphilis were limited.

Conclusions/Significance

Primary lesion was a better specimen for obtaining T. pallidum DNA than blood. There was wide geographic variation in T. pallidum subtypes. More research is needed on the relationship between clinical presentation and subtype, and further validation of ear lobe blood for obtaining T. pallidum DNA would be useful for future molecular studies of syphilis.     

Genome sequence of ground tit Pseudopodoces humilis and its adaptation to high altitude

Background

The mechanism of high-altitude adaptation has been studied in certain mammals. However, in avian species like the ground tit Pseudopodoces humilis, the adaptation mechanism remains unclear. The phylogeny of the ground tit is also controversial.

Results

Using next generation sequencing technology, we generated and assembled a draft genome sequence of the ground tit. The assembly contained 1.04 Gb of sequence that covered 95.4% of the whole genome and had higher N50 values, at the level of both scaffolds and contigs, than other sequenced avian genomes. About 1.7 million SNPs were detected, 16,998 protein-coding genes were predicted and 7% of the genome was identified as repeat sequences. Comparisons between the ground tit genome and other avian genomes revealed a conserved genome structure and confirmed the phylogeny of ground tit as not belonging to the Corvidae family. Gene family expansion and positively selected gene analysis revealed genes that were related to cardiac function. Our findings contribute to our understanding of the adaptation of this species to extreme environmental living conditions.

Conclusions

Our data and analysis contribute to the study of avian evolutionary history and provide new insights into the adaptation mechanisms to extreme conditions in animals.     

Genome analysis reveals three genomospecies in Mycobacterium abscessus

Background

Mycobacterium abscessus complex, the third most frequent mycobacterial complex responsible for community- and health care-associated infections in developed countries, comprises of M. abscessus subsp. abscessus and M. abscessus subsp. bolletii reviously referred as Mycobacterium bolletii and Mycobacterium massiliense. The diversity of this group of opportunistic pathogens is poorly described.

Results

In-depth analysis of 14 published M. abscessus complex genomes found a pan-genome of 6,153 proteins and core-genome of 3,947 (64.1%) proteins, indicating a non-conservative genome. Analysing the average percentage of amino-acid sequence identity (from 94.19% to 98.58%) discriminates three main clusters C1, C2 and C3: C1 comprises strains belonging to M. abscessus, C2 comprises strains belonging to M. massiliense and C3 comprises strains belonging to M. bolletii; and two sub-clusters in clusters C2 and C3. The phylogenomic network confirms these three clusters. The genome length (from 4.8 to 5.51-Mb) varies from 5.07-Mb in C1, 4.89-Mb in C2A, 5.01-Mb in C2B and 5.28-Mb in C3. The mean number of prophage regions (from 0 to 7) is 2 in C1; 1.33 in C2A; 3.5 in C2B and five in C3. A total of 36 genes are uniquely present in C1, 15 in C2 and 15 in C3. These genes could be used for the detection and identification of organisms in each cluster. Further, the mean number of host-interaction factors (including PE, PPE, LpqH, MCE, Yrbe and type VII secretion system ESX3 and ESX4) varies from 70 in cluster C1, 80 in cluster C2A, 74 in cluster C2B and 93 in clusters C3A and C3B. No significant differences in antibiotic resistance genes were observed between clusters, in contrast to previously reported in-vitro patterns of drug resistance. They encode both penicillin-binding proteins targeted by β-lactam antibiotics and an Ambler class A β-lactamase for which inhibitors exist.

Conclusions

Our comparative analysis indicates that M. abscessus complex comprises three genomospecies, corresponding to M. abscessus, M. bolletii, and M. massiliense. The genomics data here reported indicate differences in virulence of medical interest; and suggest targets for the refined detection and identification of M. abscessus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-359) contains supplementary material, which is available to authorized users.     

Genetic engineering of Treponema pallidum subsp. pallidum,the Syphilis Spirochete

Despite more than a century of research, genetic manipulation of Treponema pallidum subsp. pallidum (T. pallidum), the causative agent of syphilis, has not been successful. The lack of genetic engineering tools has severely limited understanding of the mechanisms behind T. pallidum success as a pathogen. A recently described method for in vitro cultivation of T. pallidum, however, has made it possible to experiment with transformation and selection protocols in this pathogen. Here, we describe an approach that successfully replaced the tprA (tp0009) pseudogene in the SS14 T. pallidum strain with a kanamycin resistance (kan^R) cassette. A suicide vector was constructed using the pUC57 plasmid backbone. In the vector, the kan^R gene was cloned downstream of the tp0574 gene promoter. The tp0574prom-kan^R cassette was then placed between two 1-kbp homology arms identical to the sequences upstream and downstream of the tprA pseudogene. To induce homologous recombination and integration of the kan^R cassette into the T. pallidum chromosome, in vitro-cultured SS14 strain spirochetes were exposed to the engineered vector in a CaCl₂-based transformation buffer and let recover for 24 hours before adding kanamycin-containing selective media. Integration of the kan^R cassette was demonstrated by qualitative PCR, droplet digital PCR (ddPCR), and whole-genome sequencing (WGS) of transformed treponemes propagated in vitro and/or in vivo. ddPCR analysis of RNA and mass spectrometry confirmed expression of the kan^R message and protein in treponemes propagated in vitro. Moreover, tprA knockout (tprA^ko-SS14) treponemes grew in kanamycin concentrations that were 64 times higher than the MIC for the wild-type SS14 (wt-SS14) strain and in infected rabbits treated with kanamycin. We demonstrated that genetic manipulation of T. pallidum is attainable. This discovery will allow the application of functional genetics techniques to study syphilis pathogenesis and improve syphilis vaccine development.     

Use of Whole Genome Sequencing to Determine the Microevolution of Mycobacterium tuberculosis during an Outbreak

Rationale

Current tools available to study the molecular epidemiology of tuberculosis do not provide information about the directionality and sequence of transmission for tuberculosis cases occurring over a short period of time, such as during an outbreak. Recently, whole genome sequencing has been used to study molecular epidemiology of Mycobacterium tuberculosis over short time periods.

Objective

To describe the microevolution of M. tuberculosis during an outbreak caused by one drug-susceptible strain.

Method and Measurements

We included 9 patients with tuberculosis diagnosed during a period of 22 months, from a population-based study of the molecular epidemiology in San Francisco. Whole genome sequencing was performed using Illumina’s sequencing by synthesis technology. A custom program written in Python was used to determine single nucleotide polymorphisms which were confirmed by PCR product Sanger sequencing.

Main results

We obtained an average of 95.7% (94.1–96.9%) coverage for each isolate and an average fold read depth of 73 (1 to 250). We found 7 single nucleotide polymorphisms among the 9 isolates. The single nucleotide polymorphisms data confirmed all except one known epidemiological link. The outbreak strain resulted in 5 bacterial variants originating from the index case A1 with 0–2 mutations per transmission event that resulted in a secondary case.

Conclusions

Whole genome sequencing analysis from a recent outbreak of tuberculosis enabled us to identify microevolutionary events observable during transmission, to determine 0–2 single nucleotide polymorphisms per transmission event that resulted in a secondary case, and to identify new epidemiologic links in the chain of transmission.     

A Complete Mitochondrial Genome Sequence from a Mesolithic Wild Aurochs (Bos primigenius)

Background

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

Methodology

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738±68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

Conclusions

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.     

GenomeFingerprinter: The Genome Fingerprint and the Universal Genome Fingerprint Analysis for Systematic Comparative Genomics

Background

No attention has been paid on comparing a set of genome sequences crossing genetic components and biological categories with far divergence over large size range. We define it as the systematic comparative genomics and aim to develop the methodology.

Results

First, we create a method, GenomeFingerprinter, to unambiguously produce a set of three-dimensional coordinates from a sequence, followed by one three-dimensional plot and six two-dimensional trajectory projections, to illustrate the genome fingerprint of a given genome sequence. Second, we develop a set of concepts and tools, and thereby establish a method called the universal genome fingerprint analysis (UGFA). Particularly, we define the total genetic component configuration (TGCC) (including chromosome, plasmid, and phage) for describing a strain as a systematic unit, the universal genome fingerprint map (UGFM) of TGCC for differentiating strains as a universal system, and the systematic comparative genomics (SCG) for comparing a set of genomes crossing genetic components and biological categories. Third, we construct a method of quantitative analysis to compare two genomes by using the outcome dataset of genome fingerprint analysis. Specifically, we define the geometric center and its geometric mean for a given genome fingerprint map, followed by the Euclidean distance, the differentiate rate, and the weighted differentiate rate to quantitatively describe the difference between two genomes of comparison. Moreover, we demonstrate the applications through case studies on various genome sequences, giving tremendous insights into the critical issues in microbial genomics and taxonomy.

Conclusions

We have created a method, GenomeFingerprinter, for rapidly computing, geometrically visualizing, intuitively comparing a set of genomes at genome fingerprint level, and hence established a method called the universal genome fingerprint analysis, as well as developed a method of quantitative analysis of the outcome dataset. These have set up the methodology of systematic comparative genomics based on the genome fingerprint analysis.