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

Searching for Mesophilic Thermotogales Bacteria: “Mesotogas” in the Wild

All cultivated Thermotogales are thermophiles or hyperthermophiles. However, optimized 16S rRNA primers successfully amplified Thermotogales sequences from temperate hydrocarbon-impacted sites, mesothermic oil reservoirs, and enrichment cultures incubated at <46°C. We conclude that distinct Thermotogales lineages commonly inhabit low-temperature environments but may be underreported, likely due to “universal” 16S rRNA gene primer bias.Thermotogales, a bacterial group in which all cultivated members are anaerobic thermophiles or hyperthermophiles (5), are rarely detected in anoxic mesothermic environments, yet their presence in corresponding enrichment cultures, bioreactors, and fermentors has been observed using metagenomic methods and 16S rRNA gene amplification (6) (see Table S1 in the supplemental material). The most commonly detected lineage is informally designated here “mesotoga M1” (see Table S1 in the supplemental material). PCR experiments indicated that mesotoga M1 sequences amplified inconsistently using “universal” 16S rRNA gene primers, perhaps explaining their poor detection in DNA isolated from environmental samples (see text and Table S2 in the supplemental material). We therefore designed three 16S rRNA PCR primer sets (Table ​(Table1)1) targeting mesotoga M1 bacteria and their closest cultivated relative, Kosmotoga olearia. Primer set A was the most successful set, detecting a wider diversity of Thermotogales sequences than set B and being more Thermotogales-specific than primer set C (Table ​(Table22).

TABLE 1.

Primers targeting mesotoga M1 bacteria constructed and used in this study

Phage Display Selection of Cyclic Peptides That Inhibit Andes Virus Infection

Specific therapy is not available for hantavirus cardiopulmonary syndrome caused by Andes virus (ANDV). Peptides capable of blocking ANDV infection in vitro were identified using antibodies against ANDV surface glycoproteins Gn and Gc to competitively elute a cyclic nonapeptide-bearing phage display library from purified ANDV particles. Phage was examined for ANDV infection inhibition in vitro, and nonapeptides were synthesized based on the most-potent phage sequences. Three peptides showed levels of viral inhibition which were significantly increased by combination treatment with anti-Gn- and anti-Gc-targeting peptides. These peptides will be valuable tools for further development of both peptide and nonpeptide therapeutic agents.Andes virus (ANDV), an NIAID category A agent linked to hantavirus cardiopulmonary syndrome (HCPS), belongs to the family Bunyaviridae and the genus Hantavirus and is carried by Oligoryzomys longicaudatus rodents (11). HCPS is characterized by pulmonary edema caused by capillary leak, with death often resulting from cardiogenic shock (9, 16). ANDV HCPS has a case fatality rate approaching 40%, and ANDV is the only hantavirus demonstrated to be capable of direct person-to-person transmission (15, 21). There is currently no specific therapy available for treatment of ANDV infection and HCPS.Peptide ligands that target a specific protein surface can have broad applications as therapeutics by blocking specific protein-protein interactions, such as preventing viral engagement of host cell receptors and thus preventing infection. Phage display libraries provide a powerful and inexpensive tool to identify such peptides. Here, we used selection of a cyclic nonapeptide-bearing phage library to identify peptides capable of binding the transmembrane surface glycoproteins of ANDV, Gn and Gc, and blocking infection in vitro.To identify peptide sequences capable of recognizing ANDV, we panned a cysteine-constrained cyclic nonapeptide-bearing phage display library (New England Biolabs) against density gradient-purified, UV-treated ANDV strain CHI-7913 (a gift from Hector Galeno, Santiago, Chile) (17, 18). To increase the specificity of the peptides identified, we eluted phage by using monoclonal antibodies (Austral Biologicals) prepared against recombinant fragments of ANDV Gn (residues 1 to 353) or Gc (residues 182 to 491) glycoproteins (antibodies 6B9/F5 and 6C5/D12, respectively). Peptide sequences were determined for phage from iterative rounds of panning, and the ability of phage to inhibit ANDV infection of Vero E6 cells was determined by immunofluorescent assay (IFA) (7). Primary IFA detection antibodies were rabbit polyclonal anti-Sin Nombre hantavirus (SNV) nucleoprotein (N) antibodies which exhibit potent cross-reactivity against other hantavirus N antigens (3). ReoPro, a commercially available Fab fragment which partially blocks infection of hantaviruses in vitro by binding the entry receptor integrin β₃ (5), was used as a positive control (80 μg/ml) along with the original antibody used for phage elution (5 μg/ml). As the maximum effectiveness of ReoPro in inhibiting hantavirus entry approaches 80%, we set this as a threshold for maximal expected efficacy for normalization. The most-potent phage identified by elution with the anti-Gn antibody 6B9/F5 bore the peptide CPSNVNNIC and inhibited hantavirus entry by greater than 60% (61%) (Table ​(Table1).1). From phage eluted with the anti-Gc antibody 6C5/D12, those bearing peptides CPMSQNPTC and CPKLHPGGC also inhibited entry by greater than 60% (66% and 72%, respectively).

TABLE 1.

Peptide-bearing phage eluted from ANDV

Impact of Viral Dose and Major Histocompatibility Complex Class IB Haplotype on Viral Outcome in Mauritian Cynomolgus Monkeys Vaccinated with Tat upon Challenge with Simian/Human Immunodeficiency Virus SHIV89.6P

The effects of the challenge dose and major histocompatibility complex (MHC) class IB alleles were analyzed in 112 Mauritian cynomolgus monkeys vaccinated (n = 67) or not vaccinated (n = 45) with Tat and challenged with simian/human immunodeficiency virus (SHIV) 89.6P_cy243. In the controls, the challenge dose (10 to 20 50% monkey infectious doses [MID₅₀]) or MHC did not affect susceptibility to infection, peak viral load, or acute CD4 T-cell loss, whereas in the chronic phase of infection, the H1 haplotype correlated with a high viral load (P = 0.0280) and CD4 loss (P = 0.0343). Vaccination reduced the rate of infection acquisition at 10 MID₅₀ (P < 0.0001), and contained acute CD4 loss at 15 MID₅₀ (P = 0.0099). Haplotypes H2 and H6 were correlated with increased susceptibility (P = 0.0199) and resistance (P = 0.0087) to infection, respectively. Vaccination also contained CD4 depletion (P = 0.0391) during chronic infection, independently of the challenge dose or haplotype.Advances in typing of the major histocompatibility complex (MHC) of Mauritian cynomolgus macaques (14, 20, 26) have provided the opportunity to address the influence of host factors on vaccine studies (13). Retrospective analysis of 22 macaques vaccinated with Tat or a Tat-expressing adenoviral vector revealed that monkeys with the H6 or H3 MHC class IB haplotype were overrepresented among aviremic or controller animals, whereas macaques with the H2 or H5 haplotype clustered in the noncontrollers (12). More recently, the H6 haplotype was reported to correlate with control of chronic infection with simian immunodeficiency virus (SIV) mac251, regardless of vaccination (18).Here, we performed a retrospective analysis of 112 Mauritian cynomolgus macaques, which included the 22 animals studied previously (12), to evaluate the impact of the challenge dose and class IB haplotype on the acquisition and severity of simian/human immunodeficiency virus (SHIV) 89.6P_cy243 infection in 45 control monkeys and 67 monkeys vaccinated with Tat from different protocols (Table ​(Table11).

TABLE 1.

Summary of treatment, challenge dose, and outcome of infection in cynomolgus monkeys

Association of Campylobacter jejuni Cj0859c Gene (fspA) Variants with Different C. jejuni Multilocus Sequence Types

Cj0859c variants fspA1 and fspA2 from 669 human, poultry, and bovine Campylobacter jejuni strains were associated with certain hosts and multilocus sequence typing (MLST) types. Among the human and poultry strains, fspA1 was significantly (P < 0.001) more common than fspA2. FspA2 amino acid sequences were the most diverse and were often truncated.Campylobacter jejuni is the leading cause of bacterial gastroenteritis worldwide and responsible for more than 90% of Campylobacter infections (7). Case-control studies have identified consumption or handling of raw and undercooked poultry meat, drinking unpasteurized milk, and swimming in natural water sources as risk factors for acquiring domestic campylobacteriosis in Finland (7, 9). Multilocus sequence typing (MLST) has been employed to study the molecular epidemiology of Campylobacter (4) and can contribute to virulotyping when combined with known virulence factors (5). FspA proteins are small, acidic, flagellum-secreted nonflagellar proteins of C. jejuni that are encoded by Cj0859c, which is expressed by a σ²⁸ promoter (8). Both FspA1 and FspA2 were shown to be immunogenic in mice and protected against disease after challenge with a homologous strain (1). However, FspA1 also protected against illness after challenge with a heterologous strain, whereas FspA2 failed to do the same at a significant level. Neither FspA1 nor FspA2 protected against colonization (1). On the other hand, FspA2 has been shown to induce apoptosis in INT407 cells, a feature not exhibited by FspA1 (8). Therefore, our aim was to study the distributions of fspA1 and fspA2 among MLST types of Finnish human, chicken, and bovine strains.In total, 367 human isolates, 183 chicken isolates, and 119 bovine isolates (n = 669) were included in the analyses (3). PCR primers for Cj0859c were used as described previously (8). Primer pgo6.13 (5′-TTGTTGCAGTTCCAGCATCGGT-3′) was designed to sequence fspA1. Fisher''s exact test or a chi-square test was used to assess the associations between sequence types (STs) and Cj0859c. The SignalP 3.0 server was used for prediction of signal peptides (2).The fspA1 and fspA2 variants were found in 62.6% and 37.4% of the strains, respectively. In 0.3% of the strains, neither isoform was found. Among the human and chicken strains, fspA1 was significantly more common, whereas fspA2 was significantly more frequent among the bovine isolates (Table ​(Table1).1). Among the MLST clonal complexes (CCs), fspA1 was associated with the ST-22, ST-45, ST-283, and ST-677 CCs and fspA2 was associated with the ST-21, ST-52, ST-61, ST-206, ST-692, and ST-1332 CCs and ST-58, ST-475, and ST-4001. Although strong CC associations of fspA1 and fspA2 were found, the ST-48 complex showed a heterogeneous distribution of fspA1 and fspA2. Most isolates carried fspA2, and ST-475 was associated with fspA2. On the contrary, ST-48 commonly carried fspA1 (Table ​(Table1).1). In our previous studies, ST-48 was found in human isolates only (6), while ST-475 was found in both human and bovine isolates (3, 6). The strict host associations and striking difference between fspA variants in human ST-48 isolates and human/bovine ST-475 isolates suggest that fspA could be important in host adaptation.

TABLE 1.

Percent distributions of fspA1 and fspA2 variants among 669 human, poultry, and bovine Campylobacter jejuni strains and their associations with hosts, STs, and CCs

Sequences from Ancestral Single-Stranded DNA Viruses in Vertebrate Genomes: the Parvoviridae and Circoviridae Are More than 40 to 50 Million Years Old

Vertebrate genomic assemblies were analyzed for endogenous sequences related to any known viruses with single-stranded DNA genomes. Numerous high-confidence examples related to the Circoviridae and two genera in the family Parvoviridae, the parvoviruses and dependoviruses, were found and were broadly distributed among 31 of the 49 vertebrate species tested. Our analyses indicate that the ages of both virus families may exceed 40 to 50 million years. Shared features of the replication strategies of these viruses may explain the high incidence of the integrations.It has long been appreciated that retroviruses can contribute significantly to the genetic makeup of host organisms. Genes related to certain other viruses with single-stranded RNA genomes, formerly considered to be most unlikely candidates for such contribution, have recently been detected throughout the vertebrate phylogenetic tree (1, 6, 13). Here, we report that viruses with single-stranded DNA (ssDNA) genomes have also contributed to the genetic makeup of many organisms, stretching back as far as the Paleocene period and possibly the late Cretaceous period of evolution.Determining the evolutionary ages of viruses can be problematic, as their mutation rates may be high and their replication may be rapid but also sporadic. To establish a lower age limit for currently circulating ssDNA viruses, we analyzed 49 published vertebrate genomic assemblies for the presence of sequences derived from the NCBI RefSeq database of 2,382 proteins from known viruses in this category, representing a total of 23 classified genera from 7 virus families. Our survey uncovered numerous high-confidence examples of endogenous sequences related to the Circoviridae and to two genera in the family Parvoviridae: the parvoviruses and dependoviruses (Fig. ​(Fig.11).Open in a separate windowFIG. 1.Phylogenetic tree of vertebrate organisms and history of ssDNA virus integrations. Times of integration of ancestral dependoviruses (yellow icosahedrons), parvoviruses (blue icosahedrons), and circoviruses (triangles) are approximate.The Dependovirus and Parvovirus genomes are typically 4 to 6 kb in length, include 2 major open reading frames (encoding replicase proteins [Rep and NS1, respectively] and capsid proteins [Cap and VP1, respectively]), and have characteristic hairpin structures at both ends (Fig. ​(Fig.2).2). For replication, these viruses depend on host enzymes that are recruited by the viral replicase proteins to the hairpin regions, where self-primed viral DNA synthesis is initiated (2). Circovirus genomes are typically ∼2-kb circles. DNA of the type species, porcine circovirus 1 (PCV-1), contains a stem-loop structure within the origin of replication (Fig. ​(Fig.2),2), and the largest open reading frame includes sequences that are homologous to the Parvovirus replicase open reading frame (9, 11). The circoviruses also depend on host enzymes for replication, and DNA synthesis is self-primed from a 3′-OH end formed by endonucleolytic cleavage of the stem-loop structure (4). The frequency of Dependovirus infection is estimated to be as high as 90% within an individual''s lifetime. None of the dependoviruses have been associated with human disease, but related viruses in the family Parvoviridae (e.g., erythrovirus B19 and possibly human bocavirus) are pathogenic for humans, and members of both the Parvoviridae and the Circoviridae can cause a variety of animal diseases (2, 4).Open in a separate windowFIG. 2.Schematics illustrating the structure and organization of Parvoviridae and Circoviridae genomes and origins of several of the longest-integrated ancestral viral sequences found in vertebrates. Integrations were aligned to the Dependovirus adeno-associated virus 2 (AAV2), the Parvovirus minute virus of mice (MVM), and the Circovirus porcine circovirus 1 (PCV-1). The inverted terminal repeat (ITR) sequences in the Dependovirus and Parvovirus genomes are depicted on an expanded scale. A linear representation of the circular genome of PCV-1 is shown with the 10-bp stem-loop structure on an expanded scale. Horizontal lines beneath the maps indicate the lengths of similar sequences that could be identified by BLAST. The numbers indicate the locations of amino acids in the viral proteins where the sequence similarities in the endogenous insertions start and end. The actual ancestral virus-derived integrated sequences may extend beyond the indicated regions.With some ancestral endogenous sequences that we identified, phylogenetic comparisons can be used to estimate age. For example, as a Dependovirus-like sequence is present at the same location in the genomes of mice and rats, the ancestral virus must have existed before their divergence, more than 20 million years ago. Some Circovirus- and Dependovirus-related integrations also predate the split between dog and panda, about 42 million years ago. However, in most other cases, we rely on an indirect method for estimating age (1). As genomic sequences evolve, they accumulate new stop codons and insertion/deletion-induced frameshifts. The rates of these events can be tied directly to the rates of neutral sequence drift and, therefore, the time of evolution. To apply this method, we first performed a BLAST search of vertebrate genomes for all known ssDNA virus proteins (BLAST options, -p tblastn -M BLOSUM62 -e 1e−4). Candidate sequences were then recorded, along with 5 kb of flanking regions, and then again aligned against the database of ssDNA viruses to find the most complete alignment (BLAST options, -t blastx -F F -w 15 -t 1500 -Z 150 -G 13 -E 1 -e 1e−2). Detected alignments were then compared with a neutral model of genome evolution, as described in the supplemental material, and the numbers of stop codons and frameshifts were converted into the expected genomic drift undergone by the sequences. The age of integration was then estimated from the known phylogeny of vertebrates (7, 10). Using these methods, we discovered that as many as 110 ssDNA virus-related sequences have been integrated into the 49 vertebrate genomes considered, during a time period ranging from the present to over 40 to 60 million years ago (Table ​(Table1;1; see also Tables S1 to S3 in the supplemental material).

TABLE 1.

Selected endogenous sequences in vertebrate genomes related to single-stranded DNA viruses

Dehalogenation Activities and Distribution of Reductive Dehalogenase Homologous Genes in Marine Subsurface Sediments

Halogenated organic compounds serve as terminal electron acceptors for anaerobic respiration in a diverse range of microorganisms. Here, we report on the widespread distribution and diversity of reductive dehalogenase homologous (rdhA) genes in marine subsurface sediments. A total of 32 putative rdhA phylotypes were detected in sediments from the southeast Pacific off Peru, the eastern equatorial Pacific, the Juan de Fuca Ridge flank off Oregon, and the northwest Pacific off Japan, collected at a maximum depth of 358 m below the seafloor. In addition, significant dehalogenation activity involving 2,4,6-tribromophenol and trichloroethene was observed in sediment slurry from the Nankai Trough Forearc Basin. These results suggest that dehalorespiration is an important energy-yielding pathway in the subseafloor microbial ecosystem.Scientific ocean drilling explorations have revealed that marine subsurface sediments harbor remarkable numbers of microbial cells that account for approximately 1/10 to 1/3 of all living biota on Earth (20, 25, 33). Thermodynamic calculations of pore-water chemistry suggest that subseafloor microbial activities are generally supported by nutrient and energy supplies from the seawater and/or underlying basaltic aquifers (6, 7). Although sulfate, nitrate, Fe(III), Mn(IV), and bicarbonate are known to be potential electron acceptors for anaerobic microbial respiration in marine subsurface sediments (5), the incidence of both the dissimilatory dehalorespiration pathway and microbial activity in halogenated organic substrates remains largely unknown.Previous molecular ecological studies using 16S rRNA gene sequences demonstrated that Chloroflexi is one of the most frequently detected phyla in subseafloor sediments of the Pacific Ocean margins (12-14). Some of the sequences within the Chloroflexi are closely related to sequences in the genus Dehalococcoides, which contains obligatory dehalorespiring bacteria that employ halogenated organic compounds as terminal electron acceptors (21, 29). The frequent detection of Dehalococcoides-related 16S rRNA genes from these environments implies the occurrence of dissimilatory dehalorespiration in marine subsurface sediments.In this study, we detected and phylogenetically analyzed the reductive dehalogenase homologous (rdhA) genes, key functional genes for dehalorespiration pathways, from frozen sediment core samples obtained by Ocean Drilling Program (ODP) Leg 201 (Peru margin and eastern equatorial Pacific) (7, 14); Integrated Ocean Drilling Program (IODP) Expedition 301 (Juan de Fuca Ridge flank) (8, 24); Chikyu Shakedown Expedition CK06-06 (Northwest Pacific off Japan) (20, 23); and IODP Expedition 315 (Nankai Trough Forearc Basin off Japan) (Table ​(Table1).1). DNA was extracted using an ISOIL bead-beating kit (Nippon Gene, Japan) and purified using a MagExtractor DNA fragment purification kit (Toyobo, Japan) according to the manufacturer''s instructions. To increase concentration, DNA was amplified by multiple displacement amplification using the phi29 polymerase supplied with a GenomiPhi kit (GE Healthcare, United Kingdom) (20). Putative rdhA genes were amplified by PCR using Ex Taq polymerase (TaKaRa, Japan) with degenerate primers RRF2 and B1R (17), dehaloF3, dehaloF4, dehaloF5, dehaloR2, dehaloR3, and dehaloR4 (32), and ceRD2S, ceRD2L, and RD7 (26) and the PCR conditions described in those studies. Amplicons of the approximate target size were gel purified and cloned into the pCR2.1 vector (Invitrogen, Japan). Sequence similarity was analyzed using FastGroupII web-based software (34), and sequences with a 95% identity were tentatively assigned to the same phylotype. Amino acid sequences were aligned by ClustalW (31), including known and putative reductive dehalogenase sequences in the genome of Dehalococcoides ethenogenes strain 195 (28), as well as several functionally characterized reductive dehalogenases from other species.

TABLE 1.

Sample locations and results of PCR amplification of rdhA

Environmental Isolates of Burkholderia pseudomallei in Ceará State,Northeastern Brazil

Melioidosis has been considered an emerging disease in Brazil since the first cases were reported to occur in the northeast region. This study investigated two municipalities in Ceará state where melioidosis cases have been confirmed to occur. Burkholderia pseudomallei was isolated in 26 (4.3%) of 600 samples in the dry and rainy seasons.Melioidosis is an endemic disease in Southeast Asia and northern Australia (2, 4) and also occurs sporadically in other parts of the world (3, 7). Human melioidosis was reported to occur in Brazil only in 2003, when a family outbreak afflicted four sisters in the rural part of the municipality of Tejuçuoca, Ceará state (14). After this episode, there was one reported case of melioidosis in 2004 in the rural area of Banabuiú, Ceará (14). And in 2005, a case of melioidosis associated with near drowning after a car accident was confirmed to occur in Aracoiaba, Ceará (11).The goal of this study was to investigate the Tejuçuoca and Banabuiú municipalities, where human cases of melioidosis have been confirmed to occur, and to gain a better understanding of the ecology of Burkholderia pseudomallei in this region.We chose as central points of the study the residences and surrounding areas of the melioidosis patients in the rural areas of Banabuiú (5°18′35″S, 38°55′14″W) and Tejuçuoca (03°59′20″S, 39°34′50′W) (Fig. ​(Fig.1).1). There are two well-defined seasons in each of these locations: one rainy (running from January to May) and one dry (from June to December). A total of 600 samples were collected at five sites in Tejuçuoca (T1, T2, T3, T4, and T5) and five in Banabuiú (B1, B2, B3, B4, and B5), distributed as follows (Fig. ​(Fig.2):2): backyards (B1 and T1), places shaded by trees (B2 and T2), water courses (B3 and T3), wet places (B4 and T4), and stock breeding areas (B5 and T5).Open in a separate windowFIG. 1.Municipalities of Banabuiú (5°18′35″S, 38°55′14″W) and Tejuçuoca (03°59′20″S, 39°34′50″W).Open in a separate windowFIG. 2.Soil sampling sites in Banabuiú and Tejuçuoca.Once a month for 12 months (a complete dry/rainy cycle), five samples were gathered at five different depths: at the surface and at 10, 20, 30 and 40 cm (Table ​(Table1).1). The samples were gathered according to the method used by Inglis et al. (9). Additionally, the sample processing and B. pseudomallei identification were carried out as previously reported (1, 8, 9).

TABLE 1.

Distribution of samples with isolates by site and soil depth