Whole-genome Single-nucleotide Polymorphism Analysis Research Articles

An Evaluation of the Lineage of Brucella Isolates in Turkey by a Whole-Genome Single-Nucleotide Polymorphism Analysis.

Brucellosis is a disease caused by the Brucella (B.) species. It is a zoonotic disease that affects farm animals and causes economic losses in many countries worldwide. Brucella has the ability to persist in the environment and infect the host at low doses. Thus, it is more important to trace brucellosis outbreaks, identify their sources of infection, and interrupt their transmission. Some countries already have initial data, but most of these data are based on a Multiple-Locus Variable-Number Tandem-Repeat Analysis (MLVA), which is completely unsuitable for studying the Brucella genome. Since brucellosis is an endemic disease in Turkey, this study aimed to examine the genome of Turkish Brucella isolates collected between 2018 and 2020, except for one isolate, which was from 2012. A total of 28 strains of B. melitensis (n = 15) and B. abortus (n = 13) were analyzed using a core-genome single-nucleotide polymorphism (cgSNP) analysis. A potential connection between the Turkish isolates and entries from Sweden, Israel, Syria, Austria, and India for B. melitensis was detected. For B. abortus, there may be potential associations with entries from China. This explains the tight ties found between Brucella strains from neighboring countries and isolates from Turkey. Therefore, it is recommended that strict measures be taken and the possible effects of uncontrolled animal introduction are emphasized.

Whole-genome sequencing reveals genomic characterization of Listeria monocytogenes from food in China.

Listeria monocytogenes is a foodborne bacterium that could persist in food and food processing environments for a long time. Understanding the population structure and genomic characterization of foodborne L. monocytogenes is essential for the prevention and control of listeriosis. A total of 322 foodborne L. monocytogenes isolates from 13 geographical locations and four food sources in China between 2000 and 2018 were selected for whole-genome sequencing. In silico subtyping divided the 322 isolates into five serogroups, 35 sequence types (STs), 26 clonal complexes (CCs) and four lineages. Serogroup IIa was the most prevalent serogroup and ST9 was the most prevalent ST of foodborne L. monocytogenes strains isolated in China. The in-depth phylogenetic analysis on CC9 revealed that ST122 clone might be original from ST9 clone. Furthermore, 23 potentially relevant clusters were identified by pair-wised whole-genome single nucleotide polymorphism analysis, indicating that persistent- and/or cross-contamination had occurred in markets in China. ST8 and ST121 were the second and third top STs of L. monocytogenes in China, which had heterogeneity with that of L. monocytogenes isolates from other countries. The antibiotic resistance genes aacA4, tetM, tetS, dfrG carried by different mobile elements were found in L. monocytogenes strains. One lineage II strain carrying Listeria Pathogenicity Island 3 was first reported. In addition, a novel type of premature stop codon in inlA gene was identified in this study. These findings revealed the genomic characteristics and evolutionary relationship of foodborne L. monocytogenes in China on a scale larger than previous studies, which further confirmed that whole-genome sequencing analysis would be a helpful tool for routine surveillance and source-tracing investigation.

Whole-genome single nucleotide polymorphism analysis for typing the pandemic pathogen Fusarium graminearum sensu stricto.

Recent improvements in microbiology and molecular epidemiology were largely stimulated by whole- genome sequencing (WGS), which provides an unprecedented resolution in discriminating highly related genetic backgrounds. WGS is becoming the method of choice in epidemiology of fungal diseases, but its application is still in a pioneer stage, mainly due to the limited number of available genomes. Fungal pathogens often belong to complexes composed of numerous cryptic species. Detecting cryptic diversity is fundamental to understand the dynamics and the evolutionary relationships underlying disease outbreaks. In this study, we explore the value of whole-genome SNP analyses in identification of the pandemic pathogen Fusarium graminearum sensu stricto (F.g.). This species is responsible for cereal diseases and negatively impacts grain production worldwide. The fungus belongs to the monophyletic fungal complex referred to as F. graminearum species complex including at least sixteen cryptic species, a few among them may be involved in cereal diseases in certain agricultural areas. We analyzed WGS data from a collection of 99 F.g. strains and 33 strains representing all known cryptic species belonging to the FGSC complex. As a first step, we performed a phylogenomic analysis to reveal species-specific clustering. A RAxML maximum likelihood tree grouped all analyzed strains of F.g. into a single clade, supporting the clustering-based identification approach. Although, phylogenetic reconstructions are essential in detecting cryptic species, a phylogenomic tree does not fulfill the criteria for rapid and cost-effective approach for identification of fungi, due to the time-consuming nature of the analysis. As an alternative, analysis of WGS information by mapping sequence data from individual strains against reference genomes may provide useful markers for the rapid identification of fungi. We provide a robust framework for typing F.g. through the web-based PhaME workflow available at EDGE bioinformatics. The method was validated through multiple comparisons of assembly genomes to F.g. reference strain PH-1. We showed that the difference between intra- and interspecies variability was at least two times higher than intraspecific variation facilitating successful typing of F.g. This is the first study which employs WGS data for typing plant pathogenic fusaria.

Genomic analysis of Shiga toxin-producing Escherichia coli O157:H7 from cattle and pork-production related environments

Three E. coli O157:H7 outbreaks have been attributed to contaminated pork in Alberta, Canada, recently. This study investigates the phylogenetic relatedness of E. coli O157:H7 from pigs, cattle, and pork-production environments for source attribution. Limited strain diversity was observed using five conventional subtyping methods, with most or all strains being in one subgroup. Whole-genome single nucleotide polymorphism analysis confirmed the recent ancestry of the isolates from all three sources. Most environmental isolates clustered closer with pig isolates than cattle isolates. Also, a direct link was observed between 2018-outbreak environmental isolates and isolates collected from a pig farm in 2018. The majority of pig isolates harbor only one Shiga toxin gene, stx2a, while 70% (35/50) of the cattle isolates have both stx1a and stx2a. The results show some E. coli O157:H7 strains could establish persistence on pig farms and as such, pigs can be a significant source of the organism.

Whole-Genome Sequencing for Investigating a Health Care-Associated Outbreak of Carbapenem-Resistant Acinetobacter baumannii.

Carbapenem-resistant Acinetobacter baumannii (CRAB) outbreaks in hospital settings challenge the treatment of patients and infection control. Understanding the relatedness of clinical isolates is important in distinguishing outbreak isolates from sporadic cases. This study investigated 11 CRAB isolates from a hospital outbreak by whole-genome sequencing (WGS), utilizing various bioinformatics tools for outbreak analysis. The results of multilocus sequence typing (MLST), single nucleotide polymorphism (SNP) analysis, and phylogenetic tree analysis by WGS through web-based tools were compared, and repetitive element polymerase chain reaction (rep-PCR) typing was performed. Through the WGS of 11 A. baumannii isolates, three clonal lineages were identified from the outbreak. The coexistence of blaOXA-23, blaOXA-66, blaADC-25, and armA with additional aminoglycoside-inactivating enzymes, predicted to confer multidrug resistance, was identified in all isolates. The MLST Oxford scheme identified three types (ST191, ST369, and ST451), and, through whole-genome MLST and whole-genome SNP analyses, different clones were found to exist within the MLST types. wgSNP showed the highest discriminatory power with the lowest similarities among the isolates. Using the various bioinformatics tools for WGS, CRAB outbreak analysis was applicable and identified three discrete clusters differentiating the separate epidemiologic relationships among the isolates.

Genomic analysis of Clostridioides difficile in two regions reveals a diversity of strains and limited transmission

Background: The epidemic NAP1/027 Clostridioides difficile strain (MLST1, ST1) that emerged in the mid-2000 is on the decline. The current distribution of C. difficile strain types and their transmission dynamics are poorly defined. We performed whole-genome sequencing (WGS) of C. difficile isolates in 2 regions to identify the predominant multilocus sequence types (MLSTs) in community- and healthcare-associated cases and potential transmission between cases using whole-genome single-nucleotide polymorphism (SNP) analysis. Methods: Isolates were collected through the CDC Emerging Infections Program population-based surveillance for C. difficile infections (CDI) for 3 months between 2016 and 2017 in 5 Minnesota counties and 1 New York county. Isolates were limited to incident cases (CDI in a county resident with no positive C. difficile test in the preceding 8 weeks). Cases were classified as healthcare associated (HA-CDI) or community associated (CA-CDI) based on healthcare exposures as previously described. WGS was performed on an Illumina Miseq. The CFSAN (FDA) pipeline was used to compute whole-genome SNPs, SPAdes was used for assembly, and MLST was assigned according to www.pubmlst.org. Results: Of 431 isolates, 269 originated from New York and 162 from Minnesota; 203 cases were classified as CA-CDI and 221 as HA-CDI. The proportion of CA-CDI cases was higher in Minnesota than in New York: 62% vs 38%. The predominant MLSTs across both sites were ST42 (9%), ST8 (8%), and ST2 (8%). MLSTs more frequently encountered in HA-CDI than CA-CDI included ST1 (note that this ST includes PCR Ribotype 027; 76% HA-CDI), ST53 (84% HA-CDI), and ST43 (80% HA-CDI). In contrast, ST110 (63% CA-CDI) and ST3 (67% CA-CDI) were more commonly isolated from CA-CDI cases. ST1 accounted for 7.6% of circulating strains and was more common in New York than Minnesota (10% vs 3%) and was concentrated among New York HA-CDI cases. Also, 412 isolates (1 per patient) were included in the final whole-genome SNP analysis. Of these, only 12 pairs were separated by 0–3 SNPs, indicating potential transmission, and most involved HA-CDI cases. ST1, ST17, and ST46 accounted for 8 of 12 pairs, with ST17 and ST46 potentially forming small clusters. Conclusions: This analysis provides a snapshot of the current genomic epidemiology of C. difficile across 2 geographically and epidemiologically distinct regions of the United States and supports other studies suggesting that the role of direct transmission in the spread of CDI may be limited.Funding: NoneDisclosures: None

Identification of Rickettsial Infections (Rickettsia sp. TH2014) in Ctenocephalides orientis Fleas (Siphonaptera: Pulicidae)

Rickettsia felis (Rickettsiales: Rickettsiaceae) is an emergent human pathogen that causes febrile illnesses in various parts of the world. This study describes the identification and growth characteristics of a R. felis-like organism (designated as Rickettsia sp. TH2014) cultured from Ctenocephalides orientis fleas in rural Malaysia. In this study, culturing of rickettsiae from filtered triturated flea lysates was performed in Aedes albopictus C6/36 cells. Cytopathic effects were observed from one of the samples 4 d post-inoculation. Electron microscopy revealed actively replicating intracytosolic coccobacillary organisms in the rickettsia-infected cells. Sequence analysis of amplified citrate synthase (gltA) gene fragment shows complete match of the rickettsia with Rickettsia sp. Rf31 in Southeast Asia, and 'Candidatus Rickettsia senegalensis' strain PU01-02 in Africa. The whole-genome sequence of Rickettsia sp. TH2014 was determined and assembled. The estimated genome size and guanine + cytosine content of the rickettsia are 1.37 Mb and 32.9%, respectively. The high values of average nucleotide identity and tetra-nucleotide signature correlation index obtained from pairwise genome comparison study suggest the identification of the rickettsia as R. felis. The whole-genome single-nucleotide polymorphism analysis demonstrates close genetic relatedness of the rickettsia with R. felis and Rickettsia asemboensis. However, based on sequence analyses of rickettsial genes (16S rDNA, gltA, ompB, and sca4), Rickettsia sp. TH2014 is found to be distinct from R. felis and R. asemboensis. The sequence analyses reveal that Rickettsia sp. TH2014 is highly similar to 'Ca. Rickettsia senegalensis' detected in fleas from Africa, Asia, and North America. Further investigation to provide insights on pathogenic potential and transmission dynamics of the rickettsia is warranted.

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis

BackgroundBrucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis.ResultsSingle sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage.ConclusionWe showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1294-x) contains supplementary material, which is available to authorized users.

Draft Genome Sequence of Bivalent Clostridium botulinum Strain IBCA10-7060, Encoding Botulinum Neurotoxin B and a New FA Mosaic Type.

Here we report the genome sequence of a Clostridium botulinum strain IBCA10-7060 producing botulinum neurotoxin serotype B and a new toxin serotype. Multilocus sequence typing analysis revealed that this strain belongs to a new sequence type, and whole-genome single nucleotide polymorphism analysis showed that this strain clustered with strains in lineage 2 from group I.

Genome-wide association study of resistance to rough dwarf disease in maize

Maize rough dwarf disease (MRDD), caused by the Maize rough dwarf virus (MRDV) is highly prevalent across summer at maize-growing areas in China. To investigate the genetic architecture underlying this viral disease, a set of 236 Chinese maize inbred lines was evaluated for resistance to MRDD in 2010 and 2011. Based on 41,101 single-nucleotide polymorphism (SNP) markers with minor allele frequencies (MAF) greater than 5 %, a genome-wide association study (GWAS) was employed to identify genomic loci for resistance to MRDD. A total of 73 SNPs were found to be associated with resistance to MRDD at a significance threshold of -log10 (P) > 4 controlling the false discovery rate (FDR) at α = 0.1. Fourteen of these SNPs were detected in both of the two environments tested. A total of 48 SNPs were identified in linkage disequilibrium (LD) blocks containing candidate resistance genes, including protein kinase genes. Using the pedigree information and whole-genome SNP analysis of five highly resistant inbred lines derived from the hybrid ‘P78599’, nine derivative fragments harbouring SNPs associated with MRDD resistance were detected. One 81.57 Mb fragment in particular located in bin 8.03, which contained six SNPs associated with MRDD resistance, and included the major quantitative trait loci (QTL) that had been identified in the previous study. These results suggest that the SNPs and fragments associated with MRDD resistance, especially those in bin 8.03, could be used for fine mapping of resistance genes and developing resistant varieties in maize.

Clinically failed eggs as a source of normal human embryo stem cells

The promise of human embryo stem cells (hESCs) for regenerative medicine is offset by the ethical and practical challenges involved in sourcing eggs and embryos for this objective. In this study we sought to isolate an hESC line from clinically failed eggs, the usage of which would not conflict with donor interests to conceive. A total of 8 blastocysts were allocated for hESC derivation from a pool of 579 eggs whose fertilization had been clinically assessed to have occurred abnormally (i.e., three pronuclei) or failed (i.e., no pronuclei) following in vitro insemination or intracytoplasmic sperm injection (ICSI). The latter were subjected to a recovery intervention consisting of either reinsemination by ICSI or parthenogenetic stimulation. One hESC line (RCM1) was obtained from a failed-to-fertilize inseminated egg recovered by parthenogenetic activation. Standard in vitro and in vivo characterization revealed this line to possess all of the properties attributed to a normal euploid hESC line. Whole-genome single-nucleotide polymorphism analysis further revealed that the line was biparental, indicating that sperm penetration had occurred, although parthenogenetic stimulation was required for activation. Our results demonstrate the viability of an alternative strategy to generate normal hESC lines from clinically failed eggs, thereby further minimizing the potential to conflict with donor reproductive interest to conceive.

Novel Mutations in NPHP4 in a Consanguineous Family With Histological Findings of Focal Segmental Glomerulosclerosis

Nephronophthisis is a form of autosomal recessive hereditary cystic kidney disease that typically progresses to end-stage renal disease by early adulthood. Conversely, focal segmental glomerulosclerosis is a histological glomerular phenotype that can be familial, primary (idiopathic), or secondary to a multitude of pathological processes affecting the kidney, including such tubulointerstitial diseases as nephronophthisis. Mutations in 6 distinct nephronophthisis genes have been described to date. We describe a consanguineous Filipino family with 2 novel sequence variants in the NPHP4 gene. Affected individuals presented with end-stage renal disease and histological features of focal segmental glomerulosclerosis on biopsy. They also had atypical radiological findings, making the clinical diagnosis of the genetic syndrome difficult. Furthermore, although ocular abnormalities and hearing loss were described previously, this is the first report of hepatic disease in patients with mutations in NPHP4. The diagnosis of nephronophthisis was made by means of mutational analysis of the NPHP4 gene after isolation of a region of homozygosity in affected individuals by using whole-genome single-nucleotide polymorphism analysis. Because establishment of the correct diagnosis has implications for therapeutic interventions, prognosis, and, in the case of heritable diseases, appropriate genetic counseling for affected individuals and their families, this report emphasizes the importance of obtaining meticulous clinical information, considering alternative diagnoses, and, when possible, performing genetic evaluation to confirm the diagnosis. We outline an approach to patients with hereditary kidney disease, focusing specifically on the molecular genetic techniques available to evaluate such families and determine a chromosomal region of interest and, subsequently, the diagnosis.