Comparative genomics unveils extensive genomic variation between populations of Listeria species in natural and food-associated environments

Comprehending bacterial genomic variation linked to distinct environments can yield novel insights into mechanisms underlying differential adaptation and transmission of microbes across environments. Gaining such insights is particularly crucial for pathogens as it benefits public health surveillance. However, the understanding of bacterial genomic variation is limited by a scarcity of investigations in genomic variation coupled with different ecological contexts. To address this limitation, we focused on Listeria, an important bacterial genus for food safety that includes the human pathogen L. monocytogenes, and analyzed a large-scale genomic dataset collected by us from natural and food-associated environments across the United States. Through comparative genomics analyses on 449 isolates from the soil and 390 isolates from agricultural water and produce processing facilities representing L. monocytogenes, L. seeligeri, L. innocua, and L. welshimeri, we find that the genomic profiles strongly differ by environments within each species. This is supported by the environment-associated subclades and differential presence of plasmids, stress islands, and accessory genes involved in cell envelope biogenesis and carbohydrate transport and metabolism. Core genomes of Listeria species are also strongly associated with environments and can accurately predict isolation sources at the lineage level in L. monocytogenes using machine learning. We find that the large environment-associated genomic variation in Listeria appears to be jointly driven by soil property, climate, land use, and accompanying bacterial species, chiefly representing Actinobacteria and Proteobacteria. Collectively, our data suggest that populations of Listeria species have genetically adapted to different environments, which may limit their transmission from natural to food-associated environments.