Whole Genome Sequence Analysis of Cryptococcus gattii from the Pacific Northwest Reveals Unexpected Diversity

A recent emergence of Cryptococcus gattii in the Pacific Northwest involves strains that fall into three primarily clonal molecular subtypes: VGIIa, VGIIb and VGIIc. Multilocus sequence typing (MLST) and variable number tandem repeat analysis appear to identify little diversity within these molecular subtypes. Given the apparent expansion of these subtypes into new geographic areas and their ability to cause disease in immunocompetent individuals, differentiation of isolates belonging to these subtypes could be very important from a public health perspective. We used whole genome sequence typing (WGST) to perform fine-scale phylogenetic analysis on 20 C. gattii isolates, 18 of which are from the VGII molecular type largely responsible for the Pacific Northwest emergence. Analysis both including and excluding (289,586 SNPs and 56,845 SNPs, respectively) molecular types VGI and VGIII isolates resulted in phylogenetic reconstructions consistent, for the most part, with MLST analysis but with far greater resolution among isolates. The WGST analysis presented here resulted in identification of over 100 SNPs among eight VGIIc isolates as well as unique genotypes for each of the VGIIa, VGIIb and VGIIc isolates. Similar levels of genetic diversity were found within each of the molecular subtype isolates, despite the fact that the VGIIb clade is thought to have emerged much earlier. The analysis presented here is the first multi-genome WGST study to focus on the C. gattii molecular subtypes involved in the Pacific Northwest emergence and describes the tools that will further our understanding of this emerging pathogen

Comparing genomic variant identification protocols for Candida auris.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies

BurkDiff: A Real-Time PCR Allelic Discrimination Assay for Burkholderia Pseudomallei and B. mallei

A real-time PCR assay, BurkDiff, was designed to target a unique conserved region in the B. pseudomallei and B. mallei genomes containing a SNP that differentiates the two species. Sensitivity and specificity were assessed by screening BurkDiff across 469 isolates of B. pseudomallei, 49 isolates of B. mallei, and 390 isolates of clinically relevant non-target species. Concordance of results with traditional speciation methods and no cross-reactivity to non-target species show BurkDiff is a robust, highly validated assay for the detection and differentiation of B. pseudomallei and B. mallei

Tracing Genetic Exchange and Biogeography of Cryptococcus neoformans var. grubii at the Global Population Level.

Cryptococcus neoformans var. grubii is the causative agent of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals, typically HIV/AIDS patients from developing countries. Despite the worldwide emergence of this ubiquitous infection, little is known about the global molecular epidemiology of this fungal pathogen. Here we sequence the genomes of 188 diverse isolates and characterized the major subdivisions, their relative diversity and the level of genetic exchange between them. While most isolates of C. neoformans var. grubii belong to one of three major lineages (VNI, VNII, and VNB), some haploid isolates show hybrid ancestry including some that appear to have recently interbred, based on the detection of large blocks of each ancestry across each chromosome. Many isolates display evidence of aneuploidy, which was detected for all chromosomes. In diploid isolates of C. neoformans var. grubii (serotype A/A) and of hybrids with C. neoformans var. neoformans (serotype A/D) such aneuploidies have resulted in loss of heterozygosity, where a chromosomal region is represented by the genotype of only one parental isolate. Phylogenetic and population genomic analyses of isolates from Brazil reveal that the previously 'African' VNB lineage occurs naturally in the South American environment. This suggests migration of the VNB lineage between Africa and South America prior to its diversification, supported by finding ancestral recombination events between isolates from different lineages and regions. The results provide evidence of substantial population structure, with all lineages showing multi-continental distributions demonstrating the highly dispersive nature of this pathogen

Molecular Investigations of a Locally Acquired Case of Melioidosis in Southern AZ, USA

Melioidosis is caused by Burkholderia pseudomallei, a Gram-negative bacillus, primarily found in soils in Southeast Asia and northern Australia. A recent case of melioidosis in non-endemic Arizona was determined to be the result of locally acquired infection, as the patient had no travel history to endemic regions and no previous history of disease. Diagnosis of the case was confirmed through multiple microbiologic and molecular techniques. To enhance the epidemiological analysis, we conducted several molecular genotyping procedures, including multi-locus sequence typing, SNP-profiling, and whole genome sequence typing. Each technique has different molecular epidemiologic advantages, all of which provided evidence that the infecting strain was most similar to those found in Southeast Asia, possibly originating in, or around, Malaysia. Advancements in new typing technologies provide genotyping resolution not previously available to public health investigators, allowing for more accurate source identification

Ambient Stable Quantitative PCR Reagents for the Detection of Yersinia pestis

Plague, caused by Yersinia pestis, is one of the oldest and most dangerous diseases in human history, and has claimed millions of lives in the three major historical pandemics. Although panic caused by the Black Death is fading, the threat of the reemergence of plague pandemics still exists, with the additional potential of misuse in biowarfare or bioterrorism. Rapid on-site detection and identification of the pathogen is of paramount significance for timely implementation of effective countermeasures. TaqMan probe-based real-time PCR assays can give quick and accurate identification; however, the need for cold delivery and storage prevents its potential on-site application. The objective of this study was to develop a stable PCR system for easy delivery and storage under room temperature, which is vital for conventional plague surveillance and for preparedness in public health emergencies. We present a solution to this particular issue, hoping that it is helpful to future applications

A standardised method for interpreting the association between mutations and phenotypic drug resistance inMycobacterium tuberculosis

A clear understanding of the genetic basis of antibiotic resistance in Mycobacterium tuberculosis is required to accelerate the development of rapid drug susceptibility testing methods based on genetic sequence. Raw genotype–phenotype correlation data were extracted as part of a comprehensive systematic review to develop a standardised analytical approach for interpreting resistance associated mutations for rifampicin, isoniazid, ofloxacin/levofloxacin, moxifloxacin, amikacin, kanamycin, capreomycin, streptomycin, ethionamide/prothionamide and pyrazinamide. Mutation frequencies in resistant and susceptible isolates were calculated, together with novel statistical measures to classify mutations as high, moderate, minimal or indeterminate confidence for predicting resistance. We identified 286 confidence-graded mutations associated with resistance. Compared to phenotypic methods, sensitivity (95% CI) for rifampicin was 90.3% (89.6–90.9%), while for isoniazid it was 78.2% (77.4–79.0%) and their specificities were 96.3% (95.7–96.8%) and 94.4% (93.1–95.5%), respectively. For second-line drugs, sensitivity varied from 67.4% (64.1–70.6%) for capreomycin to 88.2% (85.1–90.9%) for moxifloxacin, with specificity ranging from 90.0% (87.1–92.5%) for moxifloxacin to 99.5% (99.0–99.8%) for amikacin. This study provides a standardised and comprehensive approach for the interpretation of mutations as predictors of M. tuberculosis drug-resistant phenotypes. These data have implications for the clinical interpretation of molecular diagnostics and next-generation sequencing as well as efficient individualised therapy for patients with drug-resistant tuberculosis

Characterization of Ceftazidime Resistance Mechanisms in Clinical Isolates of Burkholderia pseudomallei from Australia

Burkholderia pseudomallei is a Gram-negative bacterium that causes the serious human disease, melioidosis. There is no vaccine against melioidosis and it can be fatal if not treated with a specific antibiotic regimen, which typically includes the third-generation cephalosporin, ceftazidime (CAZ). There have been several resistance mechanisms described for B. pseudomallei, of which the best described are amino acid changes that alter substrate specificity in the highly conserved class A β-lactamase, PenA. In the current study, we sequenced penA from isolates sequentially derived from two melioidosis patients with wild-type (1.5 µg/mL) and, subsequently, resistant (16 or ≥256 µg/mL) CAZ phenotypes. We identified two single-nucleotide polymorphisms (SNPs) that directly increased CAZ hydrolysis. One SNP caused an amino acid substitution (C69Y) near the active site of PenA, whereas a second novel SNP was found within the penA promoter region. In both instances, the CAZ resistance phenotype corresponded directly with the SNP genotype. Interestingly, these SNPs appeared after infection and under selection from CAZ chemotherapy. Through heterologous cloning and expression, and subsequent allelic exchange in the native bacterium, we confirmed the role of penA in generating both low-level and high-level CAZ resistance in these clinical isolates. Similar to previous studies, the amino acid substitution altered substrate specificity to other β-lactams, suggesting a potential fitness cost associated with this mutation, a finding that could be exploited to improve therapeutic outcomes in patients harboring CAZ resistant B. pseudomallei. Our study is the first to functionally characterize CAZ resistance in clinical isolates of B. pseudomallei and to provide proven and clinically relevant signatures for monitoring the development of antibiotic resistance in this important pathogen

Transmission Shifts Underlie Variability in Population Responses to Yersinia pestis Infection

Host populations for the plague bacterium, Yersinia pestis, are highly variable in their response to plague ranging from near deterministic extinction (i.e., epizootic dynamics) to a low probability of extinction despite persistent infection (i.e., enzootic dynamics). Much of the work to understand this variability has focused on specific host characteristics, such as population size and resistance, and their role in determining plague dynamics. Here, however, we advance the idea that the relative importance of alternative transmission routes may vary causing shifts from epizootic to enzootic dynamics. We present a model that incorporates host and flea ecology with multiple transmission hypotheses to study how transmission shifts determine population responses to plague. Our results suggest enzootic persistence relies on infection of an off-host flea reservoir and epizootics rely on transiently maintained flea infection loads through repeated infectious feeds by fleas. In either case, early-phase transmission by fleas (i.e., transmission immediately following an infected blood meal) has been observed in laboratory studies, and we show that it is capable of driving plague dynamics at the population level. Sensitivity analysis of model parameters revealed that host characteristics (e.g., population size and resistance) vary in importance depending on transmission dynamics, suggesting that host ecology may scale differently through different transmission routes enabling prediction of population responses in a more robust way than using either host characteristics or transmission shifts alone