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

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