Purpose: This project employed both microarray analysis (MA) and whole genomic sequencing (WGS) in an effort to provide a deeper understanding of the phylogeny of EPs found among Cronobacterspp.
Methods: One hundred thirty-nine Cronobacter strains were analyzed using a previously described pan-genomic DNA microarray which contained 37 EP probe sets representing 13 different EPs. Using the microarray annotations for these genes, gene sequences for three of these EPs, namely, kefA (a potassium efflux system), a resistance/nodulation/division (RND) EP gene, and cmeB (a tranporter) were bioinformatically captured by searching a local database of previously sequenced Cronobactergenomes. Phylogenetic analysis of these sequences was conducted using the Maximum Composite Likelihood method in Genomics CLC and MEGA7.
Results: MA revealed that the prevalence of the 37 efflux pumps followed evolutionary species lines which were previously noted. WGS analysis of sequences for kefA, the RND EP, and cmeB, supported this phylogenetic relationship found by MA.
Significance: Understanding the genomic landscape of EPs among Cronobacter spp. is essential for future studies designed to identify mechanisms used by Cronobacter spp. to survive and persist under stressful growth conditions. The results reported here demonstrate that Cronobacter species have at least 37 efflux pumps and that all of the species possess species-specific orthologues; that these genes most probably arose from a common hypothetical ancestor; and their identification is critical for reliable prediction of gene function.