Temporal and Spatial Diversity of the Tap Water Microbiota in a Norwegian Hospital

We analyzed the temporal and spatial diversity of the microbiota in a low-usage and a high-usage hospital tap. We identified a tap-specific colonization pattern, with potential human pathogens being overrepresented in the low-usage tap. We propose that founder effects and local adaptation caused the tap-specific colonization patterns. Our conclusion is that tap-specific colonization represents a potential challenge for water safety.Humans are exposed to and consume large amounts of tap water in their everyday life, with the tap water microbiota representing a potent reservoir for pathogens (8). Despite the potential impact, our knowledge about the ecological diversification processes of the tap water microbiota is limited (4, 11).The aim of the present work was to determine the temporal and spatial distribution patterns of the planktonic tap water microbiota. We compared the summer and winter microbiota from two hospital taps supplied from the same water source. We analyzed 16S rRNA gene clone libraries by using a novel alignment-independent approach for operational taxonomic unit (OTU) designation (6), while established OTU diversity and richness estimators were used for the ecological interpretations.Tap water samples (1 liter) from a high-usage kitchen and a low-usage toilet cold-water tap in Akershus University Hospital, Lørenskog, Norway, were collected in January and July 2006. The total DNA was isolated and the 16S rRNA gene PCR amplified and sequenced. Based on the sequences, we estimated the species richness and diversity, we calculated the distances between the communities, and trees were constructed to reflect the relatedness of the microbiota in the samples analyzed. Details about these analytical approaches are given in the materials and methods section in the supplemental material.Our initial analysis of species composition was done using the RDPII hierarchical classifier. We found that the majority of pathogen-related bacteria in our data set belonged to the class Gammaproteobacteria. The genera encompassed Legionella, Pseudomonas, and Vibrio (Table ​(Table1).1). We found a significant overrepresentation of pathogen-related bacteria in the toilet tap (P = 0.04), while there were no significant differences between summer and winter samples. Legionella showed the highest relative abundance for the pathogen-related bacteria. With respect to the total diversity, we found that Proteobacteria dominated the tap water microbiota (representing 86% of the taxa) (see Table S1 in the supplemental material). There was, however, a large portion (56%) of the taxa that could not be assigned to the genus level using this classifier.

TABLE 1.

Cloned sequences related to human pathogens^a

Unsuitability of Quantitative Bacteroidales 16S rRNA Gene Assays for Discerning Fecal Contamination of Drinking Water

Bacteroidales species were detected in (tap) water samples from treatment plants with three different PCR assays. 16S rRNA gene sequence analysis indicated that the sequences had an environmental rather than fecal origin. We conclude that assays for Bacteroidales 16S rRNA genes are not specific enough to discern fecal contamination of drinking water in the Netherlands.Drinking water in many countries is routinely monitored for recent fecal contamination by testing for fecal indicator organisms Escherichia coli, thermotolerant coliforms, and/or intestinal enterococci to demonstrate microbial safety (13, 21, 42). Although these indicator organisms have been used for many decades, they have some limitations: the number of E. coli/coliform/enterococcus bacteria in feces is relatively low (18, 38), and they sometimes might be able to grow in the environment (10, 11, 14, 27). Consequently, scientists have been searching for alternative indicator organisms to determine fecal contamination of water. In 1967, bacteria belonging to the genus Bacteroides were suggested as alternative indicator organisms (26). Bacteroides spp. might have some advantages over the traditional indicator organisms. The numbers of Bacteroides spp. in the intestinal tract of humans and animals are 10 to 100 times higher than the numbers of E. coli or intestinal enterococci (1, 2, 12, 26). However, the use of Bacteroides spp. as indicator organisms was hampered by the complex cultivation conditions required (1, 2). The introduction of molecular methods made it possible to detect bacterial species that belong to the order Bacteroidales, an order that includes the genus Bacteroides, without cultivation. As a result, real-time PCR methods were developed for the quantitative detection of Bacteroidales in surface and recreation water and the potential of Bacteroidales species as an indication of fecal contamination of recreational waters was demonstrated (6, 12, 16, 19, 20, 29). Bacteroidales species might be useful indicator organisms for fecal contamination of drinking water as well. However, methods to detect fecal contamination in drinking water should be more sensitive, because people ingest more drinking water and the quality assessments and standards for fecal contamination are stricter than for bathing water. Studies exploring real-time PCR for the detection of Bacteroidales genes in drinking water have not been published to our knowledge. The objective of our study was, therefore, to determine if assays for the detection of Bacteroidales 16S rRNA genes can be used to detect fecal contamination in drinking water.Unchlorinated tap water samples were obtained in November 2007 and February 2010 from one or more locations in the distribution systems of nine different drinking water treatment plants (plants A to I; Table ​Table1)1) that produced unchlorinated drinking water from confined (plants B, C, E, F, and G) and unconfined (plants A, D, H, and I) groundwater. The treatment plants are located in the central part of the Netherlands within 90 km of each other. In addition, untreated groundwater from extraction wells and/or untreated raw groundwater (mixture of groundwater from different extraction wells) was sampled in March 2008 (Table ​(Table1).1). Water samples (100 ml) were filtered over a 25-mm polycarbonate filter (0.22-μm pore size, type GTTP; Millipore, Netherlands) and a DNA fragment was added as internal control to determine the recovery efficiency of DNA isolation and PCR analysis (2a, 40). DNA was isolated using a FastDNA spin kit for soil (Qbiogene, United States) according to the supplier''s protocol. Primer sets AllBac 296f and AllBac 412r, resulting in a PCR product of 108 bp, were used in combination with TaqMan probe AllBac375Bhqr to quantitatively determine the number of Bacteroidales 16S rRNA gene copies in the water samples using a real-time PCR instrument (20). The PCR cycle after which the fluorescence signal of the amplified DNA was detected (threshold cycle [C_T]) was used to quantify the concentration of 16S rRNA gene copies. Quantification was based on comparison of the sample C_T value with the C_T values of a calibration curve graphed using known copy numbers of the Bacteroidales 16S rRNA gene, as previously described (12, 20). The correlation coefficient of the calibration curve was 0.99, and the efficiency of the PCR 95 to 105%. Finally, the Bacteroidales cell number was calculated by using the recovery rate of the internal standard and assuming five 16S rRNA gene copy numbers per cell (5). The detection limit of this gene assay was 50 Bacteroidales cells 100 ml⁻¹ (corresponding to 10 16S rRNA gene copies per reaction mixture). Furthermore, the 16S rRNA genes that were obtained from several water samples from treatment plant C with the AllBac and TotBac (12) primer sets were sequenced, and the nearest relatives were obtained from the GenBank database using BLAST searches.

TABLE 1.

Numbers of Bacteroidales cells in extraction wells, raw groundwater, and unchlorinated tap water of nine different groundwater plants in the Netherlands^a

Zinc-Independent Folate Biosynthesis: Genetic,Biochemical, and Structural Investigations Reveal New Metal Dependence for GTP Cyclohydrolase IB

GTP cyclohydrolase I (GCYH-I) is an essential Zn²⁺-dependent enzyme that catalyzes the first step of the de novo folate biosynthetic pathway in bacteria and plants, the 7-deazapurine biosynthetic pathway in Bacteria and Archaea, and the biopterin pathway in mammals. We recently reported the discovery of a new prokaryotic-specific GCYH-I (GCYH-IB) that displays no sequence identity to the canonical enzyme and is present in ∼25% of bacteria, the majority of which lack the canonical GCYH-I (renamed GCYH-IA). Genomic and genetic analyses indicate that in those organisms possessing both enzymes, e.g., Bacillus subtilis, GCYH-IA and -IB are functionally redundant, but differentially expressed. Whereas GCYH-IA is constitutively expressed, GCYH-IB is expressed only under Zn²⁺-limiting conditions. These observations are consistent with the hypothesis that GCYH-IB functions to allow folate biosynthesis during Zn²⁺ starvation. Here, we present biochemical and structural data showing that bacterial GCYH-IB, like GCYH-IA, belongs to the tunneling-fold (T-fold) superfamily. However, the GCYH-IA and -IB enzymes exhibit significant differences in global structure and active-site architecture. While GCYH-IA is a unimodular, homodecameric, Zn²⁺-dependent enzyme, GCYH-IB is a bimodular, homotetrameric enzyme activated by a variety of divalent cations. The structure of GCYH-IB and the broad metal dependence exhibited by this enzyme further underscore the mechanistic plasticity that is emerging for the T-fold superfamily. Notably, while humans possess the canonical GCYH-IA enzyme, many clinically important human pathogens possess only the GCYH-IB enzyme, suggesting that this enzyme is a potential new molecular target for antibacterial development.The Zn²⁺-dependent enzyme GTP cyclohydrolase I (GCYH-I; EC 3.5.4.16) is the first enzyme of the de novo tetrahydrofolate (THF) biosynthesis pathway (Fig. ​(Fig.1)1) (38). THF is an essential cofactor in one-carbon transfer reactions in the synthesis of purines, thymidylate, pantothenate, glycine, serine, and methionine in all kingdoms of life (38), and formylmethionyl-tRNA in bacteria (7). Recently, it has also been shown that GCYH-I is required for the biosynthesis of the 7-deazaguanosine-modified tRNA nucleosides queuosine and archaeosine produced in Bacteria and Archaea (44), respectively, as well as the 7-deazaadenosine metabolites produced in some Streptomyces species (33). GCYH-I is encoded in Escherichia coli by the folE gene (28) and catalyzes the conversion of GTP to 7,8-dihydroneopterin triphosphate (55), a complex reaction that begins with hydrolytic opening of the purine ring at C-8 of GTP to generate an N-formyl intermediate, followed by deformylation and subsequent rearrangement and cyclization of the ribosyl moiety to generate the pterin ring in THF (Fig. ​(Fig.1).1). Notably, the enzyme is dependent on an essential active-site Zn²⁺ that serves to activate a water molecule for nucleophilic attack at C-8 in the first step of the reaction (2).Open in a separate windowFIG. 1.Reaction catalyzed by GCYH-I, and metabolic fate of 7,8-dihydroneopterin triphosphate.A homologous GCYH-I is found in mammals and other higher eukaryotes, where it catalyzes the first step of the biopterin (BH₄) pathway (Fig. ​(Fig.1),1), an essential cofactor in the biosynthesis of tyrosine and neurotransmitters, such as serotonin and l-3,4-dihydroxyphenylalanine (3, 52). Recently, a distinct class of GCYH-I enzymes, GCYH-IB (encoded by the folE2 gene), was discovered in microbes (26% of sequenced Bacteria and most Archaea) (12), including several clinically important human pathogens, e.g., Neisseria and Staphylococcus species. Notably, GCYH-IB is absent in eukaryotes.The distribution of folE (gene product renamed GCYH-IA) and folE2 (GCYH-IB) in bacteria is diverse (12). The majority of organisms possess either a folE (65%; e.g., Escherichia coli) or a folE2 (14%; e.g., Neisseria gonorrhoeae) gene. A significant number (12%; e.g., B. subtilis) possess both genes (a subset of 50 bacterial species is shown in Table ​Table1),1), and 9% lack both genes, although members of the latter group are mainly intracellular or symbiotic bacteria that rely on external sources of folate. The majority of Archaea possess only a folE2 gene, and the encoded GCYH-IB appears to be necessary only for the biosynthesis of the modified tRNA nucleoside archaeosine (44) except in the few halophilic Archaea that are known to synthesize folates, such as Haloferax volcanii, where GCYH-IB is involved in both archaeosine and folate formation (13, 44).

TABLE 1.

Distribution and candidate Zur-dependent regulation of alternative GCYH-I genes in bacteria^a

Antimicrobial Activity of Simulated Solar Disinfection against Bacterial,Fungal, and Protozoan Pathogens and Its Enhancement by Riboflavin

Riboflavin significantly enhanced the efficacy of simulated solar disinfection (SODIS) at 150 watts per square meter (W m⁻²) against a variety of microorganisms, including Escherichia coli, Fusarium solani, Candida albicans, and Acanthamoeba polyphaga trophozoites (>3 to 4 log₁₀ after 2 to 6 h; P < 0.001). With A. polyphaga cysts, the kill (3.5 log₁₀ after 6 h) was obtained only in the presence of riboflavin and 250 W m⁻² irradiance.Solar disinfection (SODIS) is an established and proven technique for the generation of safer drinking water (11). Water is collected into transparent plastic polyethylene terephthalate (PET) bottles and placed in direct sunlight for 6 to 8 h prior to consumption (14). The application of SODIS has been shown to be a simple and cost-effective method for reducing the incidence of gastrointestinal infection in communities where potable water is not available (2-4). Under laboratory conditions using simulated sunlight, SODIS has been shown to inactivate pathogenic bacteria, fungi, viruses, and protozoa (6, 12, 15). Although SODIS is not fully understood, it is believed to achieve microbial killing through a combination of DNA-damaging effects of ultraviolet (UV) radiation and thermal inactivation from solar heating (21).The combination of UVA radiation and riboflavin (vitamin B₂) has recently been reported to have therapeutic application in the treatment of bacterial and fungal ocular pathogens (13, 17) and has also been proposed as a method for decontaminating donor blood products prior to transfusion (1). In the present study, we report that the addition of riboflavin significantly enhances the disinfectant efficacy of simulated SODIS against bacterial, fungal, and protozoan pathogens.Chemicals and media were obtained from Sigma (Dorset, United Kingdom), Oxoid (Basingstoke, United Kingdom), and BD (Oxford, United Kingdom). Pseudomonas aeruginosa (ATCC 9027), Staphylococcus aureus (ATCC 6538), Bacillus subtilis (ATCC 6633), Candida albicans (ATCC 10231), and Fusarium solani (ATCC 36031) were obtained from ATCC (through LGC Standards, United Kingdom). Escherichia coli (JM101) was obtained in house, and the Legionella pneumophila strain used was a recent environmental isolate.B. subtilis spores were produced from culture on a previously published defined sporulation medium (19). L. pneumophila was grown on buffered charcoal-yeast extract agar (5). All other bacteria were cultured on tryptone soy agar, and C. albicans was cultured on Sabouraud dextrose agar as described previously (9). Fusarium solani was cultured on potato dextrose agar, and conidia were prepared as reported previously (7). Acanthamoeba polyphaga (Ros) was isolated from an unpublished keratitis case at Moorfields Eye Hospital, London, United Kingdom, in 1991. Trophozoites were maintained and cysts prepared as described previously (8, 18).Assays were conducted in transparent 12-well tissue culture microtiter plates with UV-transparent lids (Helena Biosciences, United Kingdom). Test organisms (1 × 10⁶/ml) were suspended in 3 ml of one-quarter-strength Ringer''s solution or natural freshwater (as pretreated water from a reservoir in United Kingdom) with or without riboflavin (250 μM). The plates were exposed to simulated sunlight at an optical output irradiance of 150 watts per square meter (W m⁻²) delivered from an HPR125 W quartz mercury arc lamp (Philips, Guildford, United Kingdom). Optical irradiances were measured using a calibrated broadband optical power meter (Melles Griot, Netherlands). Test plates were maintained at 30°C by partial submersion in a water bath.At timed intervals for bacteria and fungi, the aliquots were plated out by using a WASP spiral plater and colonies subsequently counted by using a ProtoCOL automated colony counter (Don Whitley, West Yorkshire, United Kingdom). Acanthamoeba trophozoite and cyst viabilities were determined as described previously (6). Statistical analysis was performed using a one-way analysis of variance (ANOVA) of data from triplicate experiments via the InStat statistical software package (GraphPad, La Jolla, CA).The efficacies of simulated sunlight at an optical output irradiance of 150 W m⁻² alone (SODIS) and in the presence of 250 μM riboflavin (SODIS-R) against the test organisms are shown in Table ​Table1.1. With the exception of B. subtilis spores and A. polyphaga cysts, SODIS-R resulted in a significant increase in microbial killing compared to SODIS alone (P < 0.001). In most instances, SODIS-R achieved total inactivation by 2 h, compared to 6 h for SODIS alone (Table ​(Table1).1). For F. solani, C. albicans, ands A. polyphaga trophozoites, only SODIS-R achieved a complete organism kill after 4 to 6 h (P < 0.001). All control experiments in which the experiments were protected from the light source showed no reduction in organism viability over the time course (results not shown).

TABLE 1.

Efficacies of simulated SODIS for 6 h alone and with 250 μM riboflavin (SODIS-R)

Detection and Quantification of the Coral Pathogen Vibrio coralliilyticus by Real-Time PCR with TaqMan Fluorescent Probes

A real-time quantitative PCR-based detection assay targeting the dnaJ gene (encoding heat shock protein 40) of the coral pathogen Vibrio coralliilyticus was developed. The assay is sensitive, detecting as little as 1 CFU per ml in seawater and 10⁴ CFU per cm² of coral tissue. Moreover, inhibition by DNA and cells derived from bacteria other than V. coralliilyticus was minimal. This assay represents a novel approach to coral disease diagnosis that will advance the field of coral disease research.Vibrio coralliilyticus has recently emerged as a coral pathogen of concern on reefs throughout the Indo-Pacific. It was first implicated as the etiological agent responsible for bleaching and tissue lysis of the coral Pocillopora damicornis on Zanzibar reefs (2). More recently, V. coralliilyticus has been identified as the causative agent of white syndrome (WS) outbreaks on several Pacific reefs (14). WS is a collective term describing coral diseases characterized by a spreading band of tissue loss exposing white skeleton on Indo-Pacific scleractinian corals (16). V. coralliilyticus is an emerging model pathogen for understanding the mechanisms linking bacterial infection and coral disease (13) and therefore provides an ideal model for the development of diagnostic assays to detect coral disease. Current coral disease diagnostic methods, which are based primarily upon field-based observations of macroscopic disease signs, often detect disease only at the latest stages of infection, when control measures are least effective. The development of diagnostic tools targeting pathogens underlying coral disease pathologies may provide early indications of infection, aid the identification of disease vectors and reservoirs, and assist managers in developing strategies to prevent the spread of coral disease outbreaks. In this paper, we describe the development and validation of a TaqMan-based real-time quantitative PCR (qPCR) assay that targets a segment of the V. coralliilyticus heat shock protein 40-encoding gene (dnaJ).Nucleotide sequences of the dnaJ gene were retrieved from relevant Vibrio species, including V. coralliilyticus (LMG 20984), using the National Center for Biotechnology Information''s (NCBI) Entrez Nucleotide Database search tool (http://www.ncbi.nlm.nih.gov/). Gene sequences of strains not available in public databases (V. coralliilyticus strains LMG 21348, LMG 21349, LMG 21350, LMG 10953, LMG 20538, LMG 23696, LMG 23691, LMG 23693, LMG 23692, and LMG 23694) were obtained through extraction of total DNA using a Promega Wizard Prep DNA Purification Kit (Promega, Sydney, Australia), PCR amplification, and sequencing using primers and thermal cycling parameters described by Nhung et al. (8). A 128-bp region (nucleotides 363 to 490) containing high concentrations of single nucleotide polymorphisms (SNPs), which were conserved within V. coralliilyticus strains but differed from non-V. coralliilyticus strains, was identified, and oligonucleotide primers Vc_dnaJ_F1 (5′-CGG TTC GYG GTG TTT CAA AA-3′) and Vc_dnaJ_R1 (5′-AAC CTG ACC ATG ACC GTG ACA-3′) and a TaqMan probe, Vc_dnaJ_TMP (5′-6-FAM-CAG TGG CGC GAA G-MGBNFQ-3′; 6-FAM is 6-carboxyfluorescein and MGBNFQ is molecular groove binding nonfluorescent quencher), were designed to target this region. The qPCR assay was optimized and validated using DNA extracted from V. coralliilyticus isolates, nontarget Vibrio species, and other bacterial species grown in marine broth (MB) (Table ​(Table1),1), under the following optimal conditions: 1× TaqMan buffer A, 0.5 U of AmpliTaq Gold DNA polymerase, 200 μM deoxynucleotide triphosphates (with 400 μM dUTP replacing deoxythymidine triphosphate), 0.2 U of AmpErase uracil N-glycosylase (UNG), 3 mM MgCl₂, 0.6 μM each primer, 0.2 μM fluorophore-labeled TaqMan, 1 μl of template, and sterile MilliQ water for a total reaction volume to 20 μl. All assays were conducted on a RotoGene 300 (Corbett Research, Sydney, Australia) real-time analyzer with the following cycling parameters: 50°C for 120 s (UNG activation) and 95°C for 10 min (AmpliTaq Gold DNA polymerase activation), followed by 40 cycles of 95°C for 15 s (denaturation) and 60°C for 60 s (annealing/extension). During the annealing/extension phase of each thermal cycle, fluorescence was measured in the FAM channel (470-nm excitation and 510-nm detection).

TABLE 1.

Species, strain, and threshold cycle for all bacterial strains tested^a

Increased Crystalline Cellulose Activity via Combinations of Amino Acid Changes in the Family 9 Catalytic Domain and Family 3c Cellulose Binding Module of Thermobifida fusca Cel9A

Amino acid modifications of the Thermobifida fusca Cel9A-68 catalytic domain or carbohydrate binding module 3c (CBM3c) were combined to create enzymes with changed amino acids in both domains. Bacterial crystalline cellulose (BC) and swollen cellulose (SWC) assays of the expressed and purified enzymes showed that three combinations resulted in 150% and 200% increased activity, respectively, and also increased synergistic activity with other cellulases. Several other combinations resulted in drastically lowered activity, giving insight into the need for a balance between the binding in the catalytic cleft on either side of the cleavage site, as well as coordination between binding affinity for the catalytic domain and CBM3c. The same combinations of amino acid variants in the whole enzyme, Cel9A-90, did not increase BC or SWC activity but did have higher filter paper (FP) activity at 12% digestion.Cellulases catalyze the breakdown of cellulose into simple sugars that can be fermented to ethanol. The large amount of natural cellulose available is an exciting potential source of fuels and chemicals. However, the detailed molecular mechanisms of crystalline cellulose degradation by glycoside hydrolases are still not well understood and their low efficiency is a major barrier to cellulosic ethanol production.Thermobifida fusca is a filamentous soil bacterium that grows at 50°C in defined medium and can utilize cellulose as its sole carbon source. It is a major degrader of plant cell walls in heated organic materials such as compost piles and rotting hay and produces a set of enzymes that includes six different cellulases, three xylanases, a xyloglucanase, and two CBM33 binding proteins (12). Among them are three endocellulases, Cel9B, Cel6A, and Cel5A (7, 8), two exocellulases, Cel48A and Cel6B (6, 19), and a processive endocellulase, Cel9A (5, 7).T. fusca Cel9A-90 (Uniprot {"type":"entrez-protein","attrs":{"text":"P26221","term_id":"2506384","term_text":"P26221"}}P26221 and {"type":"entrez-protein","attrs":{"text":"YP_290232","term_id":"72162575","term_text":"YP_290232"}}YP_290232) is a multidomain enzyme consisting of a family 9 catalytic domain (CD) rigidly attached by a short linker to a family 3c cellulose binding module (CBM3c), followed by a fibronectin III-like domain and a family 2 CBM (CBM2). Cel9A-68 consists of the family 9 CD and CBM3c. The crystal structure of this species (Fig. ​(Fig.1)1) was determined by X-ray crystallography at 1.9 Å resolution (Protein Data Bank [PDB] code 4tf4) (15). Previous work has shown that E424 is the catalytic acid and D58 is the catalytic base (11, 20). H125 and Y206 were shown to play an important role in activity by forming a hydrogen bonding network with D58, an important supporting residue, D55, and Glc(−1)O1. Several enzymes with amino acid changes in subsites Glc(−1) to Glc(−4) had less than 20% activity on bacterial cellulose (BC) and markedly reduced processivity. It was proposed that these modifications disturb the coordination between product release and the subsequent binding of a cellulose chain into subsites Glc(−1) to Glc(−4) (11). Another variant enzyme with a deletion of a group of amino acids forming a block at the end of the catalytic cleft, Cel9A-68 Δ(T245-L251)R252K (DEL), showed slightly improved filter paper (FP) activity and binding to BC (20).Open in a separate windowFIG. 1.Crystal structure of Cel9A-68 (PDB code 4tf4) showing the locations of the variant residues, catalytic acid E424, catalytic base D58, hydrogen bonding network residues D55, H125, and Y206, and six glucose residues, Glc(−4) to Glc(+2). Part of the linker is visible in dark blue.The CBM3c domain is critical for hydrolysis and processivity. Cel9A-51, an enzyme with the family 9 CD and the linker but without CBM3c, had low activity on carboxymethyl cellulose (CMC), BC, and swollen cellulose (SWC) and showed no processivity (4). The role of CBM3c was investigated by mutagenesis, and one modified enzyme, R557A/E559A, had impaired activity on all of these substrates but normal binding and processivity (11). Variants with changes at five other CBM3c residues were found to slightly lower the activity of the modified enzymes, while Cel9A-68 enzymes containing either F476A, D513A, or I514H were found to have slightly increased binding and processivity (11) (see Table ​Table1).1). In the present work, CBM3c has been investigated more extensively to identify residues involved in substrate binding and processivity, understand the role of CBM3c more clearly, and study the coordination between the CD and CBM3c. An additional goal was to combine amino acid variants showing increased crystalline cellulose activity to see if this further increased activity. Finally, we have investigated whether the changes that improved the activity of Cel9A-68 also enhanced the activity of intact Cel9A-90.

TABLE 1.

Activities of Cel9A-68 CBM3c variant enzymes and CD variant enzymes used to create the double variants

Biochemical Characterization of Haloalkane Dehalogenases DrbA and DmbC,Representatives of a Novel Subfamily

This study focuses on two representatives of experimentally uncharacterized haloalkane dehalogenases from the subfamily HLD-III. We report biochemical characterization of the expression products of haloalkane dehalogenase genes drbA from Rhodopirellula baltica SH1 and dmbC from Mycobacterium bovis 5033/66. The DrbA and DmbC enzymes show highly oligomeric structures and very low activities with typical substrates of haloalkane dehalogenases.Haloalkane dehalogenases (EC 3.8.1.5.) acting on halogenated aliphatic hydrocarbons catalyze carbon-halogen bond cleavage, leading to an alcohol, a halide ion, and a proton as the reaction products (7). Haloalkane dehalogenases originating from various bacterial strains have potential for application in bioremediation technologies (4, 6, 22), construction of biosensors (2), decontamination of warfare agents (17), and synthesis of optically pure compounds (19). Recent evolutionary study of haloalkane dehalogenase sequences revealed the existence of three subfamilies, denoted HLD-I, HLD-II, and HLD-III (3). In contrast to subfamilies HLD-I and HLD-II, the subfamily HLD-III is currently lacking experimentally characterized proteins. We have therefore focused on the isolation and study of two selected representatives of the HLD-III subfamily, DrbA and DmbC.The drbA gene was amplified by PCR using the cosmid pircos.a3g10 originating from marine bacterium Rhodopirellula baltica SH1, and the dmbC gene was amplified from DNA originating from obligatory pathogen Mycobacterium bovis 5033/66. Six-histidine tails were added to the C termini of DrbA and DmbC in a cloning step, enabling single-step purification using Ni-nitrilotriacetic acid resin. Haloalkane dehalogenase DrbA was expressed under the T7 promoter and purified, with a resulting yield of 0.1 mg of protein per gram of cell mass. Haloalkane dehalogenase DmbC was obtained by expression in Mycobacterium smegmatis, with a yield of 0.07 mg of purified protein per gram of cell mass.The correct folding and secondary structures of the newly prepared enzymes were verified by circular dichroism (CD) spectroscopy. Far-UV CD spectra were recorded for DrbA and DmbC enzymes and other, related haloalkane dehalogenases. All enzymes tested exhibited CD spectra with two negative features at 208 and 222 nm and one positive peak at 195 nm, which are characteristic of α-helical content (Fig. ​(Fig.1).1). This suggested that both new enzymes, DrbA and DmbC, were folded correctly. However, DmbC exhibited more intense negative maxima which differed from other haloalkane dehalogenases in the θ₂₂₂/θ₂₀₈ ratio. This finding indicated a slight variation in the arrangement of secondary structure elements of the DmbC enzyme. Thermally induced denaturations of DrbA and DmbC were tested in parallel. Both enzymes showed changes in ellipticity during increasing temperature. The melting temperatures calculated from these curves were 45.8 ± 0.4°C for DmbC and 39.4 ± 0.1°C for DrbA. The thermostability results obtained for DrbA and DmbC were in good agreement with the range of melting temperatures determined for other, related haloalkane dehalogenases.Open in a separate windowFIG. 1.Far-UV CD spectra of DrbA, DmbC, and seven different biochemically characterized haloalkane dehalogenases. Protein concentration used for far-UV CD spectrum measurement was 0.2 mg/ml.The sizes of the purified proteins were estimated by electrophoresis under native conditions conducted using a 10% polyacrylamide gel (Fig. ​(Fig.2).2). More precise determination of the sizes of DrbA and DmbC was achieved by gel filtration chromatography performed on Sephacryl S-500 HR (GE Healthcare, Uppsala, Sweden), calibrated with blue dextran 2000, thyroglobulin (669 kDa), ferritin (440 kDa), catalase (240 kDa), conalbumin (75 kDa), and ovalbumin (43 kDa) (Fig. ​(Fig.3A).3A). Both DrbA and DmbC were eluted from the column in the fraction prior to blue dextran, indicating that both enzymes form oligomeric complexes of a size larger than 2,000 kDa (Fig. 3B and C). The haloalkane dehalogenases which have been biochemically characterized so far form monomers, except for DbjA isolated from Bradyrhizobium japonicum USDA110 (21), which shows monomeric, dimeric, and tetrameric forms according to the pH of the buffer (R. Chaloupkova, submitted for publication).Open in a separate windowFIG. 2.Native protein electrophoresis of DrbA and DmbC. Lane 1, carbonic anhydrase (29 kDa); lane 2, ovalbumin (43 kDa); lane 3, bovine albumin (67 kDa); lane 4, conalbumin (75 kDa); lane 5, catalase (240 kDa); lane 6, ferritin (440 kDa); lane 7, DrbA; lane 8, DmbC.Open in a separate windowFIG. 3.Gel filtration chromatogram of DrbA and DmbC. (A) The following calibration kit samples (0.5 ml of a concentration of 2 mg/ml protein loaded) were analyzed using 50 mM Tris-HCl with 150 mM NaCl, pH 7.5, as elution buffer: blue dextran (line 1, 9.6-ml fraction), thyroglobulin (line 2, 15.95-ml fraction), ferritin (line 3, 16.78-ml fraction), ovalbumin (line 4, 18.55-ml fraction), and RNase A (line 5, 20.08-ml fraction). (B and C) Haloalkane dehalogenase DrbA eluted in the 9.03-ml fraction (B), and haloalkane dehalogenase DmbC in the 9.31-ml fraction (C).The substrate specificities of DrbA and DmbC were investigated with a set of 30 selected chlorinated, brominated, and iodinated hydrocarbons. Standardized specific activities related to 1-chlorobutane (summarized in Table ​Table1)1) were compared with the activity profiles of other haloalkane dehalogenases (Fig. ​(Fig.4).4). DrbA and DmbC displayed similar activity patterns, with catalytic activities approximately two orders of magnitude lower than those of other known haloalkane dehalogenases (1, 5, 8-11, 13-16, 18, 20, 23). HLD-III subfamily enzymes showed a restricted specificity range and a preference for iodinated short-chain hydrocarbons. Both phenomena may be related to the composition of the catalytic pentad Asp-His-Asp+Asn-Trp, which is unique to the members of the HLD-III subfamily (3). The preference for substrates carrying an iodine substituent can be related to a pair of halide-binding residues and their spatial arrangement with the catalytic triad. These residues make up the catalytic pentad, playing a critical role in substrate binding, formation of the transition states, and the reaction intermediates of the dehalogenation reaction (12).Open in a separate windowFIG. 4.Substrate specificity profiles of DrbA, DmbC, and seven different biochemically characterized haloalkane dehalogenases. Activities were determined using a consistent set of 30 halogenated substrates (see Table ​Table1).1). Data were standardized by dividing each value by the sum of all activities determined for individual enzymes in order to mask the differences in absolute activities. Specific activities (in μmol·s⁻¹·mg⁻¹) with 1-chlorobutane are 0.0003 (DrbA), 0.0001 (DmbC), 0.0003 (DatA), 0.0133 (DbjA), 0.0010 (DbeA), 0.0128 (DhaA), 0.0231 (LinB), 0.0171 (DmbA), and 0.0117 (DhlA).

TABLE 1.

Specific activities of haloalkane dehalogenases DrbA and DmbC toward a set of 30 halogenated hydrocarbons^a

Evolutionary Strata in a Small Mating-Type-Specific Region of the Smut Fungus Microbotryum violaceum

DNA sequence analysis and genetic mapping of loci from mating-type-specific chromosomes of the smut fungus Microbotryum violaceum demonstrated that the nonrecombining mating-type-specific region in this species comprises ∼25% (∼1 Mb) of the chromosome length. Divergence between homologous mating-type-linked genes in this region varies between 0 and 8.6%, resembling the evolutionary strata of vertebrate and plant sex chromosomes.EVOLUTION of mating types or sex-determining systems often involves the suppression of recombination around the primary sex-determining or mating-type-determining locus. In animals and plants, it is often an entire or almost entire chromosome (Y or W in male or female heterogametic species, respectively) that ceases to recombine with its homologous (X or Z) chromosome (Charlesworth and Charlesworth 2000; Charlesworth 2008). Self-incompatibility loci in plants are also thought to be located in regions of suppressed recombination (Charlesworth et al. 2005; Kamau and Charlesworth 2005; Kamau et al. 2007; Li et al. 2007; Yang et al. 2007). Regardless of the phylogenetic position of a species, such nonrecombining regions are known to follow similar evolutionary trajectories. The nonrecombining region on the sex-specific chromosome expands in several steps, forming evolutionary strata—regions of different X/Y (or Z/W) divergence (Lahn and Page 1999; Handley et al. 2004; Sandstedt and Tucker 2004; Nicolas et al. 2005)—and genes in the nonrecombining regions gradually accumulate deleterious mutations that eventually render them dysfunctional (Charlesworth and Charlesworth 2005; Charlesworth 2008).Fungal mating-type systems are very diverse, with the number of mating types varying from two to several hundred (Casselton 2002). Like sex chromosomes in several animals and plants, suppressed recombination has evolved in regions near fungal mating-type loci, including in Ustilago hordei (Lee et al. 1999), Cryptococcus neoformans (Lengeler et al. 2002), and Neurospora tetrasperma (Menkis et al. 2008). These species have two mating types, but no morphologically distinct sexes. The mating-type locus (the region of suppressed recombination) of C. neoformans is small (∼100 kb) compared with known sex chromosomes and contains only ∼20 genes that, unlike many sex chromosomes (Y or W chromosomes), show no obvious signs of genetic degeneration (Lengeler et al. 2002; Fraser et al. 2004). Judging from the divergence between the homologous genes on the two mating-type-specific chromosomes, C. neoformans started to evolve sex chromosomes a long time ago because silent divergence between the two mating types in the most ancient region exceeds 100% (Fraser et al. 2004). Genes in the younger mating-type-specific region are much less diverged between the two sex chromosomes, suggesting that the evolution of the sex locus in C. neoformans might have proceeded through several steps. The nonrecombining region around the mating-type locus of N. tetrasperma is much larger than in C. neoformans (at least 6.6 Mb), and silent divergence between homologous genes on the mating-type-specific chromosomes ranges from zero to 9%, demonstrating that these mating-type-specific chromosomes evolved recently (Menkis et al. 2008).M. violaceum, which causes anther smut disease in Silene latifolia and other species in the family Caryophyllaceae, has two mating types, A1 and A2 (reviewed by Giraud et al. 2008), which are determined by the presence of mating-type-specific chromosomes (hereafter A1 and A2 chromosomes, or sex chromosomes) in the haploid stage of the life cycle (Hood 2002; Hood et al. 2004). The A1 and A2 chromosomes are distinguishable by size in pulsed-field electrophoresis, and it is possible to isolate individual chromosomes electrophoretically (Hood et al. 2004). Random fragments of A1 and A2 chromosomes have previously been isolated from mating-type-specific bands of pulsed-field separated chromosomes of M. violaceum (Hood et al. 2004). These fragments were assumed to be linked to mating type. The same method was used to isolate fragments of non-mating-type-specific chromosomes. On the basis of the analysis of their sequences, (Hood et al. 2004) proposed that mating-type-specific chromosomes in M. violaceum might be degenerate because they contained a lower proportion of protein-coding genes than other chromosomes. However, it was not determined whether the sequences isolated from the mating-type chromosomes originated from the mating-type-specific or from the recombining regions (Hood et al. 2004), and the relative sizes of these regions are not known for these M. violaceum chromosomes. We tested the mating-type specificity of 86 of these fragments and demonstrate that fewer than a quarter of these loci are located in the mating-type-specific region, suggesting that the nonrecombining region on the A1 and A2 chromosomes is quite small, while the rest of the chromosome probably recombines (like pseudoautosomal regions of sex chromosomes) and is therefore not expected to undergo genetic degeneration. Genetic mapping confirms the presence of two pseudoautosomal regions in the M. violaceum mating-type-specific chromosomes.As these chromosomes are mating type specific in the haploid stage of M. violaceum, mating-type-specific loci (or DNA fragments) can be identified by testing whether they are present exclusively in A1 or A2 haploid strains. We therefore prepared haploid A1 and A2 M. violaceum cultures from S. latifolia plants from two geographically remote locations (accessions Sl405 from Sweden and Sl127 from the French Pyrenees). Haploid sporidial cultures were isolated by a standard dilution method (Kaltz and Shykoff 1997; Oudemans and Alexander 1998). Mating types were determined by PCR amplification of each culture with primers designed for A1 and A2 pheromone receptor genes linked to A1 and A2 mating types (Yockteng et al. 2007). The primers were as follows: 5′-TGGCATCCCTCAATGTTTCC-3′ and 5′-CACCTTTTGATGAGAGGCCG-3′ for the A1 pheromone receptor (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"EF584742","term_id":"156254864","term_text":"EF584742"}}EF584742) and 5′-TGACGAGAGCATTCCTACCG-3′ and 5′-GAAGCGGAACTTGCCTTTCT-3′ for the A2 pheromone receptor (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"EF584741","term_id":"156254862","term_text":"EF584741"}}EF584741). Cultures with PCR product amplified only from an A1 or A2 pheromone receptor gene were selected for further use. The mating types of the cultures were verified by conjugating them in all combinations.The GenBank nucleotide database was searched using BLAST for sequences similar to those isolated by Hood et al. (2004). Sequences with similarity to transposable elements (TE) and other repeats were excluded. The resulting set of nonredundant sequences was used to design PCR primers for 98 fragments. Half of these were originally isolated from the A1 and half from the A2 chromosomes and are hereafter called A1-NNN or A2-NNN (where NNN is the locus number; supporting information, Table S1), which does not imply that these loci are A1 or A2 specific, but merely indicates that they were originally isolated from the A1 or A2 chromosomes. Amplification of these regions from new A1 and A2 M. violaceum cultures, independently isolated by ourselves, revealed that only 5 of the 49 loci isolated from the A1 chromosome are indeed A1 specific and only 6 of 49 isolated from the A2 chromosome are A2 specific. All other loci amplified from both A1 and A2 cultures. Figure 1 illustrates some of these results from the Swedish sample (Sl405).Open in a separate windowFigure 1.—Testing of mating-type specificity for loci isolated from A1 and A2 chromosomes. (a) PCR amplifications from haploid cultures from Sl405 using primers designed from six A1-originated loci. Loci in which a PCR product could be amplified only from A1 cultures (boxed) were classified as specific to mating type A1. (b) PCR tests of six A2-originated loci on the same set of haploids as in a. Loci in which a PCR product amplified only from A2 cultures (boxed) were classified as specific to mating type A2. Loci amplified from both A1 and A2 cultures are not mating type specific.The fragments that amplified from both A1 and A2 mating types may be in recombining regions, or they could be present in mating-type-specific regions on both A1 and A2 chromosomes. If they are in recombining regions, the A1- and A2-linked homologs should not be diverged from each other, but if they are in nonrecombining, mating-type-specific regions, the divergence of the A1- and A2-linked homologs should be roughly proportional to the time since recombination stopped in the region. We therefore sequenced and compared PCR fragments amplified from the two mating types of Sl405 or Sl127 cultures (GenBank accession nos. {"type":"entrez-nucleotide","attrs":{"text":"FI855822","term_id":"258612459","term_text":"FI855822"}}FI855822–{"type":"entrez-nucleotide","attrs":{"text":"FI856001","term_id":"258612570","term_text":"FI856001"}}FI856001). Sequencing of PCR products showed that 12 (4 A1 and 8 A2) loci have more than one copy, and they were excluded from further analysis. Sequences of 61 loci were identical between the A1 and A2 strains, and four loci demonstrated low total divergence (0.24–0.61%) between the two mating types (otintseva and D. Filatov, unpublished results). Thus, these loci might be located in the recombining part of the mating-type-specific chromosomes. Ten of 75 loci that amplified in both mating types demonstrated multiple polymorphisms fixed between the mating types rather than between the locations. Given that the strains that we used in the analysis originated from two geographically distant locations, it is highly unlikely that multiple polymorphisms distinguishing the A1 and A2 sequences arose purely by chance; thus, these loci are probably located in the nonrecombining mating-type-specific region of the M. violaceum A1 and A2 chromosomes.

TABLE 1

Loci from mating-type-specific chromosomes of M. violaceum used for PCR analysis and genetic map construction

		With nonzero A1/A2 divergenceb
Loci	Mating type specific	<1%	>1%	With zero A1/A2 divergenceb	Total
A1a	5	2 (1)	3 (3)	35 (3)	45 (7)
A2a	6	2 (0)	7 (7)	26 (3)	41 (10)
Subtotal		4 (1)	10 (10)
Total	11	14 (11)		61 (6)	86 (17)

Open in a separate window^aA1, loci originated from the A1 sex chromosome; A2, loci originated from the A2 sex chromosome.^bThe number of loci used for genetic map construction is in parentheses.To confirm the mating-type-specific or pseudoautosomal locations of the loci with and without A1/A2 divergence, we conducted genetic mapping in a family of 99 individuals, 50 of which were of mating type A1 and 49 of mating type A2. The family was generated by a cross between A1 and A2 M. violaceum strains from S. latifolia accessions Sl405 (Sweden) and Sl127 (France), respectively. The choice of strains from geographically distant locations was motivated by the hope of maximizing the number of DNA sequence differences between them that can be used as molecular genetic markers in segregation analysis. We inoculated S. latifolia seedlings with sporidial cultures of both mating types. For inoculation, petri dishes with 12-day-old seedlings of S. latifolia were flooded with 2.5 ml of inoculum suspension. Inoculum suspension consisted of equal volumes of the A1 and A2 sporidial cultures that were mixed and conjugated overnight at 14° under rotation (Biere and Honders 1996; Van Putten et al. 2003). Seedlings were potted 3 days after inoculation. Two months later, teliospores were collected from the flowers of the infected plant and grown in petri dishes on 3.6% potato dextrose agar medium. Haploid sporidia formed after meiosis were isolated and grown as separate cultures for DNA extraction. The mating types of single sporidia cultures were identified as described above. The loci analyzed in the segregation analysis were sequenced in the two parental haploid strains and in 99 (50 A1 and 49 A2) haploid strains that were generated in the cross. Single nucleotide differences between the parental strains were used as molecular genetic markers for segregation analysis in the progeny. The genetic map was constructed using MAPMAKER/EXP v3.0 (Lincoln et al. 1992) and MapDisto v1.7 (http://mapdisto.free.fr/).The resulting genetic map is shown in Figure 2. As expected, no recombination was observed between the 10 loci with diverged A1- and A2-linked copies. In addition, one marker with no A1/A2 divergence, A2-397, was also completely linked to the loci with significant A1/A2 divergence. This locus either may be very tightly linked to the nonrecombining mating-type-specific region or may have been added to that region more recently than the loci that had already accumulated some divergence between the alleles in the two mating types. The mating-type-specific pheromone receptor locus (Devier et al. 2009) and 11 mating-type-specific loci are also located in this nonrecombining region (Figure 2). Interestingly, the cluster of nonrecombining markers is flanked on both sides with markers that recombine in meiosis, demonstrating that there are pseudoautosomal regions on both ends of the mating-type-specific chromosomes.Open in a separate windowFigure 2.—Genetic map of the mating-type-determining chromosome in M. violaceum. Genetic distance (in centimorgans) and the relative positions of the markers are shown to the left and the right of the chromosome, respectively. The position of the nonrecombining region corresponds to the cluster of linked markers shown on the right of the figure. Total A1/A2 divergence is shown in parentheses. Eleven mating-type-specific markers (for which sequences are available from only one mating type), located in the nonrecombining mating-type-specific region, are not shown.Our results demonstrate that although the loci reported by Hood et al. (2004) were isolated from the A1 and A2 chromosomes, most of these loci are not located in the nonrecombining mating-type-specific regions. In fact, the nonrecombining region might be relatively small: of 86 tested fragments, only 21 appeared to be either mating type specific or linked to the mating-type locus. Assuming that these loci represent a random set of DNA fragments isolated from the A1 and A2 chromosomes, it is possible to estimate the size of the nonrecombining region using the binomial distribution: the nonrecombining region is expected to be 24.4% (95% CI: 16.7–33.6%) of the chromosome length. As the sizes of the A1 and A2 chromosomes are ∼3.4 and 4.2 Mb long (Hood 2002; Hood et al. 2004), the nonrecombining region might be ∼1 Mb long.Interestingly, total A1/A2 divergence for the 11 loci with A1- and A2-linked copies mapped to the nonrecombining region varied from 0% to 8.6% (Figure 2). In addition, 11 loci amplified from only one mating type. These genes could represent degenerated genes, some of which degenerated in A1 strains, and some in A2 strains. Alternatively, they might be highly diverged genes, such that the PCR primers amplify only one allele, and not the other. Variation in divergence may be the result of the stepwise cessation of recombination between the A1 and A2 chromosomes in M. violaceum, resembling the evolutionary strata reported for human, chicken, and white campion sex chromosomes (Lahn and Page 1999; Handley et al. 2004; Bergero et al. 2007). However, only the differences between the most and the least diverged loci are statistically significant (Devier et al. 2009), the M. violaceum mating-type region has at least three strata: one oldest stratum, including the pheromone receptor locus; a younger stratum with ∼5–9% A1/A2 divergence; and the youngest stratum with 1–4% divergence between the two mating types. There may also be an additional very recently evolved stratum containing the locus named A2-397, which is also present in all A1 strains tested, with no fixed differences between the A1 and A2 strains (

Characterization of a Thermostable Short-Chain Alcohol Dehydrogenase from the Hyperthermophilic Archaeon Thermococcus sibiricus

Short-chain alcohol dehydrogenase, encoded by the gene Tsib_0319 from the hyperthermophilic archaeon Thermococcus sibiricus, was expressed in Escherichia coli, purified and characterized as an NADPH-dependent enantioselective oxidoreductase with broad substrate specificity. The enzyme exhibits extremely high thermophilicity, thermostability, and tolerance to organic solvents and salts.Alcohol dehydrogenases (ADHs; EC 1.1.1.1.) catalyze the interconversion of alcohols to their corresponding aldehydes or ketones by using different redox-mediating cofactors. NAD(P)-dependent ADHs, due to their broad substrate specificity and enantioselectivity, have attracted particular attention as catalysts in industrial processes (5). However, mesophilic ADHs are unstable at high temperatures, sensitive to organic solvents, and often lose activity during immobilization. In this relation, there is a considerable interest in ADHs from extremophilic microorganisms; among them, Archaea are of great interest. The representatives of all groups of NAD(P)-dependent ADHs have been detected in genomes of Archaea (11, 12); however, only a few enzymes have been characterized, and the great majority of them belong to medium-chain (3, 4, 14, 16, 19) or long-chain iron-activated ADHs (1, 8, 9). Up to now, a single short-chain archaeal ADH from Pyrococcus furiosus (10, 18) and only one archaeal aldo-keto reductase also from P. furiosus (11) have been characterized.Thermococcus sibiricus is a hyperthermophilic anaerobic archaeon isolated from a high-temperature oil reservoir capable of growth on complex organic substrates (15). The complete genome sequence of T. sibiricus has been recently determined and annotated (13). Several ADHs are encoded by the T. sibiricus genome, including three short-chain ADHs (Tsib_0319, Tsib_0703, and Tsib_1998) (13). In this report, we describe the cloning and expression of the Tsib_0319 gene from T. sibiricus and the purification and the biochemical characterization of its product, the thermostable short-chain ADH (TsAdh319).The Tsib_0319 gene encodes a protein with a size of 234 amino acids and the calculated molecular mass of 26.2 kDa. TsAdh319 has an 85% degree of sequence identity with short-chain ADH from P. furiosus (AdhA; PF_0074) (18). Besides AdhA, close homologs of TsAdh319 were found among different bacterial ADHs, but not archaeal ADHs. The gene flanked by the XhoI and BamHI sites was PCR amplified using two primers (sense primer, 5′-GTTCTCGAGATGAAGGTTGCTGTGATAACAGGG-3′, and antisense primer, 5′-GCTGGATCCTCAGTATTCTGGTCTCTGGTAGACGG-3′) and cloned into the pET-15b vector. TsAdh319 was overexpressed, with an N-terminal His₆ tag in Escherichia coli Rosetta-gami (DE3) and purified to homogeneity by metallochelating chromatography (Hi-Trap chelating HP column; GE Healthcare) followed by gel filtration on Superdex 200 10/300 GL column (GE Healthcare) equilibrated in 50 mM Tris-HCl (pH 7.5) with 200 mM NaCl. The homogeneity and the correspondence to the calculated molecular mass of 28.7 kDa were verified by SDS-PAGE (7). The molecular mass of native TsAdh319 was 56 to 60 kDa, which confirmed the dimeric structure in solution.The standard ADH activity measurement was made spectrophotometrically at the optimal pH by following either the reduction of NADP (in 50 mM Gly-NaOH buffer; pH 10.5) or the oxidation of NADPH (in 0.1 M sodium phosphate buffer; pH 7.5) at 340 nm at 60°C. The enzyme exhibited a strong preference for NADP(H) and broad substrate specificity (Table ​(Table1).1). The highest oxidation rates were found with pentoses d-arabinose (2.0 U mg⁻¹) and d-xylose (2.46 U mg⁻¹), and the highest reduction rates were found with dimethylglyoxal (5.9 U mg⁻¹) and pyruvaldehyde (2.2 U mg⁻¹). The enzyme did not reduce sugars which were good substrates for the oxidation reaction. The kinetic parameters of TsAdh319 determined for the preferred substrates are shown in Table ​Table2.2. The enantioselectivity of the enzyme was estimated by measuring the conversion rates of 2-butanol enantiomers. TsAdh319 showed an evident preference, >2-fold, for (S)-2-butanol over (RS)-2-butanol. The enzyme stereoselectivity is confirmed by the preferred oxidation of d-arabinose over l-arabinose (Table ​(Table1).1). The fact that TsAdh319 is metal independent was supported by the absence of a significant effect of TsAdh319 preincubation with 10 mM Me²⁺ for 30 min before measuring the activity in the presence of 1 mM Me²⁺ or EDTA (Table ​(Table3).3). TsAdh319 also exhibited a halophilic property, so the enzyme activity increased in the presence of NaCl and KCl and the activation was maintained even at concentration of 4 M and 3 M, respectively (Table ​(Table33).

TABLE 1.

Substrate specificity of TsAdh319

A Targeted Multilocus Genotyping Assay for Lineage,Serogroup, and Epidemic Clone Typing of Listeria monocytogenes

A 30-probe assay was developed for simultaneous classification of Listeria monocytogenes isolates by lineage (I to IV), major serogroup (4b, 1/2b, 1/2a, and 1/2c), and epidemic clone (EC) type (ECI, ECIa, ECII, and ECIII). The assay was designed to facilitate rapid strain characterization and the integration of subtype data into risk-based inspection programs.Listeria monocytogenes is a facultative intracellular pathogen that can cause serious invasive illness (listeriosis) in humans and other animals. L. monocytogenes is responsible for over 25% of food-borne-disease-related deaths attributable to known pathogens and is a leading cause of food recalls due to microbial adulteration (12, 21). However, not all L. monocytogenes subtypes contribute equally to human illness, and substantial differences in the ecologies and virulence attributes of different L. monocytogenes subtypes have been identified (9, 13, 14, 23, 24, 33, 35, 36). Among the four major evolutionary lineages of L. monocytogenes, only lineages I and II are commonly isolated from contaminated food and human listeriosis patients (19, 27, 29, 33). Lineage I strains are overrepresented among human listeriosis isolates, particularly those associated with epidemic outbreaks, whereas lineage II strains are overrepresented in foods and the environment (13, 14, 24). Lineage III strains account for approximately 1% of human listeriosis cases but are common among animal listeriosis isolates and appear to be a host-adapted group that is poorly adapted to food-processing environments (6, 34-36). The ecological and virulence attributes of lineage IV are poorly understood, as this lineage is rare and was only recently described based on a small number of strains (19, 26, 29, 33).L. monocytogenes is differentiated into 13 serotypes; however, four major serogroups (4b, 1/2b, 1/2a, and 1/2c) from within lineages I and II account for more than 98% of human and food isolates (16, 31). Serogroups refer to evolutionary complexes typified by a predominant serotype but which include very rare serotypes that represent minor evolutionary variants (7, 9, 33). Phylogenetic analyses have indicated that rare serotypes may have evolved recently, or even multiple times, from one of the major serotypes (9), and numerous molecular methods fail to discriminate minor serotypes as independent groups (1, 4, 7, 9, 18, 22, 33, 38, 39). Serotyping is one of the most common methods for L. monocytogenes subtyping, and serogroup classifications are a useful component of strain characterization because ecotype divisions appear largely congruent with serogroup distinctions (16, 34). Serogroup 4b strains are of particular public health concern because contamination with these strains appears to increase the probability that a ready-to-eat (RTE) food will be implicated in listeriosis (16, 28). Serogroup 4b strains account for approximately 40% of sporadic listeriosis and also are responsible for the majority of listeriosis outbreaks despite being relatively rare contaminants of food products (9, 13, 17, 30, 34). In addition, serogroup 4b strains are associated with more severe clinical presentations and higher mortality rates than other serogroups (11, 16, 20, 31, 34). Serogroups 1/2a and 1/2b are overrepresented among food isolates but also contribute significantly to human listeriosis, whereas serogroup 1/2c rarely causes human illness and may pose a lower risk of listeriosis for humans (16). Serogroup-specific differences in association with human listeriosis are consistent with the prevalence of virulence-attenuating mutations in inlA within these serogroups (32, 34); however, a number of additional factors likely contribute to these differences.Four previously described epidemic clones (ECs; ECI, ECIa, ECII, and ECIII) of L. monocytogenes have been implicated in numerous listeriosis outbreaks and have contributed significantly to sporadic illness (15, 34). ECI, ECIa, and ECII are distinct groups within serogroup 4b that were each responsible for repeated outbreaks of listeriosis in the United States and Europe. ECIII is a lineage II clone of serotype 1/2a that persisted in the same processing facility for more than a decade prior to causing a multistate outbreak linked to contaminated turkey (15, 25). While there has been speculation that epidemic clones possess unique adaptations that explain their frequent involvement in listeriosis outbreaks (9, 34, 37), it is not clear that epidemic clones are more virulent than other strains with the same serotype. However, contamination of RTE food with EC strains would be cause for increased concern due to the previous involvement of these clones in major outbreaks of listeriosis (16).As a result of the L. monocytogenes subtype-specific differences in ecology, virulence, and association with human illness, molecular subtyping technologies have the potential to inform assessments of relative risk and to improve risk-based inspection programs. The objective of the present study was to develop a single assay for rapid and accurate classification of L. monocytogenes isolates by lineage, major serogroup, and epidemic clone in order to facilitate strain characterization and the integration of subtype data into inspection programs that are based on assessment of relative risk.A database of more than 5.3 Mb of comparative DNA sequences from 238 L. monocytogenes isolates (9, 33-35) was scanned for single nucleotide polymorphisms that could be used to differentiate lineages, major serogroups, and epidemic clones via a targeted multilocus genotyping (TMLGT) approach. The acronym TMLGT is used to distinguish this approach from previously published multilocus genotyping (MLGT) assays that were lineage specific and designed for haplotype discrimination (9, 33). To provide for simultaneous interrogation of the selected polymorphisms via TMLGT, six genomic regions (Table ​(Table1)1) were coamplified in a multiplex PCR. While the previous MLGT assays were based on three lineage-specific multiplexes and required prior identification of lineage identity, TMLGT was designed to target variation across all of the lineages simultaneously and is based on a unique set of amplicons. PCR was performed in 50-μl volumes with 1× High Fidelity PCR buffer (Invitrogen Life Technologies), 2 mM MgSO₄, 100 μM deoxynucleoside triphosphate (dNTP), 300 nM primer, 1.5 U Platinum Taq high-fidelity DNA polymerase (Invitrogen Life Technologies), and 100 ng of genomic DNA. PCR consisted of an initial denaturation of 90 s at 96°C, followed by 40 cycles of 30 s at 94°C, 30 s at 50°C, and 90 s at 68°C. Amplification products were purified using Montage PCR cleanup filter plates (Millipore) and served as a template for allele-specific primer extension (ASPE) reactions utilizing subtype-specific probes.

TABLE 1.

Primers used in multiplex amplification for the TMLGT assay