Whole-Genome Enrichment Using RNA Probes and Sequencing of Chlamydia trachomatis Directly from Clinical Samples

Contaminant DNA in bacterial sequencing experiments is a major source of false genetic variability

Background

Contaminant DNA is a well-known confounding factor in molecular biology and in genomic repositories. Strikingly, analysis workflows for whole-genome sequencing (WGS) data commonly do not account for errors potentially introduced by contamination, which could lead to the wrong assessment of allele frequency both in basic and clinical research.

Results

We used a taxonomic filter to remove contaminant reads from more than 4000 bacterial samples from 20 different studies and performed a comprehensive evaluation of the extent and impact of contaminant DNA in WGS. We found that contamination is pervasive and can introduce large biases in variant analysis. We showed that these biases can result in hundreds of false positive and negative SNPs, even for samples with slight contamination. Studies investigating complex biological traits from sequencing data can be completely biased if contamination is neglected during the bioinformatic analysis, and we demonstrate that removing contaminant reads with a taxonomic classifier permits more accurate variant calling. We used both real and simulated data to evaluate and implement reliable, contamination-aware analysis pipelines.

Conclusion

As sequencing technologies consolidate as precision tools that are increasingly adopted in the research and clinical context, our results urge for the implementation of contamination-aware analysis pipelines. Taxonomic classifiers are a powerful tool to implement such pipelines.

Whole-genome sequencing of ocular Chlamydia trachomatis isolates from Gadarif State, Sudan

BACKGROUND

Trachoma, caused by ocular Chlamydia trachomatis, is the leading infectious cause of blindness worldwide. Sudan first reported trachoma in the 1930s and has since been consistently endemic. Ocular C. trachomatis previously isolated from trachoma patients in Sudan in 1963 was antigenically identical to an isolate from Saudi Arabia (A/SA1). No contemporary ocular C. trachomatis whole genome sequences have been reported from Sudan.

METHODS

This study sequenced twenty ocular C. trachomatis isolates to improve understanding of pathogen diversity in North-East Africa and examine for genomic variation specific to Sudan, possibly related to the persistence of trachoma in surveyed communities. High quality, whole genome sequences were obtained from 12/20 isolates.

RESULTS

All isolates were serovar A and had tarP and trpA sequences typical of classical, ocular C. trachomatis isolates. The Sudanese isolates formed a closely related subclade within the T2-trachoma clade of C. trachomatis phylogeny distinct from geographically disparate ocular isolates, with little intra-population diversity. We found 333 SNPs that were conserved in Sudanese ocular isolates but rare compared to other ocular C. trachomatis populations, which were focused in two genomic loci (CTA0172-CTA0173 and CTA0482).

CONCLUSIONS

Limited intra-population diversity and geographical clustering of ocular C. trachomatis suggests minimal transmission between and slow diversification within trachoma-endemic communities. However, diversity may have been higher pre-treatment in these communities. Over-representation of Sudan-specific SNPs in three genes suggests they may have an impact on C. trachomatis growth and transmission in this population.