Improving kinship probability in analysis of ancient skeletons using identity SNPs and MPS technology.

In forensic kinship analysis and human identification cases, analysis of STRs is the gold standard. When badly preserved ancient DNA is used for kinship analysis, short identity SNPs are more promising for successful amplification. In this work, kinship analysis was performed on two skeletons from the Early Middle Ages. The surface contaminants of petrous bones were removed by chemical cleaning and UV irradiation; DNA was isolated through full demineralization and purified in an EZ1 Advanced XL machine. The PowerQuant kit was used to analyze DNA yield and degradation, and on average, 17 ng DNA/g of petrous bone was obtained. Both skeletons were typed in duplicate for STR markers using the Investigator EssplexPlus SE QS kit, and comparison of partial consensus genotypes showed shared allelic variants at most loci amplified, indicating close kinship. After statistical calculation, the full-sibling kinship probability was too low for kinship confirmation, and additional analyses were performed with PCR-MPS using the Precision ID Identity Panel. The HID Ion Chef Instrument was used to prepare the libraries and for templating and the Ion GeneStudio S5 System for sequencing. Analysis of identity SNPs produced full genetic profiles from both skeletons. For combined likelihood ratio (LR) calculation, the product rule was used, combining LR for STRs and LR for SNPs, and a combined LR of 3.3 × 107 (corresponding to a full-sibling probability of 99.999997%) was calculated. Through the SNP PCR-MPS that followed the STR analysis, full-sibling kinship between the ancient skeletons excavated from an early medieval grave was confirmed.

High allele discrimination in the typing of single nucleotide polymorphisms of miRNA.

MicroRNAs (miRNAs) belonging to the same family have similar sequences and are difficult to identify. Herein, we report the reverse transcription-hairpin-probe-polymerase chain reaction (RT-Hpro-PCR) technique, which utilises a reverse transcription (RT) primer containing a 5'-end deoxyribonucleic acid (DNA) tag, to detect miRNAs with similar sequences. This strategy follows a two-step RT-PCR method using 6-7-mer RT-primers with a ~ 10-mer tag sequence at the 5'-end and a probe with a hairpin structure (Hpro), including two C-bulges, attached. The findings demonstrate that the specificity of RT could be increased by shortening the complementary part of the RT primer containing a different base, wherein the PCR could successfully progress with the use of 5'-end DNA tag because of an increase in the length of the hybridised tagged primer. This study shows the potential of RT-Hpro-PCR to precisely detect miRNAs with similar sequences, which could help explore the roles of miRNAs in several biological processes.

Genotypic differences in CC224, CC363, CC449 and CC446 of Moraxella catarrhalis isolates based on whole genome SNP, MLST and PFGE typing.

Understanding the evolutionary path of M. catarrhalis from macrolide-susceptible to macrolide-resistant organism, is important for hindering macrolide resistance from propagation. Multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE) and whole genome SNP typing (WGST), as useful and practical typing tools, have both advantages and disadvantages. We studied the utility of these 3 typing methods, including the level of agreement, consistency and drawbacks, in characterizing M. catarrhalis clones and clonal complexes. We focused on four clonal complexes [CC224, CC363, CC449 (CCN10) and CC446 (CCN08)] and found that PFGE and WGST had a high level of agreement and a proper consistency of the same clone or very closely related clones, while MLST is less discriminatory for different clones. Furthermore, we also established an evolutionary distance cut-off value for "The same clone". Moreover, we detected macrolide-resistant M. catarrhalis in CC224, which had previously been considered as a macrolide-susceptible clonal complex. A higher number of isolates belonged to ST215 compared to ST446, implying that ST215 is more likely to be the primary founder. Our study also demonstrated that all the four clonal complexes belong to the M. catarrhalis lineage 1, which is considered to be related to increased virulence potential and serum resistance. We also observed that copB II was highly related to CC449 and LOS type B was mainly confined in CC224. In conclusion, these findings provide further insight into the evolutionary characteristics of M. catarrhalis.

Genome-wide SNP typing of ancient DNA: Determination of hair and eye color of Bronze Age humans from their skeletal remains.

OBJECTIVE: A genome-wide high-throughput single nucleotide polymorphism (SNP) typing method was tested with respect of the applicability to ancient and degraded DNA. The results were compared to mini-sequencing data achieved through single base extension (SBE) typing. The SNPs chosen for the study allow to determine the hair colors and eye colors of humans. MATERIAL AND METHODS: The DNA samples were extracted from the skeletal remains of 59 human individuals dating back to the Late Bronze Age. The 3,000 years old bones had been discovered in the Lichtenstein Cave in Lower Saxony, Germany. The simultaneous typing of 24 SNPs for each of the ancient DNA samples was carried out using the 192.24 Dynamic Array™ by Fluidigm®. RESULTS: Thirty-eight of the ancient samples (=64%) revealed full and reproducible SNP genotypes allowing hair and eye color phenotyping. In 10 samples (=17%) at least half of the SNPs were unambiguously determined, in 11 samples (=19%) the SNP typing failed. For 23 of the 59 individuals, a comparison of the SNP typing results with genotypes from an earlier performed SBE typing approach was possible. The comparison confirmed the full concordance of the results for 90% of the SNP typings. In the remaining 10% allelic dropouts were identified. DISCUSSION: The high genotyping success rate could be achieved by introducing modifications to the preamplification protocol mainly by increasing the DNA input and the amplification cycle number. The occurrence of allelic dropouts indicates that a further increase of DNA input to the preamplification step is desirable.

Whole-genome sequence-based comparison and profiling of virulence-associated genes of Bacillus cereus group isolates from diverse sources in Japan.

BACKGROUND: The complete genome sequences of 44 Bacillus cereus group isolates collected from diverse sources in Japan were analyzed to determine their genetic backgrounds and diversity levels in Japan. Multilocus sequence typing (MLST) and core-genome single-nucleotide polymorphism (SNP) typing data from whole-genome sequences were analyzed to determine genetic diversity levels. Virulence-associated gene profiles were also used to evaluate the genetic backgrounds and relationships among the isolates. RESULTS: The 44 B. cereus group isolates, including soil- and animal-derived isolates and isolates recovered from hospitalized patients and food poisoning cases, were genotyped by MLST and core-genome SNP typing. Genetic variation among the isolates was identified by the MLST and core-genome SNP phylogeny comparison against reference strains from countries outside of Japan. Exploratory principal component analysis and nonmetric multidimensional scaling (NMDS) analyses were used to assess the genetic similarities among the isolates using gene presence and absence information and isolate origins as the metadata. A significant correlation was seen between the principal components and the presence of genes encoding hemolysin BL and emetic genetic determinants in B. cereus, and the capsule proteins in B. anthracis. NMDS showed that the cluster of soil isolates overlapped with the cluster comprising animal-derived and clinical isolates. CONCLUSIONS: Molecular and epidemiological analyses of B. cereus group isolates in Japan suggest that the soil- and animal-derived bacteria from our study are not a significant risk to human health. However, because several of the clinical isolates share close genetic relationships with the environmental isolates, both molecular and epidemiological surveillance studies could be used effectively to estimate virulence in these important pathogens.

New Size-Expanded Fluorescent Thymine Analogue: Synthesis, Characterization, and Application.

Here, the synthesis, photophysical characterization, and application of a new size-expanded thymine nucleoside, diox T, is described. diox T has desirable qualities as a T surrogate, including excellent quantum yield (0.36) and high environmental sensitivity. When incorporated into single- and double-stranded DNA, diox T showed excellent photophysical characteristics including a high quantum yield (average 0.20), and unlike BgQ, demonstrated dependence on neighboring bases without significant destabilization of the duplex. Interestingly, the matched base pair of adenine (A) and diox T has the unique property that it exhibits higher fluorescence than mismatched base pairs, and diox T has self-quenching effects. As one example of the possible applications of these promising features, single nucleoside polymorphism typing is demonstrated for discrimination of A by using diox T. The results suggest that diox T can be used for a broad range of applications in chemical biology.

Single excitation three color folded DNA probe for SNP typing.

A single fluorophore based multi-color emission probing system has been developed in the current work. dA(py) was used as the fluorophore, which was incorporated at the end of folded oligonucleotide with increasing number. Depending on the number of dA(py) three different color emissions were observed, viz. blue, green and red. These interesting color patterns were further used for the SNP typing of DNA. Every folded oligonucleotide with different number of dA(py)s could discriminate the full matched target sequence from one base mismatched and two base mismatched sequences. Furthermore, this probing system has a simple structure and also exhibited large Stokes shifted signal for discriminating the target DNA.

Deciphering the evolution of the temporal and geographic distribution of French Mycobacterium bovis genotypes using a high throughput SNP-targeted amplicon sequencing method.

Mycobacterium bovis, which belongs to the Mycobacterium tuberculosis complex, is a highly clonal pathogen. However, several lineages of M. bovis have been described worldwide and nine different clusters were identified in France. Targeted amplicon sequencing using next-generation sequencing technology of eighty-eight phylogenetically informative single nucleotide polymorphisms (SNPs) were used to infer the phylogenetic relationship of 630 strains of the National Reference Laboratory isolated between 1979 and 2018 from various animal species. This study allowed classifying 618 different genotypic profiles (combination of a spoligotype and 8 loci-MIRU-VNTR profiles) into the nine previously identified clusters. A global analysis of the entire collection of the National Reference Laboratory has made it possible to represent the evolution of clonal complexes and clusters in time and space for better assessing epidemiological changes of bovine tuberculosis in France.

The complex genomic diversity of Yersinia pestis on the long-term plague foci in Qinghai-Tibet plateau.

Plague is a typical natural focus disease that circulates in different ecology of vectors and reservoir hosts. We conducted genomic population and phylogenetic analyses of the Yersinia pestis collected from the 12 natural plague foci in China with more than 20 kinds of hosts and vectors. Different ecological landscapes with specific hosts, vectors, and habitat which shape various niches for Y. pestis. The phylogeographic diversity of Y. pestis in different kinds plague foci in China showed host niches adaptation. Most natural plague foci strains are region-and focus-specific, with one predominant subpopulation; but the isolates from the Qinghai-Tibet plateau harbor a higher genetic diversity than other foci. The Y. pestis from Marmota himalayana plague foci are defined as the ancestors of different populations at the root of the evolutionary tree, suggesting several different evolutionary paths to other foci. It has the largest pan-genome and widest SNP distances with most accessory genes enriched in mobilome functions (prophages, transposons). Geological barriers play an important role in the maintenance of local Y. pestis species and block the introduction of non-native strains. This study provides new insights into the control of plague outbreaks and epidemics, deepened the understanding of the evolutionary history of MHPF (M. himalayana plague focus) in China. The population structure and identify clades among different natural foci of China renewed the space cognition of the plague.

Molecular diversity of Mycobacterium tuberculosis isolates in Treatment-experienced Patients from Andhra Pradesh, India.

INTRODUCTION: To get a comprehensive idea about the transmission and epidemiology of TB globally and locally, the use of molecular typing methods has become imperative not only for understanding genetic diversity but also the population structure of Mycobacterium tuberculosis complex (MTBC). We aimed to investigate the drug resistance pattern and genetic diversity of MTBC among previously treated patients with sputum smear-positive pulmonary tuberculosis in a South Indian population. METHODOLOGY: 104 patients with sputum smear positivity and who had previously undergone treatment were selected. Drug susceptibility testing, Spoligotyping, MIRU-VNTR, and SNP typing were performed. RESULTS: Mono-resistance to isoniazid 16 (15.38%) was the highest among all drugs. Out of 104 isolates, 24 (23%) isolates were classified as MDR strains. The distributions of most common lineages were: EAI3-Ind-20 (19.23%), EAI5-13 (12.50%), Beijing-12 (11.54%), CAS1-Delhi- 9 (8.65%), and 7 (6.73%) each of T-H37rv, Unknown and Orphan types. MIRU-VNTR-based analysis revealed that there are two major groups: CAS1-Delhi and Beijing groups. Out of 104 isolates, 82 belonged to well-defined lineages and 6 clusters, and the remaining 22 were singletons. SNP analysis showed no mutations associated with five sets of genes in 33 strains. CONCLUSIONS: The occurrence of 11.54% Beijing strains in South India is an important finding. High frequency of Isoniazid mono resistance noticed. Spoligotyping along with MIRU-VNTR and SNP typing is the best approach to the identification of strain lineages. No mutation with Antigen85C gene represents, can be used for vaccine candidates.

Stable prevalence of Coxiella burnetii in wildlife after a decade of surveillance in northern Spain.

Coxiella burnetii is an obligate intracellular zoonotic bacterium widespread in nature that causes Q fever in animals and humans. The most common sources of human infection are domestic ruminants, but wildlife can also act as reservoir. Here, spleen samples from 652 wild ungulates and 218 wild birds collected in 2011-2019 in the Basque Country (northern Spain) were analysed by real-time PCR (IS1111 gene) and the results compared with data from a past study in 2001-2006. Among wild ungulates, C. burnetii DNA was detected in 7.0% (6/86) of roe deer (Capreolus capreolus), 1.9% (9/484) of wild boar (Sus scrofa) and 2.4% (2/82) of red deer (Cervus elaphus). The prevalence in roe deer was significantly higher compared to wild boar (p = 0.006). Among wild birds, only one white stork (Ciconia ciconia) tested positive. SNP-typing of C. burnetii-positive samples showed that wild ungulates shared SNP 2, SNP 6 and SNP 8 genotypes with domestic ruminants of the region. However, the white stork harboured a C. burnetii genotype (SNP 3) never identified in the studied area before. Comparing these results with those obtained in the same area a decade before (2001-2006), no significant differences were observed in the prevalence of C. burnetii in any of the wildlife species, indicating stability in C. burnetii prevalence. Nevertheless, continuous surveillance is needed to monitor any future changes in the reservoir role of roe deer and wild boar considering the increase in density of both species observed in Europe in the last decades.

Evolutionary and genomic insights into the long-term colonization of Shigella flexneri in animals.

The enteroinvasive bacterium Shigella flexneri is known as a highly host-adapted human pathogen. There had been no known other reservoirs reported until recently. Here 34 isolates obtained from animals (yaks, dairy cows and beef cattle) from 2016 to 2017 and 268 human S. flexneri isolates from China were sequenced to determine the relationships between animal and human isolates and infer the evolutionary history of animal-associated S. flexneri. The 18 animal isolates (15 yak and 3 beef cattle isolates) in PG1 were separated into 4 lineages, and the 16 animal isolates (1 yak, 5 beef cattle and 10 dairy cow isolates) in PG3 were clustered in 8 lineages. The most recent human isolates from China belonged to PG3 whereas Chinese isolates from the 1950s-1960s belonged to PG1. PG1 S. flexneri may has been transmitted to the yaks during PG1 circulation in the human population in China and has remained in the yak population since, while PG3 S. flexneri in animals were likely recent transmissions from the human population. Increased stability of the large virulence plasmid and acquisition of abundant antimicrobial resistance determinants may have enabled PG3 to expand globally and replaced PG1 in China. Our study confirms that animals may act as a reservoir for S. flexneri. Genomic analysis revealed the evolutionary history of multiple S. flexneri lineages in animals and humans in China. However, further studies are required to determine the public health threat of S. flexneri from animals.

Sequencing-Based Genotyping of Pakistani Burkholderia mallei Strains: A Useful Way for Investigating Glanders Outbreaks.

Burkholderia (B.) mallei is a host-adapted equine pathogen that causes glanders, a re-emerging zoonotic disease, which is endemic in Pakistan and other developing countries and seriously impacts the global equine movement. Due to globalization, the geographical restriction of diseases vanishes and the lack of awareness of and experience with eradicated diseases in industrialized countries also promotes the re-introduction of infections in these regions. Owing to the high equine population, the Pakistani province Punjab is a potential hotspot where several glanders outbreaks have been seen over last two decades. For determining the genomic diversity of B. mallei in this and other equine-populated prefectures, the genomes of 19 B. mallei strains isolated between 1999 and 2020 in different locations were sequenced and their genotypes were determined. Particularly, for genetically highly homogenous pathogens like B. mallei genotyping techniques require a high discriminatory power for enabling differentiation on the strain level. Thus, core-genome single nucleotide polymorphism (cgSNP) analysis was applied for distinguishing the highly similar strains. Furthermore, a whole-genome sequence-based core genome multi locus sequence typing (cgMLST) scheme, specific to B. mallei, was developed and additionally applied to the data. It was found that B. mallei genotypes in Pakistan persisted over time and space and genotype clusters preferred connection with a time point rather than the place of isolation, probably due to frequent equine movement, which promotes the spread of glanders. The cgMLST approach proved to work in accord with SNP typing and may help to investigate future glanders outbreaks.

Genome Sequence Analysis and Characterization of Shiga Toxin 2 Production by Escherichia coli O157:H7 Strains Associated With a Laboratory Infection.

A laboratory-acquired E. coli O157:H7 infection with associated severe sequelae including hemolytic uremic syndrome occurred in an individual working in the laboratory with a mixture of nalidixic acid-resistant (NalR) O157:H7 mutant strains in a soil-biochar blend. The patient was hospitalized and treated with an intravenous combination of metronidazole and levofloxacin. The present study investigated the source of this severe laboratory acquired infection and further examined the influence of the antibiotics used during treatment on the expression and production of Shiga toxin. Genomes of two Stx2a-and eae-positive O157:H7 strains isolated from the patient's stool were sequenced along with two pairs of the wt strains and their derived NalR mutants used in the laboratory experiments. High-resolution SNP typing determined the strains' individual genetic relatedness and unambiguously identified the two laboratory-derived NalR mutant strains as the source of the researcher's life-threatening disease, rather than a conceivable ingestion of unrelated O157:H7 isolates circulating at the same time. It was further confirmed that in sublethal doses, the antibiotics increased toxin expression and production. Our results support a simultaneous co-infection with clinical strains in the laboratory, which were the causative agents of previous O157:H7 outbreaks, and further that the administration of antibiotics may have impacted the outcome of the infection.

Evidence of Extensive Circulation of Yersinia enterocolitica in Rodents and Shrews in Natural Habitats from Retrospective and Perspective Studies in South Caucasus.

Yersinia enterocolitica culture-positive rodents and shrews were reported in different territories across Georgia during 14 of 17 years of investigations conducted for the period of 1981-1997. In total, Y. enterocolitica was isolated from 2052 rodents (15 species) and 33 shrews. Most isolates were obtained from Microtus arvalis, Rattus norvegicus, Mus musculus, and Apodemus spp. During the prospective study (2017-2019), isolates of Yersinia-like bacteria were cultured from 53 rodents collected in four parts of Georgia. All the Yersinia-like isolates were confirmed as Y. enterocolitica based on the API 20E and the BD Phenix50 tests. Whole-genome (WG) sequencing of five rodents and one shrew strain of Y. enterocolitica revealed that they possessed a set of virulence genes characteristic of the potentially pathogenic strains of biogroup 1A. All isolates lacked distinguished virulence determinants for YstA, Ail, TccC, VirF, and virulence plasmid pYV but carried the genes for YstB, YmoA, HemPR-HmuVSTU, YaxAB, PhlA, PldA, ArsCBR, and a flagellar apparatus. One strain contained a gene highly homologous to heat-labile enterotoxin, a chain of E. coli, a function not previously described for Y. enterocolitica. The WG single-nucleotide polymorphism-based typing placed the isolates in four distinct phylogenetic clusters.

Genomic Characterization of Strains From a Cluster of Infant Botulism Type A in a Small Town in Colorado, United States.

Three cases of infant botulism were reported in a small Colorado town between 1981 and 1984. The first two cases occurred in 1981, 6 months apart, and the third case occurred in 1984. Clostridium botulinum type A was isolated from stool of all three case patients and from environmental samples of the patient's homes. An epidemiological investigation and follow-up study were conducted from 1981 to 1986 and concluded the cases were likely related. In this study, we sought to determine whether the C. botulinum type A clinical isolates were related to each other and to isolates obtained from environmental samples. We performed whole genome sequencing (WGS) for 17 isolates associated with this potential cluster of infant botulism. Fifteen isolates were confirmed to be C. botulinum type A(B) and contained botulinum toxin gene subtypes A1 and B5 by WGS; these strains formed a monophyletic cluster in a phylogeny and were considered closely related to each other (0-18 high-quality single-nucleotide polymorphisms), but distinct from other C. botulinum type A(B) in Colorado and elsewhere in the United States. Results of our study suggest that the three infant botulism cases could have represented a cluster due to a C. botulinum type A(B) strain present in the environment.

Evaluation of WGS-subtyping methods for epidemiological surveillance of foodborne salmonellosis.

BACKGROUND: Salmonellosis is one of the most common foodborne diseases worldwide. Although human infection by non-typhoidal Salmonella (NTS) enterica subspecies enterica is associated primarily with a self-limiting diarrhoeal illness, invasive bacterial infections (such as septicaemia, bacteraemia and meningitis) were also reported. Human outbreaks of NTS were reported in several countries all over the world including developing as well as high-income countries. Conventional laboratory methods such as pulsed field gel electrophoresis (PFGE) do not display adequate discrimination and have their limitations in epidemiological surveillance. It is therefore very crucial to use accurate, reliable and highly discriminative subtyping methods for epidemiological characterisation and outbreak investigation. METHODS: Here, we used different whole genome sequence (WGS)-based subtyping methods for retrospective investigation of two different outbreaks of Salmonella Typhimurium and Salmonella Dublin that occurred in 2013 in UK and Ireland respectively. RESULTS: Single nucleotide polymorphism (SNP)-based cluster analysis of Salmonella Typhimurium genomes revealed well supported clades, that were concordant with epidemiologically defined outbreak and confirmed the source of outbreak is due to consumption of contaminated mayonnaise. SNP-analyses of Salmonella Dublin genomes confirmed the outbreak however the source of infection could not be determined. The core genome multilocus sequence typing (cgMLST) was discriminatory and separated the outbreak strains of Salmonella Dublin from the non-outbreak strains that were concordant with the epidemiological data however cgMLST could neither discriminate between the outbreak and non-outbreak strains of Salmonella Typhimurium nor confirm that contaminated mayonnaise is the source of infection, On the other hand, other WGS-based subtyping methods including multilocus sequence typing (MLST), ribosomal MLST (rMLST), whole genome MLST (wgMLST), clustered regularly interspaced short palindromic repeats (CRISPRs), prophage sequence profiling, antibiotic resistance profile and plasmid typing methods were less discriminatory and could not confirm the source of the outbreak. CONCLUSIONS: Foodborne salmonellosis is an important concern for public health therefore, it is crucial to use accurate, reliable and highly discriminative subtyping methods for epidemiological surveillance and outbreak investigation. In this study, we showed that SNP-based analyses do not only have the ability to confirm the occurrence of the outbreak but also to provide definitive evidence of the source of the outbreak in real-time.

The Transformation of Reference Microbiology Methods and Surveillance for Salmonella With the Use of Whole Genome Sequencing in England and Wales.

The use of whole genome sequencing (WGS) as a method for supporting outbreak investigations, studying Salmonella microbial populations and improving understanding of pathogenicity has been well-described (1-3). However, performing WGS on a discrete dataset does not pose the same challenges as implementing WGS as a routine, reference microbiology service for public health surveillance. Challenges include translating WGS data into a useable format for laboratory reporting, clinical case management, Salmonella surveillance, and outbreak investigation as well as meeting the requirement to communicate that information in an understandable and universal language for clinical and public health action. Public Health England have been routinely sequencing all referred presumptive Salmonella isolates since 2014 which has transformed our approach to reference microbiology and surveillance. Here we describe an overview of the integrated methods for cross-disciplinary working, describe the challenges and provide a perspective on how WGS has impacted the laboratory and surveillance processes in England and Wales.

Exploiting the Genomic Diversity of Rice (Oryza sativa L.): SNP-Typing in 11 Early-Backcross Introgression-Breeding Populations.

This study demonstrates genotyping-by-sequencing-based single-nucleotide polymorphism (SNP)-typing in 11 early-backcross introgression populations of rice (at BC1F5), comprising a set of 564 diverse introgression lines and 12 parents. Sequencing using 10 Ion Proton runs generated a total of ∼943.4 million raw reads, out of which ∼881.6 million reads remained after trimming for low-quality bases. After alignment, 794,297 polymorphic SNPs were identified, and filtering resulted in LMD50 SNPs (low missing data, with each SNP, genotyped in at least 50% of the samples) for each sub-population. Every data point was supported by actual sequencing data without any imputation, eliminating imputation-induced errors in SNP calling. Genotyping substantiated the impacts of novel breeding strategy revealing: (a) the donor introgression patterns in ILs were characteristic with variable introgression frequency in different genomic regions, attributed mainly to stringent selection under abiotic stress and (b) considerably lower heterozygosity was observed in ILs. Functional annotation revealed 426 non-synonymous deleterious SNPs present in 102 loci with a range of 1-4 SNPs per locus and 120 novel SNPs. SNP-typing this diversity panel will further assist in the development of markers supporting genomic applications in molecular breeding programs.

Phylogeny and Classification of Yersinia pestis Through the Lens of Strains From the Plague Foci of Commonwealth of Independent States.

The established phylogeny of the etiological agent of plague, Yersinia pestis, is not perfect, as it does not take into account the strains from numerous natural foci of Commonwealth of Independent States (CIS). We have carried out PCR and SNP typing of 359 strains and whole genome sequencing of 51 strains from these plague foci and determined the phylogenetic diversity of the strains circulating here. They belong to 0.ANT3, 0.ANT5, 2.ANT3, 4.ANT branches of antique biovar, 2.MED0, 2.MED1 branches of medieval biovar and to 0.PE2, 0.PE4a. 0.PE4h, 0.PE4t branches. Based on the studies of 178 strains from 23 plague foci of CIS countries, it was determined that the population structure of 2.MED strains is subdivided into Caucasian-Caspian and Central Asian-Chinese branches. In Central-Caucasian high-mountain plague foci in the Russian Federation (RF) the most deeply diverged branch of medieval biovar, 2.MED0, has been found. With the data obtained, the current population structure of Y. pestis species has been refined. New subspecies classification is developed, comprising seven subspecies: pestis, caucasica (0.PE2), angolica (0.PE3), central asiatica (0.PE4), tibetica (0.PE7), ulegeica (0.PE5), and qinghaica (0.PE10).