Development and evaluation of a core genome multilocus sequence typing (cgMLST) scheme for Brucella spp

Opposite effects of soil pH on bacteria and fungi β diversity in forests at a continental scale

Soil microbial diversity is crucial for regulating biogeochemical cycles, including soil carbon (C) dynamics and nutrient cycling. However, how climate, plants, and soil properties influence the microbiome in forests remains unclear, especially at the continental scale, hindering us to better understand forest C-climate change feedback. Here, we investigated the spatial patterns of microbial diversity across China's forests and explored the controlling factors of microbial β diversity and network complexity. Our results showed that soil pH strongly influenced bacterial and fungal β diversity compared to climate, soil nutrient and plant properties. To further investigate the environmental preference of the microbial networks, we classified the amplicon sequence variants (ASVs) into five groups ranging from acidic to alkaline soils. Co-occurrence network analysis revealed that the topological structure of the bacterial network (e.g., edge and degree) increased with pH and was negatively correlated with β diversity but not for fungal diversity. Soil fungi exhibited higher β diversity and network complexity (i.e., degree and betweenness) than bacteria in acidic soils (pH < 5.1), and vice versa in neutral and alkaline soils (pH > 5.5). Within the pH range of 5.1-5.5, the bacterial-fungal network displayed the highest network complexity with the lowest fungal β diversity, and significant positive correlations were found between fungal β diversity and soil properties. In addition, bacterial growth in acidic soil (pH < 5.5) showed positive correlations with acid phosphatase (AP), but negative ones with β-1,4-glucosidase (BG), and vice versa in neutral and alkaline soils (pH > 5.5). Furthermore, 46 bacterial core species were identified, and their abundance had significant correlation with soil pH. These findings highlight the critical role of soil pH in driving soil microbial β diversity across China's forests and reveal the effects of pH thresholds on changes in the soil microbial network and core species.

OrthoPhyl-streamlining large-scale, orthology-based phylogenomic studies of bacteria at broad evolutionary scales

There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis.

Molecular typing methods to characterize Brucella spp. from animals: A review

Brucellosis is an infectious disease of animals that can infect humans. The disease causes significant economic losses and threatens human health. A timely and accurate disease diagnosis plays a vital role in the identification of brucellosis. In addition to traditional diagnostic methods, molecular methods allow diagnosis and typing of the causative agent of brucellosis. This review will discuss various methods, such as Bruce-ladder, Suiladder, high-resolution melt analysis, restriction fragment length polymorphism, multilocus sequence typing, multilocus variable-number tandem repeat analysis, and whole-genome sequencing single-nucleotide polymorphism, for the molecular typing of Brucella and discuss their advantages and disadvantages.

Comparative genomics of Brucella abortus and Brucella melitensis unravels the gene sharing, virulence factors and SNP diversity among the standard, vaccine and field strains

Brucella abortus and Brucella melitensis are the primary etiological agents of brucellosis in large and small ruminants, respectively. There are limited comparative genomic studies involving Brucella strains that explore the relatedness among both species. In this study, we involved strains (n=44) representing standard, vaccine and Indian field origin for pangenome, single nucleotide polymorphism (SNP) and phylogenetic analysis. Both species shared a common gene pool representing 2884 genes out of a total 3244 genes. SNP-based phylogenetic analysis indicated higher SNP diversity among B. melitensis (3824) strains in comparison to B. abortus (540) strains, and a clear demarcation was identified between standard/vaccine and field strains. The analysis for virulence genes revealed that virB3, virB7, ricA, virB5, ipx5, wbkC, wbkB, and acpXL genes were highly conserved in most of the Brucella strains. Interestingly, virB10 gene was found to have high variability among the B. abortus strains. The cgMLST analysis revealed distinct sequence types for the standard/vaccine and field strains. B. abortus strains from north-eastern India fall within similar sequence type differing from other strains. In conclusion, the analysis revealed a highly shared core genome among two Brucella species. SNP analysis revealed B. melitensis strains exhibit high diversity as compared to B. abortus strains. Strains with absence or high polymorphism of virulence genes can be exploited for the development of novel vaccine candidates effective against both B. abortus and B. melitensis.

Alteration of adeS Contributes to Tigecycline Resistance and Collateral Sensitivity to Sulbactam in Acinetobacter baumannii

The treatment of extensively drug-resistant (XDR) A. baumannii has emerged as a major problem. Tigecycline (TGC) and sulbactam (SUL) are both effective antibiotics against XDR A. baumannii. Here, we investigated the in-host evolution and mechanism of collateral sensitivity (CS) phenomenon in development of tigecycline resistance accompanied by a concomitant increase of sulbactam susceptibility. A total of four XDR A. baumannii strains were sequentially isolated from the same patient suffering from bacteremia. Core-genome multilocus sequence typing separated all the strains into two clusters. Comparative analysis of isolate pair 1 revealed that multiplication of blaOXA-23 within Tn2006 on the chromosome contributed to the change in the antimicrobial susceptibility phenotype of isolate pair 1. Additionally, we observed the emergence of CS to sulbactam in isolate pair 2, as demonstrated by an 8-fold increase in the TGC MIC with a simultaneous 4-fold decrease in the SUL MIC. Compared to the parental strain Ab-3557, YZM-0406 showed partial deletion in the two-component system sensor adeS. Reconstruction of the adeS mutant in Ab-3557 in situ suggested that TGC resistance and CS to SUL were mainly caused by the mutation of adeS. Overall, our study reported a novel CS combination of TGC and SUL in A. baumannii and further revealed a mechanism of CS attributed to the mutation of adeS. This study provides a valuable foundation for developing effective regimens and sequential combinations of tigecycline and sulbactam against XDR A. baumannii. IMPORTANCE Collateral sensitivity (CS) has become an increasingly common evolutionary trade-off during adaptive bacterial evolution. Here, we report a novel combination of tigecycline (TGC) resistance and CS to sulbactam (SUL) in A. baumannii. TGC and SUL are both effective antibiotics against XDR A. baumannii, and it is essential to reveal the mechanism of CS between TGC and SUL. In our study, the partial deletion of adeS, a two-component system sensor, was confirmed to be the key factor contributing to this CS phenomenon. This study provides a valuable foundation for developing effective regimens and sequential combinations of tigecycline and sulbactam against XDR A. baumannii.

Whole-genome sequencing for genetic diversity analysis of Iranian Brucella spp. isolated from humans and livestock

Brucellosis is one of the most common zoonoses in the Middle East. It is causing economic losses to the livestock industry and has a great public health concern. Little is known about the genetic diversity and distribution of brucellae in Iran. Therefore, forty Brucella spp. strains (B. abortus and B. melitensis) isolated from animals and humans were analyzed by whole genome sequencing (WGS) technology using single nucleotide polymorphism (SNP) analysis and core genome multilocus sequence typing (cgMLST). Brucella isolates were obtained from lymph nodes (cows and camels), milk (cows, camels and sheep), and aborted foetus samples (sheep and goats), as well as cerebrospinal fluid and blood of humans. The isolates were originating from thirteen provinces of Iran and isolated between 2015 and 2020. According to in-silico MLST, ST8 and ST2 were the most frequent sequence types in B. melitensis and B. abortus, respectively. Based on phylogeographic reconstruction using cgSNP analysis, the investigated Iranian B. melitensis strains belonged to the American and Mediterranean lineages of the B. melitensis phylogeny. Furthermore, cgSNP analysis revealed a similarity between Iranian B. abortus isolates and strains from Iraq and Egypt. Therefore, the origin of the Iranian strains can be suggested to be strains from neighboring and Middle East countries. Moreover, cgMLST analysis showed that the Iranian B. melitensis strains were closely relative to strains recovered from sheep and humans in Iraq, Afghanistan, Syria, Turkmenistan, and Pakistan. In the current panel of strains, cgMLST and cgSNP analysis provided an appropriate and accurate tool for effective traceback analyses for Brucella spp. from Iran. The results of cgSNP and cgMLST helped to understand the geographic distribution and interspecies transmission of Iranian strains and highlight the importance of specific brucellosis control measures in Iran with regard to the One-Health approach.

Genomes-based MLST, cgMLST, wgMLST and SNP analysis of Salmonella Typhimurium from animals and humans

Salmonella Typhimurium (S. Typhimurium) is an important food-borne and zoonotic pathogen that causes salmonellosis. With the development of whole genome sequencing (WGS), genome-based typing has been widely applied to bacteriology. In this study, we investigated genotyping and phylogenetic clusters of S. Typhimurium isolates from humans and animals in different provinces (including Beijing, Shandong, Guangxi, Shaanxi, Henan, and Shanghai) of China during 2009-2018 using multi locus sequence typing (MLST), core genome MLST (cgMLST), whole genome MLST (wgMLST) and single nucleotide polymorphism (SNP) based on WGS. 29 S. Typhimurium isolates from chicken (n = 22), sick pigeon (n = 2), patients (n = 4) and sick swine (n = 1) were tested. MLST analysis showed S. Typhimurium strains were divided into four STs, namely ST19 (n = 14), ST34 (n = 12), ST128 (n = 2) and ST1544 (n = 1). cgMLST and wgMLST divided 29 strains into 27 cgSTs and 29 wgST, respectively. Phylogenetic clustering showed that all isolates were divided into 4 clusters and 4 singletons. SNP analysis was used to examine MLST, cgMLST, wgMLST analysis. Finally, comparisons of MLST, cgMLST, wgMLST, and SNP were analyzed and the results showed their precision increased in order. In summary, genomic typing and phylogenetic relationships of 29 S. Typhimurium strains from different sources in China were analyzed. These findings were beneficial to investigate molecular pathogenesis, bacterial diversity, and traceability analysis of Salmonella.

Application of a core genome sequence typing (cgMLST) pipeline for surveillance of Clostridioides difficile in China

Introduction

Clostridioides difficile (C. difficile) is a nosocomial bacterial pathogen that causes antibiotic-associated diarrhea mediated by cellular exotoxins secreted into the intestine during bacterial growth. Multilocus sequence typing (MLST) and PCR ribotyping are the main molecular typing for C. difficile. Whole genome sequencing (WGS) core genome multilocus sequence typing (cgMLST) was developed for genetic evolution and outbreak investigation of C. difficile with higher precision and accuracy.

Methods

A total of 699 whole (complete and draft) genome sequences of distinct C. difficile strains were used in this study to identify core gene set (2469 core genes) and the cgMLST scheme for the phylogeny analysis of C. difficile. This cgMLST pipeline was then carried the Chinese Pathogen Identification Net (China PIN) for surveillance of C. difficile in China. Within the China PIN, 195 WGS of C. difficile and an outbreak of CDI with 12 WGS of C. difficile were used to evaluate the cgMLST pipeline.

Results

The result displayed that mostly tested C. difficile isolates could be successfully divided into 5 classic clades and the outbreak event was also successfully identified.

Discussion

The results are meaningful and offer a practicable pipeline for a national-wide surveillance of C. difficile in China.

First record of the human infection of Brucella melitensis in Kyrgyzstan: evidence from whole-genome sequencing-based analysis

BACKGROUND

Brucellosis, a zoonosis mainly transmitted by consumption of unpasteurized dairy products as well as direct contact with infected animals, is endemic in Kyrgyzstan. However, Brucella species in humans have not been investigated and the origin of the disease remains poorly known in wide parts of Сentral Asia. Thus, molecular characterization of the circulating strains is a critical first step in understanding Brucella diversity in the country.

METHODS

In this study, isolates were collected from patients with suspected brucellosis from different regions in Kyrgyzstan between 2019 and 2020. The detection and identification of Brucella was carried out by Bruce-ladder PCR. Next generation sequencing was used to sequence the 89 Brucella isolates, which were genotyped by cgSNP and cgMLST to identify epidemiological connection between Brucella isolates as well as placing them in the context of the global Brucella phylogeny.

RESULTS

The Brucella strains isolated from all regions of Kyrgyzstan were identified as B. melitensis. Based on cgSNP analysis, 18 sequence types were differentiated. The highest numbers of different sequence types were found in Batken (n = 8), Osh (n = 8) and Jalal-Abad (n = 6) oblasts. According to cgSNP and cgMLST analyses, different B. melitensis lineages circulate in Kyrgyzstan, all of them belonging to the Eastern Mediterranean group of the global Brucella phylogeny with the highest similarity to strains from Turkmenistan, Iran and Turkey.

CONCLUSION

In the present study, B. melitensis was identified as a causative agent of human brucellosis in Kyrgyzstan and different lineages could be identified. Since this study focused on isolates of human origin, the identity of Brucella species and lineages circulating among animal populations remains elusive. Implementing culture techniques and use of most recent molecular, bioinformatic and epidemiological tools are needed to set up a One Health approach to combat brucellosis in Kyrgyzstan. Further, other Сentral Asian countries need to take part in this effort as brucellosis is a transboundary disease in these regions.

Molecular characterization and phylogeny of Shiga toxin-producing Escherichia coli derived from cattle farm

Shiga toxin-producing Escherichia coli (STEC) is an important food-borne pathogen, which can cause diseases such as diarrhea, hemorrhagic enteritis, and hemolytic uremic syndrome in humans. Twelve STEC isolates were collected from beeves and feces of commercial animals in China between 2019 and 2020 for this study. In addition to the determination of serotype and Shiga toxin subtype, whole-genome sequencing (WGS) was used for determining phylogenetic relationships, antimicrobial resistance (AMR), virulence genes, and sequence type (ST) of isolates. A total of 27 AMR genes were detected, and each STEC isolate carried more than 10 AMR genes. Eight STEC isolates from ground beef and four STEC isolated from feces were screened. A total of seven serotypes were identified, and one isolate ONT:H10 was undetermined by SeroTypeFinder. Three O157:H7 strains were confirmed and the remaining five serogroups were confirmed as O26:H11, O81:H31, O105:H8, O178:H19, and O136:H12. The phylogenetic analysis showed that STEC isolates of the same serotype or ST were clustered together based on cgMLST. The comparison of the genomes of 157 STEC reference isolates worldwide with our local STEC isolates showed that STEC isolates screened in China represented various collections and could not form a separate cluster but were interspersed among the STEC reference collection, which suggested that several STEC isolates shared a common ancestor irrespective of STEC serotype isolates. cgMLST revealed that isolates of the same O serotype clustered irrespective of their H type. Further investigation is required to determine the pathogenic potential of other serotypes of STEC, particularly in regard to these rare serotypes.

Core Genome Multilocus Sequence Typing Scheme for Improved Characterization and Epidemiological Surveillance of Pathogenic Brucella

Brucellosis poses a significant burden to human and animal health worldwide. Robust and harmonized molecular epidemiological approaches and population studies that include routine disease screening are needed to efficiently track the origin and spread of Brucella strains. Core genome multilocus sequence typing (cgMLST) is a powerful genotyping system commonly used to delineate pathogen transmission routes for disease surveillance and control. Except for Brucella melitensis, cgMLST schemes for Brucella species are currently not established. Here, we describe a novel cgMLST scheme that covers multiple Brucella species. We first determined the phylogenetic breadth of the genus using 612 Brucella genomes. We selected 1,764 genes that were particularly well conserved and typeable in at least 98% of these genomes. We tested the new scheme on 600 genomes and found high agreement with the whole-genome-based single nucleotide polymorphism (SNP) analysis. Next, we applied the scheme to reanalyze the genome of Brucella strains from epidemiologically linked outbreaks. We demonstrated the applicability of the new scheme for high-resolution typing required in outbreak investigations as previously reported with whole-genome SNP methods. We also used the novel scheme to define the global population structure of the genus using 1,322 Brucella genomes. Finally, we demonstrated the possibility of tracing distribution of Brucella strains by performing cluster analysis of cgMLST profiles and found nearly identical cgMLST profiles in different countries. Our results show that sequencing depth of more than 40-fold is optimal for allele calling with this scheme. In summary, this study describes a novel Brucella-wide cgMLST scheme that is applicable in Brucella molecular epidemiology and helps in accurately tracking and thus controlling the sources of infection. The scheme is publicly accessible and should represent a valuable resource for laboratories with limited computational resources and bioinformatics expertise.

Genetic diversity and spatial distribution of Burkholderia mallei by core genome-based multilocus sequence typing analysis

Burkholderia mallei is the etiological agent of glanders, a highly contagious and often fatal disease in equids. Due to the high genetic clonality of B. mallei, high-resolution typing assays are necessary to differentiate between individual strains. Here we report on the development and validation of a robust and reproducible core genome-based Multi Locus Sequence Typing Assay (cgMLST) for B. mallei, which is based on 3328 gene targets and enables high-resolution typing at the strain level. The assay was validated using a set of 120 B. mallei genomes from public databases and 23 newly sequenced outbreak strains from in-house strain collections. In this cgMLST analysis, strains from different geographic regions were clearly distinguished by at least 70 allele differences, allowing spatial clustering while closely related and epidemiologically related strains were separated by only zero to three alleles. Neither the different sequencing technologies nor the assembly strategies had an influence on the cgMLST results. The developed cgMLST is highly robust, reproducible and can be used for outbreak investigations, source tracking and molecular characterization of new B. mallei isolates.

Functional fungal communities dominate wood decomposition and are modified by wood traits in a subtropical forest

Wood decomposition is a fundamental process of the carbon cycle in forest ecosystems and differs under varying environmental conditions. However, it remains unclear whether exposure situation and litter removal affect wood decomposition, especially in subtropical forests. Therefore, we chose wood from four dominant species and carried out an experiment with treatments consisting of placing wood in ground contact with and without litter input and above ground exposure. The experiment was performed for 2.5 consecutive years in the subtropical forest of Southwest China to reveal the potential effects of microenvironmental changes due to above ground exposure and nutrient input changes due to litter removal. In this study, neither above ground exposure nor litter removal significantly changed the fungal communities, microbial respiration rates or decomposition rates of the wood, but significant differences among tree species were observed. The abundance of Ascomycota (70.2%) was higher than that of Basidiomycota (24.3%), and there was a significant negative relationship between their abundances, suggesting competition. Moreover, negative (Ascomycota) and positive (Basidiomycota) relationships with microbial respiration and explained 21.5 and 25.5% of the variation in microbial respiration, respectively. The wood density was directly controlled by the sugar, cellulose, and lignin contents and influenced the water content in the wood. The abundances of saprotrophic and pathotrophic fungi were significantly and directly regulated by the water content of the wood. The abundance of pathotrophic fungi was unaffected by wood traits, but these fungi may limit saprotrophic fungal colonization, thereby affecting microbial respiration and decomposition processes. We confirmed that the saprotrophic fungal abundance, rather than fungal diversity, determined wood microbial respiration. These results are of great significance for the comprehensive assessment of wood decomposition and the carbon cycle in subtropical forests, although long-term fungal community dynamics and decomposition rates under different conditions require further study.

Genetic Characterization of Brucella spp.: Whole Genome Sequencing-Based Approach for the Determination of Multiple Locus Variable Number Tandem Repeat Profiles

Brucellosis is an important zoonosis that is emerging in some regions of the world, gaining increased relevance with the inclusion of the causing agent Brucella spp. in the class B bioterrorism group. Until now, multi-locus VNTR Analysis (MLVA) based on 16 loci has been considered as the gold standard for Brucella typing. However, this methodology is laborious, and, with the rampant release of Brucella genomes, the transition from the traditional MLVA to whole genome sequencing (WGS)-based typing is on course. Nevertheless, in order to avoid a disruptive transition with the loss of massive genetic data obtained throughout the last decade and considering that the transition timings will vary considerably among different countries, it is important to determine WGS-based MLVA alleles of the nowadays sequenced genomes. On this regard, we aimed to evaluate the performance of a Python script that had been previously developed for the rapid in silico extraction of the MLVA alleles, by comparing it to the PCR-based MLVA procedure over 83 strains from different Brucella species. The WGS-based MLVA approach detected 95.3% of all possible 1,328 hits (83 strains×16 loci) and showed an agreement rate with the PCR-based MLVA procedure of 96.4% for MLVA-16. According to our dataset, we suggest the use of a minimal depth of coverage of ~50x and a maximum number of ~200 contigs as guiding “boundaries” for the future application of the script. In conclusion, the evaluated script seems to be a very useful and robust tool for the in silico determination of MLVA profiles of Brucella strains, allowing retrospective and prospective molecular epidemiological studies, which are important for maintaining an active epidemiological surveillance of brucellosis.

SARS-CoV-2 surveillance in Italy through phylogenomic inferences based on Hamming distances derived from pan-SNPs, -MNPs and -InDels

BACKGROUND

Faced with the ongoing global pandemic of coronavirus disease, the 'National Reference Centre for Whole Genome Sequencing of microbial pathogens: database and bioinformatic analysis' (GENPAT) formally established at the 'Istituto Zooprofilattico Sperimentale dell'Abruzzo e del Molise' (IZSAM) in Teramo (Italy) is in charge of the SARS-CoV-2 surveillance at the genomic scale. In a context of SARS-CoV-2 surveillance requiring correct and fast assessment of epidemiological clusters from substantial amount of samples, the present study proposes an analytical workflow for identifying accurately the PANGO lineages of SARS-CoV-2 samples and building of discriminant minimum spanning trees (MST) bypassing the usual time consuming phylogenomic inferences based on multiple sequence alignment (MSA) and substitution model.

RESULTS

GENPAT constituted two collections of SARS-CoV-2 samples. The first collection consisted of SARS-CoV-2 positive swabs collected by IZSAM from the Abruzzo region (Italy), then sequenced by next generation sequencing (NGS) and analyzed in GENPAT (n = 1592), while the second collection included samples from several Italian provinces and retrieved from the reference Global Initiative on Sharing All Influenza Data (GISAID) (n = 17,201). The main results of the present work showed that (i) GENPAT and GISAID detected the same PANGO lineages, (ii) the PANGO lineages B.1.177 (i.e. historical in Italy) and B.1.1.7 (i.e. 'UK variant') are major concerns today in several Italian provinces, and the new MST-based method (iii) clusters most of the PANGO lineages together, (iv) with a higher dicriminatory power than PANGO lineages, (v) and faster that the usual phylogenomic methods based on MSA and substitution model.

CONCLUSIONS

The genome sequencing efforts of Italian provinces, combined with a structured national system of NGS data management, provided support for surveillance SARS-CoV-2 in Italy. We propose to build phylogenomic trees of SARS-CoV-2 variants through an accurate, discriminant and fast MST-based method avoiding the typical time consuming steps related to MSA and substitution model-based phylogenomic inference.

Emerging diversity and ongoing expansion of the genus Brucella

Remarkable genetic diversity and breadth of host species has been uncovered in the Brucella genus over the past decade, fundamentally changing our concept of what it means to be a Brucella. From ocean fishes and marine mammals, to pond dwelling amphibians, forest foxes, desert rodents, and cave-dwelling bats, Brucella have revealed a variety of previously unknown niches. Classical microbiological techniques have been able to help us classify many of these new strains but at times have limited our ability to see the true relationships among or within species. The closest relatives of Brucella are soil bacteria and the adaptations of Brucella spp. to live intracellularly suggest that the genus has evolved to live in vertebrate hosts. Several recently discovered species appear to have phenotypes that are intermediate between soil bacteria and core Brucella, suggesting that they may represent ancestral traits that were subsequently lost in the traditional species. Remarkably, the broad relationships among Brucella species using a variety of sequence and fragment-based approaches have been upheld when using comparative genomics with whole genomes. Nonetheless, genomes are required for fine-scale resolution of many of the relationships and for understanding the evolutionary history of the genus. We expect that the coming decades will reveal many more hosts and previously unknown diversity in a wide range of environments.

Rapid genome-wide sequence typing of African swine fever virus based on alleles

Rapid and accurate molecular typing of African swine fever virus (ASFV) during outbreaks is important to reveal diversity and sourcing of ASFV. Here we present a new way to perform rapid genome-wide multi-locus sequence typing of ASFV using an allele calling based on gene by gene approach. Using open-accessed chewBBACA software, 41 publicly available ASFV genomes were analyzed to optimize the parameters to find the alleles. Alleles as many as 127 were found for building the phylogenetic trees, which covered more than 60 % of the whole genome. Then the method was used to analyze two ASFV genomes assembled from two metagenomic sequences of a swine whole blood and a swine spleen tissue collected in Wuhan, China. It reveals that the two ASFV genomes are the closest to that of Pig/HLJ/2018 strain and DB/LN/2018 strain, which were isolated earlier in China. This proved that the ASFV in Wuhan originated from the same source causing the earlier outbreaks in Heilongjiang and Liaoning province of China. This method could identify more informative genome regions that could be used for accurate typing than other genome-wide analysis, and with less demand on computing resources. It also showed tolerance to analyze ASFV draft genomes assembled directly from metagenomic sequences. Furthermore, the ASFV-specific genetic markers found by the allele calling could be translated into clinical diagnostics or can be used broadly to identify conserved putative therapeutic candidates.

Development and evaluation of a core genome multilocus sequence typing scheme for Paenibacillus larvae, the deadly American foulbrood pathogen of honeybees

Paenibacillus larvae is the causative agent of the fatal American foulbrood disease in honeybees (Apis mellifera). Strain identification is vital for preventing the spread of the disease. To date, the most accessible and robust scheme to identify strains is the multilocus sequence typing (MLST) method. However, this approach has limited resolution, especially for epidemiological studies. As the cost of whole-genome sequencing has decreased and as it becomes increasingly available to most laboratories, an extended MLST based on the core genome (cgMLST) presents a valuable tool for high-resolution investigations. In this study, we present a standardized, robust cgMLST scheme for P. larvae typing using whole-genome sequencing. A total of 333 genomes were used to identify, validate and evaluate 2419 core genes. The cgMLST allowed fine-scale differentiation between samples that had the same profile using traditional MLST and allowed for the characterization of strains impossible by MLST. The scheme was successfully used to trace a localized Swedish outbreak, where a cluster of 38 isolates was linked to a country-wide beekeeping operation. cgMLST greatly enhances the power of a traditional typing scheme, while preserving the same stability and standardization for sharing results and methods across different laboratories.

Genomic insights into Brucella

Brucellosis is a zoonotic disease caused by certain species of Brucella. Each species has its preferred host animal, though it can infect other animals too. For a longer period, only six classical species were recognized in the genus Brucella. No vaccine is available for human brucellosis. Therefore, human brucellosis can be controlled only by controlling brucellosis in animals. The genus is now expanding with the newly isolated atypical strains from various animals, including marine mammals. Presently, 12 species of Brucella have been recognized. The first genome of Brucella was released in 2002, and today, we have more than 1500 genomes of Brucella spp. isolated worldwide. Multiple genome sequences are available for the major zoonotic species, B. abortus, B. melitensis, and B. suis. The Brucella genome has two chromosomes with the approximate sizes of 2.1 and 1.2 Mbp. The genome of Brucella is highly conserved across all the species at the nucleotide level. One of the unanswered questions is what makes host preference in different species of Brucella. Here, I summarize the recent advancements in the Brucella genomics research.

A Core Genome Multilocus Sequence Typing Scheme for Streptococcus mutans

Streptococcus mutans is one of the primary pathogens responsible for the development of dental caries. Recent whole-genome sequencing (WGS)-based core genome multilocus sequence typing (cgMLST) approaches have been employed in epidemiological studies of specific human pathogens. However, this approach has not been reported in studies of S. mutans Here, we therefore developed a cgMLST scheme for S. mutans We surveyed 199 available S. mutans genomes as a means of identifying cgMLST targets, developing a scheme that incorporated 594 targets from the S. mutans UA159 reference genome. Sixty-eight sequence types (STs) were identified in this cgMLST scheme (cgSTs) in 80 S. mutans isolates from 40 children that were sequenced in this study, compared to 35 STs identified by multilocus sequence typing (MLST). Fifty-six cgSTs (82.35%) were associated with a single isolate based on our cgMLST scheme, which is significantly higher than in the MLST scheme (11.43%). In addition, 58.06% of all MLST profiles with ≥2 isolates were further differentiated by our cgMLST scheme. Topological analyses of the maximum likelihood phylogenetic trees revealed that our cgMLST scheme was more reliable than the MLST scheme. A minimum spanning tree of 145 S. mutans isolates from 10 countries developed based upon the cgMLST scheme highlighted the diverse population structure of S. mutans This cgMLST scheme thus offers a new molecular typing method suitable for evaluating the epidemiological distribution of this pathogen and has the potential to serve as a benchmark for future global studies of the epidemiological nature of dental caries.IMPORTANCE Streptococcus mutans is regarded as a major pathogen responsible for the onset of dental caries. S. mutans can transmit among people, especially within families. In this study, we established a new epidemiological approach to S. mutans classification. This approach can effectively differentiate among closely related isolates and offers superior reliability relative to that of the traditional MLST molecular typing method. As such, it has the potential to better support effective public health strategies centered around this bacterium that are aimed at preventing and treating dental caries.

Genome-scale approach to study the genetic relatedness among Brucella melitensis strains

Brucellosis is an important zoonotic disease that affects both humans and animals. To date, laboratory surveillance is still essentially based on the traditional MLVA-16 methodology and the associated epidemiological information is frequently scarce. Our goal was to contribute to the improvement of Brucella spp. surveillance through the implementation of a whole genome sequencing (WGS) approach. We created a curated ready-to-use species-specific wgMLST scheme enrolling a panel of 2656 targets (http://doi.org/10.5281/zenodo.3575026) and used this schema to perform a retrospective analysis of the genetic relatedness among B. melitensis strains causing human infection in Portugal (a country where brucellosis is an endemic disease) from 2010 to 2018. The strains showed a phylogenetic clustering within genotype II (25 out of 36) and IV (4 out of 36), and shared clades with strains isolated from countries with which Portugal has intense food trading, tourism and similar eating habits, such as Spain, Italy and Greece. In addition, our results point to the identification of strong associations between B. melitensis strains, likely underlying missed "outbreaks" as 22 out of the 36 strains showed genetic linkage with others. In fact, the applied gene-by-gene approach grouped these strains into six genetic clusters each one containing putative epidemiological links. Nevertheless, more studies will be needed in order to define the appropriate range of cut-offs (probable non-static cut-offs) that best illustrate the association between genetic linkage and epidemiological information and may serve as alerts for the health authorities. The release of this freely available and scalable schema contributes to the required technological transition for laboratorial surveillance of brucellosis and will facilitate the assessment of ongoing and future outbreaks in order to prevent the transmission spread.