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

Cross-Subtype Neutralization Sensitivity despite Monoclonal Antibody Resistance among Early Subtype A,C, and D Envelope Variants of Human Immunodeficiency Virus Type 1

The human immunodeficiency virus type 1 (HIV-1) variants that are transmitted to newly infected individuals are the primary targets of interventions, such as vaccines and microbicides, aimed at preventing new infections. Newly acquired subtype A, B, and C variants have been the focus of neutralization studies, although many of these viruses, particularly of subtypes A and B, represent viruses circulating more than a decade ago. In order to better represent the global diversity of transmitted HIV-1 variants, an additional 31 sexually transmitted Kenyan HIV-1 env genes, representing several recent infections with subtype A, as well as subtypes A/D, C, and D, were cloned, and their neutralization profiles were characterized. Most env variants were resistant to neutralization by the monoclonal antibodies (MAbs) b12, 4E10, 2F5, and 2G12, suggesting that targeting the epitopes of these MAbs may not be effective against variants that are spreading in areas of endemicity. However, significant cross-subtype neutralization by plasma was observed, indicating that there may be other epitopes, not yet defined by the limited available MAbs, which could be recognized more broadly.Most effective viral vaccines are thought to provide protection primarily by stimulating neutralizing antibodies (NAbs) to clear cell-free virus (25, 27). Because protection by NAbs requires recognition of common viral epitopes, the extreme genetic diversity of human immunodeficiency virus type 1 (HIV-1) presents a particular challenge to NAb-based vaccine approaches. Therefore, a critical starting point for studies of immune-mediated protection against HIV-1 is a collection of newly transmitted HIV-1 variants, particularly from areas of endemicity, such as sub-Saharan Africa, in order to determine whether vaccines are appropriately targeted to common epitopes from these relevant transmitted strains.During HIV-1 transmission, a bottleneck allows only one or a few variants to be transmitted to a newly infected individual (6, 9, 16, 29, 34, 37, 39), and the sensitivity of these early transmitted strains to antibody-mediated neutralization is therefore of particular interest. Newly transmitted HIV-1 variants have demonstrated significant heterogeneity in their neutralization phenotypes both within and between subtypes (2, 3, 6-8, 11, 13-15, 22, 30, 32, 36). Panels of sexually transmitted HIV-1 envelope variants (based on the envelope gene, env) have been characterized, including subtype B variants from North America, Trinidad, and Europe, subtype C variants from South Africa and Zambia, and subtype A variants from Kenya collected between 1994 and 1996 (2, 14, 15). Here, we characterize an additional 31 envelope variants from 14 subjects with sexually transmitted HIV-1 who were infected in Kenya, where subtypes A, C, and D circulate, between 1993 and 2005 (24, 31).The env genes were cloned from samples drawn 14 to 391 (median, 65) days postinfection from individuals enrolled in a prospective cohort of high-risk women in Mombasa, Kenya (19-21). Demographic characteristics of the subjects are summarized in Table ​Table1;1; the timing of first infection was determined by both HIV-1 serology and HIV RNA testing as described previously (12). All of the subjects were presumably infected by male-to-female transmission and displayed a range of plasma viral loads at the time of env gene cloning (Table ​(Table1).1). For most individuals, full-length env genes were cloned from uncultured peripheral blood mononuclear cell (PBMC) DNA, though for two individuals, clones were obtained from DNA following short-term coculture with donor PBMCs (Table ​(Table1).1). env genes were cloned by single-copy nested PCR with primers and PCR conditions as described previously (4, 17). We tested env genes for their ability to mediate infection by transfecting env plasmid DNA into 293T cells along with an env-deficient HIV-1 subtype A proviral plasmid, Q23Δenv, to make pseudoviral particles (17). More than 80 env clones were obtained from 16 subjects; less than one-half were functional on the basis of the infectivity of pseudoviral particles in a single-round infection of TZM-bl cells (AIDS Research and Reference Reagent Program, National Institutes of Health), as observed previously for env genes cloned from proviral sequences (17); a lower fraction of functional env genes have been reported from plasma (18). We focused on the proviral sequences here because they presumably best represent the sequence closest to that of the transmitted strains. The 31 functional env variants are described in Table ​Table11.

TABLE 1.

Demographic characteristics, diversities, gp120 variable-region lengths, numbers of PNGS, and accession numbers of cloned env variants

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

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

TABLE 1.

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

Fitness Effects of Replichore Imbalance in Salmonella enterica

A fitness cost due to imbalanced replichores has been proposed to provoke chromosome rearrangements in Salmonella enterica serovars. To determine the impact of replichore imbalance on fitness, the relative fitness of isogenic Salmonella strains containing transposon-held duplications of various sizes and at various chromosomal locations was determined. Although duplication of certain genes influenced fitness, a replichore imbalance of up to 16° did not affect fitness.The bacterial chromosome is a dynamic molecule that can undergo various types of rearrangements, including inversions, translocations, and duplications. These rearrangements can alter gene order and change replichore length. Replichores are defined as the halves of the chromosome between the origin of replication and the terminus region in the vicinity of the dif site (4, 6, 10, 15). In most bacteria, both replichores are approximately the same length, with each comprising 180° around the circular chromosome. In this state, the replichores and DNA replication are balanced. Imbalance is introduced when one replichore becomes longer than the other. For example, if the replichores comprise 200° and 160° around the circular chromosome, the replichores are 20° imbalanced. Various studies over the years have investigated the effect of imbalanced DNA replication on fitness in Escherichia coli (8, 9, 16, 17). In a recent analysis, E. coli strains that contained asymmetrical interreplichore inversions were found to have growth defects when the imbalance was at least 50° (8). While these studies have demonstrated that imbalanced replichores can affect fitness, the approaches used to introduce the imbalance either do not occur or are extremely rare in nature.Duplications play important evolutionary roles because the increase in gene copy number can facilitate adaptation to certain growth conditions (20), as well as being a source for the evolution of new genes (3). Typically duplications occur at frequencies between 10⁻³ and 10⁻⁵ in unselected cultures (2) but are lost at rates of up to 1,000-fold more often if they do not provide a selective advantage (19). Duplications also result in replichore imbalance. However, the effect of duplications on fitness in relation to replichore balance has not been investigated.A major hurdle in studying the effects of duplications on fitness is that if the duplication is detrimental, haploid revertants will outgrow the parental strain containing the duplication. To circumvent this problem, we used a set of 11 isogenic Salmonella enterica strains with transposon-held duplications (5) (Fig. ​(Fig.11 and Table ​Table1).1). Culturing these strains in the presence of chloramphenicol selects for the duplication because if the duplication collapses, the transposon is lost and the cells become chloramphenicol sensitive. As the size of the duplications in these strains varies, so does the amount of introduced replichore imbalance, ranging from 5° to 23°. In addition, fitness effects due to the location of the duplication versus the size of the duplication were also discerned, as the collection includes duplications of similar sizes located in different regions of the chromosome.Open in a separate windowFIG. 1.Genetic map of S. enterica serovar Typhimurium LT2 showing genes used as endpoints in constructing transposon-held duplications of the regions between genes. Balanced replichores are indicated by the symmetry of the oriC-Ter axis. Duplications increase the length of one replichore relative to the other, imbalancing axis symmetry.

TABLE 1.

Properties of the S. enterica strains used in this study

Pulsed-Field Gel Electrophoresis Profile Changes Resulting from Spontaneous Chromosomal Deletions in Enterohemorrhagic Escherichia coli O157:H7 during Passage in Cattle

A total of 905 enterohemorrhagic Escherichia coli (EHEC) O157:H7 isolates that were recovered from experimentally infected cattle, in addition to the inoculated strain, were analyzed by pulsed-field gel electrophoresis (PFGE). Twelve PFGE profiles other than that of the inoculated strain were observed. We successfully identified five distinct chromosomal deletions that affected the PFGE profiles using whole-genome PCR scanning and DNA sequencing analysis. The changes in PFGE profiles of EHEC O157:H7 isolates after passage through the intestinal tract of cattle were partially generated by deletion of chromosomal regions.Enterohemorrhagic Escherichia coli (EHEC) O157:H7 causes hemorrhagic colitis and hemolytic-uremic syndrome in humans worldwide (18). Cattle are considered the primary reservoir for this pathogen and play a central role in transmission to humans (6). Healthy cattle transiently carry EHEC O157:H7 and shed the bacteria in their feces (5, 7). Human infections have been associated with the consumption of contaminated meat and milk, direct contact with cattle, and the consumption of vegetables, fruits, and water contaminated with cattle manure (6).Because of its high discriminatory power, pulsed-field gel electrophoresis (PFGE) has been widely employed as a molecular typing method in many epidemiological investigations to identify various outbreaks and routes of transmission of EHEC O157:H7 (1, 12, 15, 17). Simpson''s index of diversity (9) was reported to be >0.985 in previous studies (1, 15), supporting the identification of richness (the number of types among isolates) and evenness (the relative distribution of individual strains among the different types) of molecular typing using PFGE.Instability of the PFGE patterns of EHEC O157:H7 isolates has been reported. Changes in PFGE patterns were observed among strains after repeated subculturing and prolonged storage at room temperature (11). Loss of Shiga toxin genes and a large-scale inversion within the genome were identified as genetic events generating changes in PFGE patterns in vitro (10, 13). Shifts in the genotypes of EHEC O157:H7 clinical isolates from patients and cattle have been reported (3, 14). This phenomenon was also observed in EHEC O157:H7 experimental infections of cattle. Spontaneous curing of a 90-kb plasmid resulted in the loss of two restricted fragments from the PFGE profiles of EHEC O157:H7 isolates obtained from experimentally infected cattle (2). The purpose of the present study was to identify the genetic events affecting the PFGE patterns of EHEC O157:H7 after passage through the intestinal tract of cattle, especially for restriction fragments that are >90 kb long.Four 5-month-old Holstein steers were housed individually in climate-controlled biosafety level 2 containment barns in accordance with the guidelines for animal experimentation defined by the National Institute of Animal Health of Japan. The pens had individual floor drains and were cleaned twice daily with water and disinfectant. All animals were healthy and culture negative for EHEC O157:H7 strains, as determined by a previously described technique (2), prior to inoculation.EHEC O157:H7 strain Sakai-215 (12, 23), which was isolated from an outbreak in Sakai, Osaka Prefecture, in 1996 was used for inoculation. This strain harbors the genes encoding Stx1 and Stx2. A spontaneous resistant strain was selected with nalidixic acid in order to facilitate the recovery of this strain from fecal samples. All calves were inoculated using a stomach tube with an exponential-phase culture (10⁹ CFU) of the nalidixic acid-resistant Sakai-215 strain. Fecal samples were collected from the four calves daily for 45 days. Fecal culturing was performed as described previously (2). Eight non-sorbitol-fermenting colonies were selected daily from each animal and identified as EHEC O157:H7 colonies by routine diagnostic methods (25).All animals were clinically normal throughout the experimental period. The EHEC O157:H7-inoculated calves (calves 1 to 4) were culture positive for the organism 24 h after inoculation. Intermittent fecal shedding by the calves was observed until 27, 32, 26, and 39 days postinoculation for calves 1, 2, 3, and 4, respectively (Fig. ​(Fig.1).1). The numbers of EHEC O157:H7 isolates recovered from calves 1, 2, 3, and 4 were 200, 224, 200, and 281, respectively.Open in a separate windowFIG. 1.Changes in PFGE profiles of EHEC O157:H7 isolates recovered from calves 1 (A), 2 (B), 3 (C), and 4 (D). The absence of a bar indicates that no EHEC O157:H7 was detected. The open horizontal bars under the vertical bars indicate that the eight isolates obtained on a day were obtained from the enrichment culture.A total of 905 recovered isolates in addition to the inoculated strain were used for PFGE analysis. Genomic DNA from each EHEC O157:H7 isolate was prepared using the method of Persing (Mayo Clinic, Rochester, MN) described by Rice et al. (20). Agarose-embedded chromosomal DNA was cleaved with XbaI by following the manufacturer''s instructions. PFGE was performed in a 0.85% megabase agarose gel, using a CHEF DR III apparatus (Bio-Rad Laboratories). The pulse time was increased from 12 to 35 s for 18 h. The PFGE profiles of all of the EHEC O157:H7 isolates recovered from the four calves were compared with that of the inoculated strain. The number of band differences was determined by enumerating the loss and addition of fragments (22).Two hundred eighty-nine isolates had PFGE profiles different from that of the inoculated strain, and 12 distinct PFGE profiles were identified for these isolates (Table ​(Table1).1). The fact that only one to three band differences were observed for the 12 profiles suggested that these isolates were closely related (22) and were variants of the inoculated strain. In addition, the pens had individual floor drains and were cleaned twice daily with water and disinfectant, which reduced the likelihood of introduction of novel EHEC O157:H7 strains. We designated the PFGE profiles A to L. PFGE profiles A, C, and H were obtained for all four calves and accounted for 30.4% of the 905 isolates recovered. Different PFGE profiles were obtained for all animals at least 2 days postinoculation (Fig. ​(Fig.1).1). All eight isolates from calf 2 collected on day 15 postinoculation and from calf 3 collected on days 22 and 23 postinoculation had PFGE profiles different from that of the original isolate (Fig. ​(Fig.1).1). The isolates that had the same PFGE profile as the inoculated strain were detected again later.

TABLE 1.

Temporal distribution of PFGE profiles of EHEC O157:H7 isolates recovered from experimentally infected cattle

Wide Variation in the Multiplicity of HIV-1 Infection among Injection Drug Users

Recent studies indicate that sexual transmission of human immunodeficiency virus type 1 (HIV-1) generally results from productive infection by only one virus, a finding attributable to the mucosal barrier. Surprisingly, a recent study of injection drug users (IDUs) from St. Petersburg, Russia, also found most subjects to be acutely infected by a single virus. Here, we show by single-genome amplification and sequencing in a different IDU cohort that 60% of IDU subjects were infected by more than one virus, including one subject who was acutely infected by at least 16 viruses. Multivariant transmission was more common in IDUs than in heterosexuals (60% versus 19%; odds ratio, 6.14; 95% confidence interval [CI], 1.37 to 31.27; P = 0.008). These findings highlight the diversity in HIV-1 infection risks among different IDU cohorts and the challenges faced by vaccines in protecting against this mode of infection.Elucidation of virus-host interactions during and immediately following the transmission event is one of the great challenges and opportunities in human immunodeficiency virus (HIV)/AIDS prevention research (14-16, 31, 34, 45). Recent innovations involving single-genome amplification (SGA), direct amplicon sequencing, and phylogenetic inference based on a model of random virus evolution (18-20, 43) have allowed for the identification of transmitted/founder viruses that actually cross from donor to recipient, leading to productive HIV type 1 (HIV-1) infection. Our laboratory and others have made the surprising finding that HIV-1 transmission results from productive infection by a single transmitted/founder virus (or virally infected cell) in ∼80% of HIV-infected heterosexuals and in ∼60% of HIV-infected men who have sex with men (MSM) (1, 13, 18, 24). These studies thus provided a precise quantitative estimate for the long-recognized genetic bottleneck in HIV-1 transmission (6, 11-13, 17, 25, 28, 30, 35, 38, 42, 47-49) and a plausible explanation for the low acquisition rate per coital act and for graded infection risks associated with different exposure routes and behaviors (15, 36).In contrast to sexual transmission of HIV-1, virus transmission resulting from injection drug use has received relatively little attention (2, 3, 29, 42) despite the fact that injection drug use-associated transmission accounts for as many as 10% of new infections globally (26, 46). We hypothesized that SGA strategies developed for identifying transmitted/founder viruses following mucosal acquisition are applicable to deciphering transmission events following intravenous inoculation and that, due to the absence of a mucosal barrier, injection drug users (IDUs) exhibit a higher frequency of multiple-variant transmission and a wider range in numbers of transmitted viruses than do acutely infected heterosexual subjects. We obtained evidence in support of these hypotheses from the simian immunodeficiency virus (SIV)-Indian rhesus macaque infection model, where we showed that discrete low-diversity viral lineages emanating from single or multiple transmitted/founder viruses could be identified following intravenous inoculation and that the rectal mucosal barrier to infection was 2,000- to 20,000-fold greater than with intravenous inoculation (19). However, we also recognized potentially important differences between virus transmission in Indian rhesus macaques and virus transmission in humans that could complicate an IDU acquisition study. For example, in the SIV macaque model, the virus inocula can be well characterized genetically and the route and timing of virus exposure in relation to plasma sampling precisely defined, whereas in IDUs, the virus inoculum is generally undefined and the timing of virus infection only approximated based on clinical history and seroconversion testing (8). In addition, IDUs may have additional routes of potential virus acquisition due to concomitant sexual activity. Finally, there is a paucity of IDU cohorts for whom incident infection is monitored sufficiently frequently and clinical samples are collected often enough to allow for the identification and enumeration of transmitted/founder viruses. To address these special challenges, we proposed a pilot study of 10 IDU subjects designed to determine with 95% confidence if the proportion of multivariant transmissions in IDUs was more than 2-fold greater than the 20% frequency established for heterosexual transmission (1, 13, 18, 24). A secondary objective of the study was to determine whether the range in numbers of transmitted/founder viruses in IDUs exceeded the 1-to-6 range observed in heterosexuals (1, 13, 18, 24). To ensure comparability among the studies, we employed SGA-direct amplicon sequencing approaches, statistical methods, and power calculations identical to those that we had used previously to enumerate transmitted/founder viruses in heterosexual and MSM cohorts (1, 13, 18, 20, 24).We first surveyed investigators representing acute-infection cohorts in the United States, Canada, Russia, and China; only one cohort—the Montreal Primary HIV Infection Cohort (41)—had IDU clinical samples and clinical data available for study. The Montreal cohort of subjects with acute and early-stage HIV-1 infection was established in 1996 and recruits subjects from both academic and private medical centers throughout the city. Injection drug use is an important contributing factor to Montreal''s HIV burden, with IDUs comprising approximately 20% of the city''s AIDS cases and 35% of the cohort (21, 40, 41). A large proportion of Montreal''s IDUs use injection cocaine, with 50 to 69% of subjects reporting cocaine as their injection drug of choice (4, 5, 9, 22, 23).Subjects with documented serological evidence of recent HIV-1 infection and a concurrent history of injection drug use were selected for study. These individuals had few or no reported risk factors for sexual HIV-1 acquisition. Clinical history and laboratory tests of HIV-1 viremia and antibody seroconversion were used to determine the Fiebig clinical stage (8) and to estimate the date of infection (Table ​(Table1).1). One subject was determined to be in Fiebig stage III, one subject was in Fiebig stage IV, five subjects were in Fiebig stage V, and three subjects were in Fiebig stage VI. We performed SGA-direct amplicon sequencing on stored plasma samples and obtained a total of 391 3′ half-genomes (median, 25 per subject; range, 19 to 167). Nine of these sequences contained large deletions or were G-to-A hypermutated and were excluded from subsequent analysis. Sequences were aligned, visually inspected using the Highlighter tool (www.hiv.lanl.gov/content/sequence/HIGHLIGHT/highlighter.html), and analyzed by neighbor-joining (NJ) phylogenetic-tree construction. A composite NJ tree of full-length gp160 env sequences from all 10 subjects (Fig. ​(Fig.1A)1A) revealed distinct patient-specific monophyletic lineages, each with high bootstrap support and separated from the others by a mean genetic distance of 10.79% (median, 11.29%; range, 3.00 to 13.42%). Maximum within-patient env gene diversity ranged from 0.23% to 3.34% (Table ​(Table1).1). Four subjects displayed distinctly lower within-patient maximum env diversities (0.23 to 0.49%) than the other six subjects (1.48% to 3.34%). The lower maximum env diversities in the former group are consistent with infection either by a single virus or by multiple closely related viruses, while the higher diversities can be explained only by transmission of more than one virus based on empirical observations (1, 13, 18, 24) and mathematical modeling (18, 20).Open in a separate windowFIG. 1.NJ trees and Highlighter plots of HIV-1 gp160 env sequences. (A) Composite tree of 382 gp160 env sequences from all study subjects. The numerals at the nodes indicate bootstrap values for which statistical support exceeded 70%. (B) Subject ACT54869022 sequences suggest productive infection by a single virus (V1). (C) Subject HDNDRPI032 sequences suggest productive infection by as many as three viruses. (D) Subject HDNDRPI001 sequences suggest productive infection by at least five viruses with extensive interlineage recombination. Sequences are color coded to indicate viral progeny from distinct transmitted/founder viruses. Recombinant virus sequences are depicted in black. Methods for SGA, sequencing, model analysis, Highlighter plotting, and identification of transmitted/founder virus lineages are described elsewhere (18, 20, 24, 44). The horizontal scale bars represent genetic distance. nt, nucleotide.

TABLE 1.

Subject demographics and HIV-1 envelope analysis results

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

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

TABLE 1.

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