Strain-level fitness in the gut microbiome is an emergent property of glycans and a single metabolite.

The human gut microbiota resides within a diverse chemical environment challenging our ability to understand the forces shaping this ecosystem. Here, we reveal that fitness of the Bacteroidales, the dominant order of bacteria in the human gut, is an emergent property of glycans and one specific metabolite, butyrate. Distinct sugars serve as strain-variable fitness switches activating context-dependent inhibitory functions of butyrate. Differential fitness effects of butyrate within the Bacteroides are mediated by species-level variation in Acyl-CoA thioesterase activity and nucleotide polymorphisms regulating an Acyl-CoA transferase. Using in vivo multi-omic profiles, we demonstrate Bacteroides fitness in the human gut is associated together, but not independently, with Acyl-CoA transferase expression and butyrate. Our data reveal that each strain of the Bacteroides exists within a unique fitness landscape based on the interaction of chemical components unpredictable by the effect of each part alone mediated by flexibility in the core genome.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes.

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

Dynamics of the Streptomyces chromosome: chance and necessity.

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. Remarkably, these bacteria possess a large linear chromosome that is genetically compartmentalized: core genes are grouped in the central part, while the ends are populated by poorly conserved genes including antibiotic biosynthetic gene clusters. The genome is highly unstable and exhibits distinct evolutionary rates along the chromosome. Recent chromosome conformation capture (3C) and comparative genomics studies have shed new light on the interplay between genome dynamics in space and time. Here, we review insights that illustrate how the balance between chance (random genome variations) and necessity (structural and functional constraints) may have led to the emergence of spatial structuring of the Streptomyces chromosome.

Liquorilactobacillus: A Context of the Evolutionary History and Metabolic Adaptation of a Bacterial Genus from Fermentation Liquid Environments.

In the present work, we carried out a comparative genomic analysis to trace the evolutionary trajectory of the bacterial species that make up the Liquorilactobacillus genus, from the identification of genes and speciation/adaptation mechanisms in their unique characteristics to the identification of the pattern grouping these species. We present phylogenetic relationships between Liquorilactobacillus and related taxa such as Bacillus, basal lactobacilli and Ligilactobacillus, highlighting evolutionary divergences and lifestyle transitions across different taxa. The species of this genus share a core genome of 1023 genes, distributed in all COGs, which made it possible to characterize it as Liquorilactobacillus sensu lato: few amino acid auxotrophy, low genes number for resistance to antibiotics and general and specific cellular reprogramming mechanisms for environmental responses. These species were divided into four clades, with diversity being enhanced mainly by the diversity of genes involved in sugar metabolism. Clade 1 presented lower (< 70%) average amino acid identity with the other clades, with exclusive or absent genes, and greater distance in the genome compared to clades 2, 3 and 4. The data pointed to an ancestor of clades 2, 3 and 4 as being the origin of the genus Ligilactobacillus, while the species of clade 1 being closer to the ancestral Bacillus. All these traits indicated that the species of clade 1 could be soon separated in a distinct genus.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Comparative genomics of quinolone-resistant Escherichia coli from broilers and humans in Norway.

BACKGROUND: The usage of fluoroquinolones in Norwegian livestock production is very low, including in broiler production. Historically, quinolone-resistant Escherichia coli (QREC) isolated from Norwegian production animals rarely occur. However, with the introduction of a selective screening method for QREC in the Norwegian monitoring programme for antimicrobial resistance in the veterinary sector in 2014; 89.5% of broiler caecal samples and 70.7% of broiler meat samples were positive. This triggered the concern if there could be possible links between broiler and human reservoirs of QREC. We are addressing this by characterizing genomes of QREC from humans (healthy carriers and patients) and broiler isolates (meat and caecum). RESULTS: The most frequent mechanism for quinolone resistance in both broiler and human E. coli isolates were mutations in the chromosomally located gyrA and parC genes, although plasmid mediated quinolone resistance (PMQR) was also identified. There was some relatedness of the isolates within human and broiler groups, but little between these two groups. Further, some overlap was seen for isolates with the same sequence type isolated from broiler and humans, but overall, the SNP distance was high. CONCLUSION: Based on data from this study, QREC from broiler makes a limited contribution to the incidence of QREC in humans in Norway.

Genome mining, antimicrobial and plant growth-promoting potentials of halotolerant Bacillus paralicheniformis ES-1 isolated from salt mine.

Salinity severely affects crop yield by hindering nitrogen uptake and reducing plant growth. Plant growth-promoting bacteria (PGPB) are capable of providing cross-protection against biotic/abiotic stresses and facilitating plant growth. Genome-level knowledge of PGPB is necessary to translate the knowledge into a product as efficient biofertilizers and biocontrol agents. The current study aimed to isolate and characterize indigenous plant growth-promoting strains with the potential to promote plant growth under various stress conditions. In this regard, 72 bacterial strains were isolated from various saline-sodic soil/lakes; 19 exhibited multiple in vitro plant growth-promoting traits, including indole 3 acetic acid production, phosphate solubilization, siderophore synthesis, lytic enzymes production, biofilm formation, and antibacterial activities. To get an in-depth insight into genome composition and diversity, whole-genome sequence and genome mining of one promising Bacillus paralicheniformis strain ES-1 were performed. The strain ES-1 genome carries 12 biosynthetic gene clusters, at least six genomic islands, and four prophage regions. Genome mining identified plant growth-promoting conferring genes such as phosphate solubilization, nitrogen fixation, tryptophan production, siderophore, acetoin, butanediol, chitinase, hydrogen sulfate synthesis, chemotaxis, and motility. Comparative genome analysis indicates the region of genome plasticity which shapes the structure and function of B. paralicheniformis and plays a crucial role in habitat adaptation. The strain ES-1 has a relatively large accessory genome of 649 genes (~ 19%) and 180 unique genes. Overall, these results provide valuable insight into the bioactivity and genomic insight into B. paralicheniformis strain ES-1 with its potential use in sustainable agriculture.

A core-genome multilocus sequence typing scheme for the detection of genetically related Streptococcus pyogenes clusters.

The recently observed increase in invasive Streptococcus pyogenes infections causes concern in Europe. However, conventional molecular typing methods lack discriminatory power to aid investigations of outbreaks caused by S. pyogenes. Therefore, there is an urgent need for high-resolution molecular typing methods to assess genetic relatedness between S. pyogenes isolates. In the current study, we aimed to develop a novel high-resolution core-genome multilocus sequence typing (cgMLST) scheme for S. pyogenes and compared its discriminatory power to conventional molecular typing methods. The cgMLST scheme was designed with the commercial Ridom SeqSphere+ software package. To define a cluster threshold, the scheme was evaluated using publicly available data from nine defined S. pyogenes outbreaks in the United Kingdom. The cgMLST scheme was then applied to 23 isolates from a suspected S. pyogenes outbreak and 117 S. pyogenes surveillance isolates both from the Netherlands. MLST and emm-typing results were used for comparison to cgMLST results. The allelic differences between isolates from defined outbreaks ranged between 6 and 31 for isolates with the same emm-type, resulting in a proposed cluster threshold of <5 allelic differences out of 1,095 target loci. Seven out of twenty-three (30%) isolates from the suspected outbreak had an allelic difference of <2, thereby identifying a potential cluster that could not be linked to other isolates. The proposed cgMLST scheme shows a higher discriminatory ability when compared to conventional typing methods. The rapid and simple analysis workflow allows for extended detection of clusters of potential outbreak isolates and surveillance and may facilitate the sharing of sequencing results between (inter)national laboratories.

Multi-host infection and phylogenetically diverse lineages shape the recombination and gene pool dynamics of Staphylococcus aureus.

BACKGROUND: Staphylococcus aureus can infect and adapt to multiple host species. However, our understanding of the genetic and evolutionary drivers of its generalist lifestyle remains inadequate. This is particularly important when considering local populations of S. aureus, where close physical proximity between bacterial lineages and between host species may facilitate frequent and repeated interactions between them. Here, we aim to elucidate the genomic differences between human- and animal-derived S. aureus from 437 isolates sampled from disease cases in the northeast region of the United States. RESULTS: Multi-locus sequence typing revealed the existence of 75 previously recognized sequence types (ST). Our population genomic analyses revealed heterogeneity in the accessory genome content of three dominant S. aureus lineages (ST5, ST8, ST30). Genes related to antimicrobial resistance, virulence, and plasmid types were differentially distributed among isolates according to host (human versus non-human) and among the three major STs. Across the entire population, we identified a total of 1,912 recombination events that occurred in 765 genes. The frequency and impact of homologous recombination were comparable between human- and animal-derived isolates. Low-frequency STs were major donors of recombined DNA, regardless of the identity of their host. The most frequently recombined genes (clfB, aroA, sraP) function in host infection and virulence, which were also frequently shared between the rare lineages. CONCLUSIONS: Taken together, these results show that frequent but variable patterns of recombination among co-circulating S. aureus lineages, including the low-frequency lineages, that traverse host barriers shape the structure of local gene pool and the reservoir of host-associated genetic variants. Our study provides important insights to the genetic and evolutionary factors that contribute to the ability of S. aureus to colonize and cause disease in multiple host species. Our study highlights the importance of continuous surveillance of S. aureus circulating in different ecological host niches and the need to systematically sample from them. These findings will inform development of effective measures to control S. aureus colonization, infection, and transmission across the One Health continuum.

Proposal for Acinetobacter higginsii sp. nov. to accommodate organisms of human clinical origin previously classified as Acinetobacter genomic species 16.

In 1989, Bouvet and Jeanjean delineated five proteolytic genomic species (GS) of Acinetobacter, each with two to four human isolates. Three were later validly named, whereas the remaining two (GS15 and GS16) have been awaiting nomenclatural clarification. Here we present the results of the genus-wide taxonomic study of 13 human strains classified as GS16 (n=10) or GS15 (n=3). Based on core genome phylogenetic analysis, the strains formed two respective but closely related phylogroups within the Acinetobacter haemolytic clade. The intraspecies genomic average nucleotide identity based on blast (ANIb) values for GS16 and GS15 reached ≥94.9 % and ≥98.7, respectively, whereas ANIb values between them were 92.5-93.5% and those between them and the known species were ≤91.5 %. GS16 and GS15 could be differentiated from the other Acinetobacter species by their ability to lyse gelatin and sheep blood and to assimilate d,l-lactate, along with their inability to acidify d-glucose and assimilate glutarate. In contrast, GS16 and GS15 were indistinguishable from one another by metabolic/physiological features or whole-cell MALDI-TOF mass spectra. All the GS15/GS16 genomes contained genes encoding a class D ß-lactamase, Acinetobacter-derived cephalosporinase and aminoglycoside 6'-N-acetyltransferase. Searching NCBI databases revealed genome sequences of three additional isolates of GS16, but none of GS15. We conclude that our data support GS16 as representing a novel species, but leave the question of the taxonomic status of GS15 open, given its close relatedness to GS16 and the small number of available strains. We propose the name Acinetobacter higginsii sp. nov. for GS16, with the type strain NIPH 1872T (CCM 9243T=CIP 70.18T=ATCC 17988T).

Sifting through the core-genome to identify putative cross-protective antigens against Riemerella anatipestifer.

Infectious serositis of ducks, caused by Riemerella anatipestifer, is one of the main infectious diseases that harm commercial ducks. Whole-strain-based vaccines with no or few cross-protection were observed between different serotypes of R. anatipestifer, and so far, control of infection is hampered by a lack of effective vaccines, especially subunit vaccines with cross-protection. Since the concept of reverse vaccinology was introduced, it has been widely used to screen for protective antigens in important pathogens. In this study, pan-genome binding reverse vaccinology, an emerging approach to vaccine candidate screening, was used to screen for cross-protective antigens against R. anatipestifer. Thirty proteins were identified from the core-genome as potential cross-protective antigens. Three of these proteins were recombinantly expressed, and their immunoreactivity with five antisera (anti-serotypes 1, 2, 6, 10, and 11) was demonstrated by Western blotting. Our study established a method for high-throughput screening of cross-protective antigens against R. anatipestifer in silico, which will lay the foundation for the development of a cross-protective subunit vaccine controlling R. anatipestifer infection. KEY POINTS: • Pan-genome binding reverse vaccine approach was first established in R. anatipestifer to screen for subunit vaccine candidates. • Thirty potential cross-protective antigens against R. anatipestifer were identified by this method. • The reliability of the method was verified preliminarily by the results of Western blotting of three of these potential antigens.

Comparative pangenomics: analysis of 12 microbial pathogen pangenomes reveals conserved global structures of genetic and functional diversity.

BACKGROUND: With the exponential growth of publicly available genome sequences, pangenome analyses have provided increasingly complete pictures of genetic diversity for many microbial species. However, relatively few studies have scaled beyond single pangenomes to compare global genetic diversity both within and across different species. We present here several methods for "comparative pangenomics" that can be used to contextualize multi-pangenome scale genetic diversity with gene function for multiple species at multiple resolutions: pangenome shape, genes, sequence variants, and positions within variants. RESULTS: Applied to 12,676 genomes across 12 microbial pathogenic species, we observed several shared resolution-specific patterns of genetic diversity: First, pangenome openness is associated with species' phylogenetic placement. Second, relationships between gene function and frequency are conserved across species, with core genomes enriched for metabolic and ribosomal genes and accessory genomes for trafficking, secretion, and defense-associated genes. Third, genes in core genomes with the highest sequence diversity are functionally diverse. Finally, certain protein domains are consistently mutation enriched across multiple species, especially among aminoacyl-tRNA synthetases where the extent of a domain's mutation enrichment is strongly function-dependent. CONCLUSIONS: These results illustrate the value of each resolution at uncovering distinct aspects in the relationship between genetic and functional diversity across multiple species. With the continued growth of the number of sequenced genomes, these methods will reveal additional universal patterns of genetic diversity at the pangenome scale.

Salmonella enterica subsp. enterica Welikade: guideline for phylogenetic analysis of serovars rarely involved in foodborne outbreaks.

BACKGROUND: Salmonella spp. is a major foodborne pathogen with a wide variety of serovars associated with human cases and food sources. Nevertheless, in Europe a panel of ten serovars is responsible for up to 80% of confirmed human cases. Clustering studies by single nucleotide polymorphism (SNP) core-genome phylogenetic analysis of outbreaks due to these major serovars are simplified by the availability of many complete genomes in the free access databases. This is not the case for outbreaks due to less common serovars, such as Welikade, for which no reference genomes are available. In this study, we propose a method to solve this problem. We propose to perform a core genome MLST (cgMLST) analysis based on hierarchical clustering using the free-access EnteroBase to select the most suitable genome to use as a reference for SNP phylogenetic analysis. In this study, we applied this protocol to a retrospective analysis of a Salmonella enterica serovar Welikade (S. Welikade) foodborne outbreak that occurred in France in 2016. Finally, we compared the cgMLST and SNP analyses. SNP phylogenetic reconstruction was carried out considering the effect of recombination events identified by the ClonalFrameML tool. The accessory genome was also explored by phage content and virulome analyses. RESULTS: Our findings revealed high clustering concordance using cgMLST and SNP analyses. Nevertheless, SNP analysis allowed for better assessment of the genetic distance among strains. The results revealed epidemic clones of S. Welikade circulating within the poultry and dairy sectors in France, responsible for sporadic and non-sporadic human cases between 2012 and 2019. CONCLUSIONS: This study increases knowledge on this poorly described serovar and enriches public genome databases with 42 genomes from human and non-human S. Welikade strains, including the isolate collected in 1956 in Sri Lanka, which gave the name to this serovar. This is the first genomic analysis of an outbreak due to S. Welikade described to date.

Association of Phylogenomic Relatedness among Neisseria gonorrhoeae Strains with Antimicrobial Resistance, Austria, 2016-2020.

We investigated genomic determinants of antimicrobial resistance in 1,318 Neisseria gonorrhoeae strains isolated in Austria during 2016-2020. Sequence type (ST) 9363 and ST11422 isolates had high rates of azithromycin resistance, and ST7363 isolates correlated with cephalosporin resistance. These results underline the benefit of genomic surveillance for antimicrobial resistance monitoring.

CoreCruncher: Fast and Robust Construction of Core Genomes in Large Prokaryotic Data Sets.

The core genome represents the set of genes shared by all, or nearly all, strains of a given population or species of prokaryotes. Inferring the core genome is integral to many genomic analyses, however, most methods rely on the comparison of all the pairs of genomes; a step that is becoming increasingly difficult given the massive accumulation of genomic data. Here, we present CoreCruncher; a program that robustly and rapidly constructs core genomes across hundreds or thousands of genomes. CoreCruncher does not compute all pairwise genome comparisons and uses a heuristic based on the distributions of identity scores to classify sequences as orthologs or paralogs/xenologs. Although it is much faster than current methods, our results indicate that our approach is more conservative than other tools and less sensitive to the presence of paralogs and xenologs. CoreCruncher is freely available from: https://github.com/lbobay/CoreCruncher. CoreCruncher is written in Python 3.7 and can also run on Python 2.7 without modification. It requires the python library Numpy and either Usearch or Blast. Certain options require the programs muscle or mafft.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes.

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Comparison of Whole-Genome Sequence-Based Methods and PCR Ribotyping for Subtyping of Clostridioides difficile.

Clostridioides difficile is the most common cause of antibiotic-associated gastrointestinal infections. Capillary electrophoresis (CE)-PCR ribotyping is currently the gold standard for C. difficile typing but lacks the discriminatory power to study transmission and outbreaks in detail. New molecular methods have the capacity to differentiate better and provide standardized and interlaboratory exchangeable data. Using a well-characterized collection of diverse strains (N = 630; 100 unique ribotypes [RTs]), we compared the discriminatory power of core genome multilocus sequence typing (cgMLST) (SeqSphere and EnteroBase), whole-genome MLST (wgMLST) (EnteroBase), and single-nucleotide polymorphism (SNP) analysis. A unique cgMLST profile (more than six allele differences) was observed in 82 of 100 RTs, indicating that cgMLST could distinguish most, but not all, RTs. Application of cgMLST in two outbreak settings with RT078 and RT181 (known to have low intra-RT allele differences) showed no distinction between outbreak and nonoutbreak strains in contrast to wgMLST and SNP analysis. We conclude that cgMLST has the potential to be an alternative to CE-PCR ribotyping. The method is reproducible, easy to standardize, and offers higher discrimination. However, adjusted cutoff thresholds and epidemiological data are necessary to recognize outbreaks of some specific RTs. We propose to use an allelic threshold of three alleles to identify outbreaks.

Core Genome Multilocus Sequence Typing Scheme for Improved Characterization and Epidemiological Surveillance of Pathogenic Brucella.

Brucellosis poses a significant burden to human and animal health worldwide. Robust and harmonized molecular epidemiological approaches and population studies that include routine disease screening are needed to efficiently track the origin and spread of Brucella strains. Core genome multilocus sequence typing (cgMLST) is a powerful genotyping system commonly used to delineate pathogen transmission routes for disease surveillance and control. Except for Brucella melitensis, cgMLST schemes for Brucella species are currently not established. Here, we describe a novel cgMLST scheme that covers multiple Brucella species. We first determined the phylogenetic breadth of the genus using 612 Brucella genomes. We selected 1,764 genes that were particularly well conserved and typeable in at least 98% of these genomes. We tested the new scheme on 600 genomes and found high agreement with the whole-genome-based single nucleotide polymorphism (SNP) analysis. Next, we applied the scheme to reanalyze the genome of Brucella strains from epidemiologically linked outbreaks. We demonstrated the applicability of the new scheme for high-resolution typing required in outbreak investigations as previously reported with whole-genome SNP methods. We also used the novel scheme to define the global population structure of the genus using 1,322 Brucella genomes. Finally, we demonstrated the possibility of tracing distribution of Brucella strains by performing cluster analysis of cgMLST profiles and found nearly identical cgMLST profiles in different countries. Our results show that sequencing depth of more than 40-fold is optimal for allele calling with this scheme. In summary, this study describes a novel Brucella-wide cgMLST scheme that is applicable in Brucella molecular epidemiology and helps in accurately tracking and thus controlling the sources of infection. The scheme is publicly accessible and should represent a valuable resource for laboratories with limited computational resources and bioinformatics expertise.

Differentiating Closely Affiliated Dehalococcoides Lineages by a Novel Genetic Marker Identified via Computational Pangenome Analysis.

As a group, the genus Dehalococcoides dehalogenates a wide range of organohalide pollutants, but the range of organohalide compounds that can be utilized for reductive dehalogenation differs among Dehalococcoides strains. Dehalococcoides lineages cannot be reliably disambiguated in mixed communities using typical phylogenetic markers, which often confounds bioremediation efforts. Here, we describe a computational approach to identify Dehalococcoides genetic markers with improved discriminatory resolution. Screening core genes from the Dehalococcoides pangenome for degree of similarity and frequency of 100% identity found a candidate genetic marker encoding a bacterial neuraminidase repeat (BNR)-containing protein of unknown function. This gene exhibits the fewest completely identical amino acid sequences and has among the lowest average amino acid sequence identity in the core pangenome. Primers targeting BNR could effectively discriminate between 40 available BNR sequences (in silico) and 10 different Dehalococcoides isolates (in vitro). Amplicon sequencing of BNR fragments generated from 22 subsurface soil samples revealed a total of 109 amplicon sequence variants, suggesting a high diversity of Dehalococcoides distributed in the environment. Therefore, the BNR gene can serve as an alternative genetic marker to differentiate strains of Dehalococcoides in complicated microbial communities. IMPORTANCE The challenge of discriminating between phylogenetically similar but functionally distinct bacterial lineages is particularly relevant to the development of technologies seeking to exploit the metabolic or physiological characteristics of specific members of bacterial genera. A computational approach was developed to expedite screening of potential genetic markers among phylogenetically affiliated bacteria. Using this approach, a gene encoding a bacterial neuraminidase repeat (BNR)-containing protein of unknown function was selected and evaluated as a genetic marker to differentiate strains of Dehalococcoides, an environmentally relevant genus of bacteria whose members can transform and detoxify a range of halogenated organic solvents and persistent organic pollutants, in complex microbial communities to demonstrate the validity of the approach. Moreover, many apparently phylogenetically distinct, currently uncharacterized Dehalococcoides were detected in environmental samples derived from contaminated sites.

Genomic Characterization of Bacillus safensis Isolated from Mine Tailings in Peru and Evaluation of Its Cyanide-Degrading Enzyme CynD.

Understanding the biochemistry and metabolic pathways of cyanide degradation is necessary to improve the efficacy of cyanide bioremediation processes and industrial requirements. We have isolated and sequenced the genome of a cyanide-degrading Bacillus strain from water in contact with mine tailings from Lima, Peru. This strain was classified as Bacillus safensis based on 16S rRNA gene sequencing and core genome analyses and named B. safensis PER-URP-08. We searched for possible cyanide-degradation enzymes in the genome of this strain and identified a putative cyanide dihydratase (CynD) gene similar to a previously characterized CynD from Bacillus pumilus C1. Sequence analysis of CynD from B. safensis and B. pumilus allow us to identify C-terminal residues that differentiate both CynDs. We then cloned, expressed in Escherichia coli, and purified recombinant CynD from B. safensis PER-URP-08 (CynDPER-URP-08) and showed that in contrast to CynD from B. pumilus C1, this recombinant CynD remains active at up to pH 9. We also showed that oligomerization of CynDPER-URP-08 decreases as a function of increased pH. Finally, we demonstrated that transcripts of CynDPER-URP-08 in B. safensis PER-URP-08 are strongly induced in the presence of cyanide. Our results suggest that the use of B. safensis PER-URP-08 and CynDPER-URP-08 as potential tool for cyanide bioremediation warrants further investigation. IMPORTANCE Despite being of environmental concern around the world due to its toxicity, cyanide continues to be used in many important industrial processes. Thus, searching for cyanide bioremediation methods is a matter of societal concern and must be present on the political agenda of all governments. Here, we report the isolation, genome sequencing and characterization of cyanide degradation capacity of a bacterial strain isolated from an industrial mining site in Peru. We characterize a cyanide dehydratase (CynD) homolog from one of these bacteria, Bacillus safensis PER-URP-08.