Exome Sequencing of a Blastomycosis Case-Control Cohort From Manitoba and Northwestern Ontario, Canada.

BACKGROUND: Blastomycosis is a pulmonary disease caused by Blastomyces spp., a group of pathogenic dimorphic fungi endemic to a number of geographic regions, specifically Manitoba and northwestern Ontario, Canada. Immunosuppression is a major risk factor affecting disease susceptibility, yet host immunity is not well understood. Genetic immunodeficiencies can also influence disease, with variants in IL6, GATA2 and VDBP shown to influence susceptibility. Additional genetic factors in disease susceptibility and severity remain undetected. Our study seeks to identify potential genetic risk factors in a blastomycosis case-control cohort from Manitoba and northwestern Ontario, Canada. METHODS: Exomes from 18 blastomycosis cases and 9 controls were sequenced, variants were identified and filtered for accuracy and quality. We performed candidate gene prioritisation and variant aggregation to identify genetic associations and explored the full exome dataset. RESULTS: Ninety-nine genetic variants in 42 candidate genes were identified in the exome dataset. No variants associated with susceptibility were identified in a single-variant analysis although two non-synonymous variants in TYK2 were enriched among cases suggesting a possible role in susceptibility. Gene-based association analysis found variants in TLR1 enriched in controls (p = 0.024) suggesting a possible protective effect. Gene cluster analysis identified genetic variants in genes of chromatin remodelling, proteasome and intraflagellar transport significantly enriched in cases (false discovery rates < 14%). CONCLUSIONS: The findings in this study show novel associations with blastomycosis susceptibility. A better understanding of host immunity and genetic predisposition to Blastomyces infection can help to inform clinical practice for improved outcomes.

A validated in-house assay for HIV drug resistance mutation surveillance from dried blood spot specimens.

Despite increasing scale-up of antiretroviral therapy (ART) coverage, challenges related to adherence and HIV drug resistance (HIVDR) remain. The high cost of HIVDR surveillance is a persistent challenge with implementation in resource-constrained settings. Dried blood spot (DBS) specimens have been demonstrated to be a feasible alternative to plasma or serum for HIVDR genotyping and are more suitable for lower resource settings. There is a need for affordable HIVDR genotyping assays which can amplify HIV-1 sequences from DBS specimens, particularly those with low viral loads, at a low cost. Here, we present an in-house assay capable of reliably amplifying HIV-1 protease and partial reverse transcriptase genes from DBS specimens, which covers the complete World Health Organization 2009 list of drug resistance mutations under surveillance. DBS specimens were prepared using whole blood spiked with HIV-1 at concentrations of 10,000, 5000, 1000, and 500 copies/mL (n=30 for each concentration). Specimens were tested in triplicate. A two-step approach was used consisting of cDNA synthesis followed by nested PCR. The limit of detection of the assay was calculated to be approximately 5000 (95% CI: 3200-10,700) copies/mL for the protease gene and 3600 (95% CI: 2200-10,000) copies/mL for reverse transcriptase. The assay was observed to be most sensitive with higher viral load specimens (97.8% [95% CI: 92.2-99.7]) for both protease and reverse transcriptase at 10,000 copies/mL with performance decreasing with the use of specimens with lower viral loads (46.7% [36.1-57.5] and 60.0% [49.1-70.2] at 500 copies/mL for protease and reverse transcriptase, respectively). Ultimately, this assay presents a promising opportunity for use in resource-constrained settings. Future work should involve validation under field conditions including sub-optimal storage conditions and preparation of DBS with fingerprick blood in order to accurately reflect real-world collection scenarios.

Probe Capture Enrichment Methods for HIV and HCV Genome Sequencing and Drug Resistance Genotyping.

Human immunodeficiency virus (HIV) infections remain a significant public health concern worldwide. Over the years, sophisticated sequencing technologies such as next-generation sequencing (NGS) have emerged and been utilized to monitor the spread of HIV drug resistance (HIVDR), identify HIV drug resistance mutations, and characterize transmission dynamics. Similar applications also apply to the Hepatitis C virus (HCV), another bloodborne viral pathogen with significant intra-host genetic diversity. Several advantages to using NGS over conventional Sanger sequencing include increased data throughput, scalability, cost-effectiveness when batched sample testing is performed, and sensitivity for quantitative detection of minority resistant variants. However, NGS alone may fail to detect genomes from pathogens present in low copy numbers. As with all sequencing platforms, the primary determinant in achieving quality sequencing data is the quality and quantity of the initial template input. Samples containing degraded RNA/DNA and/or low copy number have been a consistent sequencing challenge. To overcome this limitation probe capture enrichment is a method that has recently been employed to target, enrich, and sequence the genome of a pathogen present in low copies, and for compromised specimens that contain poor quality nucleic acids. It involves the hybridization of sequence-specific DNA or RNA probes to a target sequence, which is followed by an enrichment step via PCR to increase the number of copies of the targeted sequences after which the samples are subjected to NGS procedures. This method has been performed on pathogens such as bacteria, fungus, and viruses and allows for the sequencing of complete genomes, with high coverage. Post NGS, data analysis can be performed through various bioinformatics pipelines which can provide information on genetic diversity, genotype, virulence, and drug resistance. This article reviews how probe capture enrichment helps to increase the likelihood of sequencing HIV and HCV samples that contain low viral loads and/or are compromised.

Evaluation of the bioMérieux VIDAS HIV Duo Quick and Anti-HCV assays for dried blood spot based serosurveillance.

Serosurveillance is central to monitoring our progress towards HIV and HCV elimination targets proposed for 2030. However, serosurveillance systems are ineffective without reliable serological assays for the detection of HIV and HCV antibodies. Assays should also be compatible with dried blood spot (DBS) samples to facilitate biological sample collection. The VIDAS HIV Duo Quick and Anti-HCV assays are sold as reagents strips and processed by the automated VIDAS benchtop immunoanalyser. While both assays have shown excellent performance in serum and plasma, performance data in DBS samples is lacking. In our study, we evaluate the performance of the VIDAS HIV Duo Quick and Anti-HCV assays in DBS (n = 725) collected during a cross-sectional serosurvey (the Transitions study). The VIDAS HIV Duo quick had a sensitivity and specificity of 94.5% (95% CI 85.1%, 98.5%) and 95.7% (95% CI 93.9%, 97.0%) respectively. Likewise, the VIDAS Anti-HCV had a sensitivity and specificity of 95.6% (95% CI 91.6%, 97.8%) and 95.6% (95% CI 93.5%, 97.0%) respectively. These assays are unlikely to be helpful in low-prevalence settings due to sub-optimal performance, but their performance could likely be improved by optimizing DBS elution protocols which was, unfortunately, not possible during our study.

External Quality Assessment for Next-Generation Sequencing-Based HIV Drug Resistance Testing: Unique Requirements and Challenges.

Over the past decade, there has been an increase in the adoption of next generation sequencing (NGS) technologies for HIV drug resistance (HIVDR) testing. NGS far outweighs conventional Sanger sequencing as it has much higher throughput, lower cost when samples are batched and, most importantly, significantly higher sensitivities for variants present at low frequencies, which may have significant clinical implications. Despite the advantages of NGS, Sanger sequencing remains the gold standard for HIVDR testing, largely due to the lack of standardization of NGS-based HIVDR testing. One important aspect of standardization includes external quality assessment (EQA) strategies and programs. Current EQA for Sanger-based HIVDR testing includes proficiency testing where samples are sent to labs and the performance of the lab conducting such assays is evaluated. The current methods for Sanger-based EQA may not apply to NGS-based tests because of the fundamental differences in their technologies and outputs. Sanger-based genotyping reports drug resistance mutations (DRMs) data as dichotomous, whereas NGS-based HIVDR genotyping also reports DRMs as numerical data (percent abundance). Here we present an overview of the need to develop EQA for NGS-based HIVDR testing and some unique challenges that may be encountered.

Dry Panels Supporting External Quality Assessment Programs for Next Generation Sequencing-Based HIV Drug Resistance Testing.

External quality assessment (EQA) is a keystone element in the validation and implementation of next generation sequencing (NGS)-based HIV drug resistance testing (DRT). Software validation and evaluation is a critical element in NGS EQA programs. While the development, sharing, and adoption of wet lab protocols is coupled with the increasing access to NGS technology worldwide, rendering it easy to produce NGS data for HIV-DRT, bioinformatic data analysis remains a bottleneck for most of the diagnostic laboratories. Several computational tools have been made available, via free or commercial sources, to automate the conversion of raw NGS data into an actionable clinical report. Although different software platforms yield equivalent results when identical raw NGS datasets are analyzed for variations at higher abundance, discrepancies arise when variations at lower frequencies are considered. This implies that validation and performance assessment of the bioinformatics tools applied in NGS HIV-DRT is critical, and the origins of the observed discrepancies should be determined. Well-characterized reference NGS datasets with ground truth on the genotype composition at all examined loci and the exact frequencies of HIV variations they may harbor, so-called dry panels, would be essential in such cases. The strategic design and construction of such panels are challenging but imperative tasks in support of EQA programs for NGS-based HIV-DRT and the validation of relevant bioinformatics tools. Here, we present criteria that can guide the design of such dry panels, which were discussed in the Second International Winnipeg Symposium themed for EQA strategies for NGS HIVDR assays.

Development and Application of Performance Assessment Criteria for Next-Generation Sequencing-Based HIV Drug Resistance Assays.

Next-generation sequencing (NGS)-based HIV drug resistance (HIVDR) assays outperform conventional Sanger sequencing in scalability, sensitivity, and quantitative detection of minority resistance variants. Thus far, HIVDR assays have been applied primarily in research but rarely in clinical settings. One main obstacle is the lack of standardized validation and performance evaluation systems that allow regulatory agencies to benchmark and accredit new assays for clinical use. By revisiting the existing principles for molecular assay validation, here we propose a new validation and performance evaluation system that helps to both qualitatively and quantitatively assess the performance of an NGS-based HIVDR assay. To accomplish this, we constructed a 70-specimen proficiency test panel that includes plasmid mixtures at known ratios, viral RNA from infectious clones, and anonymized clinical specimens. We developed assessment criteria and benchmarks for NGS-based HIVDR assays and used these to assess data from five separate MiSeq runs performed in two experienced HIVDR laboratories. This proposed platform may help to pave the way for the standardization of NGS HIVDR assay validation and performance evaluation strategies for accreditation and quality assurance purposes in both research and clinical settings.

Multi-Laboratory Comparison of Next-Generation to Sanger-Based Sequencing for HIV-1 Drug Resistance Genotyping.

Next-generation sequencing (NGS) is increasingly used for HIV-1 drug resistance genotyping. NGS methods have the potential for a more sensitive detection of low-abundance variants (LAV) compared to standard Sanger sequencing (SS) methods. A standardized threshold for reporting LAV that generates data comparable to those derived from SS is needed to allow for the comparability of data from laboratories using NGS and SS. Ten HIV-1 specimens were tested in ten laboratories using Illumina MiSeq-based methods. The consensus sequences for each specimen using LAV thresholds of 5%, 10%, 15%, and 20% were compared to each other and to the consensus of the SS sequences (protease 4-99; reverse transcriptase 38-247). The concordance among laboratories' sequences at different thresholds was evaluated by pairwise sequence comparisons. NGS sequences generated using the 20% threshold were the most similar to the SS consensus (average 99.6% identity, range 96.1-100%), compared to 15% (99.4%, 88.5-100%), 10% (99.2%, 87.4-100%), or 5% (98.5%, 86.4-100%). The average sequence identity between laboratories using thresholds of 20%, 15%, 10%, and 5% was 99.1%, 98.7%, 98.3%, and 97.3%, respectively. Using the 20% threshold, we observed an excellent agreement between NGS and SS, but significant differences at lower thresholds. Understanding how variation in NGS methods influences sequence quality is essential for NGS-based HIV-1 drug resistance genotyping.

Performance comparison of next generation sequencing analysis pipelines for HIV-1 drug resistance testing.

Next generation sequencing (NGS) is a trending new standard for genotypic HIV-1 drug resistance (HIVDR) testing. Many NGS HIVDR data analysis pipelines have been independently developed, each with variable outputs and data management protocols. Standardization of such analytical methods and comparison of available pipelines are lacking, yet may impact subsequent HIVDR interpretation and other downstream applications. Here we compared the performance of five NGS HIVDR pipelines using proficiency panel samples from NIAID Virology Quality Assurance (VQA) program. Ten VQA panel specimens were genotyped by each of six international laboratories using their own in-house NGS assays. Raw NGS data were then processed using each of the five different pipelines including HyDRA, MiCall, PASeq, Hivmmer and DEEPGEN. All pipelines detected amino acid variants (AAVs) at full range of frequencies (1~100%) and demonstrated good linearity as compared to the reference frequency values. While the sensitivity in detecting low abundance AAVs, with frequencies between 1~20%, is less a concern for all pipelines, their specificity dramatically decreased at AAV frequencies <2%, suggesting that 2% threshold may be a more reliable reporting threshold for ensured specificity in AAV calling and reporting. More variations were observed among the pipelines when low abundance AAVs are concerned, likely due to differences in their NGS read quality control strategies. Findings from this study highlight the need for standardized strategies for NGS HIVDR data analysis, especially for the detection of minority HIVDR variants.

A MiSeq-HyDRA platform for enhanced HIV drug resistance genotyping and surveillance.

Conventional HIV drug resistance (HIVDR) genotyping utilizes Sanger sequencing (SS) methods, which are limited by low data throughput and the inability of detecting low abundant drug resistant variants (LADRVs). Here we present a next generation sequencing (NGS)-based HIVDR typing platform that leverages the advantages of Illumina MiSeq and HyDRA Web. The platform consists of a fully validated sample processing protocol and HyDRA web, an open web portal that allows automated customizable NGS-based HIVDR data processing. This platform was characterized and validated using a panel of HIV-spiked plasma representing all major HIV-1 subtypes, pedigreed plasmids, HIVDR proficiency specimens and clinical specimens. All examined major HIV-1 subtypes were consistently amplified at viral loads of ≥1,000 copies/ml. The gross error rate of this platform was determined at 0.21%, and minor variations were reliably detected down to 0.50% in plasmid mixtures. All HIVDR mutations identifiable by SS were detected by the MiSeq-HyDRA protocol, while LADRVs at frequencies of 1~15% were detected by MiSeq-HyDRA only. As compared to SS approaches, the MiSeq-HyDRA platform has several notable advantages including reduced cost and labour, and increased sensitivity for LADRVs, making it suitable for routine HIVDR monitoring for both patient care and surveillance purposes.

Bioinformatic data processing pipelines in support of next-generation sequencing-based HIV drug resistance testing: the Winnipeg Consensus.

INTRODUCTION: Next-generation sequencing (NGS) has several advantages over conventional Sanger sequencing for HIV drug resistance (HIVDR) genotyping, including detection and quantitation of low-abundance variants bearing drug resistance mutations (DRMs). However, the high HIV genomic diversity, unprecedented large volume of data, complexity of analysis and potential for error pose significant challenges for data processing. Several NGS analysis pipelines have been developed and used in HIVDR research; however, the absence of uniformity in data processing strategies results in lack of consistency and comparability of outputs from different pipelines. To fill this gap, an international symposium on bioinformatic strategies for NGS-based HIVDR testing was held in February 2018 in Winnipeg, Canada, convening laboratory scientists, bioinformaticians and clinicians involved in four recently developed, publicly available NGS HIVDR pipelines. The goal of this symposium was to establish a consensus on effective bioinformatic strategies for NGS data management and its use for HIVDR reporting. DISCUSSION: Essential functionalities of an NGS HIVDR pipeline were divided into five analytic blocks: (1) NGS read quality control (QC)/quality assurance (QA); (2) NGS read alignment and reference mapping; (3) HIV variant calling and variant QC; (4) NGS HIVDR reporting; and (5) extended data applications and additional considerations for data management. The consensuses reached among the participants on all major aspects of these blocks are summarized here. They encompass not only recommended data management and analysis strategies, but also detailed bioinformatic approaches that help ensure accuracy of the derived HIVDR analysis outputs for both research and potential clinical use. CONCLUSIONS: While NGS is being adopted more broadly in HIVDR testing laboratories, data processing is often a bottleneck hindering its generalized application. The proposed standardization of NGS read QC/QA, read alignment and reference mapping, variant calling and QC, HIVDR reporting and relevant data management strategies in this "Winnipeg Consensus" may serve as a starting guideline for NGS HIVDR data processing that informs the refinement of existing pipelines and those yet to be developed. Moreover, the bioinformatic strategies presented here may apply more broadly to NGS data analysis of microbes harbouring significant genomic diversity.

Pyruvate kinase from Chlamydia trachomatis is activated by fructose-2,6-bisphosphate.

Pyruvate kinase is the final regulatory point in the catabolic Embden-Meyerhoff-Parnas pathway, which controls the carbon flux of glycolytic intermediates and regulates the level of ATP in the cell. In a previous study, we identified, cloned and sequenced pyruvate kinase from the obligate intracellular bacterium Chlamydia trachomatis and demonstrated that the enzyme was active in crude extract. Here, we report the kinetic properties of highly purified C. trachomatis pyruvate kinase. The results indicate that C. trachomatis pyruvate kinase is 53.5 kDa with a pH optima of 7.3. Kinetic studies show that C. trachomatis pyruvate kinase requires both K+ and Mg2+ ions for activity, exhibits sigmoidal kinetics with respect to phosphoenolpyruvate and Michaelis-Menten kinetics with respect to ADP. In addition, C. trachomatis pyruvate kinase is able to use alternative nucleoside diphosphates as phosphate acceptors, although it shows the greatest activity with ADP. In contrast to other bacterial pyruvate kinases that are activated by AMP, our data show that AMP, in addition to ATP and GTP, inhibits C. trachomatis pyruvate kinase. Surprisingly, unlike any other known bacterial pyruvate kinase, C. trachomatis pyruvate kinase was allosterically activated by fructose-2,6-bisphosphate, an important regulatory metabolite that has only been reported in eukaryotes.