Genome Sequencing and Comparative Genomics of Indian Isolates of Brucella melitensis

Kumaragurubaran Karthik; Subbaiyan Anbazhagan; Prasad Thomas; Murugesan Ananda Chitra; Tuticorin Maragatham Alagesan Senthilkumar; Ramaswamy Sridhar; Gopal Dhinakar Raj

doi:10.3389/fmicb.2021.698069

Genome Sequencing and Comparative Genomics of Indian Isolates of Brucella melitensis

Front Microbiol. 2021 Aug 20:12:698069. doi: 10.3389/fmicb.2021.698069. eCollection 2021.

Authors

Kumaragurubaran Karthik¹, Subbaiyan Anbazhagan², Prasad Thomas², Murugesan Ananda Chitra¹, Tuticorin Maragatham Alagesan Senthilkumar³, Ramaswamy Sridhar¹, Gopal Dhinakar Raj⁴

Affiliations

¹ Central University Laboratory, Tamil Nadu Veterinary and Animal Sciences University, Chennai, India.
² Division of Bacteriology and Mycology, Indian Veterinary Research Institute, Bareilly, India.
³ Translational Research Platform for Veterinary Biologicals, Tamil Nadu Veterinary and Animal Sciences University, Chennai, India.
⁴ Centre for Animal Health Studies, Tamil Nadu Veterinary and Animal Sciences University, Chennai, India.

Abstract

Brucella melitensis causes small ruminant brucellosis and a zoonotic pathogen prevalent worldwide. Whole genome phylogeny of all available B. melitensis genomes (n = 355) revealed that all Indian isolates (n = 16) clustered in the East Mediterranean lineage except the ADMAS-GI strain. Pangenome analysis indicated the presence of limited accessory genomes with few clades showing specific gene presence/absence pattern. A total of 43 virulence genes were predicted in all the Indian strains of B. melitensis except 2007BM-1 (ricA and wbkA are absent). Multilocus sequence typing (MLST) analysis indicated all except one Indian strain (ADMAS-GI) falling into sequence type (ST 8). In comparison with MLST, core genome phylogeny indicated two major clusters (>70% bootstrap support values) among Indian strains. Clusters with <70% bootstrap support values represent strains with diverse evolutionary origins present among animal and human hosts. Genetic relatedness among animal (sheep and goats) and human strains with 100% bootstrap values shows its zoonotic transfer potentiality. SNP-based analysis indicated similar clustering to that of core genome phylogeny. Among the Indian strains, the highest number of unique SNPs (112 SNPs) were shared by a node that involved three strains from Tamil Nadu. The node SNPs involved several peptidase genes like U32, M16 inactive domain protein, clp protease family protein, and M23 family protein and mostly represented non-synonymous (NS) substitutions. Vaccination has been followed in several parts of the world to prevent small ruminant brucellosis but not in India. Comparison of Indian strains with vaccine strains showed that M5 is genetically closer to most of the Indian strains than Rev.1 strain. The presence of most of the virulence genes among all Indian strains and conserved core genome compositions suggest the use of any circulating strain/genotypes for the development of a vaccine candidate for small ruminant brucellosis in India.

Keywords: Brucella melitensis; India; comparative genomics; pangenome; vaccine candidate; virulence factors.