# Dichloromethane Degradation Pathway from Unsequenced Hyphomicrobium sp. MC8b Rapidly Explored by Pan-Proteomics

Karim Hayoun, Emilie Geersens, Cédric C Laczny, Rashi Halder, Carmen Lázaro Sánchez, Abhijit Manna, Françoise Bringel, Michaël Ryckelynck, Paul Wilmes, Emilie Muller, et al.

► To cite this version:

Karim Hayoun, Emilie Geersens, Cédric C Laczny, Rashi Halder, Carmen Lázaro Sánchez, et al. Dichloromethane Degradation Pathway from Unsequenced Hyphomicrobium sp. MC8b Rapidly Explored by Pan-Proteomics. Microorganisms, 2020, 8 (12), 10.3390/microorganisms8121876. hal-03107382

HAL Id: hal-03107382

https://hal.science/hal-03107382

Submitted on 9 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Dichloromethane Degradation Pathway from Unsequenced *Hyphomicrobium* sp. MC8b Rapidly Explored by Pan-Proteomics

Karim Hayoun ¹ID, Emilie Geersens <fnref n="2" />,<fnref n="3" />, Cédric C. Laczny ⁴ID, Rashi Halder <fnref n="4" />,
Carmen Lázaro Sánchez <fnref n="2" />, Abhijit Manna <fnref n="2" />, Françoise Bringel ²ID, Michaël Ryckelynck ³ID,
Paul Wilmes ⁴ID, Emilie E. L. Muller ²ID, Béatrice Alpha-Bazin ¹ID, Jean Armengaud <fnref n="1" />,*,ID
and Stéphane Vuilleumier <fnref n="2" />,*,ID

Received: 11 November 2020; Accepted: 24 November 2020; Published: 27 November 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Several bacteria are able to degrade the major industrial solvent dichloromethane (DCM) by using the conserved dehalogenase DcmA, the only system for DCM degradation characterised at the sequence level so far. Using differential proteomics, we rapidly identified key determinants of DCM degradation for *Hyphomicrobium* sp. MC8b, an unsequenced facultative methylotrophic DCM-degrading strain. For this, we designed a pan-proteomics database comprising the annotated genome sequences of 13 distinct *Hyphomicrobium* strains. Compared to growth with methanol, growth with DCM induces drastic changes in the proteome of strain MC8b. Dichloromethane dehalogenase DcmA was detected by differential pan-proteomics, but only with poor sequence coverage, suggesting atypical characteristics of the DCM dehalogenation system in this strain. More peptides were assigned to DcmA by error-tolerant search, warranting subsequent sequencing of the genome of strain MC8b, which revealed a highly divergent set of *dcm* genes in this strain. This suggests that the *dcm* enzymatic system is less strongly conserved than previously believed, and that substantial molecular evolution of *dcm* genes has occurred beyond their horizontal transfer in the bacterial domain. Our study showed the power of pan-proteomics for quick characterization of new strains belonging to branches of the Tree of Life that are densely genome-sequenced.

**Keywords:** pan-proteomics; differential proteomics; *Hyphomicrobium*; dichloromethane; dehalogenation; *dcm* genes; genome sequencing; Nanopore

## 1. Introduction

Continuing advances in high-throughput sequencing technologies are revealing increasingly large diversity in the microbial world [1]. Beyond its fundamental interest, this unsuspected diversity holds

<footnote n="1">¹ Laboratoire Innovations Technologiques pour la Détection et le Diagnostic (Li2D), Service de Pharmacologie et Immunoanalyse (SPI), CEA, INRA, F-30207 Bagnols-sur-Cèze, France; Karim.HAYOUN@cea.fr (K.H.); beatrice.alpha-bazin@cea.fr (B.A.-B.)</footnote>

<footnote n="2">² Génétique Moléculaire, Génomique, Microbiologie, UMR 7156 CNRS, Université de Strasbourg, F-67000 Strasbourg, France; emilie.geersens@etu.unistra.fr (E.G.); carmen.lazaro-sanchez@etu.unistra.fr (C.L.S.); amanna@unistra.fr (A.M.); francoise.bringel@unistra.fr (F.B.); emilie.muller@unistra.fr (E.E.L.M.)</footnote>

<footnote n="3">³ Institut de Biologie Moléculaire et Cellulaire, UPR 9002 CNRS, Université de Strasbourg, F-67084 Strasbourg CEDEX, France; m.ryckelynck@unistra.fr</footnote>

<footnote n="4">⁴ Luxembourg Centre for Systems Biomedicine, University of Luxembourg, L-4362 Esch-sur-Alzette, Luxembourg; cedric.laczny@uni.lu (C.C.L.); rashi.halder@uni.lu (R.H.); paul.wilmes@uni.lu (P.W.)</footnote>

<footnote>* Correspondence: jean.armengaud@cea.fr (J.A.); vuilleumier@unistra.fr (S.V.)</footnote>
great promise for the discovery of novel enzymes for application in biocatalysis and bioremediation [2]. However, because a large part of DNA sequences retrieved from our environment is of unknown function, information derived from sequence only provides preferential access to enzyme classes that are already known, rather than to the much larger diversity of enzymes performing a given function that likely exist in nature.

Thus, complementary approaches to DNA sequencing are needed to rapidly identify new ways of performing a function of interest. Proteomics is one of the high-throughput omics techniques available today to gain information. In particular, comparative proteomics under different conditions of interest for the investigated function has the potential to detect differentially synthesised proteins, providing hints to their functional association [3]. The proteome of closely-related bacterial strains can be compared unifying all the genome sequences into a single protein database for interpreting shotgun proteomic data [4]. This pan-proteomics concept allows better quantitative proteome measurements and was recently applied to characterize the proteome of *Lactococcus lactis* strains [5], *Streptococcus agalactiae* strains [6], and mid-exponentially grown *Brucella* representatives [7], thereby promoting information attributions about unified groups of proteins within a given branch of the Tree of Life. The concept can be further extended to still uncharacterised isolates.

One microbial function that has attracted much interest is dehalogenation, also because of the wide range of potential uses for dehalogenase enzymes [8,9]. Many halogenated compounds are produced and present naturally on earth [10]. Industrial production of some of these compounds has led to the expansion and further evolution of specialized bacteria able to degrade such compounds and to use them as nutrients for growth. A number of dehalogenases have been discovered and characterised at the gene, protein, and enzymatic levels. Nevertheless, known enzymes with this function belong to a relatively small number of protein families, and may represent only a fraction of the diversity of the extant dehalogenase repertoire on our planet.

Cultivable strains that grow with a given halogenated compound, but in which the corresponding dehalogenase gene cannot be detected, hold great interest with respect to the discovery of novel dehalogenases. Dehalogenation of dichloromethane (DCM) represents a long-standing and thoroughly investigated paradigm in the field [11]. DCM is a naturally produced volatile and water-soluble toxic compound, and one of the major halogenated solvents intensively used by the industry. Thus, bacteria capable of degrading dichloromethane were among the first dehalogenating organisms to be discovered as part of initial efforts to address the issue of environmental contamination by organohalogens. Following the characterisation of DCM dehalogenases and gene identification, it was established that the DCM dehalogenase *dcmA* gene was very strongly conserved in all investigated aerobic DCM-degrading bacteria. To date, DCM dehalogenase is still the only growth-supporting dehalogenase acting on DCM that has been characterized at the molecular level [11].

The DCM-degrading strain MC8b used in pioneering studies of isotope-based characterisation of biological dehalogenation [12,13] is of high interest. Both Southern-blot and PCR-based approaches failed to detect the *dcmA* gene in strain MC8b [14]. Analysis of the ribosomal 16S rRNA gene indicated its affiliation to the Alphaproteobacterial genus *Hyphomicrobium*, but no further molecular analysis was attempted.

In this work, we used strain MC8b to demonstrate the value of a pan-proteomics-based approach circumventing whole-genome DNA sequence analysis for rapid strain interrogation of its key function of interest. Using a custom database of the predicted proteomes of available genome-sequenced *Hyphomicrobium* strains, global differential proteomics of DCM and methanol cultures showed that strain MC8b features a DCM utilisation system that is highly divergent from previously characterised DCM-degrading strains. This was confirmed by sequencing the genome of strain MC8b, which was obtained as a high quality, assembled, and closed sequence.
## 2. Materials and Methods

### 2.1. *Strain Cultivation*

Strain MC8b [12] was restreaked on solid mineral minimal medium [15] from a frozen laboratory stock, with DCM as the sole carbon and energy source for growth. For liquid cultures, the strain was routinely grown aerobically in the same medium, in gas-tight Erlenmeyer flasks fitted with Mininert caps (Supelco, Sigma, St. Louis, MO, USA), at 30 °C with agitation at 120 rpm, with either 10 mM DCM or methanol as carbon and energy source. Liquid cultures for proteomic analysis (25 mL; 5 biological replicates each with DCM or with methanol) were obtained as described above, starting from independent colonies of the strain grown with DCM on a solid medium. Upon reaching the late stationary phase, each preculture was added to a 1 L Erlenmeyer flask containing 200 mL mineral medium and 10 mM DCM or methanol and further cultivated as above. DCM-grown cultures were harvested during the exponential phase (OD~600~ 0.03–0.09) after 5 days, 130–200 mL aliquots centrifuged at 4 °C at 8000 rpm for 10 min, and cell pellets stored at –80 °C until further analysis. Methanol-grown cultures were also harvested in the exponential phase (OD~600~ 0.16–0.21), and 45 mL culture aliquots were processed in the same way as DCM cultures to yield similar amounts of cell material for the two conditions (Supplementary Table S1).

### 2.2. *Cell Lysis and Enzymatic Proteolysis*

Proteins were extracted from cell pellets as described previously [16]. Briefly, 1.7 mg of bacterial paste were resuspended in 100 µL of lithium dodecyl sulfate (LDS) 1X lysis buffer consisting of 106 mM Tris/HCl, 141 mM Tris base, 2% LDS (*w/v*), 10% glycerol (*w/v*), 0.51 mM EDTA, 0.22 mM SERVA Blue G-250, and 0.175 mM phenol red, buffered at pH 8.5 and supplemented with 5% beta-mercaptoethanol (*v/v*), boiled at 99 °C for 5 min, and then sonicated 5 min in an ultrasonic water bath (VWR ultrasonic cleaner). Resulting suspensions were transferred into 2 mL screw-cap microtubes (Sarstedt, Nümbretch, Germany) containing 200 mg of 0.1 mm silica beads (MP Biomedicals, Solon, OH, USA) and further lysed using a Precellys Evolution instrument (Bertin Technologies, Montigny-le-Bretonneux, France) operated with 3 cycles of 30 s at 7800 rpm and 30 s of pause between each cycle. After lysis, samples were centrifuged at 16,000× g for 1 min, and the supernatant was collected to a new microcentrifuge tube before incubation at 99 °C for 5 min. Enzymatic digestion was performed as previously described [17]. Proteins (15 µg, i.e., 25 µL at 0.6 µg/µL) were subjected to a short electrophoresis migration (5 min) on NuPAGE 4–12% Bis-Tris gel, at 200 V in MES/SDS 1X running buffer. The whole proteome from each sample was recovered by excision of a unique gel band. Proteins were then reduced using 25 mM dithiothreitol (Sigma-Aldrich, St. Louis, Missouri, USA) in 50 mM NH₄HCO₃ (Sigma-Aldrich) at 56 °C for 10 min and alkylated with 55 mM iodoacetamide (Sigma-Aldrich) in 50 mM NH₄HCO₃ for 10 min at room temperature in the dark. Proteins were proteolyzed with the addition of 0.4 µg of trypsin gold (Promega, Madison, WI, USA) in 0.01% of Protease Max surfactant (Promega) and 50 mM NH₄HCO₃ per sample. After 15 min incubation on ice, excess trypsin was removed, and 50 µL of 50 mM NH₄HCO₃ containing 0.01% of Protease Max surfactant were added. Proteolysis was performed at 50 °C for 60 min, followed by acidification of the resulting extracted peptides with trifluoroacetic acid (TFA) to 0.5% final concentration.

### 2.3. *Mass Spectrometry and Data Interpretation*

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis was performed using an Ultimate 3000 nano-LC system coupled to a Q-Exactive HF mass spectrometer (Thermo Scientific, Waltham, MA, USA) operated as described previously [18]. Peptide samples (2 µL, corresponding to approximately 200 ng of peptides) were loaded on a reverse-phase PepMap 100 C18 µ-precolumn (5 µm, 100 Å, 300 µm i.d. × 5 mm, Thermo Fisher, Waltham, MA, USA) and then resolved on a nanoscale PepMap 100 C18 nanoLC column (3 µm, 100 Å, 75 µm i.d. × 50 cm, Thermo Fisher) at a flow rate of 0.3 µL·min⁻¹ using a 90 min gradient (4% B for 0 to 3 min, 4–25% B from 3 to 78 min, and 25–40%
B from 78 to 93 min) of mobile phases A (0.1% HCOOH/100% H₂O) and B (0.1% HCOOH/80% CH₃CN/20% H₂O). The mass spectrometer was operated in Top20 mode, with a scan range of MS acquisition from 350 to 1800 m/z and selection and fragmentation using 10 s dynamic exclusion time for the 20 most abundant precursor ions. Only ion precursors with a 2+ or 3+ charge were selected. HCD fragmentation was performed using a 27 eV normalized collision energy. Secondary ions were isolated with a window of 1.6 m/z.

The pan-proteomics database was built by merging the predicted proteomes from the annotated genomes of 13 *Hyphomicrobium* strains: *Hyphomicrobium* sp. NDB2Meth4, *Hyphomicrobium* sp. CS1GBMeth3, *Hyphomicrobium* sp. CS1BSMeth3, *Hyphomicrobium zavarzinii* ATCC 27496, *Hyphomicrobium* sp. 99, *Hyphomicrobium* sp. 802, *Hyphomicrobium nitrativorans* NL23, *Hyphomicrobium* sp. MC1, *Hyphomicrobium sulfonivorans*, *Hyphomicrobium facile*, *Hyphomicrobium denitrificans* 1NES1, and DCM-degrading strains *Hyphomicrobium denitrificans* ATCC 51888 and *Hyphomicrobium* sp. GJ21 (downloaded from NCBI in June 2018). The resulting database comprised 47,432 polypeptide sequences for a total of 14,834,376 residues. MS/MS spectra were assigned with this database using Mascot Daemon software version 2.6.1 (Matrix Science, London, UK), set with 5 ppm peptide tolerance and 0.02 Da MS/MS fragment tolerance, 2+ and 3+ peptide charge, a maximum of two missed cleavages, carbamidomethylation of cysteine as fixed modification, oxidation of methionine as variable modification, and trypsin as proteolytic enzyme. Peptides identified at *p*-value ≤ 0.05 in homology threshold mode and proteins identified with at least two distinct peptides were parsed using IRMa 1.31.1c software [19]. The false-positive rate was estimated to be below 1% for protein identification with the MASCOT decoy option search. Shared peptides were assigned with the most detected protein as assessed by specific peptides. Spectral counts, defined as the number of MS/MS spectra assigned per protein, were counted for all validated proteins using only non-ambiguous peptides as previously described [20]. Comparison of protein abundance between culture conditions was performed using the TFold test [21] and by defining four statistical groups: blue for a fold-change ≥ 1.5 and a *p*-value ≤ 0.05, orange for a *p*-value ≤ 0.05 and fold-change ≤ 1.5, green for a fold-change ≥ 1.5 and a *p*-value ≥ 0.05, and red for a fold-change ≤ 1.5 and a *p*-value ≥ 0.05. Data interpretation using the annotated genome of *Hyphomicrobium* sp. MC8b was done similarly, with proteins validated on the basis of at least two distinct peptide sequences. In this case, no parsimony rule was used, and spectral counts were evaluated by using only non-ambiguous peptides as previously described [20].

## 2.4. Genome Sequencing, Assembly and Annotation

Total DNA was prepared from a DCM-grown culture of strain MC8b using the MasterPure Complete DNA and RNA purification kit (Epicentre, Madison, WI, USA), flash-frozen and kept at -80 °C until further use. The DNA library for long-read sequencing was prepared using the Genomic DNA ligation kit (cat. no. SQK-LSK 108, Oxford Nanopore Technologies, Oxford, UK) according to the protocol provided, with a few modifications. Briefly, the DNA repair and end preparation steps were merged, and the DNA purification step between end preparation and native barcode ligation was omitted. After adaptor ligation, the library was purified using 0.6X AMPure XP beads (cat. no. A63881, Beckman Coulter, Brea, CA, USA). The resulting library was loaded on the flowcell (FLO-MN106) and sequenced on a MinION (Oxford Nanopore Technologies). After 24 h, the freshly prepared adaptor-ligated library was loaded again on the same flowcell. For short-read sequencing, 200 ng DNA was processed for library preparation using the KAPA HyperPlus kit (cat.no. 07962410001, Roche, Basel, Switzerland) without PCR amplification, according to the manufacturer's protocol. Enzymatic fragmentation time was 20 min, aiming at 300 bp average fragment size for the 2 × 150 bp paired-end sequencing run. The library was quantified using Qubit (Invitrogen, Carlsbad, CA, USA), its quality assessed by Bioanalyzer (Agilent, Santa Clara, CA, USA), and sequenced using NextSeq500 (Illumina, San Diego, CA, USA) at the Luxembourg Center for Systems Biomedicine (LCSB) sequencing platform (University of Luxembourg).
Raw ONT sequence data were base-called using Albacore Sequencing Pipeline Software (version 2.3.3) with r94_450bps_linear.cfg to match flowcell FLO-MIN106 and kit SQK-LSK108. Base-called reads were size-selected to be at least 1 kbp in length. Illumina sequence data were preprocessed by fastp (version 0.19.5) [22], and only properly paired, preprocessed reads were conserved.

Unicycler (version 0.4.7) [23] was used to assemble the genome from Illumina and Nanopore reads with the following options: -threads 28, -no_rotate -1, FORWARD IllUMINA_READS.fq.gz -2, REVERSE IllUMINA_READS.fq.gz -l, and SIZE_SELECTED_ONT_READS.fq. Associated softwares were SPAdes (version 3.13.0), racon (version 1.3.1), bowtie2-build (version 2.3.4.3), bowtie2 (version 2.3.4.3) [24], samtools (version 1.9) [25], java (version 1.8.0_162), and pilon (version 1.23) [26]. CheckM (version 1.0.12) [27] was used to verify assembly completeness with the following parameters: lineage_wf -t 12 -f CHECKM_RESULTS_FILE.txt -x fasta DIRECTORY_OF_UNICYCLER_ASSEMBLY OUT_DIRECTORY. The dependencies of CheckM and their respective versions were HMMER (version 3.2.1) [28], prodigal (version 2.6.3) [29], pplacer (version 1.1 alpha17) [30], and Python (version 2.7.14).

## 2.5. Data

Mass spectrometry proteomics data were submitted to the ProteomeXchange Consortium via the PRIDE partner repository under dataset identifiers PXD021816 and 10.6019/PXD021816 for the pan-proteomics study and PXD021817 and 10.6019/PXD021817 for the MC8b-genome derived proteomics analysis. The assembled genome was automatically annotated at Genoscope using the MicroScope pipeline [31] and uploaded on the MicroScope web platform (https://image.genoscope.cns.fr/microscope) for genome sequence analysis, and the obtained annotations were used for proteomics analysis. The genome sequence data for this study were deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB40967 (https://www.ebi.ac.uk/ena/browse/view/PRJEB40967).

# 3. Results

Preliminary experiments confirmed that DCM-degrading strain MC8b was also capable of growing with methanol, the reference growth substrate for strains of the *Hyphomicrobium* genus [32] to which strain MC8b had been tentatively affiliated [14]. We hypothesised that synthesis of the required dehalogenase to sustain strain MC8b growth on DCM may be regulated by the presence of its growth substrate, as observed in many dehalogenating strains [33]. Thus, we performed rapid differential proteomic analysis by pan-proteomics of cultures of the strain grown with either DCM or methanol.

## 3.1. Pan-Proteomics Strategy for Characterizing DCM-Degrading Strain MC8b of Unknown Genome Sequence

The pan-proteomics approach used to rapidly characterise the proteome of *Hyphomicrobium* sp. MC8b strain of unknown genome sequence (Figure 1) consisted of using a database with the predicted proteomes of 13 taxonomically closely related strains to *Hyphomicrobium* sp. MC8b (listed in Materials and Methods). Strain MC8b was grown with DCM or with methanol as the sole carbon and energy source, in quintuple replicates for each condition (Supplementary Table S1), and comparative proteomics were performed by nanoLC-MS/MS and interpretation of the recorded spectra using the pan-proteomics database.

## 3.2. Global Changes in the Proteome *Hyphomicrobium* sp. MC8b upon Growth with Dichloromethane

In total, 573,669 spectra were recorded from quintuple replicates of strain MC8b grown with DCM or with methanol and mapped against the pan-proteomics database. We assigned 290,425 spectra to peptide sequences and detected 2118 proteins in total (Supplementary Table S2), confirming that strain MC8b is affiliated to *Hyphomicrobium* (4329 CDS on average in the sequenced *Hyphomicrobium* strains of the pan-proteomics database). Of the identified proteins, 281 showed differential abundance between
DCM and methanol conditions as assessed by label-free shotgun proteomics, with 103 proteins more abundant with DCM (Figures 1 and 2 and Supplementary Table S3), indicating that growth with DCM requires significant molecular adjustments. Differentially abundant proteins present in the two DCM-degrading strains in the pan-proteomics database, *Hyphomicrobium* strains *H. denitrificans* ATCC 51888 [34] and *H. sp.* GJ21 [35], were not over-represented (Supplementary Figure S1), suggesting that bacterial adaptation to DCM utilisation strongly involves the core genome of *Hyphomicrobium*.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Pan-proteomics workflow to decipher key protein determinants of the dichloromethane utilisation pathway from the *Hyphomicrobium* sp. MC8b strain. DCM: dichloromethane.</figcaption></figure>

### 3.3. Proteomics-Driven Identification and Sequence Prediction of Strain MC8b DCM Dehalogenase

Four abundant peptide sequences specifically matching DcmA DCM dehalogenase of known genome sequence were unexpectedly obtained, although PCR amplification of the strongly conserved *dcmA* gene had been unsuccessful [14]. DCM dehalogenase DcmA was actually the protein with the highest differential abundance between DCM and methanol conditions (Supplementary Table S3). However, these four peptides accounted for only 18% coverage of the full-length protein. This suggested that strain MC8b features a divergent DcmA sequence compared to the very conserved sequence of all experimentally characterised Alphaproteobacterial DCM-degrading strains. We explored this hypothesis by analysing obtained tandem mass spectra using a custom set of *DcmA* proteins of
*Methylobacterium extorquens* DM4, *H. sp*. GJ21, *H. denitrificans*, and *Methylophilus* sp. DM11 strains in error-tolerant mode, in order to access peptide sequences differing slightly from known DcmA sequences. In this way, eight additional peptide sequences were defined from all five replicate cultures grown with DCM, yielding additional 40% coverage of a full-length DCM dehalogenase protein (Table 1). The predicted DcmA sequence of strain MC8b strongly differed from the 96–99% identical DCM dehalogenases of Alphaproteobacterial DCM-degraders [36]. Out of the total 924 spectral counts assigned to DcmA in error-tolerant mode, only 12 were observed in the methanol condition and for two of the total twelve detected peptides only. This further indicated that the synthesis of this protein, like that of all characterised DcmA homologs, is highly dependent on the presence of DCM.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparative pan-proteomics analysis of *Hyphomicrobium* sp. MC8b. Volcano plot of protein abundances under dichloromethane versus methanol growth conditions. Proteins are distributed depending on their abundance fold-change and *p*-value. Blue proteins show significant fold-change (≥1.5) and are significantly differentially abundant (*p* value ≤ 0.05). Proteins in orange (*p*-value ≤ 0.05 but fold-change ≤ 1.5), green (fold-change ≥ 1.5 but *p*-value ≥ 0.05), and red (fold-change ≤ 1.5, *p*-value ≥ 0.05) were not further considered.</figcaption></figure>

In all DCM-degrading strains with DcmA investigated so far, the gene cluster encoding DcmA also featured near-identical genes *dcmR*, which is involved in the regulation of DCM dehalogenase gene expression, and *dcmB* and *dcmC* of unknown function [11]. At only 30% sequence identity at the protein level, the closest functionally annotated protein to DcmB is a cephalosporin biosynthesis protein involved in hydroxylation/methyl transfer [37]. DcmC lacks functionally annotated full-length relatives, an unusual cysteine protease [38] with 29% identity over 50 amino acids being its closest annotated match. The corresponding gene products DcmR, DcmB, and DcmC were all detected in a recent proteomics study of *M. extorquens* DM4 and also showed higher abundance with DCM [39]. Here, DcmR and DcmB were also detected (Supplementary Table S3), but again with poor coverage. In contrast, DcmC was not detected, unlike its homolog in strain DM4 [39].
<table><caption>**Table 1.** Peptide sequences for DcmA of strain MC8b predicted by proteomic analysis.</caption><thead><tr><th>First</th><th>Last</th><th>Sequence <sup>1</sup></th><th>#Peptides <sup>2</sup></th><th>Analysis</th></tr></thead><tbody><tr><td>16</td><td>26</td><td>LLYHPGSQPCR</td><td>1</td><td>error-tolerant mode</td></tr><tr><td>54</td><td>84</td><td>QEFRDKYNPTGQVPILVDGEFTVWENVAIAR</td><td>3</td><td>error-tolerant mode</td></tr><tr><td>85</td><td>99</td><td>YVNEKFTGAGNWFGR</td><td>2</td><td>error-tolerant mode</td></tr><tr><td>105</td><td>119</td><td>AHINQFLQWYAYTLR</td><td>1</td><td>error-tolerant mode</td></tr><tr><td>120</td><td>143</td><td>LGGGAFHWNIFGCLIYGEKPYSPK</td><td>1</td><td>no mismatch</td></tr><tr><td>144</td><td>152</td><td>FSAEQNKGR</td><td>2</td><td>error-tolerant mode</td></tr><tr><td>153</td><td>168</td><td>ALLYEAMGTLENYWLR</td><td>3</td><td>error-tolerant mode</td></tr><tr><td>197</td><td>208</td><td>IIPDQVWQGFPK</td><td>2</td><td>error-tolerant mode</td></tr><tr><td>209</td><td>215</td><td>IAAWFKK</td><td>1</td><td>no mismatch</td></tr><tr><td>215</td><td>223</td><td>KLSARPHAK</td><td>2</td><td>error-tolerant mode</td></tr><tr><td>236</td><td>248</td><td>IIRGELTASMFKR</td><td>4</td><td>no mismatch</td></tr><tr><td>248</td><td>254</td><td>RKTAVLK</td><td>1</td><td>no mismatch</td></tr></tbody></table>

### 3.4. The Genome of Strain MC8b Features the Most Divergent Set of Dcm Genes Known So Far

Insights about the atypical DCM dehalogenation system of strain MC8b obtained by rapid genome-independent pan-proteomics analysis warranted sequencing of its genome. Sequencing of total DNA from DCM-grown cultures of strain MC8b using a combination of Oxford Nanopore and Illumina technologies yielded a circularised high-quality sequence for this strain. No plasmids were detected. Genome analysis was performed on the MicroScope platform at Genoscope [31]. At 4274 kb in size, the genome sequence of strain MC8b is typical of the *Hyphomicrobium* genus, with a single rRNA operon and a GC content of 59.65%. The genome was predicted to be complete by checkM [27] (100% completeness and 0.2% contamination with 1 marker duplicated). In total, 4574 CDS were predicted from the genome sequence [31]. Comparative analysis with other *Hyphomicrobium* genomes showed that strain MC8b is most closely related to the chloromethane-degrading strain *Hyphomicrobium* sp. MC1 [40] and to *H. sp.* 802, with about 80% closely homologous protein-encoding genes in synteny between the genome of strain MC8b and these two genomes.

We then checked whether the protein sequence of strain MC8b DCM dehalogenase predicted by pan-proteomics was confirmed by genome sequencing. In addition, we investigated how efficient analysis of our ad hoc pan-proteomics database of 13 predicted proteomes from genome-sequenced *Hyphomicrobium* strains had been in defining the DCM-specific proteome of strain MC8b, comparing it to an analysis performed with the predicted proteome of strain MC8b derived from its genome sequence.

The eight DcmA peptides predicted by pan-proteomics analysis in error-tolerant mode, together with the four peptides identical to known DcmA sequences (Table 1), were validated by the genome sequence of strain MC8b, confirming that its DCM dehalogenase significantly differed from other known Alphaproteobacterial sequences, with only the sequence of Betaproteobacterium *Methylophilus* sp. strain DM11 being more divergent [36] (Figure 3). Interestingly, the peptide correctly identified as YVNEKFAGTGNWFGR in the genome-based proteomic analysis had been predicted as the isomeric sequence YVNEKFTGAGNWFGR in the original analysis with the pan-proteomics database (Table 1). This is because the peptide sequence derived from the DNA sequence involves a two-residue difference from sequences in the pan-proteomics database, above the single residue mismatch threshold of error-tolerant proteomic analysis that is computationally practicable today.

<footnote n="1">¹ Residues in green highlight amino acid variations from reference DcmA sequences detected in error-tolerant mode.</footnote>

<footnote n="2">² Number (#) of different peptides detected by LC-MS/MS that correspond to all or part of the indicated sequence.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Sequence alignment of selected DCM dehalogenases covering the known sequence diversity of DCM dehalogenases, highlighting the highly divergent sequence of *Hyphomicrobium* sp. MC8b DcmA sequence detected and partially predicted using pan-proteomics analysis, and verified by sequencing of the MC8b genome. All peptides detected using the pan-proteomics *Hyphomicrobium* database analysed in normal or error-tolerant mode are shown in red and green, respectively. Additional peptides detected by analysis of obtained spectra using the predicted proteome of strain MC8b based on the genome sequence of the strain are shown by black lines. A star denotes a sequence difference in the DcmA sequence of strain MC8b compared to that of strain *H. sp.* GJ21 and/or *H. denitrificans* ATCC 51888 (HDEN). Green stars indicate sequence variations in the MC8b DcmA sequence that prevented detection of peptides by analysis of the *Hyphomicrobium* database in the absence of error tolerance, and black stars indicate sequence variations extending beyond the one-mismatch threshold for peptide detection in error-tolerant mode, respectively.</figcaption></figure>
Analysis of the genome context of gene *dcmA* gave further indications that the *dcmA*-dependent DCM dehalogenation system of strain MC8b differed from that encoded by the strongly conserved *dcmRABC* gene cluster of other DCM-degrading Alphaproteobacteria (Figure 4). DcmA, DcmB, and DcmR differ significantly from their homologs in other DCM-degrading strains. Further, gene *dcmC* was only present as a relic in strain MC8b (Figure 4). The corresponding short (57 aa)orf downstream of DcmB showed only 68% identity over 17 residues with the 184-residue reference DcmC of *Hyphomicrobium* strains ATCC 51888 and GJ21 and of strain DM4, and remained undetected by proteomics. The strong sequence variation observed in strain MC8b in the genomic context of *dcmA*, including the absence of transposases, usually flanking *dcm* genes in DCM-degrading strains (Figure 4), also explained the failure to detect the *dcm* gene cluster until now and confirmed that strain MC8b features the most divergent DcmA-based system discovered so far.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Comparison of *dcm* genes of *Hyphomicrobium* strains MC8b, ATCC 51888, and GJ21, and of the reference DCM-degrading Alphaproteobacterial strain *Methylobacterium extorquens* DM4. Homologous genes are shown in the same colour, with percentage identities to strain MC8b given at the protein level. Transposase genes are shown in grey, with homologous genes showing the same fill pattern. Numbers under the arrows for *Hyphomicrobium* strains ATCC 51888 and GJ21 and *Methylobacterium extorquens* strain DM4 refer to corresponding gene identifiers in Genbank.</figcaption></figure>

### 3.5. Further Insights from Proteomic Analysis Underlines the Power of the Pan-Proteomics Approach

In order to assess the value of results previously obtained by the pan-proteomics approach without a genome sequence for strain MC8b, proteomic data were then matched to the predicted proteome of strain MC8b derived from its genome sequence. In comparison to the 2118 proteins identified using the pan-proteomics database, 2101 were confirmed with the predicted proteome of strain MC8b. With regard to differentially abundant proteins (Figure 5), 328 proteins were identified with the specific MC8b proteome database (Supplementary Table S4), as compared to 281 with the *Hyphomicrobium* pan-proteomics database (Supplementary Table S3). About a third (126) additional proteins were identified as more abundant with DCM using the strain-specific proteome (Figure 5, Supplementary Table S4). Along the same lines, we also checked whether the distribution of COG categories [41] for proteins with significant differential abundance between DCM and methanol conditions differed from that encoded in the genome. Proteins associated with DNA repair (L) and cell wall structure and biogenesis (M) were over-represented in the proteome (Supplementary Table S5). Conversely, proteins associated with transcription and translation were under-represented. This was not only observed in the analysis of the MC8b proteome derived from its genome sequence but already in the initial analysis with the ad hoc database constructed from the theoretical proteomes of the 13 *Hyphomicrobium* strains of known genome sequence (Supplementary Table S5). This confirmed the power of the pan-proteomics approach for rapid analysis of strains with functions of interest, within a well-investigated taxonomical framework, and in the absence of specific genome sequence information.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** *Hyphomicrobium* sp. MC8b differential proteomics using the genome sequence of the strain. Volcano plot of protein abundances in dichloromethane versus methanol growth conditions. Colour code is as in Figure 2, with open blue symbols highlighting proteins with significant differential abundance that were additionally identified from the proteome of strain MC8b predicted from the genome sequence.</figcaption></figure>

**4. Discussion**

The pan-proteomics concept [4,7] can be applied to a relatively large number of prokaryotic isolates, as the genomes of many types of microorganisms have now been sequenced. While the small number of available genomes for eukaryotes still represents a limitation for the application of pan-proteomics, the similar concept of "homology-driven proteomics" [42,43] has proved helpful for the analysis of proteomes of unsequenced animals and plants [44,45]. The power of the pan-proteomics approach was supported by the results obtained in the present study. First, it was unexpected that strain *Hyphomicrobium* sp. MC8b contains DCM dehalogenase DcmA. The atypical DcmA-based system of strain MC8b was evidenced by unprecedented differences in DcmA and DcmB sequences, while DcmC remained undetected. Pan-proteomics results were then confirmed by genome sequencing. Notably, the gene encoding the differentially abundant putative nitrilase (Figure 4, Table S3), whose homolog is located nearby *dcmRABC* genes in the genome of *M. extorquens* DM4 (Figure 4), was found next to the *dcmC* relic in the MC8b genome, suggesting that it may be associated with DCM metabolism as well.

Indeed, the possibility of analysing differential abundance in the condition of interest compared to a reference condition represents a strong asset to identify proteins involved in a function of interest by pan-proteomics. In the field of bacterial dehalogenation, in particular, synthesis of the proteins involved in dehalogenation of a given organohalide often depends on its presence [33]. Identification of proteins of interest for a particular function will be more challenging when only one cultivation condition is available. For *Candidatus* Dichloromethanomonas elyunquensis growing with DCM under strictly anoxic conditions, for example, potential dehalogenase candidates were tentatively proposed based on DCM-grown cultures as the only available growth condition for this system [46].

Differential pan-proteomics may also provide clues on associated key proteins and potentially also corresponding metabolic pathways and adaptations related to the function of interest. For instance,
proteomic analysis of the reference DCM-degrading strain *M. extorquens* DM4 [39] suggested that DCM metabolism by a DcmA-dependent system triggers adaptations related to DNA genotoxicity, acid and chloride production, and membrane integrity, thereby confirming previous work based on mutagenesis studies [11,47]. These findings were confirmed here at the level of general gene functional classes (COGs) (Supplementary Table S5). The observed shift in COG distribution of proteins more abundant in the DCM condition compared to that in the theoretical proteome predicted from the genome sequence highlighted general functions associated with DCM metabolism. This shift was already detected by pan-proteomics analysis, i.e., without knowledge of the strain-specific theoretical proteome (Supplementary Table S5).

Nevertheless, the modest overlap in proteins with differential abundance with the previous study on the reference DCM-degrading strain *M. extorquens* DM4 is noteworthy. Only nine proteins, including DCM dehalogenase DcmA and DcmB protein of unknown function, were identified as differentially abundant in both *Hyphomicrobium* strain MC8b and *M. extorquens* strain DM4 (Supplementary Table S3). On the one hand, two of the proteins detected as differentially abundant in both DM4 and MC8b proteomics studies had already been associated with DCM metabolism. The squalene hopene cyclase *shc* gene was identified as essential for growth with DCM [47,48], while a transglycosylase/transpeptidase homolog (HYPMC8B_3734, annotated as penicillin-binding protein; METDI4661 in strain DM4) showed DCM-dependent synthesis [47]. These findings confirm the likely importance of these two proteins as part of a specific ensemble of proteins associated with bacterial growth on DCM and involving DNA repair and envelope processes (Supplementary Table S5). On the other hand, the small number of differentially abundant proteins shared by strains MC8b and DM4 growing with DCM contrasts with the fact that the two strains share 1136 homologous proteins with over 50% identity at the protein level, with a similar number of detected proteins in the two studies (2453 proteins for strain MC8b versus 2878 for strain DM4 [39]). To us, this suggests that adaptation to dehalogenation of DCM involves specific changes in expression of the taxonomically defined core genome following acquisition of genes for DCM utilisation [49,50], in keeping with the broad functional categories associated with transformation of DCM (Supplementary Table S5), and as suggested by transcriptional studies [51]. In other words, genes involved in adaptation to DCM are not limited to a specific set of genes of DCM-degrading strains, as evidenced by differentially abundant proteins shared by many and sometimes all 13 *Hyphomicrobium* strains of the pan-proteomics database (Supplementary Figure S1).

Finally, and with regard to the process of horizontal transfer of *dcm* genes itself, the lack of IS elements flanking *dcm* genes in strain MC8b, as well as the unusual arrangement and sequence of its *dcm* gene cluster (Figure 4), raises new questions on the evolution and subsequent dissemination of the capacity to grow with DCM in the bacterial world. Clearly, the *dcm* gene cluster of DCM-degrading strains with DcmA DCM dehalogenase may not be as conserved or as essential in defining DCM-degrading strains as generally believed until now.

**5. Conclusions**

Pan-proteomics allowed us to uncover a hitherto undetected dehalogenase enzymatic system in strain MC8b. Moreover, differential proteomics allowed to identify DCM-induced proteins that may also be associated with dehalogenation. Thus, a custom pan-proteomics database allows one to propose proteins involved in key functions of interest in the absence of genome sequence. The power of this approach will increase with the number of sequenced genomes in databases, and make possible fast fingerprinting of strain metabolism under conditions where genome or transcriptome sequencing is unfeasible.

The results obtained for strain MC8b also put the evolution of DCM dehalogenases in a new light. Indeed, the strong differences in *dcm* genes observed for strain MC8b strongly suggest that horizontal gene transfer, while important today in environments contaminated with DCM, may represent a relatively recent feature of the evolution of DCM dehalogenases capable of supporting bacterial growth with DCM. This, to us, is an incentive to renew explorations of enzymatic systems of DCM
dehalogenation using state-of-the-art function-based approaches. More generally, pan-proteomics may help to rapidly discover in unsequenced strains still unchartered solutions developed by the microbial world to degrade organohalides over the eons.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2076-2607/8/12/1876/s1.
Table S1: Cultures used for comparative proteomics analysis; Table S2: Overview of obtained data from LC MS/MS analysis of the proteome of strain MC8b grown with DCM or methanol; Table S3: Proteins detected by pan-proteomics. See separate Excel file; Table S4: Proteins detected from the annotated genome sequence of *Hyphomicrobium* sp. MC8b. See separate Excel file; Table S5: Distribution (in percent) of COG categories for differentially abundant proteins of *Hyphomicrobium* sp. MC8b analysed with the pan-proteomics database or the MC8bgenome-derived database, and for the total predicted proteome of strain MC8b; Figure S1: Proteomes in the pan-proteomics database.

**Author Contributions:** S.V. and J.A. conceived the study; J.A., M.R., P.W. and S.V. arranged for resources, materials, and access to equipment; K.H., E.G., C.L.S., A.M. and R.H. performed experiments; K.H., E.G., C.C.L., F.B., E.E.L.M., B.A.-B., J.A. and S.V. analysed and interpreted results; K.H., B.A.-B., J.A. and S.V. drafted the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the French Agence Nationale de la Recherche (ANR), grant number ANR-17-CE07-0009 to the dehalofluidX project.

**Acknowledgments:** Neil C. Sturchio is gratefully acknowledged for his generous gift of a liquid culture of strain MC8b to S.V. at the tail end of the last millennium.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**Abbreviations**

<table><tr><td>CDS</td><td>coding DNA sequence</td></tr><tr><td>COG</td><td>cluster of orthologous groups</td></tr><tr><td>DCM</td><td>dichloromethane</td></tr><tr><td>LC-MS/MS</td><td>liquid chromatography-tandem mass spectrometry</td></tr><tr><td>PSM</td><td>peptide spectrum match</td></tr></table>

**References**

1. Thompson, L.R.; Sanders, J.G.; McDonald, D.; Amir, A.; Ladau, J.; Locey, K.J.; Prill, R.J.; Tripathi, A.; Gibbons, S.M.; Ackermann, G.; et al. A communal catalogue reveals Earth's multiscale microbial diversity. *Nature* **2017**, *551*, 457–463. [CrossRef]

2. Bull, A.T. *Microbial Diversity and Bioprospecting*; American Society of Microbiology: Washington, DC, USA, 2004. [CrossRef]

3. Gouveia, D.; Grenga, L.; Pible, O.; Armengaud, J. Quick microbial molecular phenotyping by differential shotgun proteomics. *Environ. Microbiol.* **2020**, *22*, 2996–3004. [CrossRef]

4. Broadbent, J.A.; Broszczak, D.A.; Tennakoon, I.U.K.; Huygens, F. Pan-proteomics, a concept for unifying quantitative proteome measurements when comparing closely-related bacterial strains. *Expert Rev. Proteom.* **2016**, *13*, 355–365. [CrossRef]

5. Silva, W.M.; Sousa, C.S.; Oliveira, L.C.; Soares, S.C.; Souza, G.F.M.H.; Tavares, G.C.; Resende, C.P.; Folador, E.L.; Pereira, F.L.; Figueiredo, H.; et al. Comparative proteomic analysis of four biotechnological strains *Lactococcus lactis* through label-free quantitative proteomics. *Microb. Biotechnol.* **2019**, *12*, 265–274. [CrossRef]

6. Tavares, G.C.; Pereira, F.L.; Barony, G.M.; Rezende, C.P.; da Silva, W.M.; de Souza, G.H.M.F.; Verano-Braga, T.; de Carvalho Azevedo, V.A.; Leal, C.A.G.; Figueiredo, H.C.P. Delineation of the pan-proteome of fish-pathogenic *Streptococcus agalactiae* strains using a label-free shotgun approach. *BMC Genom.* **2019**, *20*, 11. [CrossRef]

7. Murugaiyan, J.; Eravci, M.; Weise, C.; Roesler, U.; Sprague, L.D.; Neubauer, H.; Wareth, G. Pan-proteomic analysis and elucidation of protein abundance among the closely related *Brucella* species, *Brucella abortus* and *Brucella melitensis*. *Biomolecules* **2020**, *10*, 836. [CrossRef]

8. Atashgahi, S.; Liebensteiner, M.G.; Janssen, D.B.; Smidt, H.; Stams, A.J.M.; Sipkema, D. Microbial synthesis and transformation of inorganic and organic chlorine compounds. *Front. Microbiol.* **2018**, *9*, 3079. [CrossRef]
9. Janssen, D.B. Biocatalysis by dehalogenating enzymes. *Adv. Appl. Microbiol.* **2007**, *61*, 233–252. [CrossRef]

10. Gribble, G.W. Newly discovered naturally occurring organohalogens. *Arkivoc* **2018**, *2018*, 372–410. [CrossRef]

11. Muller, E.E.L.; Bringel, F.; Vuilleumier, S. Dichloromethane-degrading bacteria in the genomic age. *Res. Microbiol.* **2011**, *162*, 869–876. [CrossRef]

12. Ergas, S.J.; Kinney, K.; Fuller, M.E.; Scow, K.M. Characterization of compost biofiltration system degrading dichloromethane. *Biotechnol. Bioeng.* **1994**, *44*, 1048–1054. [CrossRef]

13. Heraty, L.J.; Fuller, M.E.; Huang, L.; Abrajano, T.; Sturchio, N.C. Isotopic fractionation of carbon and chlorine by microbial degradation of dichloromethane. *Org. Geochem.* **1999**, *30*, 793–799. [CrossRef]

14. Nikolausz, M.; Kappelmeyer, U.; Nijenhuis, I.; Ziller, K.; Kästner, M. Molecular characterization of dichloromethane-degrading *Hyphomicrobium* strains using 16S rDNA and DCM dehalogenase gene sequences. *Syst. Appl. Microbiol.* **2005**, *28*, 582–587. [CrossRef]

15. Roselli, S.; Nadalig, T.; Vuilleumier, S.; Bringel, F. The 380 kb pCMU01 plasmid encodes chloromethane utilization genes and redundant genes for vitamin B12- and tetrahydrofolate-dependent chloromethane metabolism in *Methylobacterium extorquens* CM4: A proteomic and bioinformatics study. *PLoS ONE* **2013**, *8*, e56598. [CrossRef]

16. Mappa, C.; Pible, O.; Armengaud, J.; Alpha-Bazin, B. Assessing the ratio of *Bacillus* spores and vegetative cells by shotgun proteomics. *Environ. Sci. Pollut. Res.* **2018**, 1–9. [CrossRef]

17. Hartmann, E.M.; Allain, F.; Gaillard, J.C.; Pible, O.; Armengaud, J. Taking the shortcut for high-throughput shotgun proteomic analysis of bacteria. *Methods Mol. Biol.* **2014**, *1197*, 275–285. [CrossRef]

18. Klein, G.; Mathé, C.; Biola-Clier, M.; Devineau, S.; Drouineau, E.; Hatem, E.; Marichal, L.; Alonso, B.; Gaillard, J.C.; Lagniel, G.; et al. RNA-binding proteins are a major target of silica nanoparticles in cell extracts. *Nanotoxicology* **2016**, *10*, 1555–1564. [CrossRef]

19. Dupierris, V.; Masselon, C.; Court, M.; Kieffer-Jaquinod, S.; Bruley, C. A toolbox for validation of mass spectrometry peptides identification and generation of database: IRMa. *Bioinformatics* **2009**, *25*, 1980–1981. [CrossRef]

20. Cogne, Y.; Almunia, C.; Gouveia, D.; Pible, O.; François, A.; Degli-Esposti, D.; Geffard, O.; Armengaud, J.; Chaumot, A. Comparative proteomics in the wild: Accounting for intrapopulation variability improves describing proteome response in a *Gammarus pulex* field population exposed to cadmium. *Aquat. Toxicol.* **2019**, *214*, 105244. [CrossRef]

21. Carvalho, P.C.; Yates, J.R., III; Barbosa, V.C. Improving the TFold test for differential shotgun proteomics. *Bioinformatics* **2012**, *28*, 1652–1654. [CrossRef]

22. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. *Bioinformatics* **2018**, *34*, i884–i890. [CrossRef] [PubMed]

23. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. *PLoS Comput. Biol.* **2017**, *13*, e1005595. [CrossRef] [PubMed]

24. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. *Nat. Methods* **2012**, *9*, 357–359. [CrossRef] [PubMed]

25. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The sequence alignment/map format and SAMtools. *Bioinformatics* **2009**, *25*, 2078–2079. [CrossRef]

26. Walker, B.J.; Abeel, T.; Shea, T.; Priest, M.; Abouelliel, A.; Sakthikumar, S.; Cuomo, C.A.; Zeng, Q.; Wortman, J.; Young, S.K.; et al. Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement. *PLoS ONE* **2014**, *9*, e112963. [CrossRef]

27. Parks, D.H.; Imelfort, M.; Skennerton, C.T.; Hugenholtz, P.; Tyson, G.W. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. *Genome Res.* **2015**, *25*, 1043–1055. [CrossRef]

28. Eddy, S.R. Profile hidden Markov models. *Bioinformatics* **1998**, *14*, 755–763. [CrossRef]

29. Hyatt, D.; Chen, G.-L.; LoCascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. *BMC Bioinform.* **2010**, *11*, 119. [CrossRef]

30. Matsen, F.A.; Kodner, R.B.; Armbrust, E.V. pplacer: Linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree. *BMC Bioinform.* **2010**, *11*, 538. [CrossRef]
31. Vallenet, D.; Calteau, A.; Dubois, M.; Amours, P.; Bazin, A.; Beuvin, M.; Burlot, L.; Bussell, X.; Fouteau, S.; Gautreau, G.; et al. MicroScope: An integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis. *Nucleic Acids Res.* **2020**, *48*, D579–D589. [CrossRef]

32. Gliesche, C.; Fesefeldt, A.; Hirsch, P. Hyphomicrobium. In *Bergey's Manual of Systematics of Archaea and Bacteria*; Wiley: Hoboken, NJ, USA, 2015. [CrossRef]

33. Maucourt, B.; Vuilleumier, S.; Bringel, F. Transcriptional regulation of organohalide pollutant utilisation in bacteria. *FEMS Microbiol. Rev.* **2020**, *44*, 189–207. [CrossRef] [PubMed]

34. Brown, P.J.B.; Kysela, D.T.; Buechlein, A.; Hemmerich, C.; Brun, Y.V. Genome sequences of eight morphologically diverse Alphaproteobacteria. *J. Bacteriol.* **2011**, *193*, 4567–4568. [CrossRef] [PubMed]

35. Bringel, F.; Postema, C.P.; Mangenot, S.; Bibi-Triki, S.; Chaignaud, P.; Haque, M.F.U.; Gruffaz, C.; Hermon, L.; Louhichi, Y.; Maucourt, B.; et al. Genome sequence of the dichloromethane-degrading bacterium *Hyphomicrobium* sp. strain GJ21. *Genome Announc.* **2017**, *5*, 622. [CrossRef] [PubMed]

36. Vuilleumier, S.; Ivoš, N.; Dean, M.; Leisinger, T. Sequence variation in dichloromethane dehalogenases/glutathione S-transferases. *Microbiology* **2001**, *147*, 611–619. [CrossRef] [PubMed]

37. Coque, J.J.R.; Enguita, F.J.; Martín, J.F.; Liras, P. A 2-protein component 7-alpha-cephem-methoxylase encoded by 2 genes of the cephamycin-C cluster converts cephalosporin-C to 7-methoxycephalosporin-C. *J. Bacteriol.* **1995**, *177*, 2230–2235. [CrossRef]

38. Mottram, J.C.; North, M.J.; Barry, J.D.; Coombs, G.H. A cysteine proteinase cDNA from *Trypanosoma brucei* predicts an enzyme with an unusual C-terminal extension. *FEBS Lett.* **1989**, *258*, 211–215. [CrossRef]

39. Bibi-Triki, S.; Husson, G.; Maucourt, B.; Vuilleumier, S.; Carapito, C.; Bringel, F. N-terminome and proteogenomic analysis of the *Methylobacterium extorquens* DM4 reference strain for dichloromethane utilization. *J. Proteom.* **2018**, *179*, 131–139. [CrossRef]

40. Vuilleumier, S.; Nadalig, T.; Haque, M.F.U.; Magdelenat, G.; Lajus, A.; Roselli, S.; Muller, E.E.L.; Gruffaz, C.; Barbe, V.; Médigue, C.; et al. Complete genome sequence of the chloromethane-degrading *Hyphomicrobium* sp. strain MC1. *J. Bacteriol.* **2011**, *193*, 5035–5036. [CrossRef]

41. Galperin, M.Y.; Kristensen, D.M.; Makarova, K.S.; Wolf, Y.I.; Koonin, E.V. Microbial genome analysis: The COG approach. *Brief. Bioinform.* **2019**, *20*, 1063–1070. [CrossRef]

42. Junqueira, M.; Spirin, V.; Balbuena, T.S.; Thomas, H.; Adzhubei, I.; Sunyaev, S.; Shevchenko, A. Protein identification pipeline for the homology-driven proteomics. *J. Proteom.* **2008**, *71*, 346–356. [CrossRef]

43. Waridel, P.; Frank, A.; Thomas, H.; Surendranath, V.; Sunyaev, S.; Pevzner, P.; Shevchenko, A. Sequence similarity-driven proteomics in organisms with unknown genomes by LC-MS/MS and automated *de novo* sequencing. *Proteomics* **2007**, *7*, 2318–2329. [CrossRef] [PubMed]

44. Liska, A.J.; Popov, A.V.; Sunyaev, S.; Coughlin, P.; Habermann, B.; Shevchenko, A.; Bork, P.; Karsenti, E.; Shevchenko, A. Homology-based functional proteomics by mass spectrometry: Application to the *Xenopus* microtubule-associated proteome. *Proteomics* **2004**, *4*, 2707–2721. [CrossRef] [PubMed]

45. Saha, B.; Bhattacharya, S.G. Charting novel allergens from date palm pollen (*Phoenix sylvestris*) using homology driven proteomics. *J. Proteom.* **2017**, *165*, 1–10. [CrossRef] [PubMed]

46. Kleindienst, S.; Chourey, K.; Chen, G.; Murdoch, R.W.; Higgins, S.A.; Iyer, R.; Campagna, S.R.; Mack, E.E.; Seger, E.S.; Hettich, R.L.; et al. Proteogenomics reveals novel reductive dehalogenases and methyltransferases expressed during anaerobic dichloromethane metabolism. *Appl. Environ. Microbiol.* **2019**, *85*, 2768. [CrossRef]

47. Muller, E.E.L.; Hourcade, E.; Louhichi-Jelail, Y.; Hammann, P.; Vuilleumier, S.; Bringel, F. Functional genomics of dichloromethane utilization in *Methylobacterium extorquens* DM4. *Environ. Microbiol.* **2011**, *13*, 2518–2535. [CrossRef]

48. Bradley, A.S.; Swanson, P.K.; Muller, E.E.L.; Bringel, F.; Caroll, S.M.; Pearson, A.; Vuilleumier, S.; Marx, C.J. Hopanoid-free *Methylobacterium extorquens* DM4 overproduces carotenoids and has widespread growth impairment. *PLoS ONE* **2017**, *12*, e0173323. [CrossRef]

49. Michener, J.K.; Neves, A.A.C.; Vuilleumier, S.; Bringel, F.; Marx, C.J. Effective use of a horizontally-transferred pathway for dichloromethane catabolism requires post-transfer refinement. *eLife* **2014**, *3*, 04279. [CrossRef]
50. Michener, J.K.; Vuilleumier, S.; Bringel, F.; Marx, C.J. Transfer of a catabolic pathway for chloromethane in *Methylobacterium* strains highlights different limitations for growth with chloromethane or with dichloromethane. *Front. Microbiol.* **2016**, *7*, 1116. [CrossRef]

51. Chaignaud, P.; Maucourt, B.; Weiman, M.; Alberti, A.; Kolb, S.; Cruveiller, S.; Vuilleumier, S.; Bringel, F. Genomic and transcriptomic analysis of growth-supporting dehalogenation of chlorinated methanes in *Methylobacterium*. *Front. Microbiol.* **2017**, *8*, 1600. [CrossRef]

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).