Comparison of long- and short-read metagenomic assembly for low-abundance species and resistance genes

Recent technological and computational advances have made metagenomic assembly a viable approach to achieving high-resolution views of complex microbial communities. In previous benchmarking, short-read (SR) metagenomic assemblers had the highest accuracy, long-read (LR) assemblers generated the most contiguous sequences and hybrid (HY) assemblers balanced length and accuracy. However, no assessments have specifically compared the performance of these assemblers on low-abundance species, which include clinically relevant organisms in the gut. We generated semi-synthetic LR and SR datasets by spiking small and increasing amounts of Escherichia coli isolate reads into fecal metagenomes and, using different assemblers, examined E. coli contigs and the presence of antibiotic resistance genes (ARGs). For ARG assembly, although SR assemblers recovered more ARGs with high accuracy, even at low coverages, LR assemblies allowed for the placement of ARGs within longer, E. coli-specific contigs, thus pinpointing their taxonomic origin. HY assemblies identified resistance genes with high accuracy and had lower contiguity than LR assemblies. Each assembler type's strengths were maintained even when our isolate was spiked in with a competing strain, which fragmented and reduced the accuracy of all assemblies. For strain characterization and determining gene context, LR assembly is optimal, while for base-accurate gene identification, SR assemblers outperform other options. HY assembly offers contiguity and base accuracy, but requires generating data on multiple platforms, and may suffer high misassembly rates when strain diversity exists. Our results highlight the trade-offs associated with each approach for recovering low-abundance taxa, and that the optimal approach is goal-dependent.

Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are an urgent global health threat. Inferring the dynamics of local CRE dissemination is currently limited by our inability to confidently trace the spread of resistance determinants to unrelated bacterial hosts. Whole-genome sequence comparison is useful for identifying CRE clonal transmission and outbreaks, but high-frequency horizontal gene transfer (HGT) of carbapenem resistance genes and subsequent genome rearrangement complicate tracing the local persistence and mobilization of these genes across organisms.

METHODS

To overcome this limitation, we developed a new approach to identify recent HGT of large, near-identical plasmid segments across species boundaries, which also allowed us to overcome technical challenges with genome assembly. We applied this to complete and near-complete genome assemblies to examine the local spread of CRE in a systematic, prospective collection of all CRE, as well as time- and species-matched carbapenem-susceptible Enterobacterales, isolated from patients from four US hospitals over nearly 5 years.

RESULTS

Our CRE collection comprised a diverse range of species, lineages, and carbapenem resistance mechanisms, many of which were encoded on a variety of promiscuous plasmid types. We found and quantified rearrangement, persistence, and repeated transfer of plasmid segments, including those harboring carbapenemases, between organisms over multiple years. Some plasmid segments were found to be strongly associated with specific locales, thus representing geographic signatures that make it possible to trace recent and localized HGT events. Functional analysis of these signatures revealed genes commonly found in plasmids of nosocomial pathogens, such as functions required for plasmid retention and spread, as well survival against a variety of antibiotic and antiseptics common to the hospital environment.

CONCLUSIONS

Collectively, the framework we developed provides a clearer, high-resolution picture of the epidemiology of antibiotic resistance importation, spread, and persistence in patients and healthcare networks.

StrainGE: a toolkit to track and characterize low-abundance strains in complex microbial communities

Human-associated microbial communities comprise not only complex mixtures of bacterial species, but also mixtures of conspecific strains, the implications of which are mostly unknown since strain level dynamics are underexplored due to the difficulties of studying them. We introduce the Strain Genome Explorer (StrainGE) toolkit, which deconvolves strain mixtures and characterizes component strains at the nucleotide level from short-read metagenomic sequencing with higher sensitivity and resolution than other tools. StrainGE is able to identify strains at 0.1x coverage and detect variants for multiple conspecific strains within a sample from coverages as low as 0.5x.

Limited effects of long-term daily cranberry consumption on the gut microbiome in a placebo-controlled study of women with recurrent urinary tract infections

BACKGROUND

Urinary tract infections (UTIs) affect 15 million women each year in the United States, with > 20% experiencing frequent recurrent UTIs. A recent placebo-controlled clinical trial found a 39% reduction in UTI symptoms among recurrent UTI sufferers who consumed a daily cranberry beverage for 24 weeks. Using metagenomic sequencing of stool from a subset of these trial participants, we assessed the impact of cranberry consumption on the gut microbiota, a reservoir for UTI-causing pathogens such as Escherichia coli, which causes > 80% of UTIs.

RESULTS

The overall taxonomic composition, community diversity, carriage of functional pathways and gene families, and relative abundances of the vast majority of observed bacterial taxa, including E. coli, were not changed significantly by cranberry consumption. However, one unnamed Flavonifractor species (OTU41), which represented ≤1% of the overall metagenome, was significantly less abundant in cranberry consumers compared to placebo at trial completion. Given Flavonifractor's association with negative human health effects, we sought to determine OTU41 characteristic genes that may explain its differential abundance and/or relationship to key host functions. Using comparative genomic and metagenomic techniques, we identified genes in OTU41 related to transport and metabolism of various compounds, including tryptophan and cobalamin, which have been shown to play roles in host-microbe interactions.

CONCLUSION

While our results indicated that cranberry juice consumption had little impact on global measures of the microbiome, we found one unnamed Flavonifractor species differed significantly between study arms. This suggests further studies are needed to assess the role of cranberry consumption and Flavonifractor in health and wellbeing in the context of recurrent UTI.

TRIAL REGISTRATION

Clinical trial registration number: ClinicalTrials.gov NCT01776021 .

QuantTB - a method to classify mixed Mycobacterium tuberculosis infections within whole genome sequencing data

BACKGROUND

Mixed infections of Mycobacterium tuberculosis and antibiotic heteroresistance continue to complicate tuberculosis (TB) diagnosis and treatment. Detection of mixed infections has been limited to molecular genotyping techniques, which lack the sensitivity and resolution to accurately estimate the multiplicity of TB infections. In contrast, whole genome sequencing offers sensitive views of the genetic differences between strains of M. tuberculosis within a sample. Although metagenomic tools exist to classify strains in a metagenomic sample, most tools have been developed for more divergent species, and therefore cannot provide the sensitivity required to disentangle strains within closely related bacterial species such as M. tuberculosis. Here we present QuantTB, a method to identify and quantify individual M. tuberculosis strains in whole genome sequencing data. QuantTB uses SNP markers to determine the combination of strains that best explain the allelic variation observed in a sample. QuantTB outputs a list of identified strains, their corresponding relative abundances, and a list of drugs for which resistance-conferring mutations (or heteroresistance) have been predicted within the sample.

RESULTS

We show that QuantTB has a high degree of resolution and is capable of differentiating communities differing by less than 25 SNPs and identifying strains down to 1× coverage. Using simulated data, we found QuantTB outperformed other metagenomic strain identification tools at detecting strains and quantifying strain multiplicity. In a real-world scenario, using a dataset of 50 paired clinical isolates from a study of patients with either reinfections or relapses, we found that QuantTB could detect mixed infections and reinfections at rates concordant with a manually curated approach.

CONCLUSION

QuantTB can determine infection multiplicity, identify hetero-resistance patterns, enable differentiation between relapse and re-infection, and clarify transmission events across seemingly unrelated patients - even in low-coverage (1×) samples. QuantTB outperforms existing tools and promises to serve as a valuable resource for both clinicians and researchers working with clinical TB samples.

1202. Multimodal Sequencing of a Clonal Case Cluster of Carbapenem-Resistant Citrobacter Reveals Unexpectedly Rapid Dynamics of KPC3-Containing Plasmids

Background

Carbapenem-resistant Enterobacteriaceae (CRE) are a major public health threat. We report four clonally related Citrobacter freundii isolates harboring the bla_KPC-3 carbapenemase in April–May 2017 that are nearly identical to a strain from 2014 at the same institution. Despite differing by ≤5 single nucleotide polymorphisms (SNPs), these isolates exhibited dramatic differences in carbapenemase plasmid architecture.

Methods

We sequenced four carbapenem-resistant C. freundii isolates from 2017 and compared them with an ongoing CRE surveillance project at our institution. SNPs were identified from Illumina MiSeq data aligned to a reference genome using the variant caller Pilon. Plasmids were assembled from Illumina and Oxford Nanopore sequencing data using Unicycler.

Results

The four 2017 isolates differed from one another by 0–5 chromosomal SNPs; two were identical. With one exception, these isolates differed by >38,000 SNPs from 25 C. freundii isolates sequenced from 2013 to 2017 at the same institution for CRE surveillance. The exception was a 2014 isolate that differed by 13–16 SNPs from each 2017 isolate, with 13 SNPs common to all four. Each C. freundii isolate harbored wild-type bla_KPC-3. Despite the close relationship among the 2017 cluster, the plasmids harboring the bla_KPC-3 genes differed dramatically: the carbapenemase occurred in one of the two different plasmids, with rearrangements between these plasmids across isolates. The related 2014 isolate harbored both plasmids, each with a separate copy of bla_KPC-3. No transmission chains were found between any of the affected patients.

Conclusion

WGS confirmed clonality among four contemporaneous bla_KPC-3-containing C. freundii isolates, and marked similarity with a 2014 isolate, within an institution. That only 13–16 SNPs varied between the 2014 and 2017 isolates suggests durable persistence of the bla_KPC-3 gene within this lineage in a hospital ecosystem. The plasmids harboring these carbapenemase genes proved remarkably plastic, with plasmid loss and rearrangements occurring on the same time scale as two to three chromosomal point mutations. Combining short and long-read sequencing in a case cluster uniquely revealed unexpectedly rapid dynamics of carbapenemase plasmids, providing critical insight into their manner of spread.

Disclosures

M. J. Ferraro, SeLux Diagnostics: Scientific Advisor and Shareholder, Consulting fee. D. C. Hooper, SeLux Diagnostics: Scientific Advisor, Consulting fee.

Genomic and functional analyses of Mycobacterium tuberculosis strains implicate ald in D-cycloserine resistance

A more complete understanding of the genetic basis of drug resistance in Mycobacterium tuberculosis is critical for prompt diagnosis and optimal treatment, particularly for toxic second-line drugs such as D-cycloserine. Here we used the whole-genome sequences from 498 strains of M. tuberculosis to identify new resistance-conferring genotypes. By combining association and correlated evolution tests with strategies for amplifying signal from rare variants, we found that loss-of-function mutations in ald (Rv2780), encoding L-alanine dehydrogenase, were associated with unexplained drug resistance. Convergent evolution of this loss of function was observed exclusively among multidrug-resistant strains. Drug susceptibility testing established that ald loss of function conferred resistance to D-cycloserine, and susceptibility to the drug was partially restored by complementation of ald. Clinical strains with mutations in ald and alr exhibited increased resistance to D-cycloserine when cultured in vitro. Incorporation of D-cycloserine resistance in novel molecular diagnostics could allow for targeted use of this toxic drug among patients with susceptible infections.

Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal

BACKGROUND

The continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.

METHODS AND FINDINGS

We performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937-1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974-1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988-1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.

CONCLUSIONS

In the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.

Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3–5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.

Sequence-based discovery of Bradyrhizobium enterica in cord colitis syndrome

BACKGROUND

Immunosuppression is associated with a variety of idiopathic clinical syndromes that may have infectious causes. It has been hypothesized that the cord colitis syndrome, a complication of umbilical-cord hematopoietic stem-cell transplantation, is infectious in origin.

METHODS

We performed shotgun DNA sequencing on four archived, paraffin-embedded endoscopic colon-biopsy specimens obtained from two patients with cord colitis. Computational subtraction of human and known microbial sequences and assembly of residual sequences into a bacterial draft genome were performed. We used polymerase-chain-reaction (PCR) assays and fluorescence in situ hybridization to determine whether the corresponding bacterium was present in additional patients and controls.

RESULTS

DNA sequencing of the biopsy specimens revealed more than 2.5 million sequencing reads that did not match known organisms. These sequences were computationally assembled into a 7.65-Mb draft genome showing a high degree of homology with genomes of bacteria in the bradyrhizobium genus. The corresponding newly discovered bacterium was provisionally named Bradyrhizobium enterica. PCR identified B. enterica nucleotide sequences in biopsy specimens from all three additional patients with cord colitis whose samples were tested, whereas B. enterica sequences were absent in samples obtained from healthy controls and patients with colon cancer or graft-versus-host disease.

CONCLUSIONS

We assembled a novel bacterial draft genome from the direct sequencing of tissue specimens from patients with cord colitis. Association of these sequences with cord colitis suggests that B. enterica may be an opportunistic human pathogen. (Funded by the National Cancer Institute and others.)

Finished bacterial genomes from shotgun sequence data

Exceptionally accurate genome reference sequences have proven to be of great value to microbial researchers. Thus, to date, about 1800 bacterial genome assemblies have been “finished” at great expense with the aid of manual laboratory and computational processes that typically iterate over a period of months or even years. By applying a new laboratory design and new assembly algorithm to 16 samples, we demonstrate that assemblies exceeding finished quality can be obtained from whole-genome shotgun data and automated computation. Cost and time requirements are thus dramatically reduced.

Asymmetric preference of serine proteases toward phosphonate and phosphinate esters

We have previously reported the asymmetric synthesis of (alpha-aminoalkyl) diphenylphosphonate and phosphinate derivatives designed as inhibitors of chymotrypsin- and elastase-like proteases. This paper reports the first kinetic evaluation of individual epimers of the (alpha-aminoalkyl) diphenylphosphonates as inactivators of chymotrypsin, cathepsin G and neutrophil elastase (HNE). Results show that the (R)-epimers consistently function as more potent irreversible inactivators of their respective target proteases than the corresponding (S)-epimers. Additionally, phosphinate analogues were found to be consistently superior to their diphenylphosphonate counterparts. For example, Cbz. Phe(P)(OPh)-(CH(2))(2)-CO(2)Et inactivates cathepsin G approximately 45-fold more rapidly (k(i)/K(i) = 1.2 x 10(5) M(-1). min(-1)) than the analogous Cbz.Phe(P)(OPh)(2) (2.6 x 10(3) M(-1). min(-1)). Similarly, Cbz.Val(P)(OPh)-(CH(2))(2)-CO(2)Et was found to inactivate HNE some 3-fold more efficiently than Cbz.Val(P)(OPh)(2) (6.5 x 10(3) and 2.0 x 10(3) M(-1). min(-1), respectively).

Synthesis and proteinase inhibitory properties of diphenyl phosphonate analogues of aspartic and glutamic acids

The synthesis of diphenyl phosphonate analogues of aspartic and glutamic acid, and their inhibitory activity against S. aureus V8 protease and granzyme B, is described. The study has revealed difficulties with protecting group compatibility in the synthesis of these analogues. Two analogues, Acetyl. AspP (OPh)2 and Acetyl.GluP (OPh)2 were found to function as irreversible inactivators of V8 proteinase, yet exhibit no activity against granzyme B.

Absence of food effects on the pharmacokinetics of terfenadine

Administration of terfenadine (Seldane) immediately after a high fat breakfast slightly affects the rate but not the extent of absorption relative to fasting administration. Mean peak levels of the active metabolite were increased by 13 per cent but delayed by 0.9 h while AUC was virtually the same as when terfenadine was administered while fasting. Changes in rate of absorption may be due to delayed gastric emptying and more rapid terfenadine solubilization. In any case, these rate differences are unlikely to be clinically important in the absence of differences in extent of absorption.

Pharmacokinetics of terfenadine in healthy elderly subjects

The pharmacokinetics of the terfenadine active metabolite, metabolite I, was examined in ten healthy elderly adults and ten younger adults after single-dose oral administration of 120-mg terfenadine. All subjects successfully completed the study without reporting sedation or other adverse events. Absorption was rapid in both the young and elderly. The mean Cmax was the same for both groups, 501 ng/mL, and occurred at 2.3 hours in the young subjects and 2.5 hours in elderly subjects. However, the apparent clearance was reduced by about 25% in the elderly. After correcting clearance for bodyweight, this difference was not statistically significant.

Mutational analysis of the Drosophila miniature-dusky (m-dy) locus: effects on cell size and circadian rhythms

A mutational analysis has been performed to explore the function of the Drosophila melanogaster miniature-dusky (m-dy) locus. Mutations at this locus affect wing development, fertility and behavior. The genetic characterization of 13 different mutations suggests that m and dy variants are alleles of a single complex gene. All of these mutations alter wing size, apparently by reducing the volume of individual epidermal cells of the developing wing. In m mutants, epidermal cell boundaries persist in the mature wing, whereas they normally degenerate 1-2 hr after eclosion in wild-type or dy flies. This has permitted the direct visualization of cell size differences among several m mutants. Mutations at the m-dy locus also affect behavioral processes. Three out of nine dy alleles (dyn1, dyn3 and dyn4) lengthen the circadian period of the activity and eclosion rhythms by approximately 1.5 hr. In contrast, m mutants have normal circadian periods, but an abnormally large percentage of individuals express aperiodic bouts of activity. These behavior genetic studies also indicate that an existing "rhythm" mutation known as Andante is an allele of the m-dy locus. The differential effects of certain m-dy mutations on wing and behavioral phenotypes suggest that separable domains of function exist within this locus.

Somatostatin and water excretion in man: an intrarenal action

Intravenous infusion of somatostatin in six water loaded normal human subjects produced a prompt reduction in urine flow accompanied by a rise in urine osmolality and a decrease in free water clearance. Plasma AVP levels did not change. Somatostatin would seem to exert an antidiuretic effect directly on the kidney.

Assay of microalbuminuria using gel electroimmunoassay

The possibility of using electroimmunoassay as a technique for assessing microalbuminuria in the diabetic population has been studied. The method was found to be precise (given inter-batch coefficients of variation of 5.7% and 5.8% at levels of 16 mg/l and 23 mg/l respectively), showed adequate sensitivity, and produced results which correlated well (r = 0.971) with these obtained from a routine radioimmunoassay procedure. It is concluded that electroimmunoassay provides a reliable alternative to more sophisticated techniques in the non-specialized laboratory dealing with only moderate numbers of specimens.

Bioavailability of terfenadine in man

Fourteen normal male subjects were given either 60mg or 180mg of terfenadine suspension in a randomized two-way crossover study. Peak plasma concentrations of 1.544 +/- 0.726 (mean +/- S.D.) ng ml-1 were obtained in 0.786 h following the 60 mg dose and displayed an AUC or 11.864 +/- 3.369 ng h ml-1. Whereas peak plasma concentrations of 4.519 +/- 2.002 ng ml-1 in 1.071 +/- 0.514 h were obtained following the 180 mg dose. The AUC following the 180 mg dose was 44.341 +/- 22.041 ng h ml-1. When 60 mg of 14C terfenadine was given to six additional subjects, the peak plasma concentrations of 351 +/- 43 ng equivalents per ml were obtained in 1.67 +/- 0.41 h and the AUC was 2297.71 +/- 310.85 ng-equivalents h ml-1. This indicates that approximately 99.5 per cent of the terfenadine related material that is absorbed undergoes biotransformation. Urinary excretion of 14C accounted for 39.89 +/- 5.29 per cent of the dose while 60.58 +/- 2.44 per cent of the dose was recovered in the feces in twelve days. Thin-layer chromatographic (TLC) examination of fecal extracts showed only a trace of material chromatographing with terfenadine. This may indicate that the 14C present in the feces is not due to lack of absorption.