Whole-genome Sequencing Single Nucleotide Polymorphism Research Articles

Whole genome sequencing analysis demonstrates therapy-induced echinocandin resistance in Candida auris isolates.

Candida auris is an emerging, multidrug-resistant yeast, causing outbreaks in healthcare facilities. Echinocandins are the antifungal drugs of choice to treat candidiasis, as they cause few side effects and resistance is rarely found. Previously, immunocompromised patients from Kuwait with C. auris colonisation or infection were treated with echinocandins, and within days to months, resistance was reported in urine isolates. To determine whether the development of echinocandin resistance was due to independent introductions of resistant strains or resulted from intra-patient resistance development, whole genome sequencing (WGS) single-nucleotide polymorphism (SNP) analysis was performed on susceptible (n = 26) and echinocandin-resistant (n = 6) isolates from seven patients. WGS SNP analysis identified three distinct clusters differing 17-127 SNPs from two patients, and the remaining isolates from five patients, respectively. Sequential isolates within patients had a maximum of 11 SNP differences over a time period of 1-10 months. The majority of isolates with reduced susceptibility displayed unique FKS1 substitutions including a novel FKS1M690V substitution, and nearly all were genetically related, ranging from only three to six SNP differences compared to susceptible isolates from the same patient. Resistant isolates from three patients shared the common FKS1S639F substitution; however, WGS analysis did not suggest a common source. These findings strongly indicate that echinocandin resistance is induced during antifungal treatment. Future studies should determine whether such echinocandin-resistant strains are capable of long-term colonisation, cause subsequent breakthrough candidiasis, have a propensity to cross-infect other patients, or remain viable for longer time periods in the hospital environment.

Genotyping and susceptibility testing uncovers large azole-resistant Candida tropicalis clade in Alexandria, Egypt

Candida tropicalis is an emerging medically relevant Candida species. The yeast primarily causes opportunistic infections in intensive care units and is highly prevalent in tropical countries. The genetic diversity within this species is high, and nosocomial transmission has been reported. C. tropicalis genotyping of isolates from low- and middle-income countries is underrepresented when compared with that from high-income countries. Also, in Egypt, only limited genotyping has been conducted for C. tropicalis isolates, while antifungal resistance seems to increase, especially against azoles. Antifungal susceptibility testing was performed on 64 C. tropicalis isolates from ICU patients collected from multiple hospitals in Alexandria, Egypt. Genotyping by means of short tandem repeat (STR) and whole genome sequencing (WGS) single nucleotide polymorphism (SNP) analysis was performed. Using antifungal susceptibility testing, fluconazole resistance was observed in 24 isolates (38%), of which 23 harboured an ERG11 G464S substitution, previously shown to cause resistance in Candida albicans. STR genotyping showed that these 23 isolates were related, forming a distinct resistant clade. WGS SNP analysis subsequently confirmed this genetic relationship, although isolates within this clade differed in at least 429 SNPs, suggesting that these were independently introduced. Overall, STR and WGS SNP analysis of this collection indicates limited C. tropicalis nosocomial transmission in Alexandria, while the presence of this large azole-resistant C. tropicalis clade within this city hampers the treatment of intensive care unit patients.

Original and introduced lineages co-driving the persistence of Brucella abortus circulating in West Africa.

Brucellosis, a serious public health issue affecting animals and humans, is neglected in West Africa (WA). In the present study, bio-typing, multi-locus sequence typing (MLST), multiple-locus variable-number tandem repeat analysis (MLVA), and whole genome sequencing single-nucleotide polymorphism (WGS-SNP) analysis were used to characterize the Brucella abortus (B. abortus) strains from WA. All of the 309 strains analyzed in this study were extracted and downloaded from the international MLVA bank and were from 10 hosts (cattle, humans, ovine, buffalo, dromedaries, horse, sheep, zebu, dog, and cat) distributed in 17 countries in WA. Based on the bio-typing, three biovars, dominated by B. abortus bv.3, were observed and reported across seven decades (1958-2019). With MLST, 129 B. abortus strains from the present study were sorted into 14 STs, with ST34 as the predicted founder. These 14 STs clustered into the global MLST data into three clone complexes (C I-C III) with the majority of strains clustering in C I, while C II forms an independent branch, and C III harbors three STs shared by different continents. These data revealed that most cases were caused by strains from native lineages. According to the MLVA-11 comparison, 309 strains were divided into 22 MLVA-11 genotypes, 15 of which were unique to WA and the remaining seven had a global distribution. MLVA-16 analysis showed that there were no epidemiological links among these strains. Based on the MLVA data, B. abortus strains from WA have high genetic diversity, and predominated genotypes were descended from a native lineage. While the MLVA-16 globally highlights that the dominant native and few introduced lineages (from Brazil, the USA, South Korea, Argentina, India, Italy, Portugal, the UK, Costa Rica, and China) co-driving the B. abortus ongoing prevalence in WA. The high-resolution SNP analysis implied the existence of introduced B. abortus lineages, which may be reasonably explained by the movement and trade of dominant hosts (cattle) and/or their products. Our results indicated that B. abortus strains in WA consist of native and introduced strains that necessitate control such as vaccination, testing, slaughtering, and movement control by the relevant country authorities to reduce brucellosis in livestock.

Development and Application of a Short Tandem Repeat Multiplex Typing Assay for Candida tropicalis.

Candida tropicalis is a clinically important yeast that causes candidemia in humans with a high mortality rate. The yeast primarily infects immunocompromised patients, and causes outbreaks in health care facilities. Antifungal resistant isolates have been reported. We developed a short tandem repeat (STR) typing scheme for C. tropicalis to enable fast, cost-effective, and high-resolution genotyping. For the development of the typing scheme, 6 novel STR markers were selected, combined into 2 multiplex PCRs. In total, 117 C. tropicalis isolates were typed, resulting in the identification of 104 different genotypes. Subsequently, the outcome of STR typing of 10 isolates was compared to single nucleotide polymorphism (SNP) calling from whole-genome sequencing (WGS). Isolates with more than 111 SNPs were differentiated by the typing assay. Two isolates, which were identical according to SNP analysis, were separated by STR typing in 1 marker. To test specificity, the STR typing was applied to 15 related yeast species, and we found no amplification of these targets. For reproducibility testing, 2 isolates were independently typed five times, which showed identical results in each experiment. In summary, we developed a reliable and multiplex STR genotyping for C. tropicalis, which was found to correlate well to SNP calling by WGS. WGS analysis from and extensive collection of isolates is required to establish the precise resolution of this STR assay. IMPORTANCE Candida tropicalis frequently causes candidemia in immunocompromised patients. C. tropicalis infections have a high mortality rate, and the yeast is able to cause outbreaks in health care facilities. Further, antifungal resistant isolates are on the rise. Genotyping is necessary to investigate potential outbreaks. Here, we developed and applied a STR genotyping scheme in order to rapidly genotype isolates with a high-resolution. WGS SNP outcomes were highly comparable with STR typing results. Altogether, we developed a rapid, high-resolution, and specific STR genotyping scheme for C. tropicalis.

The design and application of a 50 K SNP chip for a threatened Aotearoa New Zealand passerine, the hihi.

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or "SNP chips", enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species.

Investigation of Genetic Relatedness of Brucella Strains in Countries Along the Silk Road.

In this study, MLVA (multiple-locus variable-number tandem repeat analysis) genotype data of Brucella strains from 11 countries along the Silk Road were downloaded from the MLVAbank. MLVA data of strains were applied to the constructed Minimum Spanning Tree to explore the species/biovars distribution, geographic origins, and genetic relationships of the strains analyzed. Moreover, whole-genome sequencing–single-nucleotide polymorphism (WGS-SNP) phylogenetic analysis of the genome of Brucella melitensis strains from GenBank was performed to discriminate the relatedness of strains further and investigate the transmission pattern of B. melitensis brucellosis. A total of 1,503 Brucella strains were analyzed in this study: 431 Brucella abortus strains (29.8%), 1,009 B. melitensis strains (65.7%), and 63 Brucella suis strains (4.5%). B. melitensis biovar 3 was the dominant species and was shown to be widespread in all of the examined regions, suggesting that the prevention and surveillance of the B. melitensis population are a main challenge in these countries. A wide host spectrum was observed for this Brucella population; many animal reservoirs are a potential reason for the continuous brucellosis circulation in these countries. Although the B. abortus strains from the examined regions had common geographic origins, only a few shared genotypes were observed in different countries. These data revealed that the majority of B. abortus strains were spreading within the national borders. However, the B. melitensis strains from Italy originated from a Western Mediterranean lineage; strains from the other 10 countries originated from Eastern Mediterranean lineage, and this lineage was shared by strains from three to nine different countries, suggesting that the introduction and reintroduction of the disease in the 10 countries might have occurred in the past. Furthermore, the most shared MLVA-16 genotypes were formed in the B. melitensis strains from China, Kazakhstan, and Turkey, suggesting that the introduction and trade in sheep and goats have occurred frequently in these countries. WGS-SNP analysis showed that the B. melitensis in this study originated from the Malta (Italy) region. According to their territorial affiliation between four clade strains from these countries in genotype B, the absence of a clear differentiation suggests that strains continuously expand and spread in countries along with Silk Road. Active exchange and trade of animals (sheep and goats) among these countries are reasonable explanations. B. suis strains from different nations showed unique geographic origins and epidemiological characteristics. Therefore, there is an urgent need for the control of transfer and trade of infected sheep (goats) in countries along the Silk Road, namely, the strengthening of the entry–exit quarantine of sheep and goats and improvements in the diagnosis of animal brucellosis.

Evaluation of a Combined Multilocus Sequence Typing and Whole-Genome Sequencing Two-Step Algorithm for Routine Typing of Clostridioides difficile.

Multilocus sequence typing (MLST) is a low-resolution but rapid genotyping method for Clostridioides difficile Whole-genome sequencing (WGS) has emerged as the new gold standard for C. difficile typing, but cost and lack of standardization still limit broad utilization. In this study, we evaluated the potential to combine the portability of MLST with the increased resolution of WGS for a cost-saving approach to routine C. difficile typing. C. difficile strains from two New York City hospitals (hospital A and hospital B) were selected. WGS single-nucleotide polymorphism (wgSNP) was performed using established methods. Sequence types (ST) were determined using PubMLST, while wgSNP analysis was performed using the Bionumerics software. An additional analysis of a subset of data (hospital A) was made comparing the Bionumerics software to the CosmosID pipeline. Cost and turnaround time to results were compared for the algorithmic approach of MLST followed by wgSNP versus direct wgSNP. Among the 202 C. difficile isolates typed, 91% (n = 185/203) clustered within the representative ST, showing a high agreement between MLST and wgSNP. While clustering was similar between the Bionumerics and CosmosID pipelines, large differences in the overall number of SNPs were noted. A two-step algorithm for routine typing results in significantly lower cost than routine use of WGS. Our results suggest that using MLST as a first step in routine typing of C. difficile followed by WGS for MLST concordant strains is a less technically demanding, cost-saving approach for performing C. difficile typing than WGS alone without loss of discriminatory power.

Extensive Homoplasy but No Evidence of Convergent Evolution of Repeat Numbers at MIRU Loci in Modern Mycobacterium tuberculosis Lineages.

More human deaths have been attributable to Mycobacterium tuberculosis than any other pathogen, and the epidemic is sustained by ongoing transmission. Various typing schemes have been developed to identify strain-specific differences and track transmission dynamics in affected communities, with recent introduction of whole genome sequencing providing the most accurate assessment. Mycobacterial interspersed repetitive unit (MIRU) typing is a family of variable number tandem repeat schemes that have been widely used to study the molecular epidemiology of M. tuberculosis. MIRU typing was used in most well-resourced settings to perform routine molecular epidemiology. Instances of MIRU homoplasy have been observed in comparison with sequence-based phylogenies, limiting its discriminatory value. A fundamental question is whether the observed homoplasy arises purely through stochastic processes, or whether there is evidence of natural selection. We compared repeat numbers at 24 MIRU loci with a whole genome sequence-based phylogeny of 245 isolates representing three modern M. tuberculosis lineages. This analysis demonstrated extensive homoplasy of repeat numbers, but did not detect any evidence of natural selection of repeat numbers, at least since the ancestral branching of the three modern lineages of M. tuberculosis. In addition, we observed good sensitivity but poor specificity and positive predictive values of MIRU-24 to detect clusters of recent transmission, as defined by whole-genome single nucleotide polymorphism analysis. These findings provide mechanistic insight, and support a transition away from VNTR-based typing toward sequence-based typing schemes for both research and public health purposes.

Insights into the long-term persistence of Legionella in facilities from whole-genome sequencing

We investigated the value of whole-genome sequencing (WGS) and single nucleotide polymorphism (SNP) analyses in determining the relationships among and evolutionary rates of Legionella species with long-term persistence in three healthcare facilities. We examined retrospective clinical and environmental isolates of Legionella micdadei and Legionella pneumophila serogroup 1 isolates with identical PFGE DNA fingerprints sampled over the course of up to 18 years. WGS analyses demonstrated that heterogeneous populations of Legionella were present within each facility despite displaying the same PFGE profiles. Additionally, clustering of some clinical isolates with those from a separate but related institution exposed a source of infection not previously detected, underscoring the importance of considering phylogenetic relationships when assessing epidemiological links. The data supported an average substitution rate of 0.80 SNPs per genome per year for L. micdadei but a reliable estimate for L. pneumophila serogroup 1 could not be obtained due to complicating factors such as non-chronological links among isolates and inadequate sampling depths. While the substitution rate for L. micdadei is consistent with previous estimates for L. pneumophila, the lack of a temporal signal in our sequence data for L. pneuomphila serogroup 1 isolates suggests either insufficient change to provide an estimate or variable evolutionary rates, which could reflect the presence of both actively dividing and viable but non-culturable Legionella spp. in the built environment. This study highlights the increased discriminatory power of WGS SNP analysis as compared to PFGE, emphasizes the need for extended sampling, and provides insight into the evolution of Legionella from longitudinal investigations.

Natural Genetic Variation in the Caenorhabditis elegans Response to Pseudomonas aeruginosa.

Caenorhabditis elegans responds to pathogenic microorganisms by activating its innate immune system, which consists of physical barriers, behavioral responses, and microbial killing mechanisms. We examined whether natural variation plays a role in the response of C. elegans to Pseudomonas aeruginosa using two C. elegans strains that carry the same allele of npr-1, a gene that encodes a G-protein-coupled receptor related to mammalian neuropeptide Y receptors, but that differ in their genetic backgrounds. Strains carrying an allele for the NPR-1 215F isoform have been shown to exhibit lack of pathogen avoidance behavior and deficient immune response toward P. aeruginosa relative to the wild-type (N2) strain. We found that the wild isolate from Germany RC301, which carries the allele for NPR-1 215F, shows an enhanced resistance to P. aeruginosa infection when compared with strain DA650, which also carries NPR-1 215F but in an N2 background. Using a whole-genome sequencing single-nucleotide polymorphism (WGS-SNP) mapping strategy, we determined that the resistance to P. aeruginosa infection maps to a region on chromosome V. Furthermore, we demonstrated that the mechanism for the enhanced resistance to P. aeruginosa infection relies exclusively on strong P. aeruginosa avoidance behavior, and does not involve the main immune, stress, and lifespan extension pathways in C. elegans. Our findings underscore the importance of pathogen-specific behavioral immune defense in the wild, which seems to be favored over the more energy-costly mechanism of activation of physiological cellular defenses.