HAPHPIPE: Haplotype Reconstruction and Phylodynamics for Deep Sequencing of Intrahost Viral Populations.

Deep sequencing of viral populations using next-generation sequencing (NGS) offers opportunities to understand and investigate evolution, transmission dynamics, and population genetics. Currently, the standard practice for processing NGS data to study viral populations is to summarize all the observed sequences from a sample as a single consensus sequence, thus discarding valuable information about the intrahost viral molecular epidemiology. Furthermore, existing analytical pipelines may only analyze genomic regions involved in drug resistance, thus are not suited for full viral genome analysis. Here, we present HAPHPIPE, a HAplotype and PHylodynamics PIPEline for genome-wide assembly of viral consensus sequences and haplotypes. The HAPHPIPE protocol includes modules for quality trimming, error correction, de novo assembly, alignment, and haplotype reconstruction. The resulting consensus sequences, haplotypes, and alignments can be further analyzed using a variety of phylogenetic and population genetic software. HAPHPIPE is designed to provide users with a single pipeline to rapidly analyze sequences from viral populations generated from NGS platforms and provide quality output properly formatted for downstream evolutionary analyses.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses.

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

Insights Into the Resistome of Bovine Clinical Mastitis Microbiome, a Key Factor in Disease Complication.

Bovine clinical mastitis (CM) is one of the most prevalent diseases caused by a wide range of resident microbes. The emergence of antimicrobial resistance in CM bacteria is well-known, however, the genomic resistance composition (the resistome) at the microbiome-level is not well characterized. In this study, we applied whole metagenome sequencing (WMS) to characterize the resistome of the CM microbiome, focusing on antibiotics and metals resistance, biofilm formation (BF), and quorum sensing (QS) along with in vitro resistance assays of six selected pathogens isolated from the same CM samples. The WMS generated an average of 21.13 million reads (post-processing) from 25 CM samples that mapped to 519 bacterial strains, of which 30.06% were previously unreported. We found a significant (P = 0.001) association between the resistomes and microbiome composition with no association with cattle breed, despite significant differences in microbiome diversity among breeds. The in vitro investigation determined that 76.2% of six selected pathogens considered "biofilm formers" actually formed biofilms and were also highly resistant to tetracycline, doxycycline, nalidixic acid, ampicillin, and chloramphenicol and remained sensitive to metals (Cr, Co, Ni, Cu, Zn) at varying concentrations. We also found bacterial flagellar movement and chemotaxis, regulation and cell signaling, and oxidative stress to be significantly associated with the pathophysiology of CM. Thus, identifying CM microbiomes, and analyzing their resistomes and genomic potentials will help improve the optimization of therapeutic schemes involving antibiotics and/or metals usage in the prevention and control of bovine CM.

Drug Resistance Prediction Using Deep Learning Techniques on HIV-1 Sequence Data.

The fast replication rate and lack of repair mechanisms of human immunodeficiency virus (HIV) contribute to its high mutation frequency, with some mutations resulting in the evolution of resistance to antiretroviral therapies (ART). As such, studying HIV drug resistance allows for real-time evaluation of evolutionary mechanisms. Characterizing the biological process of drug resistance is also critically important for sustained effectiveness of ART. Investigating the link between "black box" deep learning methods applied to this problem and evolutionary principles governing drug resistance has been overlooked to date. Here, we utilized publicly available HIV-1 sequence data and drug resistance assay results for 18 ART drugs to evaluate the performance of three architectures (multilayer perceptron, bidirectional recurrent neural network, and convolutional neural network) for drug resistance prediction, jointly with biological analysis. We identified convolutional neural networks as the best performing architecture and displayed a correspondence between the importance of biologically relevant features in the classifier and overall performance. Our results suggest that the high classification performance of deep learning models is indeed dependent on drug resistance mutations (DRMs). These models heavily weighted several features that are not known DRM locations, indicating the utility of model interpretability to address causal relationships in viral genotype-phenotype data.

Author Correction: A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Evaluation of haplotype callers for next-generation sequencing of viruses.

Currently, the standard practice for assembling next-generation sequencing (NGS) reads of viral genomes is to summarize thousands of individual short reads into a single consensus sequence, thus confounding useful intra-host diversity information for molecular phylodynamic inference. It is hypothesized that a few viral strains may dominate the intra-host genetic diversity with a variety of lower frequency strains comprising the rest of the population. Several software tools currently exist to convert NGS sequence variants into haplotypes. Previous benchmarks of viral haplotype reconstruction programs used simulation scenarios that are useful from a mathematical perspective but do not reflect viral evolution and epidemiology. Here, we tested twelve NGS haplotype reconstruction methods using viral populations simulated under realistic evolutionary dynamics. We simulated coalescent-based populations that spanned known levels of viral genetic diversity, including mutation rates, sample size and effective population size, to test the limits of the haplotype reconstruction methods and to ensure coverage of predicted intra-host viral diversity levels (especially HIV-1). All twelve investigated haplotype callers showed variable performance and produced drastically different results that were mainly driven by differences in mutation rate and, to a lesser extent, in effective population size. Most methods were able to accurately reconstruct haplotypes when genetic diversity was low. However, under higher levels of diversity (e.g., those seen intra-host HIV-1 infections), haplotype reconstruction quality was highly variable and, on average, poor. All haplotype reconstruction tools, except QuasiRecomb and ShoRAH, greatly underestimated intra-host diversity and the true number of haplotypes. PredictHaplo outperformed, in regard to highest precision, recall, and lowest UniFrac distance values, the other haplotype reconstruction tools followed by CliqueSNV, which, given more computational time, may have outperformed PredictHaplo. Here, we present an extensive comparison of available viral haplotype reconstruction tools and provide insights for future improvements in haplotype reconstruction tools using both short-read and long-read technologies.

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing.

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

Corrigendum: A 28-Year History of HIV-1 Drug Resistance and Transmission in Washington, DC.

[This corrects the article DOI: 10.3389/fmicb.2019.00369.].

Fecal Transplant in Children With Clostridioides difficile Gives Sustained Reduction in Antimicrobial Resistance and Potential Pathogen Burden.

BACKGROUND: Fecal microbiota transplantation (FMT) treats Clostridioides difficile infection (CDI). Little is known regarding the changes in antimicrobial resistance (AMR) genes and potential pathogen burden that occur in pediatric recipients of FMT. The aim of this study was to investigate changes in AMR genes, potential pathogens, species, and functional pathways with FMT in children. METHODS: Nine children with recurrent CDI underwent FMT. Stool was collected from donor and recipient pre-FMT and longitudinally post-FMT for up to 24 weeks. Shotgun metagenomic sequencing was performed. Reads were analyzed using PathoScope 2.0. RESULTS: All children had resolution of CDI. AMR genes decreased post-FMT (P < .001), with a sustained decrease in multidrug resistance genes (P < .001). Tetracycline resistance genes increased post-FMT (P < .001). Very low levels of potential pathogens were identified in donors and recipients, with an overall decrease post-FMT (P < .001). Prevotella sp. 109 expanded in all recipients post-FMT, and no recipients had any clinical infection. Alpha diversity was lower in recipients vs donors pre-FMT (P < .001), with an increase post-FMT (P ≤ .002) that was sustained. Beta diversity differed significantly in pre- vs post-FMT recipient samples (P < .001). Bacterial species Faecalibacterium prausnitzii and Bacteroides ovatus showed higher abundance in donors than recipients (P = .008 and P = .040, respectively), with expansion post-FMT. Biosynthetic pathways predominated in the donor and increased in the recipient post-FMT. CONCLUSIONS: FMT for CDI in children decreases AMR genes and potential pathogens and changes microbiota composition and function. However, acquisition of certain AMR genes post-FMT combined with low levels of potential pathogens found in donors suggests that further study is warranted regarding screening donors using metagenomics sequencing before FMT.

Gut microbiome differences between wild and captive black rhinoceros - implications for rhino health.

A number of recent studies have shown the importance of the mammalian gut microbiome in host health. In the context of endangered species, a few studies have examined the relationship between the gut microbiome in wild versus captive populations due to digestive and other health issues. Unfortunately, the results seem to vary across taxa in terms of captive animals having higher, lower, or equivalent microbiome diversity relative to their wild counterparts. Here, we focus on the black rhinoceros as captive animals suffer from a number of potentially dietary related health effects. We compared gut microbiomes of wild and captive black rhinos to test for differences in taxonomic diversity (alpha and beta) and in functional diversity of the microbiome. We incorporated a more powerful metagenomic shotgun sequencing approach rather than a targeted amplification of the 16S gene for taxonomic assignment of the microbiome. Our results showed no significant differences in the alpha diversity levels between wild and captive black rhinos, but significant differences in beta diversity. We found that bacterial taxa traditionally associated with ruminant guts of domesticated animals had higher relative abundances in captive rhinos. Our metagenomic sequencing results suggest that unknown gut microbes of wild rhinos are being replaced by those found in conventional human-domesticated livestock. Wild rhinos have significantly different functional bacterial communities compared to their captive counterparts. Functional profiling results showed greater abundance of glycolysis and amino acid synthesis pathways in captive rhino microbiomes, representing an animal receiving sub-optimal nutrition with a readily available source of glucose but possibly an imbalance of necessary macro and micronutrients. Given the differences observed between wild and captive rhino gut microbiomes, we make a number of recommendations for potentially modifying captive gut microbiome to better reflect their wild counterparts and thereby hopefully improve overall rhino health in captivity.

A 28-Year History of HIV-1 Drug Resistance and Transmission in Washington, DC.

Washington, DC consistently has one of the highest annual rates of new HIV-1 diagnoses in the United States over the last 10 years. To guide intervention and prevention strategies to combat DC HIV infection, it is helpful to understand HIV transmission dynamics in a historical context. Toward this aim, we conducted a retrospective study (years 1987-2015) of 3,349 HIV pol sequences (1,026 bp) from 1,996 individuals living in the DC area belonging to three different cohorts. We coupled HIV sequence data with clinical information (sex, risk factor, race/ethnicity, viral load, subtype, anti-retroviral regimen) to identify circulating drug resistant mutations (DRM) and transmission clusters and assess their persistence over time. Of the transmission clusters identified in the DC area, 78.0 and 31.7% involved MSM and heterosexuals, respectively. The longest spread of time for a single cluster was 5 years (2007-2012) using a distance-based network inference approach and 27 years (1987-2014) using a maximum likelihood phylogenetic approach. We found eight subtypes and nine recombinants. Genetic diversity increased steadily over time with a slight peak in 2009 and remained constant thereafter until 2015. Nucleotide diversity also increased over time while relative genetic diversity (BEAST) remained relatively steady over the last 28 years with slight increases since 2000 in subtypes B and C. Sequences from individuals on drug therapy contained the highest total number of DRMs (1,104-1,600) and unique DRMs (63-97) and the highest proportion (>20%) of resistant individuals. Heterosexuals (43.94%), MSM (40.13%), and unknown (44.26%) risk factors showed similar prevalence of DRMs, while injection drug users had a lower prevalence (33.33%). Finally, there was a 60% spike in the number of codons with DRMs between 2007 and 2010. Past patterns of HIV transmission and DRM accumulation over time described here will help to predict future efficacy of ART drugs based on DRMs persisting over time and identify risk groups of interest for prevention and intervention efforts within the DC population. Our results show how longitudinal data can help to understand the temporal dynamics of HIV-1 at the local level.

Benchmark Evaluation of True Single Molecular Sequencing to Determine Cystic Fibrosis Airway Microbiome Diversity.

Cystic fibrosis (CF) is an autosomal recessive disease associated with recurrent lung infections that can lead to morbidity and mortality. The impact of antibiotics for treatment of acute pulmonary exacerbations on the CF airway microbiome remains unclear with prior studies giving conflicting results and being limited by their use of 16S ribosomal RNA sequencing. Our primary objective was to validate the use of true single molecular sequencing (tSMS) and PathoScope in the analysis of the CF airway microbiome. Three control samples were created with differing amounts of Burkholderia cepacia, Pseudomonas aeruginosa, and Prevotella melaninogenica, three common bacteria found in cystic fibrosis lungs. Paired sputa were also obtained from three study participants with CF before and >6 days after initiation of antibiotics. Antibiotic resistant B. cepacia and P. aeruginosa were identified in concurrently obtained respiratory cultures. Direct sequencing was performed using tSMS, and filtered reads were aligned to reference genomes from NCBI using PathoScope and Kraken and unique clade-specific marker genes using MetaPhlAn. A total of 180-518 K of 6-12 million filtered reads were aligned for each sample. Detection of known pathogens in control samples was most successful using PathoScope. In the CF sputa, alpha diversity measures varied based on the alignment method used, but similar trends were found between pre- and post-antibiotic samples. PathoScope outperformed Kraken and MetaPhlAn in our validation study of artificial bacterial community controls and also has advantages over Kraken and MetaPhlAn of being able to determine bacterial strains and the presence of fungal organisms. PathoScope can be confidently used when evaluating metagenomic data to determine CF airway microbiome diversity.