Next Article in Journal
EPA and DHA Inhibit Myogenesis and Downregulate the Expression of Muscle-related Genes in C2C12 Myoblasts
Next Article in Special Issue
Metal-Resistance in Bacteria: Why Care?
Previous Article in Journal
A High-Quality De novo Genome Assembly from a Single Mosquito Using PacBio Sequencing
Previous Article in Special Issue
Using a Chemical Genetic Screen to Enhance Our Understanding of the Antimicrobial Properties of Gallium against Escherichia coli
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 1
10.3390/genes10010063
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic and Transcriptomic Changes That Mediate Increased Platinum Resistance in Cupriavidus metallidurans

by
Md Muntasir Ali
1,2,
Ann Provoost
1,
Laurens Maertens
1,3,
Natalie Leys
1,
Pieter Monsieurs
1,
Daniel Charlier
2 and
Rob Van Houdt
1,*
1
Microbiology Unit, Belgian Nuclear Research Centre (SCK•CEN), 2400 Mol, Belgium
2
Research Group of Microbiology, Department of Bioengineering Sciences, Vrije Universiteit Brussel, 1050 Brussel, Belgium
3
Research Unit in Biology of Microorganisms (URBM), Faculty of Sciences, UNamur, 5000 Namur, Belgium
*
Author to whom correspondence should be addressed.
Genes 2019, 10(1), 63; https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010063
Submission received: 7 December 2018 / Revised: 11 January 2019 / Accepted: 15 January 2019 / Published: 18 January 2019
(This article belongs to the Special Issue Genomics of Bacterial Metal Resistance)
Browse Figures
Versions Notes

Abstract

:
The extensive anthropogenic use of platinum, a rare element found in low natural abundance in the Earth’s continental crust and one of the critical raw materials in the EU innovation partnership framework, has resulted in increased concentrations in surface environments. To minimize its spread and increase its recovery from the environment, biological recovery via different microbial systems is explored. In contrast, studies focusing on the effects of prolonged exposure to Pt are limited. In this study, we used the metal-resistant Cupriavidus metallidurans NA4 strain to explore the adaptation of environmental bacteria to platinum exposure. We used a combined Nanopore–Illumina sequencing approach to fully resolve all six replicons of the C. metallidurans NA4 genome, and compared them with the C. metallidurans CH34 genome, revealing an important role in metal resistance for its chromid rather than its megaplasmids. In addition, we identified the genomic and transcriptomic changes in a laboratory-evolved strain, displaying resistance to 160 µM Pt4+. The latter carried 20 mutations, including a large 69.9 kb deletion in its plasmid pNA4_D (89.6 kb in size), and 226 differentially-expressed genes compared to its parental strain. Many membrane-related processes were affected, including up-regulation of cytochrome c and a lytic transglycosylase, down-regulation of flagellar and pili-related genes, and loss of the pNA4_D conjugative machinery, pointing towards a significant role in the adaptation to platinum.

1. Introduction

Platinum (Pt) is a rare element that is found in low natural abundance (0.4 parts per billion) in the Earth’s continental crust [1,2]. It is extensively used in industry, vehicle exhaust catalysts (VECs), and anticancer drugs [3], with cisplatin being one of the potent anti-cancer drugs in use [4]. Anthropogenic uses and emissions of platinum have resulted in increased concentrations (0.5–1.4 ton year−1) in surface environments, which could negatively impact natural habitats, especially because of the solubility of some forms of platinum [3,5]. It can enter waters, sediments, and soils and eventually reach the food chain [3]. Therefore, effective measures must be taken to minimize its spread and increase its recovery from the environment.
Platinum is also one of the critical raw materials in the EU innovation partnership framework, which were selected because of their high economic importance and high supply risk [6]. This makes platinum an important candidate for biological recovery from waste streams and other environmental niches. As bacterial communities have been naturally associated with platinum-group mineral grains [7], efforts have been made to explore the usability of microorganisms for the recovery of platinum group metals from the environment [8]. For instance, biological recovery of platinum has been shown for halophilic microbial communities, indicating that Pt from waste streams can be transformed into Pt-rich biomass, which in turn can be used as input for the refinery of precious metals [9]. In addition, Pt biosorption has also been studied for axenic bacterial cultures, including sulfate-reducing bacteria such as Desulfovibrio desulfuricans, Desulfovibrio fructosivorans, and Desulfovibrio vulgaris [10], as well as Shewanella oneidensis, Cupriavidus metallidurans, Geobacter metallireducens, Pseudomonas stutzeri, and Bacillus toyonensis [11,12]. The sulfate-reducing Desulfovibrio spp. and metal-ion reducing Shewanella algae have the capability to reduce Pt to zero state and form Pt nanoparticles in their periplasmic space [13,14]. It has been hypothesized that Desulfovibrio spp. can use Pt, as well as palladium, as terminal electron acceptors in their energy production pathway via cytochrome c3 [15,16,17]. This promotes nanoparticle formation on the cell surface, preventing re-entry and acting as catalysts for further metal reduction [17]. Cupriavidus sp. also showed similar nanoparticle formation in the presence of palladium, according to an equivalent strategy [18].
Studying the interaction between platinum and bacteria showed that platinum inhibits cell division and enhances filamentous growth of Escherichia coli [19], Caulobacter crescentus and Hyphomicrobium sp. [20]. It inhibits DNA synthesis and DNA repair functions were shown to be essential for growth in the presence of platinum, as E. coli mutants deficient in DNA repair functions are unable to grow in the presence of platinum [21,22]. Similar to other DNA synthesis inhibitors, Pt also induces prophages from lysogenic E. coli strains [23]. The mutagenic ability of different platinum compounds has also been demonstrated in Salmonella enterica subsp. enterica serovar Typhimurium strains [24].
It is clear that most studies have analyzed the biological immobilization of Pt and its possible applications. However, only a limited number of studies focused on the effect of prolonged exposure to Pt, as would be the case for environmental bacteria in Pt-contaminated waters, soils, and sediments. For instance, Maboeta et al. showed that enzymatic activities and viable biomass were impacted in a platinum tailing disposal facility associated with mining activities [25].
In this study, we used C. metallidurans NA4 as a model to explore the adaptation of environmental bacteria to platinum exposure. It has been extensively studied for its resistance to a variety of metal (oxyan)ions [26,27,28]. We used a combinatorial sequencing approach to fully resolve the C. metallidurans NA4 genome consisting of six replicons [29], compared its genome with that of C. metallidurans CH34, and identified the genomic and transcriptomic changes in a laboratory-evolved strain, displaying increased resistance to platinum.

2. Materials and Methods

2.1. Strains, Media, and Culture Conditions

C. metallidurans NA4 was routinely cultured at 30 °C in Lysogeny broth (LB) or Tris-buffered mineral medium (MM284) supplemented with 0.2% (w/v) gluconate [30]. For culturing on solid medium, 1.5% agar (Thermo Scientific, Oxoid, Hampshire, UK) was added; liquid cultures were grown in the dark on a rotary shaker at 150 rpm. Metal salts used included PtCl4, Na2PdCl4, ZnSO4.7H2O, NiCl2.6H2O, CuSO4.5H2O and AgNO3. (Sigma-Aldrich, Overijse, Belgium).

2.2. Determination of the Minimal Inhibitory Concentration and Generation of Pt-Resistant Mutants

The minimal inhibitory concentration (MIC) of Pt4+, Pd2+, Zn2+, Ni2+, Cu2+, and Ag+ was determined using the broth dilution method in a 96-well plate containing a concentration gradient of the corresponding metals [31]. To select for C. metallidurans mutants displaying increased platinum resistance, a serial passage experiment was performed by continuous exposure to subinhibitory concentrations of Pt4+ using the gradient MIC method [32].

2.3. Plasmid Isolation and Restriction Digestion

The extraction of megaplasmids was based on the method proposed by Andrup et al. [33]. Extracted plasmid DNA was separated by horizontal gel electrophoresis (23 cm-long 0.5% Certified Megabase agarose gel (Bio-Rad, Temse, Belgium) in 1X Tris-Borate-EDTA buffer, 100 V, 20 h) in a precooled (4 °C) electrophoresis chamber. After GelRed staining (30 min + overnight destaining at 4 °C in ultrapure water), DNA was visualized and images captured under UV light transillumination (Fusion Fx, Vilber Lourmat, Collégien, France). To confirm the presence and size of the smaller plasmid in NA4Pt (Pt4+ resistant mutant of NA4), plasmid DNA was isolated with the Wizard® Plus SV Miniprep DNA Purification System (Promega, Leiden, The Netherlands). The isolated DNA was used for restriction digestion with PagI (Fisher Scientific, Merelbeke, Belgium). The products were separated on a 0.6% agarose (Molecular Biology Grade, Eurogentec, Belgium) gel to visualize the individual fragments together with the GeneRuler 1 kb plus ladder (Fisher Scientific, Merelbeke, Belgium).

2.4. Motility, Scanning Electron Microscopy (SEM), and Flow Cytometry

For testing motility, C. metallidurans NA4 and NA4Pt were grown in LB media until the OD600 reached 0.6. Five µL of the culture was then stab inoculated onto a LB plate containing 0.3% agar. The radius of the growth pattern was measured after 24 h.
For SEM, cells were grown in MM284 in normal growth conditions, centrifuged (5000 rpm for 8 min), washed in Milli-Q water, and fixed with 3% glutaraldehyde solution at 4 °C (3 h). Cells were sputter coated (22 nm) with gold and examined under SEM at an accelerating voltage of 10 kV.
For flow cytometry, bacterial cell suspensions (OD600 = 0.6) were diluted 1000 times in 0.2 μm filtered Tris-buffered mineral medium (MM284), Next, SYBR green (Sigma Aldrich) dye was added and incubated at 37 °C for 20 min. Stained bacterial suspensions were analyzed on the Accuri C6 flow cytometer (BD, Erembodegem, Belgium).

2.5. Genome Sequencing

Total DNA from C. metallidurans NA4 and NA4Pt was isolated using the QIAamp DNA mini kit (Qiagen, Venlo, the Netherlands). The parental strain NA4 [29] was resequenced using a combination of Illumina and Nanopore sequencing. Illumina sequencing was performed on the Illumina HiSeq 2500 platform using 2 × 75 bp paired-end sequencing (Baseclear, Leiden, Netherlands). Nanopore sequencing was performed in-house using the MinION device with an R9 flow cell and the Rapid Sequencing kit. Strain NA4Pt was sequenced using the Illumina Miseq platform (40× coverage; MicrobesNG, Birmingham, UK).

2.6. Genome Assembly

Genome assembly was performed using the pre-assembled contigs based on the 454 sequencing data as “trusted contigs” combined with the illumina and nanopore sequencing data as input for the SPAdes algorithm (version 3.11.1, default parameter settings) [34,35]. Subsequent genome polishing was performed by consecutive runs of an in-house Perl script, where the original reads were realigned against the resulting assembly using the Burrows-Wheeler Aligner (BWA). Based on this output, the genome assembly was updated accordingly until no further single-nucleotide polymorphisms (SNPs) and indels were detected.
The C. metallidurans NA4Pt genome was compared to the parental strain at two levels. For the small SNPs and indels, the output of two algorithms was combined: BreSeq 0.32.0 and the Genome Analysis Toolkit (GATK) [36,37]. BreSeq was run using the default parameter settings. Before running GATK, we first converted the raw BWA output to a sorted and index Binary Alignment Map (BAM) file using the view, sort and index command of the SAMtools package version 0.1.18 [38]. SNP prediction was performed on this BAM file by following the pipeline described in Van der Auwera et al. [37], with default parameters, and using ploidy = 1 when running the HaplotyperCaller command. For larger structural variations (insertions and deletion), an in-house developed Python script, specifically focused on the identification of structural variations caused by mobile genetic elements, was used. This program exploits the paired-end information and insert size distributions to predict these variations.

2.7. Transcriptomic Analysis Using RNA-Seq

Gene expression in NA4Pt was compared with the parental strain NA4 under non-selective growth conditions. Three independent C. metallidurans NA4 and NA4Pt cultures were allowed to grow until an OD600 of 0.6 was reached. Each culture was subdivided in 2 mL portions and cells were harvested by centrifugation for 2 min at 10,000× g. Bacterial pellets were flash frozen by immersion into liquid nitrogen and kept frozen at −80 °C at all times. Total RNA was extracted using the Promega SV Total RNA Isolation System kit (Promega, Leiden, The Netherlands). RNA sequencing (directional mRNA library, RiboZero rRNA depletion and 2 × 125 bp paired-end sequencing) was performed by Eurofins genomics (Ebersberg, Germany).

2.8. RNA-Seq Data Analysis

Obtained RNA-Seq reads were aligned using BWA software and the default parameters [39]. Raw counts per gene were calculated based on the latest genome annotation of C. metallidurans NA4, as available on the MaGe platform. Reads were allowed to map 50 bp upstream of the start codon or 50 bp downstream of the stop codon. Reads that were mapped to ribosomal or transfer RNA were removed from the raw count data to prevent bias in detecting differential expression. Differential expression was calculated using the edgeR package (version 3.2.4) [40] in BioConductor (release 3.0, R version 3.1.2), resulting in a fold change value and a corresponding p value corrected for multiple testing for each gene. Genes were found to be differentially expressed if they show an absolute log2 fold change higher than 0.80 and a false discovery rate (FDR) value lower than 0.05.

2.9. Functional Analysis

Homologous genes and synteny groups were computed via the MaGe platform [41]. Homologous genes were based on the bidirectional best hit criterion and a blastP alignment threshold (at least 35% amino-acid identity on 80% of the length of the smallest protein). Synteny, orthologous gene sets that have the same local organization are based on the bidirectional best hit criterion or a blastP alignment threshold (at least 30% amino-acid identity on 80% of the length of the smallest protein), and co-localization (with the maximum number of consecutive genes not involved in a synteny group being five).
Distribution of insertion sequence (IS) elements was determined by identification and annotation of IS elements with ISsaga [42] and manual curation.
MOB typing of plasmids was performed using the mob_typer script, part of the MOB-suite [43]. All analyses were performed using standard parameters. The most recent reference database was downloaded on the 19th of November, 2018.
Plasmids were aligned to pNA4_D by the AliTV perl interface in the AliTV package [44] using standard parameters. Output json files were visualized with the AliTV web service (http://alitvteam.github.io/AliTV/d3/AliTV.html). Subsystem categories were assigned via the online implementation of RASTtk (Rapid Annotation using Subsystem Technology) (http://rast.nmpdr.org/) [45]. Circos plots from genomic data were constructed with Circa 1.2.1 (http://omgenomics.com/circa/).

3. Results and Discussion

3.1. Genome Analysis of C. metallidurans NA4

The C. metallidurans NA4 genome was previously sequenced via the 454 GS-FLX sequencing platform [29] and assembled into 109 contigs. Since all Cupriavidus genomes have a multipartite organization composed of at least a chromosome, a chromid and one or more megaplasmids [27], we tested if combining sequencing platforms could produce completely closed replicon sequences. C. metallidurans NA4, in particular, contains six replicons [29,46]. Complete assembly of multipartite genomes is necessary to fully understand the genomic structure and invaluable to correctly assess genomic events. An additional Illumina and Nanopore sequencing was performed. Integration of the pre-assembled contigs (based on the 454 pyrosquencing data) and the Illumina and Nanopore sequencing data, using the SPAdes platform, resulted in complete assembly and closure of the six C. metallidurans NA4 replicons (Table 1). The long reads from Nanopore sequencing resolved the complex repeat regions, while Illumina sequencing provided more accurate sequencing results.
The size, number of coding sequences (CDSs) and GC content of the C. metallidurans NA4 chromosome and chromid were very similar to that of type strain C. metallidurans CH34 [26] (Figure 1). To compare the different replicons of CH34 and NA4, homologous genes were computed and visualized by Circa (Figure 2). Both the chromosome (CHR1) and chromid (CHR2) were well conserved. However, this comparison revealed an extensive number of homologous genes between plasmid pMOL30 from CH34 and the chromid of NA4 (Figure 2). To zoom in on this, large (>15 genes) synteny groups between pMOL30 and the chromid of NA4 were visualized separately (Figure 3), as well as homologous genes shared between pMOL28 and pMOL30, and pNA4_A and pNA4_B (Figure 4).
The latter indicated that the main large metal resistance clusters on pMOL30 are syntenic with gene clusters on NA4’s chromid instead of NA4’s plasmids. Previously, we showed that pMOL28 and pMOL30 contain large genomic islands that harbor all plasmid-borne genes involved in the response to heavy metals [27,47,48]. Plasmid pMOL30 carries two genomic islands, CMGI-30a and -30b, that convey resistance to cadmium, zinc, cobalt, lead, and mercury, and copper and silver, respectively. A 25-gene cluster within CMGI-30a, containing the czc cluster (related to cadmium, zinc and cobalt resistance) and genes involved in membrane-related functions, and highly conserved (>99% nucleotide similarity) in all C. metallidurans strains [48], was found on NA4’s chromid. Interestingly, this cluster is always flanked by a tyrosine-based site-specific recombinase (TBSSR) associated with a conserved protein of unknown function making up a bipartite module (BIM) [48,49]. Our observation adds evidence to the mobility of this cluster. The cluster on NA4’s chromid related to the pbr (lead resistance) and mer (mercury resistance) cluster (Tn4380) of pMOL30 was flanked at both sides by (remnants of) TBSSRs (only at one side for pMOL30). Three clusters on NA4’s chromid were homologous to pMOL30’s cop cluster (copper resistance) (Figure 3, Supplementary Figure S1). The latter was part of genomic island CMGI-30b of pMOL30 that contains a 33-gene copper-related cluster almost completely induced by Cu2+ (coding for the efflux PI-type ATPase CopF, the heavy metal efflux (HME) resistance nodulation cell-division (RND) system SilCBA, the periplasmic detoxification system CopABCDI, and accessory and membrane-related functions) as well as the nre/ncc cluster. One NA4 cluster comprised almost the complete CMGI30-b island except copV and copT. A second cluster contained the copper-related cluster as well as some gene fragments and was flanked by TBSSRs. The third cluster comprised the cop cluster except copV and copT (Supplementary Figure S1). In addition, a cluster carrying resistance to chromate and cobalt/nickel, similar to that on genomic island CMGI-28a of pMOL28, was also found on NA4’s chromid (next to plasmid pNA4_B). The cluster was delimited by ISRme1 and a TBBSR, and was located between the pbr and mer cluster, and the second cop cluster. It did not contain the additional five genes that are carried by the pMOL28 and pNA4_B cluster [27,50,51].
Although most metal resistance determinants are conserved between C. metallidurans NA4 and CH34, which results in similar growth profiles in the presence of metals [52], they are harbored by different replicons. We showed that, unlike the characteristic megaplasmid pMOL30 in C. metallidurans CH34 that is specialized in heavy metal resistance, pNA4_A does not fulfill this role in C. metallidurans NA4. Plasmid pNA4_B does harbor the pMOL28-like metal resistance determinants and, in addition, a second nccYXHCBAN locus coding for an RND-driven efflux system homologous to that of C. metallidurans 31A and KT02, which mediates resistance to 40 mM Ni2+.

3.1.1. Characterization of (Mega)Plasmids

As our sequencing efforts resulted in the complete closure of the four plasmid replicons, we were able to characterize them in more detail. Numerous proteins are involved in the horizontal transmission of plasmids via conjugation and establishment in the recipient cell. Conjugative plasmids contain an origin of transfer (oriT), a DNA relaxase, a Type IV coupling protein (T4CP), and a membrane-associated mating pair formation (MPF) complex, which is a form of Type IV secretion system (T4SS). Transmissible (mobilizable) plasmids require an oriT and a relaxase that can be provided in trans [53,54,55,56,57]. Different classification or typing schemes for plasmids have been developed, but the principal ones are replicon and MOB typing, relying on plasmid replication and mobility, respectively [55,56,58]. The phylogenetic relationship among relaxases has been thoroughly studied and resulted in the classification of conjugative systems into six MOB families: MOBF, MOBH, MOBQ, MOBC, MOBP, and MOBV [43,56]. The archetype plasmids defining the families are R388 (MOBF), R27 (MOBH), RP4 (MOBP), RSF1010 (MOBQ), pMV158 (MOBV), and CloDF13 (MOBC). We used MOB-suite to type and characterize pNA4_A, pNA4_B, pNA4_C, and pNA4_D, as well as pMOL28 and pMOL30 (Table 2).
Our analysis indicated that pNA4_A and pMOL30, and pNA4_B and pMOL28 were similar based on replicon, relaxase and mating pair formation (MPF) family. Plasmid pNA4_A was classified as being non-mobilizable, whereas pMOL30 as conjugative. Although low-frequency transfer of pMOL30 has been observed, this transfer could be mediated via other conjugative systems. For instance, plasmid RP4 can enhance transfer frequency to 10−3 by cointegrate formation via transposition of Tn4380 [47]. Plasmid pNA4_A carried 416 CDSs, most of them code for unknown proteins (71.6%). Identifiable functions were, next to plasmid replication, maintenance and conjugation, related to metal resistance (partial czc cluster and mercury transposon; see Figure 4) and alkaline phosphatase. Plasmid pNA4_B carried 274 CDSs, most of them also code for unknown proteins (59.5%). Identifiable functions were, next to plasmid replication, maintenance and conjugation, related to metal resistance (cluster carrying resistance to chromate and cobalt/nickel, similar to that on genomic island CMGI-28a of pMOL28; see Figure 4).
Plasmid pNA4_C was classified as non-mobilizable and carried 209 CDSs, most of them also code for unknown proteins (74.6%). No accessory plasmid functions could be identified. Plasmid pNA4_D was classified as conjugative and was very similar to plasmids from Acidovorax carolinensis P3 (plasmid pACP3.3) and P4 (plasmid pACP4.4), Acidovorax sp. JS42 (plasmid pAOVO01), Alicycliphilus denitrificans K601 (plasmid pALIDE201), and Pandoraea pnomenusa MCB032 (unnamed plasmid) (Supplementary Figure S2). Characterized proteins encoded by these plasmids were mainly related to conjugational transfer and replication.

3.1.2. Insertion Sequence Elements Distribution

Insertion sequences (IS) are simple mobile genetic elements that play an important role in genome plasticity and activity of these IS elements are often correlated with the adaptive potential to promote genetic variability under different environmental challenges [59,60]. An initial assessment of the number and identity of IS elements in NA4 was previously performed based on the draft genome assembly [52]. However, this could lead to an underestimation of the number of IS elements because possible identical IS elements will only be represented as one contig [61]. Therefore, we reanalyzed and determined the correct number of IS elements in NA4 (Table 3). In total, 21 intact IS elements were identified. C. metallidurans NA4 carried much less IS elements in comparison with type strain C. metallidurans CH34, which carried 57 intact IS elements.

3.2. Analysis of C. metallidurans NA4Pt

3.2.1. Determination of Minimal Inhibitory Concentration and Generation of a Pt4+ Resistant Mutant

We used C. metallidurans NA4, which is able to survive in oligotrophic conditions for many months [46], to scrutinize adaptation to Pt4+. The minimal inhibitory concentration (MIC) of Pt4+ for C. metallidurans NA4 in Tris-buffered mineral medium was 70 µM (Table 1), which was similar to that of C. metallidurans CH34. This already indicated that NA4 has a high level of resistance to Pt4+ when compared to other strains such as Klebsiella pneumoniae (20 µM), Acinetobacter baumannii (30 µM), and Enterococcus faecium (60 µM) [63]. Furthermore, the MIC determinations in that study were performed in rich broth medium, which could affect the metal bioavailability and lead to an overestimation of the MIC [63]. Next, C. metallidurans NA4 was exposed to a subinhibitory concentration of 62.5 µM Pt4+ during 30 days (eight serial passages). After passage on non-selective medium, a mutant (designated NA4Pt) that displayed a higher resistance to Pt4+ (MIC of 160 µM) was obtained. No differences were observed between NA4 and NA4Pt when grown in non-selective conditions (Figure 5). In addition, no differences in MIC of Pd2+, Ag+, Zn2+, Ni2+, and Cu2+ were observed (Table 4). From this data it is evident that the NA4Pt mutant is specifically resistant to Pt4+ and does not have a higher resistance to Pd2+, in contrast to what was described for Desulfovibrio sp. preference during biosorption [14].

3.2.2. Sequence Analysis of C. metallidurans NA4Pt

The laboratory-evolved mutant NA4Pt was sequenced to identify genomic changes such as insertions, deletions and SNPs (Supplementary Table S1). These genomic changes were not observed in other adaptive laboratory evolution experiments with NA4 (unpublished results). In the chromosome, an insertion (+56 bp) in the upstream region of a metal-dependent hydrolase and a point mutation in pseudouridine synthase (rluB) (resulting in R64H substitution) were observed. In the chromid, several point mutations (mostly synonymous mutations) and two deletions (a single bp and a 119-bp region), both located in the copB gene coding for an outer membrane protein involved in copper resistance, were found. The latter did not affect copper resistance of NA4Pt (Table 4). The biggest change, a large deletion of 69.9 kbp, was observed in plasmid pNA4_D. The 19.7 kbp remaining fragment (positions 52,551 to 72,345) included genes encoding a DNA primase, a C-5 cytosine-specific DNA methylase, a single-stranded DNA-binding (ssb) protein, and the replication initiator protein RepA (Supplementary Table S2), which is responsible for plasmid replication in bacteria [64,65]. This suggested that the remaining part could be maintained as a smaller plasmid instead of being integrated in one of the other replicons. The absence of the native pNA4_D and presence of the smaller plasmid were confirmed by plasmid DNA extraction and PagI digestion (Figure 6).
The large deletion, which could probably be mediated by the presence of a 221 bp direct repeat (Supplementary Figure S3), resulted in the loss of 89 CDSs. (Mega)plasmids are nonessential and dispensable for cell viability in most environments [66], and large deletions have also been observed in other laboratory-evolved strains. For instance, a large deletion occurred in the megaplasmid pAtC58 from Agrobacterium tumefaciens laboratory-evolved strains, resulting in increased virulence gene expression and reduced fitness cost [67]. Prolonged cultivation of the Gram-positive actinobacterium Rhodococcus opacus 1CP containing megaplasmid p1CP (740 kb) under non-selective conditions led to the isolation of mutants 1CP.01 and 1CP.02, harboring the shortened plasmid variants p1CP.01 (500 kb) and p1CP.02 (400 kb) [68]. Methylobacterium extorquens AM1 lost 10% of its megaplasmid in an evolution experiment, which were beneficial in the applied selective environment, but disadvantageous in alternative environments [69]. Gene loss is a very common evolutionary process in bacteria and provides an increased fitness under one or several growth conditions [70]. The large deletion in pNA4_D resulted in loss of the conjugative machinery (T4SS), which is known to impose a burden [71,72]. For instance, the growth rate of an experimentally evolved E. coli increased by IS26-mediated loss of the T4SS on its plasmid pKP33 [73]. Several mechanisms are put forward to explain the reduced fitness costs, such as ribosome occupancy, reduction of energy demands for DNA replication, transcription and translation, and negative interactions between chromosomal pathways and (mega)plasmid-encoded proteins [66,67,74]. In addition, selective processes favoring adaptation to specific stressors/environments can be a driving force behind gene loss [66,70,75,76]. Therefore, although no growth differences were observed between NA4 and NA4Pt, this loss could have an effect in challenging environments with increased selection pressure (e.g., high platinum concentration).

3.2.3. Transcriptome Analysis

Resistance is the result of natural selection for resistance-conferring mutations (i.e., random mutations that allow growth under selection, outcompeting the parent, and subsequent isolation of the adapted mutant) [77]. Therefore, and similar to other studies [78,79,80], the global shift in transcriptome, resulting from the altered genotype of the evolved strain (NA4Pt) as compared with the parental strain (NA4), was examined by RNA-Seq in non-selective conditions (average total number of reads was 4,324,281 ± 463,666). Up- and down-regulated genes were selected based on log2 fold change (<−0.8 and >0.8) and significance (p < 0.05), which resulted in 111 up- and 115 down-regulated genes (Figure 7 and Supplementary Tables S2 and S3).
The functional relevance of differentially-expressed genes was explored by using functional categories from the eggNOG classification system [81]. Genes of different categories were found to be differently expressed (Figure 8). However, none of the categories were significantly over-represented (Fisher’s exact test).
Overall, genes involved in defense mechanisms, intracellular trafficking, signal transduction and membrane-related genes were more up-regulated than down-regulated. On the other hand, genes involved in carbohydrate transport, nucleotide transport, cell motility and cell-cycle-related genes were more down-regulated than up-regulated. Up-regulated defense mechanism genes included mainly RND-driven efflux systems, which are abundant in the NA4 genome [52]. Next to systems putatively involved in the efflux of chemicals (acridines), the expression of two HME-RND genes was increased. The cusC gene is part of the CusCBA efflux pump responsible for copper and silver resistance [82]. In C. metallidurans, cusDCBAF genes are also up-regulated by silver ions [28], and induction of the CusC protein synthesis was also observed in the presence of 1 µM silver or 0.85 mM copper ions [83]. The cnrC gene is part of cnrCBAT efflux system that mediates nickel and cobalt resistance [84]. However, the resistance of NA4Pt to Cu2+, Zn2+, Ag+ and Ni2+ was unaffected compared to that of the parental strain NA4 (Table 4). No other metal resistance genes were up- or down-regulated.
Up-regulated genes of the signal transducing pathways belong to different two-component-system-related proteins (Supplementary Table S2). Namely, ompR-family regulators, members of the largest response regulator family involved in many signal transduction processes [85], and cheY, which is involved in chemotaxis and modulates motility by regulating flagellar motor switch proteins [86]. In contrast to cheY, cell-motility-related genes were down-regulated (flgC, flgD, flgG, pilW, and pilX). FlgC and FlgG proteins form the rod part of the flagellar basal body [87] and FlgG polymerizes to form the distal rod on top of the proximal rod, acting as a hook cap [88]. PilX and PilW are involved in biogenesis of Type IV fimbriae, which are surface filaments mediating attachment to host epithelial cells and flagella-independent twitching motility [89]. These differences did not affect motility as no significant differences were observed between NA4Pt and NA4 in motility assays (Supplementary Figure S4).
The class of membrane-related up-regulated genes contains a highly up-regulated (36.3-fold) lytic transglycosylase (LT), which usually plays an important role in shaping the periplasmic space and is tightly regulated, as over activity can have deleterious effects [90]. LT is responsible for creating space within the peptidoglycan layer for cell division and the insertion of cell-envelope spanning structures such as flagella and secretion systems [90]. In C. metallidurans, different LTs are induced in the presence of Cu2+, Cd2+, Pb2+ and Zn2+ [91]. Therefore, it may be hypothesized that the LTs might play a role in coping with toxic metal ion concentrations. Next to LTs, membrane-bound cytochrome c, which enables electron transfer as well as the catalysis of various redox reactions [92], is also up-regulated. In D. vulgaris, this periplasmic protein is directly involved in H2-mediated metal reduction [16]. Desulfovibrio sp. cytochrome c3 mediates electron transfer to palladium and platinum complexes, thereby reducing them to zero-valent nanoparticles [10,17]. The latter supports a role of this overexpressed gene in the Pt-resistant phenotype. Furthermore, nm-scale colloidal platinum was found in C. metallidurans primarily along the cell envelope, where energy generation/electron transport occurs [12].
To scrutinize the impact of differentially-expressed membrane-related genes (LT, flagellar- and pili-related) on morphology, scanning electron microscopy (SEM) was performed and showed that NA4Pt cells appeared to be more elongated (Figure 9). This correlates with the induction of filaments during platinum metal stress [19]. Nevertheless, the presence of a subpopulation of elongated cells could not be confirmed by flow cytometry (Supplementary Figure S5).
Up-regulated intracellular trafficking genes were mainly related to a Type II secretion system (T2SS), which secretes proteins to the extracellular environment [93]. In C. metallidurans, the T2SS secretes alkaline phosphatase in the extracellular environment [94].
Finally, none of the CDSs affected by the genomic changes in NA4Pt (either directly or via alteration in their promoter region) were differentially expressed (Supplementary Tables S1–S3). The latter points towards a significant role of the large deletion in pNA4_D on the overall expression pattern, as well as adaptation to platinum.

4. Conclusions

We resolved the multireplicon genome of C. metallidurans NA4 by combining Illumina and Nanopore sequencing and; thereby, paved the way for genetic, genomic and evolution-related studies. Further analysis of the replicons showed a distinctive pattern in the location of metal resistance determinants, with the chromid playing a pivotal role in contrast to type strain CH34, for which the megaplasmids pMOL28 and pMOL30 are the main actors. An NA4 derivative (NA4Pt) that showed increased resistance to platinum was generated via an adaptive laboratory evolution experiment. On the basis of our observations, one might speculate that the increased resistance to Pt4+ in C. metallidurans NA4Pt is not mediated by a dedicated mechanism, but by pleiotropic alterations in membrane-related processes, such as pili, peptidoglycan turnover and electron transfer, probably elicited by the large deletion in its 98-kb plasmid.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/10/1/63/s1, Figure S1: C. metallidurans chromid clusters homologous to the cop cluster on pMOL30; Figure S2: Whole plasmid alignments of pNA4_D from C. metallidurans NA4 and pACP3.3 and pACP4.4 from A. carolinensis P3 and P4, pAOVO01 from Acidovorax sp. JS42, pALIDE201 from A. denitrificans K601, and an unnamed plasmid from P. pnomenusa MCB032; Figure S3: Alignment of 221 bp direct repeat flanking large deletion in pNA4_D; Figure S4: Motility assay in soft (0.3%) agar. Diameter measurement of the haloes formed by C. metallidurans NA4 and NA4Pt after 24 h; Figure S5: Flow cytometry histogram of C. metallidurans NA4 and NA4Pt based on forward-scattered light (FSC) and side-scattered light (SSC); Table S1: List of genomic changes in the laboratory-evolved C. metallidurans NA4Pt; Table S2: RNA-Seq-derived up-regulated genes (log2 fold > 0.8; p < 0.05) in C. metallidurans NA4Pt, compared to its parental strain in non-selective conditions; Table S3: RNA-Seq-derived down-regulated genes (log2 fold < −0.8; p < 0.05) in C. metallidurans NA4Pt, compared to its parental strain in non-selective conditions.

Author Contributions

Conceptualization, R.V.H.; methodology, M.M.A., A.P., and R.V.H.; software, M.M.A, L.M., P.M., and R.V.H.; validation, M.M.A., P.M., D.C., and R.V.H.; formal analysis, M.M. A., L.M., A.P., P.M., and R.V.H.; investigation, M.M.A., A.P., and R.V.H.; resources, M.M.A. and R.V.H.; writing—original draft preparation, M.M.A. and R.V.H.; writing—review and editing, M.M.A., P.M., D.C., and R.V.H.; visualization, L.M. and R.V.H.; supervision, N.L., D.C., and R.V.H.; project administration, N.L. and R.V.H.; funding acquisition, N.L. and R.V.H.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Rauch, S.; Morrison, G.M. Environmental relevance of the platinum-group elements. Elements 2008, 4, 259–263. [Google Scholar] [CrossRef]
  2. Wedepohl, K.H. The Composition of the Continental-Crust. Geochim. Cosmochim. Acta 1995, 59, 1217–1232. [Google Scholar] [CrossRef]
  3. Ravindra, K.; Bencs, L.; Van Grieken, R. Platinum group elements in the environment and their health risk. Sci. Total Environ. 2004. [Google Scholar] [CrossRef]
  4. Florea, A.M.; Busselberg, D. Cisplatin as an anti-tumor drug: Cellular mechanisms of activity, drug resistance and induced side effects. Cancers 2011, 3, 1351–1371. [Google Scholar] [CrossRef] [PubMed]
  5. Maes, S.; Props, R.; Fitts, J.P.; Smet, R.D.; Vilchez-Vargas, R.; Vital, M.; Pieper, D.H.; Vanhaecke, F.; Boon, N.; Hennebel, T. Platinum Recovery from Synthetic Extreme Environments by Halophilic Bacteria. Environ. Sci. Technol. 2016, 50, 2619–2626. [Google Scholar] [CrossRef] [PubMed]
  6. Hennebel, T.; Boon, N.; Maes, S.; Lenz, M. Biotechnologies for critical raw material recovery from primary and secondary sources: R&D priorities and future perspectives. New Biotechnol. 2015, 32, 121–127. [Google Scholar] [CrossRef]
  7. Reith, F.; Zammit, C.M.; Shar, S.S.; Etschmann, B.; Bottrill, R.; Southam, G.; Ta, C.; Kilburn, M.; Oberthür, T.; Ball, A.S.; et al. Biological role in the transformation of platinum-group mineral grains. Nat. Geosci. 2016, 9, 294–298. [Google Scholar] [CrossRef]
  8. Zereini, F.; Wiseman, C.L.S. Platinum Metals in the Environment; Springer-Verlag: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  9. Maes, S.; Claus, M.; Verbeken, K.; Wallaert, E.; De Smet, R.; Vanhaecke, F.; Boon, N.; Hennebel, T. Platinum recovery from industrial process streams by halophilic bacteria: Influence of salt species and platinum speciation. Water Res. 2016, 105, 436–443. [Google Scholar] [CrossRef]
  10. de Vargas, I.; Macaskie, L.E.; Guibal, E. Biosorption of palladium and platinum by sulfate-reducing bacteria. J. Chem. Technol. Biotechnol. 2004, 79, 49–56. [Google Scholar] [CrossRef]
  11. Maes, S.; Props, R.; Fitts, J.P.; De Smet, R.; Vanhaecke, F.; Boon, N.; Hennebel, T. Biological Recovery of Platinum Complexes from Diluted Aqueous Streams by Axenic Cultures. PLoS ONE 2017. [Google Scholar] [CrossRef]
  12. Campbell, G.; MacLean, L.; Reith, F.; Brewe, D.; Gordon, R.A.; Southam, G. Immobilisation of Platinum by Cupriavidus metallidurans. Minerals 2018, 8, 10. [Google Scholar] [CrossRef]
  13. Konishi, Y.; Ohno, K.; Saitoh, N.; Nomura, T.; Nagamine, S.; Hishida, H.; Takahashi, Y.; Uruga, T. Bioreductive deposition of platinum nanoparticles on the bacterium Shewanella algae. J. Biotechnol. 2007, 128, 648–653. [Google Scholar] [CrossRef] [PubMed]
  14. Yong, P.; Rowson, N.A.; Farr, J.P.G.; Harris, I.R.; Macaskie, L.E. Bioaccumulation of palladium by Desulfovibrio desulfuricans. J. Chem. Technol. Biotechnol. 2002, 77, 593–601. [Google Scholar] [CrossRef]
  15. Rapp-Giles, B.J.; Casalot, L.; English, R.S.; Ringbauer, J.A.; Dolla, A.; Wall, J.D. Cytochrome c3 Mutants of Desulfovibrio desulfuricans. Appl. Environ. Microbiol. 2000, 66, 671–677. [Google Scholar] [CrossRef] [PubMed]
  16. Elias, D.A.; Suflita, J.M.; McInerney, M.J.; Krumholz, L.R. Periplasmic cytochrome c3 of Desulfovibrio vulgaris is directly involved in H2-mediated metal but not sulfate reduction. Appl. Environ. Microbiol. 2004, 70, 413–420. [Google Scholar] [CrossRef] [PubMed]
  17. Capeness, M.J.; Edmundson, M.C.; Horsfall, L.E. Nickel and platinum group metal nanoparticle production by Desulfovibrio alaskensis G20. New Biotechnol. 2015, 32, 727–731. [Google Scholar] [CrossRef] [PubMed]
  18. Gauthier, D.; Sobjerg, L.S.; Jensen, K.M.; Lindhardt, A.T.; Bunge, M.; Finster, K.; Meyer, R.L.; Skrydstrup, T. Environmentally benign recovery and reactivation of palladium from industrial waste by using gram-negative bacteria. ChemSusChem 2010, 3, 1036–1039. [Google Scholar] [CrossRef]
  19. Rosenberg, B.; Renshaw, E.; Vancamp, L.; Hartwick, J.; Drobnik, J. Platinum-induced filamentous growth in Escherichia coli. J. Bacteriol. 1967, 93, 716–721. [Google Scholar]
  20. Moore, R.L.; Brubaker, R.R. Effect of cis-platinum(II)diamminodichloride on cell division of Hyphomicrobium and Caulobacter. J. Bacteriol. 1976, 125, 317–323. [Google Scholar]
  21. Beck, D.J.; Brubaker, R.R. Effect of cis-platinum(II)diamminodichloride on wild type and deoxyribonucleic acid repair deficient mutants of Escherichia coli. J. Bacteriol. 1973, 116, 1247–1252. [Google Scholar]
  22. Drobnik, J.; Urbankova, M.; Krekulova, A. The effect of cis-dichlorodiammineplatinum(II) on Escherichia coli B. The role of fil, exr and hcr markers. Mutat. Res. 1973, 17, 13–20. [Google Scholar] [CrossRef]
  23. Reslova, S. The induction of lysogenic strains of Escherichia coli by cis-dichlorodiammineplatinum (II). Chem. Biol. Interact. 1971, 4, 66–70. [Google Scholar] [CrossRef]
  24. Beck, D.J.; Fisch, J.E. Mutagenicity of platinum coordination complexes in Salmonella typhimurium. Mutat. Res. 1980, 77, 45–54. [Google Scholar] [CrossRef]
  25. Maboeta, M.S.; Claassens, S.; Van Rensburg, L.; Van Rensburg, P.J.J. The effects of platinum mining on the environment from a soil microbial perspective. Water Air Soil Pollut. 2006, 175, 149–161. [Google Scholar] [CrossRef]
  26. Janssen, P.J.; Van Houdt, R.; Moors, H.; Monsieurs, P.; Morin, N.; Michaux, A.; Benotmane, M.A.; Leys, N.; Vallaeys, T.; Lapidus, A.; et al. The complete genome sequence of Cupriavidus metallidurans strain CH34, a master survivalist in harsh and anthropogenic environments. PLoS ONE 2010, 5, e10433. [Google Scholar] [CrossRef] [PubMed]
  27. Mergeay, M.; Van Houdt, R. Metal Response in Cupriavidus metallidurans, Volume I: FROM Habitats to Genes and Proteins; Springer International Publishing: Basel, Switzerland, 2015; p. 89. [Google Scholar]
  28. Monsieurs, P.; Moors, H.; Van Houdt, R.; Janssen, P.J.; Janssen, A.; Coninx, I.; Mergeay, M.; Leys, N. Heavy metal resistance in Cupriavidus metallidurans CH34 is governed by an intricate transcriptional network. Biometals 2011, 24, 1133–1151. [Google Scholar] [CrossRef]
  29. Monsieurs, P.; Mijnendonckx, K.; Provoost, A.; Venkateswaran, K.; Ott, C.M.; Leys, N.; Van Houdt, R. Genome Sequences of Cupriavidus metallidurans Strains NA1, NA4, and NE12, Isolated from Space Equipment. Genome Announc. 2014. [Google Scholar] [CrossRef]
  30. Mergeay, M.; Nies, D.; Schlegel, H.G.; Gerits, J.; Charles, P.; Van Gijsegem, F. Alcaligenes eutrophus CH34 is a facultative chemolithotroph with plasmid-bound resistance to heavy metals. J. Bacteriol. 1985, 162, 328–334. [Google Scholar]
  31. Wiegand, I.; Hilpert, K.; Hancock, R.E. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat. Protoc. 2008, 3, 163–175. [Google Scholar] [CrossRef]
  32. Randall, C.P.; Oyama, L.B.; Bostock, J.M.; Chopra, I.; O’Neill, A.J. The silver cation (Ag+): Antistaphylococcal activity, mode of action and resistance studies. J. Antimicrob. Chemother. 2013, 68, 131–138. [Google Scholar] [CrossRef]
  33. Andrup, L.; Barfod, K.K.; Jensen, G.B.; Smidt, L. Detection of large plasmids from the Bacillus cereus group. Plasmid 2008, 59, 139–143. [Google Scholar] [CrossRef] [PubMed]
  34. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed]
  35. Antipov, D.; Korobeynikov, A.; McLean, J.S.; Pevzner, P.A. hybridSPAdes: An algorithm for hybrid assembly of short and long reads. Bioinformatics 2016, 32, 1009–1015. [Google Scholar] [CrossRef] [PubMed]
  36. Deatherage, D.E.; Barrick, J.E. Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods Mol. Biol. 2014, 1151, 165–188. [Google Scholar] [CrossRef] [PubMed]
  37. Van der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ data to high confidence variant calls: The Genome Analysis Toolkit best practices pipeline. Curr. Protoc. Bioinform. 2013. [Google Scholar] [CrossRef]
  38. 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. [Google Scholar] [CrossRef] [Green Version]
  39. Li, H.; Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 2010, 26, 589–595. [Google Scholar] [CrossRef]
  40. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef]
  41. Vallenet, D.; Belda, E.; Calteau, A.; Cruveiller, S.; Engelen, S.; Lajus, A.; Le Fevre, F.; Longin, C.; Mornico, D.; Roche, D.; et al. MicroScope—An integrated microbial resource for the curation and comparative analysis of genomic and metabolic data. Nucleic Acids Res. 2013, 41, D636–D647. [Google Scholar] [CrossRef]
  42. Varani, A.M.; Siguier, P.; Gourbeyre, E.; Charneau, V.; Chandler, M. ISsaga is an ensemble of web-based methods for high throughput identification and semi-automatic annotation of insertion sequences in prokaryotic genomes. Genome Biol. 2011, 12, R30. [Google Scholar] [CrossRef]
  43. Robertson, J.; Nash, J.H.E. MOB-suite: Software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb. Genom. 2018. [Google Scholar] [CrossRef] [PubMed]
  44. Ankenbrand, M.J.; Hohlfeld, S.; Hackl, T.; Forster, F. AliTV-interactive visualization of whole genome comparisons. PEERJ Comput. Sci. 2017. [Google Scholar] [CrossRef]
  45. Brettin, T.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Olsen, G.J.; Olson, R.; Overbeek, R.; Parrello, B.; Pusch, G.D.; et al. RASTtk: A modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 2015, 5, 8365. [Google Scholar] [CrossRef] [PubMed]
  46. Mijnendonckx, K.; Provoost, A.; Ott, C.M.; Venkateswaran, K.; Mahillon, J.; Leys, N.; Van Houdt, R. Characterization of the survival ability of Cupriavidus metallidurans and Ralstonia pickettii from space-related environments. Microb. Ecol. 2013, 65, 347–360. [Google Scholar] [CrossRef] [PubMed]
  47. Mergeay, M.; Monchy, S.; Janssen, P.; Houdt, R.V.; Leys, N. Megaplasmids in Cupriavidus Genus and Metal Resistance. In Microbial Megaplasmids; Schwartz, E., Ed.; Springer Berlin Heidelberg: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  48. Van Houdt, R.; Monsieurs, P.; Mijnendonckx, K.; Provoost, A.; Janssen, A.; Mergeay, M.; Leys, N. Variation in genomic islands contribute to genome plasticity in Cupriavidus metallidurans. BMC Genom. 2012, 13, 111. [Google Scholar] [CrossRef] [PubMed]
  49. Van Houdt, R.; Monchy, S.; Leys, N.; Mergeay, M. New mobile genetic elements in Cupriavidus metallidurans CH34, their possible roles and occurrence in other bacteria. Antonie Van Leeuwenhoek 2009, 96, 205–226. [Google Scholar] [CrossRef] [PubMed]
  50. Henne, K.L.; Nakatsu, C.H.; Thompson, D.K.; Konopka, A.E. High-level chromate resistance in Arthrobacter sp. strain FB24 requires previously uncharacterized accessory genes. BMC Microbiol. 2009, 9, 199. [Google Scholar] [CrossRef]
  51. Henne, K.L.; Turse, J.E.; Nicora, C.D.; Lipton, M.S.; Tollaksen, S.L.; Lindberg, C.; Babnigg, G.; Giometti, C.S.; Nakatsu, C.H.; Thompson, D.K.; et al. Global proteomic analysis of the chromate response in Arthrobacter sp. strain FB24. J. Proteome Res. 2009, 8, 1704–1716. [Google Scholar] [CrossRef]
  52. Van Houdt, R.; Provoost, A.; Van Assche, A.; Leys, N.; Lievens, B.; Mijnendonckx, K.; Monsieurs, P. Cupriavidus metallidurans Strains with Different Mobilomes and from Distinct Environments Have Comparable Phenomes. Genes 2018, 9, 507. [Google Scholar] [CrossRef]
  53. Shintani, M.; Sanchez, Z.K.; Kimbara, K. Genomics of microbial plasmids: Classification and identification based on replication and transfer systems and host taxonomy. Front. Microbiol. 2015, 6, 242. [Google Scholar] [CrossRef]
  54. Ramsay, J.P.; Kwong, S.M.; Murphy, R.J.; Yui Eto, K.; Price, K.J.; Nguyen, Q.T.; O’Brien, F.G.; Grubb, W.B.; Coombs, G.W.; Firth, N. An updated view of plasmid conjugation and mobilization in Staphylococcus. Mob. Genet. Elements 2016, 6, e1208317. [Google Scholar] [CrossRef] [PubMed]
  55. Carattoli, A.; Bertini, A.; Villa, L.; Falbo, V.; Hopkins, K.L.; Threlfall, E.J. Identification of plasmids by PCR-based replicon typing. J. Microbiol. Methods 2005, 63, 219–228. [Google Scholar] [CrossRef]
  56. Garcillan-Barcia, M.P.; Francia, M.V.; de la Cruz, F. The diversity of conjugative relaxases and its application in plasmid classification. FEMS Microbiol. Rev. 2009, 33, 657–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Ilangovan, A.; Connery, S.; Waksman, G. Structural biology of the Gram-negative bacterial conjugation systems. Trends Microbiol. 2015, 23, 301–310. [Google Scholar] [CrossRef]
  58. Orlek, A.; Stoesser, N.; Anjum, M.F.; Doumith, M.; Ellington, M.J.; Peto, T.; Crook, D.; Woodford, N.; Walker, A.S.; Phan, H.; et al. Plasmid Classification in an Era of Whole-Genome Sequencing: Application in Studies of Antibiotic Resistance Epidemiology. Front. Microbiol. 2017, 8, 182. [Google Scholar] [CrossRef] [PubMed]
  59. Vandecraen, J.; Chandler, M.; Aertsen, A.; Van Houdt, R. The impact of insertion sequences on bacterial genome plasticity and adaptability. Crit. Rev. Microbiol. 2017, 43, 709–730. [Google Scholar] [CrossRef]
  60. Siguier, P.; Gourbeyre, E.; Chandler, M. Bacterial insertion sequences: Their genomic impact and diversity. FEMS Microbiol. Rev. 2014, 38, 865–891. [Google Scholar] [CrossRef]
  61. Ricker, N.; Qian, H.; Fulthorpe, R.R. The limitations of draft assemblies for understanding prokaryotic adaptation and evolution. Genomics 2012, 100, 167–175. [Google Scholar] [CrossRef] [Green Version]
  62. Siguier, P.; Perochon, J.; Lestrade, L.; Mahillon, J.; Chandler, M. ISfinder: The reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34, D32–D36. [Google Scholar] [CrossRef]
  63. Vaidya, M.Y.; McBain, A.J.; Butler, J.A.; Banks, C.E.; Whitehead, K.A. Antimicrobial Efficacy and Synergy of Metal Ions against Enterococcus faecium, Klebsiella pneumoniae and Acinetobacter baumannii in Planktonic and Biofilm Phenotypes. Sci. Rep. 2017, 7, 5911. [Google Scholar] [CrossRef]
  64. Chattoraj, D.K.; Snyder, K.M.; Abeles, A.L. P1 Plasmid Replication—Multiple Functions of Repa Protein at the Origin. Proc. Natl. Acad. Sci. USA 1985, 82, 2588–2592. [Google Scholar] [CrossRef] [PubMed]
  65. Dı́az-López, T.; Lages-Gonzalo, M.; Serrano-López, A.; Alfonso, C.; Rivas, G.; Dı́az-Orejas, R.; Giraldo, R. Structural changes in RepA, a plasmid replication initiator, upon binding to origin DNA. J. Biol. Chem. 2003, 278, 18606–18616. [Google Scholar] [CrossRef] [PubMed]
  66. diCenzo, G.C.; Finan, T.M. The Divided Bacterial Genome: Structure, Function, and Evolution. Microbiol. Mol. Biol. Rev. 2017, 81, e00019-17. [Google Scholar] [CrossRef]
  67. Morton, E.R.; Merritt, P.M.; Bever, J.D.; Fuqua, C. Large Deletions in the pAtC58 Megaplasmid of Agrobacterium tumefaciens Can Confer Reduced Carriage Cost and Increased Expression of Virulence Genes. Genome Biol. Evol. 2013, 5, 1353–1364. [Google Scholar] [CrossRef] [PubMed]
  68. Konig, C.; Eulberg, D.; Groning, J.; Lakner, S.; Seibert, V.; Kaschabek, S.R.; Schlomann, M. A linear megaplasmid, p1CP, carrying the genes for chlorocatechol catabolism of Rhodococcus opacus 1CP. Microbiology 2004, 150, 3075–3087. [Google Scholar] [CrossRef] [PubMed]
  69. Lee, M.C.; Marx, C.J. Repeated, selection-driven genome reduction of accessory genes in experimental populations. PLoS Genet. 2012, 8, e1002651. [Google Scholar] [CrossRef]
  70. Koskiniemi, S.; Sun, S.; Berg, O.G.; Andersson, D.I. Selection-driven gene loss in bacteria. PLoS Genet. 2012, 8, e1002787. [Google Scholar] [CrossRef]
  71. Fernandez-Lopez, R.; Del Campo, I.; Revilla, C.; Cuevas, A.; de la Cruz, F. Negative feedback and transcriptional overshooting in a regulatory network for horizontal gene transfer. PLoS Genet. 2014, 10, e1004171. [Google Scholar] [CrossRef]
  72. Zahrl, D.; Wagner, M.; Bischof, K.; Koraimann, G. Expression and assembly of a functional type IV secretion system elicit extracytoplasmic and cytoplasmic stress responses in Escherichia coli. J. Bacteriol. 2006, 188, 6611–6621. [Google Scholar] [CrossRef]
  73. Porse, A.; Schonning, K.; Munck, C.; Sommer, M.O. Survival and Evolution of a Large Multidrug Resistance Plasmid in New Clinical Bacterial Hosts. Mol. Biol. Evol. 2016, 33, 2860–2873. [Google Scholar] [CrossRef]
  74. Romanchuk, A.; Jones, C.D.; Karkare, K.; Moore, A.; Smith, B.A.; Jones, C.; Dougherty, K.; Baltrus, D.A. Bigger is not always better: Transmission and fitness burden of approximately 1MB Pseudomonas syringae megaplasmid pMPPla107. Plasmid 2014, 73, 16–25. [Google Scholar] [CrossRef] [PubMed]
  75. Moran, N.A. Microbial minimalism: Genome reduction in bacterial pathogens. Cell 2002, 108, 583–586. [Google Scholar] [CrossRef]
  76. Dufresne, A.; Garczarek, L.; Partensky, F. Accelerated evolution associated with genome reduction in a free-living prokaryote. Genome Biol. 2005, 6, R14. [Google Scholar] [CrossRef] [PubMed]
  77. Lenski, R.E. What is adaptation by natural selection? Perspectives of an experimental microbiologist. PLoS Genet. 2017, 13, e1006668. [Google Scholar] [CrossRef] [PubMed]
  78. LaCroix, R.A.; Sandberg, T.E.; O’Brien, E.J.; Utrilla, J.; Ebrahim, A.; Guzman, G.I.; Szubin, R.; Palsson, B.O.; Feist, A.M. Use of Adaptive Laboratory Evolution To Discover Key Mutations Enabling Rapid Growth of Escherichia coli K-12 MG1655 on Glucose Minimal Medium. Appl. Environ. Microbiol. 2015, 81, 17–30. [Google Scholar] [CrossRef] [PubMed]
  79. McCloskey, D.; Xu, S.; Sandberg, T.E.; Brunk, E.; Hefner, Y.; Szubin, R.; Feist, A.M.; Palsson, B.O. Growth Adaptation of gnd and sdhCB Escherichia coli Deletion Strains Diverges From a Similar Initial Perturbation of the Transcriptome. Front. Microbiol. 2018, 9, 1793. [Google Scholar] [CrossRef] [PubMed]
  80. Sandberg, T.E.; Pedersen, M.; LaCroix, R.A.; Ebrahim, A.; Bonde, M.; Herrgard, M.J.; Palsson, B.O.; Sommer, M.; Feist, A.M. Evolution of Escherichia coli to 42 degrees C and subsequent genetic engineering reveals adaptive mechanisms and novel mutations. Mol. Biol. Evol. 2014, 31, 2647–2662. [Google Scholar] [CrossRef]
  81. Powell, S.; Szklarczyk, D.; Trachana, K.; Roth, A.; Kuhn, M.; Muller, J.; Arnold, R.; Rattei, T.; Letunic, I.; Doerks, T.; et al. eggNOG v3.0: Orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic Acids Res. 2012, 40, D284–D289. [Google Scholar] [CrossRef]
  82. Franke, S.; Grass, G.; Rensing, C.; Nies, D.H. Molecular analysis of the copper-transporting efflux system CusCFBA of Escherichia coli. J. Bacteriol. 2003, 185, 3804–3812. [Google Scholar] [CrossRef]
  83. Mergeay, M.; Monchy, S.; Vallaeys, T.; Auquier, V.; Benotmane, A.; Bertin, P.; Taghavi, S.; Dunn, J.; van der Lelie, D.; Wattiez, R. Ralstonia metallidurans, a bacterium specifically adapted to toxic metals: Towards a catalogue of metal-responsive genes. FEMS Microbiol. Rev. 2003, 27, 385–410. [Google Scholar] [CrossRef]
  84. Liesegang, H.; Lemke, K.; Siddiqui, R.A.; Schlegel, H.G. Characterization of the inducible nickel and cobalt resistance determinant cnr from pMOL28 of Alcaligenes eutrophus CH34. J. Bacteriol. 1993, 175, 767–778. [Google Scholar] [CrossRef] [PubMed]
  85. Nguyen, M.P.; Yoon, J.M.; Cho, M.H.; Lee, S.W. Prokaryotic 2-component systems and the OmpR/PhoB superfamily. Can. J. Microbiol. 2015, 61, 799–810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Szurmant, H.; Ordal, G.W. Diversity in chemotaxis mechanisms among the bacteria and archaea. Microbiol. Mol. Biol. Rev. 2004, 68, 301–319. [Google Scholar] [CrossRef] [PubMed]
  87. Fujii, T.; Kato, T.; Hiraoka, K.D.; Miyata, T.; Minamino, T.; Chevance, F.F.; Hughes, K.T.; Namba, K. Identical folds used for distinct mechanical functions of the bacterial flagellar rod and hook. Nat. Commun. 2017, 8, 14276. [Google Scholar] [CrossRef] [PubMed]
  88. Cohen, E.J.; Hughes, K.T. Rod-to-hook transition for extracellular flagellum assembly is catalyzed by the L-ring-dependent rod scaffold removal. J. Bacteriol. 2014, 196, 2387–2395. [Google Scholar] [CrossRef]
  89. Alm, R.A.; Hallinan, J.P.; Watson, A.A.; Mattick, J.S. Fimbrial biogenesis genes of Pseudomonas aeruginosa: PILW and PILX increase the similarity of type 4 fimbriae to the GSP protein-secretion systems and pilY1 encodes a gonococcal PilC homologue. Mol. Microbiol. 1996, 22, 161–173. [Google Scholar]
  90. Scheurwater, E.; Reid, C.W.; Clarke, A.J. Lytic transglycosylases: Bacterial space-making autolysins. Int. J. Biochem. Cell Biol. 2008, 40, 586–591. [Google Scholar] [CrossRef] [Green Version]
  91. Monchy, S.; Benotmane, M.A.; Janssen, P.; Vallaeys, T.; Taghavi, S.; van der Lelie, D.; Mergeay, M. Plasmids pMOL28 and pMOL30 of Cupriavidus metallidurans are specialized in the maximal viable response to heavy metals. J. Bacteriol. 2007, 189, 7417–7425. [Google Scholar] [CrossRef]
  92. Bertini, I.; Cavallaro, G.; Rosato, A. Cytochrome c: Occurrence and functions. Chem. Rev. 2006, 106, 90–115. [Google Scholar] [CrossRef]
  93. Cianciotto, N.P.; White, R.C. Expanding Role of Type II Secretion in Bacterial Pathogenesis and Beyond. Infect. Immun. 2017. [Google Scholar] [CrossRef]
  94. Xu, H.; Denny, T.P. Native and Foreign Proteins Secreted by the Cupriavidus metallidurans Type II System and an Alternative Mechanism. J. Microbiol. Biotechnol. 2017, 27, 791–807. [Google Scholar] [CrossRef] [PubMed]
Genes 10 00063 g001 550
Figure 1. Comparison of the size (bp), number of coding sequences (CDSs), and GC content (% GC) of the C. metallidurans NA4 (red) and CH34 (blue) genome replicons (CHR1: chromosome; CHR2: chromid; p: plasmids).
Figure 1. Comparison of the size (bp), number of coding sequences (CDSs), and GC content (% GC) of the C. metallidurans NA4 (red) and CH34 (blue) genome replicons (CHR1: chromosome; CHR2: chromid; p: plasmids).
Genes 10 00063 g001
Genes 10 00063 g002 550
Figure 2. Circa plot of the C. metallidurans CH34 and NA4 genomes (CHR1: chromosome; CHR2: chromid) and plasmids (28: pMOL28; 30: pMOL30; pA: pNA4_A; pB: pNA4_B; pC: pNA4_C; pD: pNA4_D). Connections correspond to homologous genes based on the bidirectional best hit criterion and a blastP alignment threshold (at least 35% amino-acid identity on 80% of the length of the smallest protein).
Figure 2. Circa plot of the C. metallidurans CH34 and NA4 genomes (CHR1: chromosome; CHR2: chromid) and plasmids (28: pMOL28; 30: pMOL30; pA: pNA4_A; pB: pNA4_B; pC: pNA4_C; pD: pNA4_D). Connections correspond to homologous genes based on the bidirectional best hit criterion and a blastP alignment threshold (at least 35% amino-acid identity on 80% of the length of the smallest protein).
Genes 10 00063 g002
Genes 10 00063 g003 550
Figure 3. Circa plot of pMOL30 from C. metallidurans CH34 and sections of the C. metallidurans NA4 chromid (CHR2). Ribbons correspond to synteny groups, orthologous gene sets that have the same local organization, based on the bidirectional best hit criterion or a blastP alignment threshold (at least 30% amino-acid identity on 80% of the length of the smallest protein) and co-localization (with the maximum number of consecutive genes not involved in a synteny group being five). The pMOL30 genes are colored according to their function.
Figure 3. Circa plot of pMOL30 from C. metallidurans CH34 and sections of the C. metallidurans NA4 chromid (CHR2). Ribbons correspond to synteny groups, orthologous gene sets that have the same local organization, based on the bidirectional best hit criterion or a blastP alignment threshold (at least 30% amino-acid identity on 80% of the length of the smallest protein) and co-localization (with the maximum number of consecutive genes not involved in a synteny group being five). The pMOL30 genes are colored according to their function.
Genes 10 00063 g003
Genes 10 00063 g004 550
Figure 4. Circa plot of pMOL30 and pMOL28 from C. metallidurans CH34 and pNA4_A and pNA4_B from C. metallidurans NA4. Connections correspond to homologous genes based on the bidirectional best hit criterion and a blastP alignment threshold (at least 35% amino-acid identity on 80% of the length of the smallest protein). Colored connections correspond to genes involved in metal resistance. The pMOL28 and pMOL30 genes are colored according to their function.
Figure 4. Circa plot of pMOL30 and pMOL28 from C. metallidurans CH34 and pNA4_A and pNA4_B from C. metallidurans NA4. Connections correspond to homologous genes based on the bidirectional best hit criterion and a blastP alignment threshold (at least 35% amino-acid identity on 80% of the length of the smallest protein). Colored connections correspond to genes involved in metal resistance. The pMOL28 and pMOL30 genes are colored according to their function.
Genes 10 00063 g004
Genes 10 00063 g005 550
Figure 5. Growth of C. metallidurans NA4 and NA4Pt in the presence of different PtCl4 concentrations (in µM).
Figure 5. Growth of C. metallidurans NA4 and NA4Pt in the presence of different PtCl4 concentrations (in µM).
Genes 10 00063 g005
Genes 10 00063 g006 550
Figure 6. Agarose gel electrophoresis of C. metallidurans NA4 (C) and NA4Pt (Pt) plasmid DNA. (a) Megaplasmids of NA4 vs. NA4Pt; (b) PagI digest of pNA4_D in NA4 vs. NA4Pt (only small plasmids were extracted; therefore, no discrete bands are visible for NA4); and (c) theoretical PagI digest of NA4Pt pNA4_D.
Figure 6. Agarose gel electrophoresis of C. metallidurans NA4 (C) and NA4Pt (Pt) plasmid DNA. (a) Megaplasmids of NA4 vs. NA4Pt; (b) PagI digest of pNA4_D in NA4 vs. NA4Pt (only small plasmids were extracted; therefore, no discrete bands are visible for NA4); and (c) theoretical PagI digest of NA4Pt pNA4_D.
Genes 10 00063 g006
Genes 10 00063 g007 550
Figure 7. Scatter plot of RNA-Seq-derived gene expression of C. metallidurans NA4Pt compared to its parental strain in non-selective conditions. Dots (blue p < 0.05) represent Log2 ratios with red, black, and green lines corresponding to −0.8, 8, and 0.8, respectively.
Figure 7. Scatter plot of RNA-Seq-derived gene expression of C. metallidurans NA4Pt compared to its parental strain in non-selective conditions. Dots (blue p < 0.05) represent Log2 ratios with red, black, and green lines corresponding to −0.8, 8, and 0.8, respectively.
Genes 10 00063 g007
Genes 10 00063 g008 550
Figure 8. EggNOG classification of differentially expressed genes in C. metallidurans NA4Pt, based on RNA-Seq-derived gene expression of C. metallidurans NA4Pt compared to its parental strain in non-selective conditions. Percentages were calculated by normalizing the number of up- and down-regulated genes in each category by the total number of genes in the NA4 genome grouped in the corresponding category. (D: Cell cycle control, cell division, and chromosome partitioning; M: Cell wall/membrane/envelope biogenesis; N: Cell motility; O: Posttranslational modification, protein turnover, and chaperones; T: Signal transduction mechanisms; U: Intracellular trafficking, secretion, and vesicular transport; V: Defense mechanisms; A: RNA processing and modification; B: Chromatin structure and dynamics; J: Translation, ribosomal structure, and biogenesis; K: Transcription; L: Replication, recombination, and repair; C: Energy production and conversion; E: Amino acid transport and metabolism; F: Nucleotide transport and metabolism; G: Carbohydrate transport and metabolism; H: Coenzyme transport and metabolism; I: Lipid transport and metabolism; P: Inorganic ion transport and metabolism; Q: Secondary metabolites biosynthesis, transport, and catabolism; S: Poorly characterized).
Figure 8. EggNOG classification of differentially expressed genes in C. metallidurans NA4Pt, based on RNA-Seq-derived gene expression of C. metallidurans NA4Pt compared to its parental strain in non-selective conditions. Percentages were calculated by normalizing the number of up- and down-regulated genes in each category by the total number of genes in the NA4 genome grouped in the corresponding category. (D: Cell cycle control, cell division, and chromosome partitioning; M: Cell wall/membrane/envelope biogenesis; N: Cell motility; O: Posttranslational modification, protein turnover, and chaperones; T: Signal transduction mechanisms; U: Intracellular trafficking, secretion, and vesicular transport; V: Defense mechanisms; A: RNA processing and modification; B: Chromatin structure and dynamics; J: Translation, ribosomal structure, and biogenesis; K: Transcription; L: Replication, recombination, and repair; C: Energy production and conversion; E: Amino acid transport and metabolism; F: Nucleotide transport and metabolism; G: Carbohydrate transport and metabolism; H: Coenzyme transport and metabolism; I: Lipid transport and metabolism; P: Inorganic ion transport and metabolism; Q: Secondary metabolites biosynthesis, transport, and catabolism; S: Poorly characterized).
Genes 10 00063 g008
Genes 10 00063 g009 550
Figure 9. Scanning electron microscopy (SEM) image of C. metallidurans NA4 and NA4Pt.
Figure 9. Scanning electron microscopy (SEM) image of C. metallidurans NA4 and NA4Pt.
Genes 10 00063 g009
Table
Table 1. Overall characterization of the replicons constituting the C. metallidurans NA4 genome.
Table 1. Overall characterization of the replicons constituting the C. metallidurans NA4 genome.
RepliconLength (bp)% GC 1# CDSs 2
Chromosome3,838,19563.843818
Chromid2,776,39563.592774
Plasmid pNA4_A294,57559.45416
Plasmid pNA4_B227,79660.57274
Plasmid pNA4_C155,04158.90209
Plasmid pNA4_D89,60661.56120
1 Guanine-cytosine content; 2 Coding sequences.
Table
Table 2. Plasmid characterization.
Table 2. Plasmid characterization.
RepliconRep TypeRelaxase Type 1MPF Type 2Predicted MobilityMASH Nearest Neighbor 3
pMOL28591MOBHMPFTConjugative-
pMOL30332(MOBP)MPFFConjugative-
pNA4_A332-MPFFNon-mobilizablepMOL30
pNA4_B591MOBHMPFTConjugativepMOL28
pNA4_C---Non-mobilizablepHS87a (Pseudomonas aeruginosa)
pNA4_D1864MOBFMPFFConjugativepACP3.3 (Acidovorax sp. P3)
1 Archetype relaxase: plasmid R388 (MOBF), plasmid R27 (MOBH), and plasmid RP4 (MOBP). 2 Archetype mating pair formation (MPF) complex: plasmid F (MPFF) and plasmid Ti (MPFT). 3 MASH Nearest Neighbor are not included for pMOL28 and pMOL30, as the closest database match was a self-hit.
Table
Table 3. Distribution of insertion sequence elements in C. metallidurans NA4.
Table 3. Distribution of insertion sequence elements in C. metallidurans NA4.
ElementFamilySize (%) 1CHR1CHR2pNA4_ApNA4_BpNA4_CpNA4_D
IS1071Tn33204 (99.9%) 1
2991 (93.4%) 2
ISRme4IS212469 (100%)43
ISRme9IS212674 (94.8%) 1
ISRme10IS301063 (100%) 1
ISRme3IS31288 (100%)12 2
ISPst3IS212605 (97.8%)12
ISPa45IS41637 (100%)1
1 Size (bp) of the element and % nucleotide sequence similarity to the insertion sequence (IS) element as defined in ISFinder [62].
Table
Table 4. Minimal inhibitory concentration of selected metals for C. metallidurans NA4 and NA4Pt.
Table 4. Minimal inhibitory concentration of selected metals for C. metallidurans NA4 and NA4Pt.
Pt4+ (µM)Pd2+ (µM)Ag+ (µM)Zn2+ (mM)Ni2+ (mM)Cu2+ (mM)
NA47012.5112406
NA4Pt16012.5112406

Share and Cite

MDPI and ACS Style

Ali, M.M.; Provoost, A.; Maertens, L.; Leys, N.; Monsieurs, P.; Charlier, D.; Van Houdt, R. Genomic and Transcriptomic Changes That Mediate Increased Platinum Resistance in Cupriavidus metallidurans. Genes 2019, 10, 63. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010063

AMA Style

Ali MM, Provoost A, Maertens L, Leys N, Monsieurs P, Charlier D, Van Houdt R. Genomic and Transcriptomic Changes That Mediate Increased Platinum Resistance in Cupriavidus metallidurans. Genes. 2019; 10(1):63. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010063

Chicago/Turabian Style

Ali, Md Muntasir, Ann Provoost, Laurens Maertens, Natalie Leys, Pieter Monsieurs, Daniel Charlier, and Rob Van Houdt. 2019. "Genomic and Transcriptomic Changes That Mediate Increased Platinum Resistance in Cupriavidus metallidurans" Genes 10, no. 1: 63. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10010063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop