Core Genome Multilocus Sequence Typing: a Standardized Approach for Molecular Typing of Mycoplasma gallisepticum

Mycoplasma gallisepticum and Mycoplasma synoviae in commercial poultry: current control strategies and future challenges

Mycoplasma gallisepticum (Mg) and Mycoplasma synoviae (Ms) are regarded as the most important avian mycoplasma species for today's chicken and turkey farming industry from clinical and economical perspectives. Control strategies for Mg and Ms have become more efficient due to investments in mycoplasma research over the last 70 years. These investments have contributed to the further implementation of serological and molecular testing, the development of vaccines, and the improvement of antimicrobial treatment strategies. However, the increasing spotlight on welfare, the pressure on prudent use of antimicrobials, and the expected global increase in poultry production, are going to have an impact on the future control of avian mycoplasmas in commercial poultry. In this paper a group of avian mycoplasma experts discuss the future challenges in mycoplasma control considering the background of these expected changes and the relevance for future avian mycoplasma research.

Ribotyping Staphylococcus epidermidis Using Probabilistic Sequence Analysis and Levenshtein Distance Algorithm

Staphylococcus epidermidis (S. epidermidis) live in different human locations and natural environments. For ribotyping S. epidermidis sub-species, 2507 PCR-amplified reads of 16S rRNA genes of S. epidermidis in a public dataset were used for probabilistic sequence analysis. A sequence probability logo (sequence pLogo) as a reference sequence of 16S rRNA genes of S. epidermidis was constructed. Through implementation of Levenshtein Distance algorithm, two 20-base pairs (bp) motifs, commonly present in 2507 PCR-amplified reads, were identified. The top 38 S. epidermidis isolates, which carried 16S rRNA nucleotide domains that were made of different sequences but have high similarity scores to two 20-bp motifs, were found from 11 human, 8 animal, 9 plant and 10 environmental samples, indicating that these two 20-bp motifs were broadly present in diverse S. epidermidis isolates. Thirty-one PCR-amplified reads of 16S rRNA genes, which were currently not in the dataset, were utilized to verify the feasibility of using two 20-bp motifs for ribotyping S. epidermidis sub-species. S. epidermidis S1, S3, but not S2, isolates on the human scalp carried a 20-bp sequence domain with high similarities to a 20-bp motif in the sequence pLogo. The phylogenetic tree showed that S. epidermidis S1, S2 and S3 were not from a single common ancestor. Two newly identified 20-bp motifs here, thus, provided reference nucleotide residues for ribotyping S. epidermidis.

Research Progress in the Development of Vaccines against Mycoplasma gallisepticum and Mycoplasma synoviae

Mycoplasma gallisepticum (MG) and Mycoplasma synoviae (MS) are the primary agents responsible for mycoplasma disease in poultry. MG has been identified as a significant cause of chronic respiratory disease in chickens, while MS has been linked to the development of tenosynovitis, joint swelling and other symptoms in chickens, leading to considerable economic losses for the poultry industry. Unfortunately, there is no specific drug for treatment and vaccination is the most important way to control the disease. There are some different types of vaccines, including live vaccines, inactivated vaccines, sub-unit vaccines and vector vaccines. This paper provides a comprehensive review of the development of vaccines for MG and MS.

gcPathogen: a comprehensive genomic resource of human pathogens for public health

Here, we present the manually curated Global Catalogue of Pathogens (gcPathogen), an extensive genomic resource designed to facilitate rapid and accurate pathogen analysis, epidemiological exploration and monitoring of antibiotic resistance features and virulence factors. The catalogue seamlessly integrates and analyzes genomic data and associated metadata for human pathogens isolated from infected patients, animal hosts, food and the environment. The pathogen list is supported by evidence from medical or government pathogenic lists and publications. The current version of gcPathogen boasts an impressive collection of 1 164 974 assemblies comprising 986 044 strains from 497 bacterial taxa, 4794 assemblies encompassing 4319 strains from 265 fungal taxa, 89 965 assemblies featuring 13 687 strains from 222 viral taxa, and 646 assemblies including 387 strains from 159 parasitic taxa. Through this database, researchers gain access to a comprehensive 'one-stop shop' that facilitates global, long-term public health surveillance while enabling in-depth analysis of genomes, sequence types, antibiotic resistance genes, virulence factors and mobile genetic elements across different countries, diseases and hosts. To access and explore the data and statistics, an interactive web interface has been developed, which can be accessed at https://nmdc.cn/gcpathogen/. This user-friendly platform allows seamless querying and exploration of the extensive information housed within the gcPathogen database.

Development of molecular assays for the analysis of genetic relationships of Mycoplasma iowae

Mycoplasma iowae is a worldwide spread and economically important avian pathogen that mostly infects turkeys. Currently, multi-locus sequence typing (MLST) serves as the gold standard method for strain identification in M. iowae. However, additional robust genotyping methods are required to effectively monitor M. iowae infections and conduct epidemiological investigations. The first aim of this study was to develop genotyping assays with high resolution, that specifically target M. iowae, namely a multiple-locus variable number of tandem-repeats analysis (MLVA) and a core genome multi-locus sequence typing (cgMLST) schema. The second aim was the determination of relationships among a diverse selection of M. iowae strains and clinical isolates with a previous and the newly developed assays. The MLVA was designed based on the analyses of tandem-repeat (TR) regions in the six serotype reference strains (I, J, K, N, Q and R). The cgMLST schema was developed based on the coding sequences (CDSs) common in 95% of the examined 99 isolates. The samples were submitted for a previously published MLST assay for comparison with the developed methods. Out of 94 TR regions identified, 17 alleles were selected for further evaluation by PCR. Finally, seven alleles were chosen to establish the MLVA assay. Additionally, whole genome sequence analyses identified a total of 676 CDSs shared by 95% of the isolates, all of which were included into the developed cgMLST schema. The MLVA discriminated 19 distinct genotypes (GT), while with the cgMLST assay 79 sequence types (ST) could be determined with Simpson's diversity indices of 0.810 (MLVA) and 0.989 (cgMLST). The applied assays consistently identified the same main clusters among the diverse selection of isolates, thereby demonstrating their suitability for various genetic analyses and their ability to yield congruent results.

A core-genome multilocus sequence typing scheme for the detection of genetically related Streptococcus pyogenes clusters

The recently observed increase in invasive Streptococcus pyogenes infections causes concern in Europe. However, conventional molecular typing methods lack discriminatory power to aid investigations of outbreaks caused by S. pyogenes. Therefore, there is an urgent need for high-resolution molecular typing methods to assess genetic relatedness between S. pyogenes isolates. In the current study, we aimed to develop a novel high-resolution core-genome multilocus sequence typing (cgMLST) scheme for S. pyogenes and compared its discriminatory power to conventional molecular typing methods. The cgMLST scheme was designed with the commercial Ridom SeqSphere+ software package. To define a cluster threshold, the scheme was evaluated using publicly available data from nine defined S. pyogenes outbreaks in the United Kingdom. The cgMLST scheme was then applied to 23 isolates from a suspected S. pyogenes outbreak and 117 S. pyogenes surveillance isolates both from the Netherlands. MLST and emm-typing results were used for comparison to cgMLST results. The allelic differences between isolates from defined outbreaks ranged between 6 and 31 for isolates with the same emm-type, resulting in a proposed cluster threshold of <5 allelic differences out of 1,095 target loci. Seven out of twenty-three (30%) isolates from the suspected outbreak had an allelic difference of <2, thereby identifying a potential cluster that could not be linked to other isolates. The proposed cgMLST scheme shows a higher discriminatory ability when compared to conventional typing methods. The rapid and simple analysis workflow allows for extended detection of clusters of potential outbreak isolates and surveillance and may facilitate the sharing of sequencing results between (inter)national laboratories.

Epidemiological investigations and multilocus sequence typing of Mycoplasma gallisepticum collected in China

Mycoplasma gallisepticum (MG) is one of the important pathogens in poultry industry and has led to major economic losses. Understanding the epidemiology is crucial to improve the control and eradication program of MG. This study collected 1,250 chicken samples, including trachea and lung, from China in 2022 to investigate the epidemiology of MG. Among the collected samples, 938 samples were positive for MG infection, resulting in an average positive rate of 75.04%. Additionally, 570 samples were positive for both MG and Mycoplasma synoviae (MS) coinfection, with an average positive rate of 45.60%. A total of 183 MG infection positive samples in this study were selected for genotyping, and the multilocus sequence typing (MLST) method based on 7 housekeeping genes was used. As a result, 183 samples belonged to 11 sequence types (STs), with ST-78 being the most prevalent. After BURST analysis, all 183 sequences were divided into group 3. Besides, 119 reference sequences from database and 183 sequences of this study were selected to construct the phylogenetic tree using the neighbor-joining method. The results revealed that the sequences from China, total 196 sequences, were classified into 4 branches. The findings suggest that the MG strains in China exhibit diverse genotypes, which may be related to international trade and the use of live vaccines. Furthermore, we detected the drug susceptibility of 10 isolated strains randomly, which may be helpful to guide the clinical use of drugs to control MG infection.

Drug sensitivity and genome-wide analysis of two strains of Mycoplasma gallisepticum with different biofilm intensity

Mycoplasma gallisepticum (MG) is one of the major causative agents of chronic respiratory diseases in poultry. The biofilms of MG are highly correlated to its chronic infection. However data on genes involved in biofilm formation ability are still scarse. MG strains with distinct biofilm intensity were screened by crystal violet staining morphotyped and characterized for the drug sensitivity. Two MG strains NX-01 and NX-02 showed contrasted ability to biofilm formation. The biofilm formation ability of NX-01 strain was significantly higher than that of NX-02 strain (p < 0.01). The drug sensitivity test showed that the stronger the ability of MG stain to form biofilms, the weaker its sensitivity to 17 antibiotic drugs. Moreover, putative key genes related to biofilm formation were screened by genome-wide analysis. A total of 13 genes and proteins related to biofilm formation, including ManB, oppA, oppD, PDH, eno, RelA, msbA, deoA, gapA, rpoS, Adhesin P1 precursor, S-adenosine methionine synthetase, and methionyl tRNA synthetase were identified. There were five major discrepancies between the two isolated MG strains and the five NCBI-published MG strains. These findings provide potential targets for inhibiting the formation of biofilm of MG, and lay a foundation for treating chronic infection.

Space-Time Patterns of Poultry Pathogens in the USA: A Case Study of Ornithobacterium rhinotracheale and Pasteurella multocida in Turkey Populations

Respiratory infections caused by Ornithobacterium rhinotrachealis (ORT) and Pasteurella multocida (PM) bacteria are significant threats to the poultry industry by causing economic losses and welfare issues. Due to characterization difficulties and underutilization of epidemiological tools, description of the spatio-temporal spread of these diseases in the field is limited. The objectives of this retrospective observational cross-sectional study were to (a) investigate the existence of space-time clusters (hotspots); and (b) investigate the association between genetic similarity and spatial proximity for both pathogens using molecular typing and a recently developed Core-Genome Multilocus Sequencing Typing (cgMLST) scheme. ORT (n = 103) and PM (n = 69) isolates from confirmed disease outbreaks from one commercial company between 2013 and 2021 were obtained from a veterinary diagnostic laboratory, characterized using a cgMLST scheme and visualized using a minimum spanning tree. Spatio-temporal cluster analysis using SaTScanTM and a Spearman's rank correlation were performed to investigate clustering and any association between allelic diversity and geospatial distance. The cgMLST sequencing revealed three allelic clusters for ORT and thirteen clusters for PM. The spatio-temporal analysis revealed two significant clusters for PM, one with a 259.3 km cluster containing six cases between May and July 2018 and a 9 km cluster containing five cases between February 2019 and February 2021. No spatio-temporal clusters were found for ORT. A weak negative correlation between allelic diversity and geospatial distance was observed for ORT (r = -0.04, p < 0.01) and a weak positive correlation was observed for PM (r = 0.11, p < 0.01). This study revealed regional spatio-temporal clusters for PM in commercial turkey sites between 2018 and 2021 and provided additional insight into bacterial strain subgroups and the geographical spread of ORT and PM over time.

Hydroxytyrosol mitigates Mycoplasma gallisepticum-induced pulmonary injury through downregulation of the NF-κB/NLRP3/IL-1β signaling pathway in chicken

In this study, the anti-inflammatory and antiapoptotic effects of hydroxytyrosol (HT) in Mycoplasma gallisepticum (MG)-infected chicken were investigated, and the underlying molecular mechanisms were explored. The results revealed severe ultrastructural pathological changes after MG infection in the lung tissue of chicken, including inflammatory cell infiltration, thickening of the lung chamber wall, visible cell swelling, mitochondrial cristae rupture, and ribosome shedding. MG possibly activated the nuclear factor κB (NF-κB)/nucleotide-binding oligomerization domain-like receptor protein 3 (NLRP3)/interleukin (IL)-1β signaling pathway in the lung. However, HT treatment significantly ameliorated MG-induced pathological damage of the lung. HT reduced the magnitude of pulmonary injury after MG infection by reducing apoptosis and releasing the proinflammatory factors. Compared with the MG-infected group, the HT-treated group exhibited significant inhibition of the expression of NF-κB/NLRP3/IL-1β signaling-pathway-related genes; for example, the expressions of NF-κB, NLRP3, caspase-1, IL-1β, IL-2, IL-6, IL-18, and TNF-α significantly decreased (P < 0.01 or <0.05). In conclusion, HT effectively inhibited MG-induced inflammatory response and apoptosis and protected the lung by blocking the activation of NF-κB/NLRP3/IL-1β signaling pathway and reducing the damage caused by MG infection in chicken. This study revealed that HT may be a suitable and effective anti-inflammatory drug against MG infection in chicken.

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.

A core genome multilocus sequence typing scheme for Mycoplasma hyorhinis

Mycoplasma (M.) hyorhinis is a commensal and pathobiont residing in the upper respiratory tract in swine and with the ability to spread systemically, mainly causing polyserositis and polyarthritis in nursery pigs. Since little is known on the epidemiology of M. hyorhinis infection, whole genome sequences of 73 strains isolated from pigs in Austria (n = 71) and Germany (n = 2), that have been isolated from clinically affected pigs during routine diagnostics, and publicly available genomes of eight M. hyorhinis strains were analyzed in the presented study. For this purpose, a core genome multi locus sequence typing (cgMLST) scheme encompassing 453 target genes was developed using the Ridom© SeqSphere + software. Results were compared to two previously described conventional MLST schemes and to a core genome single nucleotide polymorphism (cgSNP) analysis approach. Core genome MLST showed high diversity among the M. hyorhinis strains studied and while certain isolates from one farm or a single animal formed cgMLST clusters (≤ 8 allele differences), no isolates with identical allele profiles were identified. In addition, cgMLST had superior discriminatory power (Simpson's ID = 0.995) over conventional MLST (Simpson's ID = 0.952 and 0.985), while demonstrating a lack of congruence between conventional MLST and genome-wide relationship. Core genome SNP results were highly congruent with cgMLST results but lacked in resolution when comparing closely related isolates. Thus, cgMLST is the most suitable method for epidemiological investigations such as outbreak analysis, and to gain insights into M. hyorhinis population structure.

Complete Sequence-Based Genotyping of mgc2/pvpA Genes in Chicken-Derived Mycoplasma gallisepticum Isolates of Iran

Mycoplasma gallisepticum (MG) is a major pathogen of the poultry industry throughout the world. MG causes chronic respiratory disease in chickens and infectious sinusitis in turkeys. Despite constant improvements in the biosecurity of the poultry industry in Iran, MG infection still occurs and causes significant economic issues. To evaluate genetic variability, 10 Iranian MG isolates along with 17 available sequences were characterized by gene-targeted sequencing (GTS) analysis of complete mgc2/pvpA genes. According to the findings, 21 different sequence types within the sample set of 27 strains were typed by this method. The discriminatory power of this typing assay was established to be 0.97. Although no insertions and deletions of nucleotides were observed in the mgc2 gene among the Iranian strains, different lengths of pvpA genes with 1086, 1095, and 1101 nucleotides were detected within direct repeats (DRs) 1 and 2. Generally, eight tetrapeptides Pro-Arg-Pro-Met/Gln/Asn were found in the DRs of PvpA. Analysis of the carboxyl ends of PvpA proteins exhibited various repeats of prolines. In the phylogenetic tree of partial and complete mgc2/pvpA genes, all Iranian MG isolates were clustered into two distinct groups. Because this typing assay could provide a higher discriminatory power than the previously reported GTS scheme of partial mgc2/ pvpA genes, these results can be considered a blueprint for future national control and diagnostic strategies. Furthermore, consistent surveillance with larger datasets will be needed to clarify the epidemiologic characteristics of MG outbreaks in different poultry hosts.

Application of Four Genotyping Methods to Mycoplasma bovis Isolates Derived from Western Canadian Feedlot Cattle

Mycoplasma bovis is a significant pathogen of feedlot cattle, responsible for chronic pneumonia and polyarthritis syndrome (CPPS). M. bovis isolates (n = 129) were used to compare four methods of phylogenetic analysis and to determine if the isolates' genotypes were associated with phenotypes. Metadata included the health status of the animal from which an isolate was derived (healthy, diseased, or dead), anatomical location (nasopharynx, lung, or joint), feedlot, and production year (2006 to 2018). Four in silico phylogenetic typing methods were used: multilocus sequence typing (MLST), core genome MLST (cgMLST), core genome single nucleotide variant (cgSNV) analysis, and whole-genome SNV (wgSNV) analysis. Using Simpson's diversity index (D) as a proxy for resolution, MLST had the lowest resolution (D = 0.932); cgSNV (D = 0.984) and cgMLST (D = 0.987) generated comparable results; and wgSNV (D = 1.000) provided the highest resolution. Visual inspection of the minimum spanning trees found that the memberships of the clonal complexes and clades had similar structural appearances. Although MLST had the lowest resolution, this methodology was intuitive and easy to apply, and the PubMLST database facilitates the comparison of sequence types across studies. The cg methods had higher resolution than MLST, and the graphical interface software was user-friendly for nonbioinformaticians, but the proprietary software is relatively expensive. The wgSNV approach was the most robust for processing poor-quality sequence data while offering the highest resolution; however, application of its software requires specialized training. None of the four methods could associate genotypes with phenotypes.

Clinical expression, epidemiology and monitoring of Mycoplasma gallisepticum and Mycoplasma synoviae: an update

Mycoplasma gallisepticum (MG) and Mycoplasma synoviae (MS) are of clinical and economic importance for the global poultry industry. Many countries and integrations are involved in monitoring programs to control both mycoplasma species. This review provides an extensive historic overview of the last seven decades on the development of the knowledge regarding the factors that influence the clinical expression of the disease, the epidemiology and monitoring of both MG and MS. This includes the detection of new virulent strains, studies unravelling the transmission routes, survival characteristics and the role of other avian hosts. Also the role of molecular typing tests in unravelling epidemiology, and factors that complicate the interpretation of test results such as heterologous mycoplasma infections, use of heterologous oil-emulsion vaccines, use of antibiotic treatments, occurrence of MG and MS strains with low virulence, and last but not least the use of live and/or inactivated MG and MS vaccines are discussed.

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.

Molecular Differentiation of Mycoplasma gallisepticum Outbreaks: A Last Decade Study on Italian Farms Using GTS and MLST

Mycoplasma gallisepticum (MG) infects many avian species and leads to significant economic losses in the poultry industry. Transmission of this pathogen occurs both horizontally and vertically, and strategies to avoid the spread of MG rely on vaccination and the application of biosecurity measures to maintain breeder groups as pathogen-free. Two live attenuated MG vaccine strains are licensed in Italy: 6/85 and ts-11. After their introduction, the implementation of adequate genotyping tools became necessary to distinguish between field and vaccine strains and to guarantee proper infection monitoring activity. In this study, 40 Italian MG isolates collected between 2010–2019 from both vaccinated and unvaccinated farms were genotyped using gene-targeted sequencing (GTS) of the cythadesin gene mgc2 and multilocus sequence typing (MLST) based on six housekeeping genes. The discriminatory power of GTS typing ensures 6/85-like strain identification, but the technique does not allow the identification ts-11 strains; conversely, MLST differentiates both vaccine strains, describing more detailed interrelation structures. Our study describes MG genetic scenario within a mixed farming context. In conclusion, the use of adequate typing methods is essential to understand the evolutionary dynamics of MG strains in a particular area and to conduct epidemiological investigations in the avian population.

Investigating extradomiciliary transmission of tuberculosis: An exploratory approach using social network patterns of TB cases and controls and the genotyping of Mycobacterium tuberculosis

Extradomiciliary contacts have been overlooked in the study of TB transmission due to difficulties in identifying actual contacts in large populations. Complex network analysis provides a framework to model the structure of contacts, specially extradomiciliary ones. We conducted a study of incident sputum-positive TB cases and healthy controls occurring in a moderate TB burden city. Cases and controls were interviewed to obtain data regarding the usual locations of residence, work, study, and leisure. Mycobacterium tuberculosis isolated from sputum was genotyped. The collected data were used to build networks based on a framework of putative social interactions indicating possible TB transmission. A user-friendly open source environment (GraphTube) was setup to extract information from the collected data. Networks based on the likelihood of patient-patient, patient-healthy, and healthy-healthy contacts were setup, depending on a constraint of geographical distance of places attended by the volunteers. Using a threshold for the geographical distance of 300 m, the differences between TB cases and controls are revealed. Several clusters formed by social network nodes with high genotypic similarity were characterized. The developed framework provided consistent results and can be used to support the targeted search of potentially infected individuals and to help to understand the TB transmission.

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.

The core genome multi-locus sequence typing of Mycoplasma anserisalpingitidis

BACKGROUND

Mycoplasma anserisalpingitidis is a waterfowl pathogen that mainly infects geese, can cause significant economic losses and is present worldwide. With the advance of whole genome sequencing technologies, new methods are available for the researchers; one emerging methodology is the core genome Multi-Locus Sequence Typing (cgMLST). The core genome contains a high percentage of the coding DNA sequence (CDS) set of the studied strains. The cgMLST schemas are powerful genotyping tools allowing for the investigation of potential epidemics, and precise and reliable classification of the strains. Although whole genome sequences of M. anserisalpingitidis strains are available, to date, no cgMLST schema has been published for this species.

RESULTS

In this study, Illumina short reads of 81 M. anserisalpingitidis strains were used, including samples from Hungary, Poland, Sweden, and China. Draft genomes were assembled with the SPAdes software and analysed with the online available chewBBACA program. User made modifications in the program enabled analysis of mycoplasmas and provided similar results as the conventional SeqSphere+ software. The threshold of the presence of CDS in the strains was set to 93% due to the quality of the draft genomes, resulting in the most accurate and robust schema. Three hundred thirty-one CDSs constituted our cgMLST schema (representing 42,77% of the whole CDS set of M. anserisalpingitidis ATCC BAA-2147), and a Neighbor joining tree was created using the allelic profiles. The correlation was observed between the strains' cgMLST profile and geographical origin; however, strains from the same integration but different locations also showed close relationship. Strains isolated from different tissue samples of the same animal revealed highly similar cgMLST profiles.

CONCLUSIONS

The Neighbor joining tree from the cgMLST schema closely resembled the real-life spatial and temporal relationships of the strains. The incongruences between background data and the cgMLST profile in the strains from the same integration can be because of the higher probability of contacts between the flocks. This schema can help with the epidemiological investigation and can be used as a basis for further studies.

Development of a Multilocus Sequence Typing Assay for Mycoplasma gallisepticum

Mycoplasma gallisepticum (MG) is the most pathogenic avian mycoplasma species. It affects commercial and noncommercial poultry and wild birds. Current MG sequence typing methods rely on the partial sequence of one or more surface antigen genes. Multilocus sequence typing (MLST), a widely used typing method for many human and animal pathogens, relies on conserved housekeeping genes. Recently, MLST assays have been developed for Mycoplasma synoviae (MS) and Mycoplasma iowae. Additionally, a whole genome-based core genome MLST (cgMLST) assay has been developed for MG and MS. However, cgMLST can be implemented only on pure isolates and cannot be applied to clinical samples. Here, we have developed a seven-locus-based MLST scheme for MG that can be applied directly on clinical samples without the need for isolation. The seven loci were selected from 425 genes recently used for the cgMLST assay. A total of 101 diverse MG samples, including isolates and clinical samples, were typed with the newly developed seven-locus MLST. The phylogeny and discriminatory power of the seven-locus MLST were evaluated and compared with the cgMLST and gene-targeted sequencing methods currently used for MG sequence typing. The seven-locus MLST provided optimum discriminatory power and congruent phylogeny to cgMLST. Additionally, a database for MG MLST was created and is currently available for public use online. This assay will increase accessibility to MG sequence typing and provide a stable and expandable nomenclature compatible with cgMLST. The seven-locus MLST assay represents an important tool for epidemiologic investigation of MG that can contribute to better control and eradication efforts.

Development of mismatch amplification mutation assay for the rapid differentiation of Mycoplasma gallisepticum K vaccine strain from field isolates

Mycoplasma gallisepticum causes respiratory diseases and reproduction disorders in turkeys and chickens. The infection has considerable economic impact due to reduced meat and egg production. Because elimination programmes are not feasible in a large number of poultry farms, vaccination remains the only effective measure of disease control. Differentiating vaccine strains from field isolates is necessary in the control of vaccination programmes and diagnostics. The aim of this study was to develop a polymerase chain reaction based mismatch amplification mutation assay (MAMA) for the discrimination of K vaccine strain (K 5831, Vaxxinova Japan K.K.). After determining the whole genome sequence of the K strain, primers were designed to detect seven different vaccine-specific single nucleotide polymorphisms. After evaluating preliminary results, the MAMA-K-fruA test detecting a single guanine-adenine substitution within the fruA gene (G88A) was found to be the most applicable assay to distinguish the K vaccine strain from field isolates. The detected K strain-specific single nucleotide polymorphism showed genetic stability after serial passage in vitro, but this stability test should still be evaluated in vivo as well, investigating a large number of K strain re-isolates. The MAMA-K-fruA assay was tested on a total of 280 culture and field samples. The designed assay had 10² and 10³ template copy number/µl sensitivity in melt-curve analysis based and agarose-gel based assays, respectively, and showed no cross reaction with other avian Mycoplasma species. The new MAMA provides a time- and cost-effective molecular tool for the control of vaccination programmes and for diagnostics.

Development of Molecular Methods for Rapid Differentiation of Mycoplasma gallisepticum Vaccine Strains from Field Isolates

Mycoplasma gallisepticum is among the most economically significant mycoplasmas causing production losses in poultry. Seven melt-curve and agarose gel-based mismatch amplification mutation assays (MAMAs) and one PCR are provided in the present study to distinguish the M. gallisepticum vaccine strains and field isolates based on mutations in the crmA, gapA, lpd, plpA, potC, glpK, and hlp2 genes. A total of 239 samples (M. gallisepticum vaccine and type strains, pure cultures, and clinical samples) originating from 16 countries and from at least eight avian species were submitted to the presented assays for validation or in blind tests. A comparison of the data from 126 samples (including sequences available at GenBank) examined by the developed assays and a recently developed multilocus sequence typing assay showed congruent typing results. The sensitivity of the melt-MAMA assays varied between 10¹ and 10⁴ M. gallisepticum template copies/reaction, while that of the agarose-MAMAs ranged from 10³ to 10⁵ template copies/reaction, and no cross-reactions occurred with other Mycoplasma species colonizing birds. The presented assays are also suitable for discriminating multiple strains in a single sample. The developed assays enable the differentiation of live vaccine strains by targeting two or three markers/vaccine strain; however, considering the high variability of the species, the combined use of all assays is recommended. The suggested combination provides a reliable tool for routine diagnostics due to the sensitivity and specificity of the assays, and they can be performed directly on clinical samples and in laboratories with basic PCR equipment.

Virulence gene profiles and phylogeny of Shiga toxin-positive Escherichia coli strains isolated from FDA regulated foods during 2010-2017

Illnesses caused by Shiga toxin-producing Escherichia coli (STECs) can be life threatening, such as hemolytic uremic syndrome (HUS). The STECs most frequently identified by USDA's Microbiological Data Program (MDP) carried toxin gene subtypes stx1a and/or stx2a. Here we described the genome sequences of 331 STECs isolated from foods regulated by the FDA 2010-2017, and determined their genomic identity, serotype, sequence type, virulence potential, and prevalence of antimicrobial resistance. Isolates were selected from the MDP archive, routine food testing by FDA field labs (ORA), and food testing by a contract company. Only 276 (83%) strains were confirmed as STECs by in silico analysis. Foods from which STECs were recovered included cilantro (6%), spinach (25%), lettuce (11%), and flour (9%). Phylogenetic analysis using core genome MLST revealed these STEC genomes were highly variable, with some clustering associated with ST types and serotypes. We detected 95 different sequence types (ST); several ST were previously associated with HUS: ST21 and ST29 (O26:H11), ST11 (O157:H7), ST33 (O91:H14), ST17 (O103:H2), and ST16 (O111:H-). in silico virulome analyses showed ~ 51% of these strains were potentially pathogenic [besides stx gene they also carried eae (25%) or 26% saa (26%)]. Virulence gene prevalence was also determined: stx1 only (19%); stx2 only (66%); and stx1/sxt2 (15%). Our data form a new WGS dataset that can be used to support food safety investigations and monitor the recurrence/emergence of E. coli in foods.

Genotyping Mycoplasma gallisepticum by multilocus sequence typing

Mycoplasma gallisepticum causes chronic respiratory disease and reproductive disorders in many bird species, resulting in considerable economic losses to the poultry industry. Maintenance of M. gallisepticum-free flocks is the most adequate method to control infection. To this end, monitoring systems and vaccination programs with live vaccine strains are applied worldwide. There is strong demand for efficient epidemiological investigation tools to distinguish M. gallisepticum strains in order to control disease. Up to now, multilocus sequence typing (MLST) has been regarded as gold standard for genotyping bacteria due to its good reproducibility and high discriminatory power. The aim of this study was to develop an MLST assay which can determine phylogenetic distances between M. gallisepticum strains. After analysing more than 30 housekeeping genes, six loci (atpG, dnaA, fusA, rpoB, ruvB, uvrA) were selected for the MLST assay due to their genomic location and high diversity. Examination of 130 M. gallisepticum strains with this MLST method yielded 57 unique sequence types (STs) with a 0.96 Simpson's index of diversity. Considering the large number of STs and high diversity index, this MLST method was found to be appropriate to discriminate M. gallisepticum strains. In addition, the developed method was shown to be suitable for epidemiological investigations, as it confirmed linkage between related strains from outbreaks in different farms. Besides, MLST also suggested high impact of extensive international trade on the spread of different M. gallisepticum strains. Furthermore this method can be used for differentiation among vaccine and field strains.

Multiple Recombination Events Drive the Current Genetic Structure of Xanthomonas perforans in Florida

Prior to the identification of Xanthomonas perforans associated with bacterial spot of tomato in 1991, X. euvesicatoria was the only known species in Florida. Currently, X. perforans is the Xanthomonas sp. associated with tomato in Florida. Changes in pathogenic race and sequence alleles over time signify shifts in the dominant X. perforans genotype in Florida. We previously reported recombination of X. perforans strains with closely related Xanthomonas species as a potential driving factor for X. perforans evolution. However, the extent of recombination across the X. perforans genomes was unknown. We used a core genome multilocus sequence analysis approach to identify conserved genes and evaluated recombination-associated evolution of these genes in X. perforans. A total of 1,356 genes were determined to be "core" genes conserved among the 58 X. perforans genomes used in the study. Our approach identified three genetic groups of X. perforans in Florida based on the principal component analysis (PCA) using core genes. Nucleotide variation in 241 genes defined these groups, that are referred as Phylogenetic-group Defining (PgD) genes. Furthermore, alleles of many of these PgD genes showed 100% sequence identity with X. euvesicatoria, suggesting that variation likely has been introduced by recombination at multiple locations throughout the bacterial chromosome. Site-specific recombinase genes along with plasmid mobilization and phage associated genes were observed at different frequencies in the three phylogenetic groups and were associated with clusters of recombinant genes. Our analysis of core genes revealed the extent, source, and mechanisms of recombination events that shaped the current population and genomic structure of X. perforans in Florida.

A Core Genome Multilocus Sequence Typing Scheme for Enterococcus faecalis

Among enterococci, Enterococcus faecalis occurs ubiquitously, with the highest incidence of human and animal infections. The high genetic plasticity of E. faecalis complicates both molecular investigations and phylogenetic analyses. Whole-genome sequencing (WGS) enables unraveling of epidemiological linkages and putative transmission events between humans, animals, and food. Core genome multilocus sequence typing (cgMLST) aims to combine the discriminatory power of classical multilocus sequence typing (MLST) with the extensive genetic data obtained by WGS. By sequencing a representative collection of 146 E. faecalis strains isolated from hospital outbreaks, food, animals, and colonization of healthy human individuals, we established a novel cgMLST scheme with 1,972 gene targets within the Ridom SeqSphere⁺ software. To test the E. faecalis cgMLST scheme and assess the typing performance, different collections comprising environmental and bacteremia isolates, as well as all publicly available genome sequences from the NCBI and SRA databases, were analyzed. In more than 98.6% of the tested genomes, >95% good cgMLST target genes were detected (mean, 99.2% target genes). Our genotyping results not only corroborate the known epidemiological background of the isolates but exceed previous typing resolution. In conclusion, we have created a powerful typing scheme, hence providing an international standardized nomenclature that is suitable for surveillance approaches in various sectors, linking public health, veterinary public health, and food safety in a true One Health fashion.

Development of Mycoplasma synoviae (MS) core genome multilocus sequence typing (cgMLST) scheme

Mycoplasma synoviae (MS) is a poultry pathogen with reported increased prevalence and virulence in recent years. MS strain identification is essential for prevention, control efforts and epidemiological outbreak investigations. Multiple multilocus based sequence typing schemes have been developed for MS, yet the resolution of these schemes could be limited for outbreak investigation. The cost of whole genome sequencing became close to that of sequencing the seven MLST targets; however, there is no standardized method for typing MS strains based on whole genome sequences. In this paper, we propose a core genome multilocus sequence typing (cgMLST) scheme as a standardized and reproducible method for typing MS based whole genome sequences. A diverse set of 25 MS whole genome sequences were used to identify 302 core genome genes as cgMLST targets (35.5% of MS genome) and 44 whole genome sequences of MS isolates from six countries in four continents were used for typing applying this scheme. cgMLST based phylogenetic trees displayed a high degree of agreement with core genome SNP based analysis and available epidemiological information. cgMLST allowed evaluation of two conventional MLST schemes of MS. The high discriminatory power of cgMLST allowed differentiation between samples of the same conventional MLST type. cgMLST represents a standardized, accurate, highly discriminatory, and reproducible method for differentiation between MS isolates. Like conventional MLST, it provides stable and expandable nomenclature, allowing for comparing and sharing the typing results between different laboratories worldwide.