Clade-Specific Evolution Mediated by HLA-B*57/5801 in Human Immunodeficiency Virus Type 1 Clade A1 p24

HLA-B*57-mediated selection pressure leads to a typical escape pathway in human immunodeficiency virus type 1 (HIV-1) CD8 epitopes such as TW10. Whether this T242N pathway is shared by all clades remains unknown. We therefore assessed the nature of HLA-B*57 selection in a large, observational Kenyan cohort where clades A1 and D predominate. While T242N was ubiquitous in clade D HLA-B*57⁺ subjects, this mutation was rare (15%) in clade A1. Instead, P243T and I247L were selected by clade A1-infected HLA-B*57 subjects but not by HLA-B*5801⁺ subjects. Our data suggest that clade A1 consensus proline at Gag residue 243 might represent an inherent block to T242N escape in clade A1. We confirmed immunologically that P243T and I247L likely represent escape mutations. HLA-B*57 evolution also differed between clades in the KF11 and IW9 epitopes. A better understanding of clade-specific evolution is important for the development of HIV vaccines in regions with multiple clades.Human immunodeficiency virus type 1 (HIV-1) displays extreme genetic diversity, with nine clades (subtypes) described in group M, and frequent genomic recombination among and within the clades (7, 44). HIV is also capable of rapid evolution, which can lead to mutational escape from immune control (43). Escape from CD8⁺ T-cell responses occurs frequently in HIV-1 infection through mutations that affect epitope processing, HLA class I binding, and/or T-cell receptor recognition (23). In early HIV-1 infection, the majority of amino acid substitutions are associated HLA class I alleles (1). The timing and consequences of mutational escape from CD8⁺ T-cell responses vary considerably (8, 22).HLA-B*57, and to a lesser extent HLA-B*5801, has been associated with slower progression to AIDS in several studies (18, 27, 39), and HLA-B*5701 was associated with a lower viremia set point in a genome-wide association study (16). Several attributes of HLA-B*57-restricted CD8⁺ T-cell responses may contribute to their protectiveness, including dominant responses in acute infection (2), recognition of protective epitopes in HIV-1 p24 (33), better recognition of epitope variation (45), and retention of proliferative capability in chronic infection (24).HLA-B*57/5801 also exert powerful selection pressure on HIV to avoid CD8⁺ T-cell recognition. This was first demonstrated in the HLA-B*57-restricted TW10 epitope (TSTLQEQIGW [Gag_240-249]), which accounts for >30% of overall HIV-specific CD8⁺ T-cell responses in acutely infected HLA-B*57⁺ subjects (3). Escape in this epitope usually occurs early in infection, which coincidently is when HLA-B*57 is most protective (18). In clade B and C infections, >75 to 100% of HLA-B*57/5801⁺ subjects develop the T242N escape mutation, while HLA-B*57/5701-negative subjects rarely display polymorphism at this residue (5, 9, 10, 15, 32, 35, 38, 41). When T242N is transmitted to HLA-B*57/5801-negative subjects, it rapidly reverts to the consensus, suggesting that T242N is associated with a fitness defect (32, 35).While CD8⁺ T-cell cross-clade recognition has been tested extensively (6, 11, 19, 36, 48), few studies have addressed the possibility of clade-specific escape from CD8⁺ T-cell responses. This may be especially relevant where clade consensus sequences differ in immunologically relevant epitopes. Here we demonstrate in a large Kenyan cohort substantial differences in HLA-B*57/B*5801-mediated selection among HIV clades.Participants were enrolled from a Nairobi, Kenya-based cohort, and the relevant ethical review boards approved the study. HLA typing was performed as described previously (34). CD4 counts were measured longitudinally at biannual visits. Multiple and other clade infections were excluded. The HIV-1 p24 gene was amplified from proviral HIV DNA or RNA using a nested PCR approach and sequenced, and viral subtyping was carried out as described previously (42). Previously described HLA-B*57 epitopes IW9 (ISPRTLNAW), KF11 (KAFSPEVIPMF), and TW10 (TSTLQEQIGW) and selected variants were tested in immunological assays and described where relevant. Gamma interferon enzyme-linked immunospot (ELISPOT) assays were performed as described previously (37) using blood samples from HLA-B*57⁺ and -B*5801⁺ subjects. All peptides were tested at concentrations of 10 μg, 1 μg, 0.1 μg, and 0.01 μg/ml. Responses were considered positive if they were more than two times higher than that of the negative control and were measured at ≥100 spot-forming units ml⁻¹. Fisher''s exact test and chi-square analyses were used to determine differences among groups in categorical analyses. Mann-Whitney U tests were used to compare response magnitudes and disease progression between groups.We confirmed the protective effects of HLA-B*57 in clade A1 infection (mean of 9.9 years versus 7.8 years until CD4 counts were <200, P = 0.041) (Fig. ​(Fig.1).1). Slow progressors were overrepresented in HLA-B*57⁺ clade A1⁺ subjects (52.2%) compared to both HLA-B*5801⁺ clade A1⁺ (13.3%, P = 0.02) and HLA-B*57/5801-negative clade A1⁺ (27.8%, P = 0.028) subjects (Fig. ​(Fig.1b).1b). In contrast to what has been shown for other clades (2, 27), protection was not observed for clade A1-infected HLA-B*5801⁺ subjects (mean of 6.5 years versus 7.8 years until CD4 counts were <200, P > 0.3) (Fig. ​(Fig.11).Open in a separate windowFIG. 1.HLA-B*57, but not HLA-B*5801, is associated with a lower rate of disease progression in clade A1-infected subjects than that of the overall cohort. (a) Number of years from cohort entry until sequential CD4 counts fell below 200/μl. (b) Slow progressors (>10 years with CD4 counts of >200) were also more common in HLA-B*57⁺ clade A1-infected subjects than in those expressing HLA-B*5801 or neither.Stratification of TW10 (Gag_240-249) proviral sequences on the basis of HLA allele and clade revealed several differences in selection between clades A1 and D (Fig. ​(Fig.2a).2a). We observed the expected T242N substitution in 100% of HLA-B*57⁺ clade D-infected subjects (7/7), compared to only 14.7% variability at Gag residue 242 in HLA-B*57/5801-negative subjects (13/88, P = 3.26 × 10⁻⁹) (Fig. ​(Fig.2b).2b). In contrast, T242N was found infrequently in clade A1-infected HLA-B*57⁺ subjects (15%, 5/33, P = 0.0004). Instead, variants containing the mutations P243T and I247L were more frequently observed (both observed in 11/33 subjects). Overall, variation at residues 243 and 247 was more common in HLA-B*57⁺ subjects (51% and 15%, respectively; P = 2.92 × 10⁻⁶) than in HLA-B*57/5801-negative clade A1⁺ subjects (33% and 9%, respectively; P = 0.0008). Selection at both residues 243 and 247 was observed only in 2/33 HLA-B*57⁺ subjects, suggesting that these mutations are independent. Selection at residue 248, observed in clade B infection (32), was not evident in either clade A1 or D. While I247X selection has been described in other clades at low frequencies and in elite controllers (21, 40), HLA-B*57-mediated selection at Gag residue 243 has not yet been described. In summary, the T242N mutation, which is typical of other clades, does not appear to be the primary escape mutant in clade A1.Open in a separate windowFIG. 2.HLA-B*57-mediated selection in TW10 differs between clade A1 and clade D. (a) TW10 sequences were stratified by HLA-B*57, HLA-B*5801, or other alleles (HLA-B*57/5801⁻) and compared between clades A1 and D, based on the clade B consensus TW10 sequence. Each subject is represented by one sequence, and the numbers of subjects with a given sequence are shown in parentheses. A summary of variation from the TW10 consensus at Gag residues 242, 243, 247, and 248 is shown for HLA-B*57⁺ (b) and -B*5801⁺ (c) subjects. (b and c) Clade D is shown at the top, and clade A1 is shown at the bottom. Variation is shown in dark gray, and consensus is shown in light gray. (d) The proportions of clade A1-infected subjects with selection at Gag residue 242 only, and those with selection at residues 242 and 243 in combination, are shown. *, P < 0.05; **, P < 0.005; ***, P < 0.0005.Previous studies have suggested that HLA-B*5801 places selection pressure on TW10, similar to that of HLA-B*57 (35). Similar to clades B and C, selection of T242N was evident in HLA-B*5801⁺ clade D-infected subjects (TW10 variation in 8/11 HLA-B*5801⁺ subjects versus 13/88 HLA-B*57/5801-negative subjects; P = 0.0069) (Fig. ​(Fig.2c).2c). Limited T242N selection was observed in clade A1-infected HLA-B*5801⁺ subjects, and in contrast to HLA-B*57, there were no HLA-B*5801-associated substitutions at residues 243 and 247 in clade A1 (P values of 0.75 and 0.29, respectively) (Fig. ​(Fig.2c).2c). In summary, these data suggest that in addition to HLA-B*5801 not being associated with protection in clade A1 (Fig. ​(Fig.1),1), HLA-B*5801 does not select the HLA-B*57-associated clade A1 TW10 escape mutations.Inclusion of all clade A1 sequences with T242X substitutions (regardless of the HLA allele) reveals that in every case (10/10), there is an accompanying residue 243 mutation. Polymorphisms at these sites correlate very strongly (P = 3.71 × 10⁻⁸) (Fig. ​(Fig.2d).2d). Together, these data suggest that residue 242 polymorphism in clade A1 is incompatible with proline at residue 243, which is the clade A1 consensus.We next assessed the immunological implications of novel clade A1 variants in HLA-B*57⁺ (n = 12) and -B*5801⁺ (n = 6) subjects infected primarily by clade A1. Clade A1-infected subjects commonly made anamnestic, low-avidity responses to TW10. The majority of HLA-B*57⁺ subjects who recognized clade A1 TW10 did not respond to P243T or I247L in ELISPOT assays (Fig. ​(Fig.3),3), supporting the hypothesis that these represent escape mutations. Those who did recognize P243T and I247L had lower magnitude responses than those who recognized clade A1 TW10 at the 10-μg/ml peptide concentration (P of 0.0005 for both) (Fig. ​(Fig.3a).3a). Similarly, these variants were not well recognized by CD8⁺ T cells from HLA-B*5801⁺ subjects (Fig. ​(Fig.3b).3b). For two HLA-B*57⁺ subjects, P243T and I247L responses had lower avidity than clade A1 TW10 responses (Fig. ​(Fig.3c).3c). These data support the hypothesis that P243T and I247L likely represent escape mutations.Open in a separate windowFIG. 3.Peptides with novel TW10 clade A1-selected mutations are poorly recognized in ex vivo gamma interferon ELISPOT avidity assays, suggestive of escape mutations. ELISPOT responses to TW10 and variants at 10 μg/ml peptide by HLA-B*57⁺ (a) and HLA-B*5801⁺ (b) subjects are shown. (c) The functional avidity of TW10 and variants for two HLA-B*57⁺ subjects is shown, suggesting that P243T and I247L are less recognized than TW10, particularly at lower peptide concentrations. Sequence names are described in the text. (d) Elispot responses to A1 TW10 and T242N correlated at 10 μg/ml. SFU/m, spot-forming units/million PBMCs.Recognition of the clade B/D consensus (TSTLQEQIGW) was diminished compared to that of clade A1 TW10. However, despite the presumed absence of this variant in these subjects'' autologous sequences, the clade B/D escape variant (TSNLQEQIGW [T242N]) was recognized at a magnitude similar to that of the consensus clade A1 TW10 (r = 0.71, P = 0.0099) (Fig. ​(Fig.3d).3d). No responses to T242N/G248A were observed (not shown), as described previously (32). These data suggest that clade A1 and B TW10, and their escape variants, are immunologically distinct from one another.We next assessed whether clade-specific selection was evident in other immunodominant HLA-B*57 p24 epitopes that are commonly targeted in chronic clade B infection (2). In clade D IW9 (ISPRTLNAW [Gag_147-155]), variants containing the escape variant I147L (14) were more common in HLA-B*57⁺ subjects than in HLA-B*57/5801-negative subjects (86% and 30%, respectively; P = 0.0055) (Table ​(Table1).1). However, this variant was not selected in clade A1 (variation in 30% versus 21% subjects; P value was not significant), where leucine is the consensus. Interestingly, ELISPOT data indicated substantial cross-reactivity between 147L and 147I in clade A1-infected subjects (10 μg/ml, r = 0.987, P < 0.0001) (data not shown), suggesting that infection with an escape variant from one clade (clade D) does not necessarily preclude recognition of this epitope in another one (clade A1). Other amino acids (F, M, and P) were common in HLA-B*57⁺ subjects at residue 147 (>30% versus 3% in HLA-B*57/5801-negative subjects, P = 3.85 × 10⁻⁶). Although the immunological consequences are unknown, these HLA-associated substitutions could represent novel escape variants.

TABLE 1.

p24 sequences in HLA-B*57⁺ subjects infected by clades A1 and D

Hepatitis C Virus (HCV) Sequence Variation Induces an HCV-Specific T-Cell Phenotype Analogous to Spontaneous Resolution

Hepatitis C virus (HCV)-specific CD8⁺ T cells in persistent HCV infection are low in frequency and paradoxically show a phenotype associated with controlled infections, expressing the memory marker CD127. We addressed to what extent this phenotype is dependent on the presence of cognate antigen. We analyzed virus-specific responses in acute and chronic HCV infections and sequenced autologous virus. We show that CD127 expression is associated with decreased antigenic stimulation after either viral clearance or viral variation. Our data indicate that most CD8 T-cell responses in chronic HCV infection do not target the circulating virus and that the appearance of HCV-specific CD127⁺ T cells is driven by viral variation.Hepatitis C virus (HCV) persists in the majority of acutely infected individuals, potentially leading to chronic hepatitis, liver cirrhosis, and hepatocellular carcinoma. The cellular immune response has been shown to play a significant role in viral control and protection from liver disease. Phenotypic and functional studies of virus-specific T cells have attempted to define the determinants of a successful versus an unsuccessful T-cell response in viral infections (10). So far these studies have failed to identify consistent distinguishing features between a T-cell response that results in self-limiting versus chronic HCV infection; similarly, the impact of viral persistence on HCV-specific memory T-cell formation is poorly understood.Interleukin-7 (IL-7) receptor alpha chain (CD127) is a key molecule associated with the maintenance of memory T-cell populations. Expression of CD127 on CD8 T cells is typically only observed when the respective antigen is controlled and in the presence of significant CD4⁺ T-cell help (9). Accordingly, cells specific for persistent viruses (e.g., HIV, cytomegalovirus [CMV], and Epstein-Barr virus [EBV]) have been shown to express low levels of CD127 (6, 12, 14) and to be dependent on antigen restimulation for their maintenance. In contrast, T cells specific for acute resolving virus infections, such as influenza virus, respiratory syncytial virus (RSV), hepatitis B virus (HBV), and vaccinia virus typically acquire expression of CD127 rapidly with the control of viremia (5, 12, 14). Results for HCV have been inconclusive. The expected increase in CD127 levels in acute resolving but not acute persisting infection has been found, while a substantial proportion of cells with high CD127 expression have been observed in long-established chronic infection (2). We tried to reconcile these observations by studying both subjects with acute and chronic HCV infection and identified the presence of antigen as the determinant of CD127 expression.Using HLA-peptide multimers we analyzed CD8⁺ HCV-specific T-cell responses and CD127 expression levels in acute and chronic HCV infection. We assessed a cohort of 18 chronically infected subjects as well as 9 individuals with previously resolved infection. In addition, we longitudinally studied 9 acutely infected subjects (5 individuals who resolved infection spontaneously and 4 individuals who remain chronically infected) (Tables ​(Tables11 and ​and2).2). Informed consent in writing was obtained from each patient, and the study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki, as reflected in a priori approval from the local institutional review boards. HLA-multimeric complexes were obtained commercially from Proimmune (Oxford, United Kingdom) and Beckman Coulter (CA). The staining and analysis procedure was as described previously (10). Peripheral blood mononuclear cells (PBMCs) were stained with the following antibodies: CD3 from Caltag; CD8, CD27, CCR7, CD127, and CD38 from BD Pharmingen; and PD-1 (kindly provided by Gordon Freeman). Primer sets were designed for different genotypes based on alignments of all available sequences from the public HCV database (http://hcvpub.ibcp.fr). Sequence analysis was performed as previously described (8).

TABLE 1.

Patient information and autologous sequence analysis for patients with chronic and resolved HCV infection

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

Transfer of a Phage T4 Gene into Enterobacteriaceae,Determined at the Single-Cell Level

The transfer range of phage genes was investigated at the single-cell level by using an in situ DNA amplification technique. After absorption of phages, a phage T4 gene was maintained in the genomes of non-plaque-forming bacteria at frequencies of 10⁻² gene copies per cell. The gene transfer decreased the mutation frequencies in nonhost recipients.Recently, whole-genome analyses have revealed that many bacterial genomes contain foreign genes, especially phage genes (9). The phage genes in bacterial genomes include genes for virulence or fitness factors such as extracellular toxins, superantigens, lipopolysaccharide-modifying enzymes, and proteins conferring serum resistance, etc. (1). These findings suggest that the horizontal transfer of phage genes has contributed significantly to the acquisition of new genetic traits and to the genetic diversity of bacteria (1, 9, 10). To truly appreciate the mechanisms behind phage-associated evolution, it is important to understand the frequency and range of transfer of phage genes.Most phage genomes consist of many genes derived from different origins (5, 8). Some genes are similar to those of other phages with phylogenetically different hosts or are found in the genomes of bacteria that are not the phage hosts. The mosaic nature of phage genomes has been known for some time, and a body of molecular genetic studies of phages to explain the mechanisms that drive this feature have been attempted previously (1, 5). More importantly, the horizontal transfer of phage genes has emerged as a major factor in the evolution of the phage genome. Since recombination between phage and phage/prophage can occur when these elements coexist in the same cell, coinfection with multiple phage species may result in the production of hybrid phage genomes (5). The pathways by which phages exchange genetic material vary dramatically in concert with host ranges. However, conventional plaque assays have shown that the host ranges of the phages studied are narrow. We hypothesized that phage genes can be transferred to more diverse species than previously thought.In order to accurately quantify DNA movement, gene targeting that does not require cultivation or gene expression is necessary (7). In situ DNA amplification methods allow the visualization of specific DNA sequences inside bacterial cells. In this study, we employed cycling primed in situ amplification-fluorescent in situ hybridization (CPRINS-FISH) to examine the possible range and frequency of the transfer of phage genes. CPRINS uses one primer and results in linear amplification of the target DNA inside cells, and multiply labeled fluorescent probe sets are applied for detection of the amplicons to improve the specificity and sensitivity of CPRINS (3). Previously, CPRINS-FISH did clarify the movement of DNA of a specific gene among Escherichia coli cells at the single-cell level (4).Enterobacterial phages P1 and T4 infect E. coli and have been well studied. P1 can exist as circular DNA within the bacterial cell as if it were a plasmid. Phage T4 is capable of undergoing only a lytic life cycle and not the lysogenic life cycle. Conventional methods using plaque assays have shown that the host of P1 and T4 is E. coli, but orthologous phage genes have been found in bacteria other than E. coli (6, 8). In the present study, strains of Enterobacteriaceae were allowed to grow on agar medium after the phage was adsorbed, and the maintenance of the transferred phage gene in the bacterial genomes was examined at the community level by quantitative real-time PCR and at the single-cell level by CPRINS-FISH.The following bacterial strains were used for maintenance experiments: Citrobacter freundii IFO 12681, Enterobacter aerogenes BM 2688, E. coli NBRC 12713, a Proteus mirabilis clinical isolate, Salmonella enterica serovar Enteritidis IID 640, and Yersinia enterocolitica IID 981. The bacterial strains were grown in Luria-Bertani (LB) medium (1% tryptone, 0.5% yeast extract, 0.5% NaCl; Nacalai Tesque Inc., Kyoto, Japan) at 37°C overnight.Stationary-phase cultures of 500 μl were incubated with 500 μl of SM buffer (50 mmol liter⁻¹ Tris-HCl [pH 7.5], 100 mmol liter⁻¹ NaCl, 8 mmol liter⁻¹ MgSO₄, 0.01% gelatin) containing the phage P1kc NBRC 20008 (2) or T4GT7 (11) at 37°C for 10 min at a multiplicity of infection of 1:1 (ratio of PFU of the phage to CFU of the recipient bacterium). The concentration of bacterial cells was adjusted to 10⁹ cells ml⁻¹. After 10 min of incubation, the diluted cell suspension (10⁵ cells) was filtered through a polycarbonate filter (Advantec, Tokyo, Japan) with a pore size of 0.2 μm and a diameter of 25 mm. Cells trapped on the filter were cultured on LB agar medium at 37°C for 24 h. The filter was transferred into a microtube, and cells on the filter were suspended in 1 ml of sterile deionized water. The numbers of cells in the suspension and cells remaining on the filter were determined by using an epifluorescence microscope (see below) after staining of the samples with 1 μg ml⁻¹ of 4′,6-diamidino-2-phenylindole (DAPI; Sigma-Aldrich Japan, Tokyo). The level of recovery of cells from the filter into sterile deionized water was about 99%. The cultured cells were subjected to real-time PCR and CPRINS-FISH.For real-time PCR, bacterial DNA was extracted using a QIAamp DNA isolation kit (Qiagen, Tokyo, Japan). The cell suspension was mixed with 10 mg ml⁻¹ of lysozyme solution and incubated at 37°C for 1 h. DNA extraction was then performed according to the kit manufacturer''s instructions. Table ​Table11 lists the oligonucleotide primers for PCR and CPRINS and the polynucleotide probes used in the present study. Tail fiber genes from phages P1kc and T4GT7 were quantified by real-time PCR with a LightCycler system (Roche Diagnostics, Tokyo, Japan). LightCycler FastStart DNA master SYBR green I (Roche Diagnostics) was used with 5 mmol liter⁻¹ Mg²⁺ and 0.5 μmol liter⁻¹ (each) primers targeting the tail fiber genes of P1kc (P1-tail931f and P1-tail1148r) and T4GT7 (T4-tail2770f and T4-tail2983r). After a hot start for 10 min at 95°C, 40 cycles of PCR were run with denaturation at 94°C for 15 s, annealing at 60°C for 10 s, extension at 72°C for 10 s, and fluorescence detection at 83°C for 5 s. The known amounts of PCR products from the phage DNA (10¹ to 10⁷ copies per reaction) were used for the standard curves to quantify the target DNA. To confirm the specificity of the reaction after real-time PCR, the PCR mixture was collected in a glass capillary and subjected to agarose gel electrophoresis in addition to a melting-curve analysis with the LightCycler system. The maintenance frequencies determined by real-time PCR were recorded as the copy number of the phage tail fiber gene per bacterial genome detected by staining with PicoGreen (Invitrogen, Tokyo, Japan) after cultivation of cells on LB agar medium for 24 h as described above. The frequencies were determined in triplicate for each sample. The increase in the phage gene copy number was determined by comparing the copy numbers in cells on the filter before and after cultivation. The phage gene copy number in cells on the filter was determined by the following formula: (total number of cells determined by DAPI staining) × (phage tail fiber gene copy number determined by real-time PCR)/(bacterial genome copy number determined by PicoGreen staining).

TABLE 1.

Probes and primers designed in this study

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

Different Minimal Signal Peptide Lengths Recognized by the Archaeal Prepilin-Like Peptidases FlaK and PibD

In Archaea, the preflagellin peptidase (a type IV prepilin-like peptidase designated FlaK in Methanococcus voltae and Methanococcus maripaludis) is the enzyme that cleaves the N-terminal signal peptide from preflagellins. In methanogens and several other archaeal species, the typical flagellin signal peptide length is 11 to 12 amino acids, while in other archaea preflagellins possess extremely short signal peptides. A systematic approach to address the signal peptide length requirement for preflagellin processing is presented in this study. M. voltae preflagellin FlaB2 proteins with signal peptides 3 to 12 amino acids in length were generated and used as a substrate in an in vitro assay utilizing M. voltae membranes as an enzyme source. Processing by FlaK was observed in FlaB2 proteins containing signal peptides shortened to 5 amino acids; signal peptides 4 or 3 amino acids in length were unprocessed. In the case of Sulfolobus solfataricus, where the preflagellin peptidase PibD has broader substrate specificity, some predicted substrates have predicted signal peptides as short as 3 amino acids. Interestingly, the shorter signal peptides of the various mutant FlaB2 proteins not processed by FlaK were processed by PibD, suggesting that some archaeal preflagellin peptidases are likely adapted toward cleaving shorter signal peptides. The functional complementation of signal peptidase activity by FlaK and PibD in an M. maripaludis ΔflaK mutant indicated that processing of preflagellins was detected by complementation with either FlaK or PibD, yet only FlaK-complemented cells were flagellated. This suggested that a block in an assembly step subsequent to signal peptide removal occurred in the PibD complementation.The bacterial type IV prepilin peptidase (TFPP) is a well-characterized enzyme belonging to a family of novel aspartic acid proteases (20). It is responsible for the cleavage of N-terminal signal peptides from prepilins and pseudopilins, prior to their incorporation into the type IV pilus structure (22, 30, 31). The prepilin peptidase is also responsible for the processing of prepilin-like proteins needed for type II secretion (22). In Archaea, the existence of bacterial TFPP-like enzymes has also been reported, and they have been most extensively studied in relation to the assembly of the archaeal flagellum. In the euryarchaeotes Methanococcus maripaludis and Methanococcus voltae, the preflagellin peptidase FlaK was demonstrated to be responsible for cleaving the N-terminal signal peptide from the preflagellin prior to its incorporation into the growing flagellar filament, a step essential to flagellar assembly (6, 7, 26). In Sulfolobus solfataricus, an acidophilic crenarchaeote, the equivalent enzyme, PibD, was also shown to process preflagellins (4). Site-directed mutagenesis of FlaK and PibD demonstrated that both aspartic acid residues that aligned with aspartic acid residues essential for bacterial TFPP activity were also essential in the archaeal enzymes (6, 32), indicating that the two archaeal peptidases belong with the bacterial TFPPs in this novel family of aspartic acid proteases (20). More recently, an additional archaeal TFPP was found to be required for cleavage of the prepilin substrates (33) that are assembled into the unique pili of M. maripaludis (37).The substrate specificity of the archaeal preflagellin peptidase remains an open question. Like prepilin peptidases, FlaK in M. voltae has stringent requirements for the amino acids surrounding the cleavage site of the substrate, especially the −1 glycine, −2 and −3 lysines, and the +3 glycine (numbers given relative to the cleavage site) (35); the last position was conserved in all archaeal flagellins (25). Upon N-terminal sequence alignment of all available archaeal flagellin amino acid sequences at the predicted cleavage site, it was found that most archaeal preflagellin signal peptides are quite conserved in length, with the typical flagellin signal peptide being 11 to 12 amino acids in length (Table ​(Table1).1). It is speculated that while a certain amount of flexibility might exist, some optimum and minimum length probably exists that is crucial for the juxtaposition of the signal peptide and signal peptidase with respect to each other and the membrane (18). A recent study examining possible type IV pilin-like substrates in archaea using the FlaFind program indicated that such substrates may be more widespread than initially thought (33). Since in Methanococcus the pilins are processed by a second TFPP (EppA) (33), it is very possible that the preflagellins might be the only substrates of FlaK in these archaea.

TABLE 1.

N-terminal amino acid alignment of selected archaeal flagellin sequences^a

Detection and Identification of tdh- and trh-Positive Vibrio parahaemolyticus Strains from Four Species of Cultured Bivalve Molluscs on the Spanish Mediterranean Coast

Presented here is the first report describing the detection of potentially diarrheal Vibrio parahaemolyticus strains isolated from cultured bivalves on the Mediterranean coast, providing data on the presence of both tdh- and trh-positive isolates. Potentially diarrheal V. parahaemolyticus strains were isolated from four species of bivalves collected from both bays of the Ebro delta, Spain.Gastroenteritis caused by Vibrio parahaemolyticus has been reported worldwide, though only sporadic cases have been reported in Europe (7, 14). The bacterium can be naturally present in seafood, but pathogenic isolates capable of inducing gastroenteritis in humans are rare in environmental samples (2 to 3%) (15) and are often not detected (10, 19, 20).The virulence of V. parahaemolyticus is based on the presence of a thermostable direct hemolysin (tdh) and/or the thermostable direct hemolysin-related gene (trh) (1, 5). Both are associated with gastrointestinal illnesses (2, 9).Spain is not only the second-largest producer in the world of live bivalve molluscs but also one of the largest consumers of bivalve molluscs, and Catalonia is the second-most important bivalve producer of the Spanish Autonomous Regions. Currently, the cultivation of bivalves in this area is concentrated in the delta region of the Ebro River. The risk of potentially pathogenic Vibrio spp. in products placed on the market is not assessed by existing legislative indices of food safety in the European Union, which emphasizes the need for a better knowledge of the prevalence of diarrheal vibrios in seafood products. The aim of this study was to investigate the distribution and pathogenic potential of V. parahaemolyticus in bivalve species exploited in the bays of the Ebro delta.Thirty animals of each species of Mytilus galloprovincialis, Crassostrea gigas, Ruditapes decussatus, and Ruditapes philippinarum were collected. They were sampled from six sites of the culture area, three in each bay of the Ebro River delta, at the beginning (40°37′112"N, 0°37′092"E [Alfacs]; 40°46′723"N, 0°43′943"E [Fangar]), middle (40°37′125"N, 0°38′570"E [Alfacs]; 40°46′666"N, 0°45′855"E [Fangar]), and end (40°37′309"N, 0°39′934"E [Alfacs]; 40°46′338"N, 0°44′941"E [Fangar]) of the culture polygon. Clams were sampled from only one site per bay as follows: in the Alfacs Bay from a natural bed of R. decussatus (40°37′44"N, 0°38′0"E) and in the Fangar Bay from an aquaculture bed of R. philippinarum (40°47′3"N, 0°43′8"E). In total, 367 samples were analyzed in 2006 (180 oysters, 127 mussels, 30 carpet shell clams, and 30 Manila clams) and 417 samples were analyzed in 2008 (178 oysters, 179 mussels, 30 carpet shell clams, and 30 Manila clams).All animals were individually processed and homogenized, and 1 ml of the homogenate was inoculated into 9 ml of alkaline peptone water (Scharlau, Spain). Following a 6-h incubation at 37°C, one loopful of the contents of each tube of alkaline peptone water was streaked onto CHROMagar vibrio plates (CHROMagar, France) and incubated for 18 h at 37°C. Mauve-purple colonies were purified, and each purified isolate was cryopreserved at −80°C (135 isolates in 2006 and 96 in 2008). From the initial homogenate portion, 100 μl was inoculated onto marine agar (Scharlau, Spain) and onto thiosulfate citrate-bile salts-sucrose agar (Scharlau, Spain) for total heterotrophic marine bacteria counts and total vibrio counts, respectively (Table ​(Table11).

TABLE 1.

Vibrio parahaemolyticus isolates, serotypes, and origins and total number of vibrios/heterotrophic bacteria contained in the bivalve^a

Physicochemical Basis for the Inhibitory Effects of Organic and Inorganic Salts on the Growth of Pectobacterium carotovorum subsp. carotovorum and Pectobacterium atrosepticum

Twenty-one salts were tested for their effects on the growth of Pectobacterium carotovorum subsp. carotovorum and Pectobacterium atrosepticum. In liquid medium, 11 salts (0.2 M) exhibited strong inhibition of bacterial growth. The inhibitory action of salts relates to the water-ionizing capacity and the lipophilicity of their constituent ions.Different biochemical mechanisms have been put forth to explain the antimicrobial activity of organic and inorganic salts, including inhibition of several steps of the energy metabolism (benzoate, bicarbonate, propionate, sorbate, and sulfite salts) (2, 3, 11, 16, 17, 19, 25) and complexation to DNA and RNA (aluminum and sulfites) (12, 13, 15, 20, 27, 28). However, little is known about the physicochemical basis for the general antimicrobial action of salts. The objective of this work was to gain an understanding of the relationship between the inhibitory action of salts on bacterial growth and their physicochemical properties by using the bacteria Pectobacterium carotovorum subsp. carotovorum (formerly Erwinia carotovora subsp. carotovora) and Pectobacterium atrosepticum (formerly Erwinia carotovora subsp. atroseptica). These bacteria are responsible for soft rot, a disease of economic importance affecting numerous stored vegetable crops (14, 22).Pectobacterium carotovorum subsp. carotovorum (strain Ecc 1367) and P. atrosepticum (strain Eca 709), provided by the Laboratoire de Diagnostic en Phytoprotection (MAPAQ, Québec, Canada), were grown in 250-ml flasks containing 50 ml of 20% tryptic soy broth (Difco Laboratories, Becton Dickinson, Sparks, MD) amended with salts (200 mM) or unamended (control), by incubation at 24°C with agitation (150 rpm; Lab-Line Instruments Inc., Melrose Park, IL) for 24 h. The pHs of the media were not adjusted but varied with the type of salts, unless stated otherwise. Flasks were inoculated with 100 μl of each bacterial suspension (1 × 10⁷ CFU/ml). Bacterial growth was determined by turbidimetry at 600 nm with a UV/visible spectrophotometer (Ultrospec 2000; Pharmacia Biotech Ltd, Cambridge, United Kingdom), using appropriate blanks. Results were expressed as the percentage of growth inhibition compared with the growth of the control. A completely randomized experimental design with three replicates was used, the experimental unit being a flask. Analysis of variance was carried out with the GLM (general linear model) procedure of SAS (SAS Institute, Cary, NC) software. When they were significant (P < 0.05), treatment means were compared using Fisher''s protected least-significant-difference test.Among the 21 salts tested, sodium carbonate, sodium metabisulfite, trisodium phosphate, aluminum lactate, aluminum chloride, sodium bicarbonate, sodium propionate, ammonium acetate, aluminum dihydroxy acetate, potassium sorbate, and sodium benzoate exhibited strong inhibition (≥97%) of the growth of both P. carotovorum subsp. carotovorum and P. atrosepticum (Table ​(Table1).1). Calcium chloride, sodium formate, sodium acetate, ammonium hydrogen phosphate, and sodium hydrogen phosphate exhibited a moderately inhibitory effect; sodium lactate and tartrate had no effect. On the other hand, ammonium chloride, potassium chloride, and sodium chloride stimulated the growth of P. atrosepticum.

TABLE 1.

Effect of salts on the growth of P. atrosepticum and P. carotovorum subsp. carotovorum

Escherichia coli O157:H7 and Other E. coli Strains Share Physiological Properties Associated with Intestinal Colonization

Escherichia coli isolates (72 commensal and 10 O157:H7 isolates) were compared with regard to physiological and growth parameters related to their ability to survive and persist in the gastrointestinal tract and found to be similar. We propose that nonhuman hosts in E. coli O157:H7 strains function similarly to other E. coli strains in regard to attributes relevant to gastrointestinal colonization.Escherichia coli is well known for its ecological versatility (15). A life cycle which includes both gastrointestinal and environmental stages has been stressed by both Savageau (15) and Adamowicz et al. (1). The gastrointestinal stage would be subjected to acid and detergent stress. The environmental stage is implicit in E. coli having transport systems for fungal siderophores (4) as well as pyrroloquinoline quinone-dependent periplasmic glucose utilization (1) because their presence indicates evolution in a location containing fungal siderophores and pyrroloquinoline quinone (1).Since its recognition as a food-borne pathogen, there have been numerous outbreaks of food-borne infection due to E. coli O157:H7, in both ground beef and vegetable crops (6, 13). Cattle are widely considered to be the primary reservoir of E. coli O157:H7 (14), but E. coli O157:H7 does not appear to cause disease in cattle. To what extent is E. coli O157:H7 physiologically unique compared to the other naturally occurring E. coli strains? We feel that the uniqueness of E. coli O157:H7 should be evaluated against a backdrop of other wild-type E. coli strains, and in this regard, we chose the 72-strain ECOR reference collection originally described by Ochman and Selander (10). These strains were chosen from a collection of 2,600 E. coli isolates to provide diversity with regard to host species, geographical distribution, and electromorph profiles at 11 enzyme loci (10).In our study we compared the 72 strains of the ECOR collection against 10 strains of E. coli O157:H7 and six strains of E. coli which had been in laboratory use for many years (Table ​(Table1).1). The in vitro comparisons were made with regard to factors potentially relevant to the bacteria''s ability to colonize animal guts, i.e., acid tolerance, detergent tolerance, and the presence of the Entner-Doudoroff (ED) pathway (Table ​(Table2).2). Our longstanding interest in the ED pathway (11) derives in part from work by Paul Cohen''s group (16, 17) showing that the ED pathway is important for E. coli colonization of the mouse large intestine. Growth was assessed by replica plating 88 strains of E. coli under 40 conditions (Table ​(Table2).2). These included two LB controls (aerobic and anaerobic), 14 for detergent stress (sodium dodecyl sulfate [SDS], hexadecyltrimethylammonium bromide [CTAB], and benzalkonium chloride, both aerobic and anaerobic), 16 for acid stress (pH 6.5, 6.0, 5.0, 4.6, 4.3, 4.2, 4.1, and 4.0), four for the ability to grow in a defined minimal medium (M63 glucose salts with and without thiamine), and four for the presence or absence of a functional ED pathway (M63 with gluconate or glucuronate). All tests were done with duplicate plates in two or three separate trials. The data are available in Tables S1 to S14 in the supplemental material, and they are summarized in Table ​Table22.

TABLE 1.

E. coli strains used in this study

Single Nucleotide Polymorphism-Based Diagnostic System for Crop-Associated Sclerotinia Species

A molecular diagnostic system using single nucleotide polymorphisms (SNPs) was developed to identify four Sclerotinia species: S. sclerotiorum (Lib.) de Bary, S. minor Jagger, S. trifoliorum Erikss., and the undescribed species Sclerotinia species 1. DNAs of samples are hybridized with each of five 15-bp oligonucleotide probes containing an SNP site midsequence unique to each species. For additional verification, hybridizations were performed using diagnostic single nucleotide substitutions at a 17-bp sequence of the calmodulin locus. The accuracy of these procedures was compared to that of a restriction fragment length polymorphism (RFLP) method based on Southern hybridizations of EcoRI-digested genomic DNA probed with the ribosomal DNA-containing plasmid probe pMF2, previously shown to differentiate S. sclerotiorum, S. minor, and S. trifoliorum. The efficiency of the SNP-based assay as a diagnostic test was evaluated in a blind screening of 48 Sclerotinia isolates from agricultural and wild hosts. One isolate of Botrytis cinerea was used as a negative control. The SNP-based assay accurately identified 96% of Sclerotinia isolates and could be performed faster than RFLP profiling using pMF2. This method shows promise for accurate, high-throughput species identification.Sclerotinia is distinguished morphologically from other genera in the Sclerotiniaceae (Ascomycota, Pezizomycotina, Leotiomycetes) by the production of tuberoid sclerotia that do not incorporate host tissue, by the production of microconidia that function as spermatia but not as a disseminative asexual state, and by the development of a layer of textura globulosa composing the outer tissue of apothecia (8). Two hundred forty-six species of Sclerotinia have been reported, most distinguished morphotaxonomically (Index Fungorum [www.indexfungorum.org]). These include the four species of agricultural importance now recognized plus many that are imperfectly known, seldom collected, or apparently endemic to relatively small geographic areas (2, 5, 6, 7, 8, 9, 17).The main species of phythopathological interest in the genus Sclerotinia are S. sclerotiorum (Lib.) de Bary, S. minor Jagger, S. trifoliorum Erikss., and the undescribed species Sclerotinia species 1. Sclerotinia species 1 is an important cause of disease in vegetables in Alaska (16) and has been found in association with wild Taraxacum sp., Caltha palustris, and Aconitum septentrionalis in Norway (7). It is morphologically indistinguishable from S. sclerotiorum, but it was shown to be a distinct species based on distinctive polymorphisms in sequences from internal transcribed spacer 2 (ITS2) of the nuclear ribosomal repeat (7). The other three species have been delimited using morphological, cytological, biochemical, and molecular characters (3, 8, 9, 10, 12, 15). Interestingly, given that the ITS is sufficiently polymorphic in many fungal genera to resolve species, in Sclerotinia, only species 1 and S. trifoliorum are distinguished by characteristic ITS sequence polymorphisms; S. sclerotiorum and S. minor cannot be distinguished based on ITS sequence (2, 7).Sclerotinia sclerotiorum is a necrotrophic pathogen with a broad host range (1). S. minor has a more restricted host range but causes disease in a variety of important crops such as lettuce, peanut, and sunflower crops (11). S. trifoliorum has a much narrower host range, limited to the Fabaceae (3, 8, 9). Sclerotial and ascospore characteristics also serve to differentiate among the three species. Sclerotinia minor has small sclerotia that develop throughout the colony in vitro and aggregate to form crusts on the host, while the sclerotia of S. sclerotiorum and S. trifoliorum are large and form at the colony periphery in vitro, remaining separate on the host (8, 9). The failure of an isolate to produce sclerotia or apothecia in vitro is not unusual, especially after serial cultivation (8). The presence of dimorphic, tetranucleate ascospores characterizes S. trifoliorum, while S. sclerotiorum and S. minor both have uniformly sized ascospores that are binucleate and tetranucleate, respectively (9, 14).With the apparent exception of Sclerotinia species 1, morphological characteristics are sufficient to delimit Sclerotinia species given that workers have all manifestations of the life cycle in hand. In cultures freshly isolated from infected plants, investigators usually have mycelia and sclerotia but not apothecia. Restriction fragment length polymorphisms (RFLPs) in ribosomal DNA (rDNA) are diagnostic for Sclerotinia species (3, 10), but the assay requires cloned probes (usually accessed from other laboratories) hybridized to Southern blots from vertical gels, an impractical procedure for large samples. We have analyzed sequence data from previous phylogenetic studies (2) and have identified diagnostic variation for the rapid identification of the four Sclerotinia species. The single nucleotide polymorphism (SNP) assay that we report here is amenable to a high throughput of samples and requires only PCR amplification with a standard set of primers and oligonucleotide hybridizations to Southern blots in a dot format.The SNP assay was performed using two independent sets of species-specific oligonucleotide probes, all with SNP sites shown to differentiate the four Sclerotinia species (Fig. ​(Fig.1).1). A panel of 49 anonymously coded isolates (Table ​(Table1)1) was screened using these species-specific SNP probes, as outlined in Fig. ​Fig.1.1. The assay was validated by comparison to Southern hybridizations of EcoRI-digested genomic DNA hybridized with pMF2, a plasmid probe containing the portion of the rDNA repeat with the 18S, 5.8S, and 26S rRNA cistrons of Neurospora crassa (4, 10).Open in a separate windowFIG. 1.Protocol for the SNP-based identification of Sclerotinia species, with diagnostic SNP sites underlined and in boldface type for each hybridization probe.

TABLE 1.

Isolates and hybridization results for all SNP-based oligonucleotide probes^f

Distribution of Sulfonamide Resistance Genes in Escherichia coli and Salmonella Isolates from Swine and Chickens at Abattoirs in Ontario and Québec,Canada

Sulfonamide-resistant Escherichia coli and Salmonella isolates from pigs and chickens in Ontario and Québec were screened for sul1, sul2, and sul3 by PCR. Each sul gene was distributed differently across populations, with a significant difference between distribution in commensal E. coli and Salmonella isolates and sul3 restricted mainly to porcine E. coli isolates.Resistance to sulfonamides is frequent in bacteria from farm animals (7, 8, 9, 10) and is usually caused by the acquisition of the genes sul1, sul2, and sul3 (20, 22). The objectives of this study were (i) to assess the distribution of these genes in Escherichia coli and Salmonella enterica isolates in swine and chickens from two major provinces in Canada, (ii) to assess whether differences occur in the distribution of these genes among bacterial species found within two different animal host species, and (iii) to assess whether significant differences in the distribution of these genes are present between the commensal E. coli strains used as indicators for surveillance of antimicrobial resistance and the zoonotic Salmonella pathogens found in the same ecological niche. In contrast to previous studies, a multivariable logistic regression model was used to analyze the data, control for confounding factors, and assess the interaction effect between animal and bacterial species in terms of the probability of an isolate carrying a specific sul gene. The distribution of sulfonamide resistance genes among sulfonamide-resistant E. coli (393 isolates from chickens and 311 from swine) and Salmonella (13 isolates from chickens and 221 from swine) isolates was assessed. These isolates were collected by the Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) between 2003 and 2005 from ceca of apparently healthy animals at abattoirs in Ontario (n = 435) and Québec (n = 503). The methods used by CIPARS are presented in detail elsewhere (8-10). The isolates were screened with a previously published multiplex PCR for sul1, sul2, and sul3 (16). The sul1 and sul2 genes were found in E. coli and Salmonella isolates from both animal species. The sul3 gene was detected in both E. coli and Salmonella isolates from swine but only in E. coli isolates from chickens (Table ​(Table1).1). Three percent of the isolates had no detectable sul gene, 12.5% possessed two genes, and two isolates carried three genes (Table ​(Table1).1). Similar (2, 3, 14, 19) or higher (4, 11, 17) values for multiple genes have been reported by others. The overall higher prevalence of sul1 in Salmonella isolates and of sul2 and sul3 in E. coli isolates was in agreement with the results of previous studies (2-4, 11, 12, 14, 21).

TABLE 1.

Frequency of the three resistance genes sul1, sul2, and sul3 in sulfonamide-resistant E. coli and Salmonella isolates from chickens and swine in Ontario and Québec between 2003 and 2005

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

Bar-Coded Pyrosequencing of 16S rRNA Gene Amplicons Reveals Changes in Ileal Porcine Bacterial Communities Due to High Dietary Zinc Intake

Feeding high levels of zinc oxide to piglets significantly increased the relative abundance of ileal Weissella spp., Leuconostoc spp., and Streptococcus spp., reduced the occurrence of Sarcina spp. and Neisseria spp., and led to numerical increases of all Gram-negative facultative anaerobic genera. High dietary zinc oxide intake has a major impact on the porcine ileal bacterial composition.Zinc oxide (ZnO) is used as a feed additive for diarrhea prophylaxis in piglets (23). However, the mode of action of ZnO is not fully understood. Besides its effects on the host (10, 30, 31), high dietary zinc levels may affect the diversity of intestinal microbial communities (2, 11, 20). The prevention of postweaning diarrhea in piglets due to high dietary ZnO intake may not be directly related to a reduction of pathogenic E. coli (8) but, rather, to the diversity of the coliform community (15). Studies on the impact of high ZnO levels on the porcine ileal bacterial community are scarce but nevertheless important, as bacterial diarrhea is initiated in the small intestine (9, 17). The small intestine is a very complex habitat with many different factors shaping the bacterial community. Studies on the ecophysiology (22) and maturation of the porcine ileal microbiota (13, 27) indicate a drastic impact directly after weaning and a gradual decline of modifications during the following 2 weeks. Thus, the time point for analysis chosen in this study (14 days postweaning) does reflect a more stable period of the ileal porcine microbiota. In this study, we used bar-coded pyrosequencing of 16S rRNA genes to gain further insight into the mode of action of pharmacological levels of ZnO in the gastrointestinal tract of young pigs.Total DNA was extracted from the ileal digesta of 40- to 42-day-old piglets using a commercial kit (Qiagen stool kit; Qiagen, Hilden, Germany) and PCR amplified with unique bar-coded primer sets targeting the V1-to-V3 and the V6-to-V8 hypervariable regions (see the supplemental material for detailed methods). The rationale behind this approach was derived from the fact that no single “universal” primer pair can completely cover a complex bacterial habitat (4, 24, 32, 33). Furthermore, these studies also show that in silico information on the coverage of selected primer sets diverges from empirical results, and hence, two hypervariable regions were chosen in this study to maximize the detection of phylogenetically diverse bacterial groups.Equimolar dilutions of all samples were combined into one master sample. Pyrosequencing was performed by Agowa (Berlin, Germany) on a Roche genome sequencer FLX system using a Titanium series PicoTiterPlate. The resulting data files were uploaded to the MG-RAST server (http://metagenomics.nmpdr.org/) (19) and processed with its SEED software tool using the RDP database (5) as the reference database. After automated sequence analysis, all sequences with less than five identical reads per sample were deleted in order to increase the confidence of sequence reads and reduce bias from possible sequencing errors (12, 16). Thus, 0.43% of all sequences were not considered (1,882 of 433,302 sequences). These sequences were assigned to a total of 238 genera, of which most only occurred in a few samples (see the supplemental material). Furthermore, all unclassified sequences were removed (8.7%; 41,467 of 474,769 sequences). Due to the use of the RDP reference database, the SEED software incorrectly assigned the majority of unclassified sequences as unclassified Deferribacterales (83%; 34,393 sequences), which were actually identified as 16S soybean or wheat chloroplasts by BLAST or as cyanobacterial chloroplasts by the RDP II seqmatch tool.The pyrosequencing results for the two primer combinations were merged by taking only sequences from the primer combination that yielded the higher number of reads for a specific sequence assignment in a sample. The remaining reads were used to calculate the relative contribution of assigned sequences to total sequence reads in a sample.The Firmicutes phylum dominated the small intestinal bacterial communities in both the control group and the group with high dietary ZnO intake, with 98.3% and 97.0% of total sequence reads, respectively. No significant influence of high dietary ZnO intake was found for the main phyla Proteobacteria (0.92% versus 1.84%), Actinobacteria (0.61% versus 0.75%), Bacteroidetes (0.15% versus 0.17%), and Fusobacteria (0.09% versus 0.12%).On the order level, a total of 20 bacterial orders were detected (data not shown). Lactobacillales dominated bacterial communities in the control and high-dietary-ZnO-intake groups, with 83.37% and 93.24% of total reads. Lactic acid bacteria are well known to dominate the bacterial community in the ileum of piglets (11, 22). No significant difference between the control group and the group with high dietary ZnO intake was observed on the order level, although high dietary ZnO intake led to a strong numerical decrease for Clostridiales (14.4 ± 24.0% [mean ± standard deviation] versus 2.8 ± 1.7%), as well as to numerical increases for Pseudomonadales (0.3 ± 0.3% versus 0.6 ± 0.6%) and Enterobacteriales (0.2 ± 0.2% versus 0.5 ± 0.6%).On the genus level, a total of 103 genera were detected. Table ​Table11 summarizes the main 31 genera which exceeded 0.05% of total reads (see the supplemental material for a complete list). Lactobacilli clearly dominated the bacterial communities in both trial groups, but they also were numerically lower due to high dietary ZnO intake.

TABLE 1.

Bacterial genera in the ileum of piglets fed diets supplemented with 200 or 3,000 ppm ZnO

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