Comparative Analysis of Myxococcus Predation on Soil Bacteria

Predator-prey relationships among prokaryotes have received little attention but are likely to be important determinants of the composition, structure, and dynamics of microbial communities. Many species of the soil-dwelling myxobacteria are predators of other microbes, but their predation range is poorly characterized. To better understand the predatory capabilities of myxobacteria in nature, we analyzed the predation performance of numerous Myxococcus isolates across 12 diverse species of bacteria. All predator isolates could utilize most potential prey species to effectively fuel colony expansion, although one species hindered predator swarming relative to a control treatment with no growth substrate. Predator strains varied significantly in their relative performance across prey types, but most variation in predatory performance was determined by prey type, with Gram-negative prey species supporting more Myxococcus growth than Gram-positive species. There was evidence for specialized predator performance in some predator-prey combinations. Such specialization may reduce resource competition among sympatric strains in natural habitats. The broad prey range of the Myxococcus genus coupled with its ubiquity in the soil suggests that myxobacteria are likely to have very important ecological and evolutionary effects on many species of soil prokaryotes.Predation plays a major role in shaping both the ecology and evolution of biological communities. The population and evolutionary dynamics of predators and their prey are often tightly coupled and can greatly influence the dynamics of other organisms as well (1). Predation has been invoked as a major cause of diversity in ecosystems (11, 12). For example, predators may mediate coexistence between superior and inferior competitors (2, 13), and differential trajectories of predator-prey coevolution can lead to divergence between separate populations (70).Predation has been investigated extensively in higher organisms but relatively little among prokaryotes. Predation between prokaryotes is one of the most ancient forms of predation (27), and it has been proposed that this process may have been the origin of eukaryotic cells (16). Prokaryotes are key players in primary biomass production (44) and global nutrient cycling (22), and predation of some prokaryotes by others is likely to significantly affect these processes. Most studies of predatory prokaryotes have focused on Bdellovibrionaceae species (e.g., see references 51, 55, and 67). These small deltaproteobacteria prey on other Gram-negative cells, using flagella to swim rapidly until they collide with a prey cell. After collision, the predator cells then enter the periplasmic space of the prey cell, consume the host cell from within, elongate, and divide into new cells that are released upon host cell lysis (41). Although often described as predatory, the Bdellovibrionaceae may also be considered to be parasitic, as they typically depend (apart from host-independent strains that have been observed [60]) on the infection and death of their host for their reproduction (47).In this study, we examined predation among the myxobacteria, which are also deltaproteobacteria but constitute a monophyletic clade divergent from the Bdellovibrionaceae (17). Myxobacteria are found in most terrestrial soils and in many aquatic environments as well (17, 53, 74). Many myxobacteria, including the model species Myxococcus xanthus, exhibit several complex social traits, including fruiting body formation and spore formation (14, 18, 34, 62, 71), cooperative swarming with two motility systems (64, 87), and group (or “wolf pack”) predation on both bacteria and fungi (4, 5, 8, 9, 15, 50). Using representatives of the genus Myxococcus, we tested for both intra- and interspecific variation in myxobacterial predatory performance across a broad range of prey types. Moreover, we examined whether prey vary substantially in the degree to which they support predatory growth by the myxobacteria and whether patterns of variation in predator performance are constant or variable across prey environments. The latter outcome may reflect adaptive specialization and help to maintain diversity in natural populations (57, 59).Although closely related to the Bdellovibrionaceae (both are deltaproteobacteria), myxobacteria employ a highly divergent mode of predation. Myxobacteria use gliding motility (64) to search the soil matrix for prey and produce a wide range of antibiotics and lytic compounds that kill and decompose prey cells and break down complex polymers, thereby releasing substrates for growth (66). Myxobacterial predation is cooperative both in its “searching” component (6, 31, 82; for details on cooperative swarming, see reference 64) and in its “handling” component (10, 29, 31, 32), in which secreted enzymes turn prey cells into consumable growth substrates (56, 83). There is evidence that M. xanthus employs chemotaxis-like genes in its attack on prey cells (5) and that predation is stimulated by close contact with prey cells (48).Recent studies have revealed great genetic and phenotypic diversity within natural populations of M. xanthus, on both global (79) and local (down to centimeter) scales (78). Phenotypic diversity includes variation in social compatibility (24, 81), the density and nutrient thresholds triggering development (33, 38), developmental timing (38), motility rates and patterns (80), and secondary metabolite production (40). Although natural populations are spatially structured and both genetic diversity and population differentiation decrease with spatial scale (79), substantial genetic diversity is present even among centimeter-scale isolates (78). No study has yet systematically investigated quantitative natural variation in myxobacterial predation phenotypes across a large number of predator genotypes.Given the previous discovery of large variation in all examined phenotypes, even among genetically extremely similar strains, we anticipated extensive predatory variation as well. Using a phylogenetically broad range of prey, we compared and contrasted the predatory performance of 16 natural M. xanthus isolates, sampled from global to local scales, as well as the commonly studied laboratory reference strain DK1622 and representatives of three additional Myxococcus species: M. flavescens (86), M. macrosporus (42), and M. virescens (63) (Table ​(Table1).1). In particular, we measured myxobacterial swarm expansion rates on prey lawns spread on buffered agar (31, 50) and on control plates with no nutrients or with prehydrolyzed growth substrate.

TABLE 1.

List of myxobacteria used, with geographical origin

Evidence for a New Avian Paramyxovirus Serotype 10 Detected in Rockhopper Penguins from the Falkland Islands

The biological, serological, and genomic characterization of a paramyxovirus recently isolated from rockhopper penguins (Eudyptes chrysocome) suggested that this virus represented a new avian paramyxovirus (APMV) group, APMV10. This penguin virus resembled other APMVs by electron microscopy; however, its viral hemagglutination (HA) activity was not inhibited by antisera against any of the nine defined APMV serotypes. In addition, antiserum generated against this penguin virus did not inhibit the HA of representative viruses of the other APMV serotypes. Sequence data produced using random priming methods revealed a genomic structure typical of APMV. Phylogenetic evaluation of coding regions revealed that amino acid sequences of all six proteins were most closely related to APMV2 and APMV8. The calculation of evolutionary distances among proteins and distances at the nucleotide level confirmed that APMV2, APMV8, and the penguin virus all were sufficiently divergent from each other to be considered different serotypes. We propose that this isolate, named APMV10/penguin/Falkland Islands/324/2007, be the prototype virus for APMV10. Because of the known problems associated with serology, such as antiserum cross-reactivity and one-way immunogenicity, in addition to the reliance on the immune response to a single protein, the hemagglutinin-neuraminidase, as the sole base for viral classification, we suggest the need for new classification guidelines that incorporate genome sequence comparisons.Viruses from the Paramyxoviridae family have caused disease in humans and animals for centuries. Over the last 40 years, many paramyxoviruses isolated from animals and people have been newly described (16, 17, 22, 29, 31, 32, 36, 42, 44, 46, 49, 58, 59, 62-64). Viruses from this family are pleomorphic, enveloped, single-stranded, nonsegmented, negative-sense RNA viruses that demonstrate serological cross-reactivity with other paramyxoviruses related to them (30, 46). The subfamily Paramyxovirinae is divided into five genera: Respirovirus, Morbillivirus, Rubulavirus, Henipavirus, and Avulavirus (30). The Avulavirus genus contains nine distinct avian paramyxovirus (APMV) serotypes (Table ​(Table1),1), and information on the discovery of each has been reported elsewhere (4, 6, 7, 9, 12, 34, 41, 50, 51, 60, 68).

TABLE 1.

Characteristics of prototype viruses APMV1 to APMV9 and the penguin virus

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

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

TABLE 1.

Primers used in multiplex amplification for the TMLGT assay

Structural Characterization of the Dual Glycan Binding Adeno-Associated Virus Serotype 6

The three-dimensional structure of adeno-associated virus (AAV) serotype 6 (AAV6) was determined using cryo-electron microscopy and image reconstruction and using X-ray crystallography to 9.7- and 3.0-Å resolution, respectively. The AAV6 capsid contains a highly conserved, eight-stranded (βB to βI) β-barrel core and large loop regions between the strands which form the capsid surface, as observed in other AAV structures. The loops show conformational variation compared to other AAVs, consistent with previous reports that amino acids in these loop regions are involved in differentiating AAV receptor binding, transduction efficiency, and antigenicity properties. Toward structure-function annotation of AAV6 with respect to its unique dual glycan receptor (heparan sulfate and sialic acid) utilization for cellular recognition, and its enhanced lung epithelial transduction compared to other AAVs, the capsid structure was compared to that of AAV1, which binds sialic acid and differs from AAV6 in only 6 out of 736 amino acids. Five of these residues are located at or close to the icosahedral 3-fold axis of the capsid, thereby identifying this region as imparting important functions, such as receptor attachment and transduction phenotype. Two of the five observed amino acids are located in the capsid interior, suggesting that differential AAV infection properties are also controlled by postentry intracellular events. Density ordered inside the capsid, under the 3-fold axis in a previously reported, conserved AAV DNA binding pocket, was modeled as a nucleotide and a base, further implicating this capsid region in AAV genome recognition and/or stabilization.Adeno-associated viruses (AAVs) are nonpathogenic single-stranded DNA (ssDNA) parvoviruses that belong to the Dependovirus genus and require helper viruses, such as Adenovirus or Herpesvirus, for lytic infection (4, 8, 22, 67). These viruses package a genome of ∼4.7 kb inside an icosahedral capsid (∼260 Å in diameter) with a triangulation number equal to 1 assembled from a total of 60 copies of their overlapping capsid viral protein (VP) 1 (VP1), VP2, and VP3 in a predicted ratio of 1:1:8/10 (10). The VPs are encoded from a cap open reading frame (ORF). VP3 is 61 kDa and constitutes 90% of the capsid''s protein composition. The less abundant VPs, VP1 (87 kDa) and VP2 (73 kDa), share the same C-terminal amino acid sequence with VP3 but have additional N-terminal sequences. A rep ORF codes for four overlapping proteins required for replication and DNA packaging.To date, more than 100 AAV isolates have been identified (21). Among the human and nonhuman primate AAVs isolated, 12 serotypes (AAV serotype 1 [AAV1] to AAV12) have been described and are classified into six phylogenetic clades on the basis of their VP sequences and antigenic reactivities, with AAV4 and AAV5 considered to be clonal isolates (21). AAV1 and AAV6, which represent clade A, differ by only 6 out of 736 VP1 amino acids (5 amino acids within VP3) and are antigenically cross-reactive. Other clade representatives include AAV2 (clade B), AAV2-AAV3 hybrid (clade C), AAV7 (clade D), AAV8 (clade E), and AAV9 (clade F) (21).The AAVs are under development as clinical gene delivery vectors (e.g., see references 5, 9, 12, 13, 24, 25, 53, and 61), with AAV2, the prototype member of the genus, being the most extensively studied serotype for this application. AAV2 has been successfully used to treat several disorders, but its broad tissue tropism makes it less effective for tissue-specific applications and the prevalence of preexisting neutralizing antibodies in the human population (11, 43) limits its utilization, especially when readministration is required to achieve a therapeutic outcome. Efforts have thus focused on characterizing the capsid-associated tissue tropism and transduction properties conferred by the capsid of representative serotypes of other clades (21). Outcomes of these studies include the observation that AAV1 and AAV6, for example, transduce liver, muscle, and airway epithelial cells more efficiently (e.g., up to 200-fold) than AAV2 (27, 28, 30). In addition, the six residues (Table ​(Table1)1) that differ between the VPs of AAV1 and AAV6 (a natural recombinant of AAV1 and AAV2 [56]) confer functional disparity between these two viruses. For example, AAV6 shows ∼3-fold higher lung cell epithelium transduction than AAV1 (27), and AAV1 and AAV6 bind terminally sialylated proteoglycans as their primary receptor, whereas AAV6 additionally binds to heparan sulfate (HS) proteoglycans with moderate affinity (70, 71). Therefore, a comparison of the AAV1 and AAV6 serotypes and, in particular, their capsid structures can help pinpoint the capsid regions that confer differences in cellular recognition and tissue transduction.

TABLE 1.

Amino acid differences between AAV1 and AAV6 and their reported mutants

Improved Molecular Detection of Angiostrongylus cantonensis in Mollusks and Other Environmental Samples with a Species-Specific Internal Transcribed Spacer 1-Based TaqMan Assay

Angiostrongylus cantonensis is the most common cause of human eosinophilic meningitis. Humans become infected by ingesting food items contaminated with third-stage larvae that develop in mollusks. We report the development of a real-time PCR assay for the species-specific identification of A. cantonensis in mollusk tissue.Angiostrongylus cantonensis is the most common agent associated with eosinophilic meningitis in humans. Young adult worms develop in the brains of rodents and are carried to pulmonary arteries to reach sexual maturity. Eggs are laid in lung tissues, and first-stage (L1) larvae break into air spaces, migrate to the trachea, are swallowed, and are passed with rodent feces. The L1 larvae must infect mollusks to develop into third-stage (L3) larvae; L3 is the infective stage for rodents and other mammals. Humans become infected by ingesting raw produce contaminated with L3 larvae or infected raw or undercooked mollusks or paratenic hosts. The immature worms remain in the human brain, creating tissue damage and inflammation (2, 19, 21).A. cantonensis is endemic in Southeast Asia, parts of the Caribbean, and the Pacific Islands, including Hawaii (7, 12, 15-17). The worm has been detected in host animals in Louisiana (5, 14) and in one human patient from New Orleans (18), but it is currently unclear to what extent the nematode has spread into other U.S. states (8, 9). Ascertaining the geographic presence of the parasite is important to manage and prevent new cases of eosinophilic meningitis associated with ingestion of infective larvae (12, 18).Detection of A. cantonensis in mollusks can be performed by releasing the larvae from the tissue with pepsin digestion (11). However, that procedure requires access to living mollusks, which complicates analysis of large numbers of samples. After a recent outbreak of angiostrongyliasis in Hawaii (12), we developed a conventional PCR assay and applied it to survey the Hawaiian mollusk population using frozen tissue (20). That PCR assay, as well as morphological identification using pepsin digestion, can only identify the larvae on the superfamily level, so additional molecular work is required for species-specific classification. Here we describe a new real-time PCR assay that allows for a direct detection of A. cantonensis at the species level.The 18S rRNA gene is too conserved among nematode species to allow species-specific detection. The first and second internal transcribed spacers (ITS1 and ITS2) are comparatively more variable than the rRNA coding regions and have thus been used for differentiation of closely related species (1, 4, 6, 10, 22, 23). We PCR amplified and sequenced ITS1 from A. costaricensis (two laboratory strains from Costa Rica and Brazil), A. vasorum (from naturally infected hosts in United Kingdom), and A. cantonensis from three geographical regions (one laboratory strain from Japan plus nine environmental isolates from Hawaii and New Orleans, LA) to assess the variability of this potential PCR target. The oligonucleotide primers used were AngioF1674 (5′-GTCGTAACAAGGTATCTGTAGGTG-3′) and 58SR4 (5′-TAGCTGCGTTTTTCATCGATA-3′). The reaction mixtures contained 0.4 μM each primer and AmpliTaq Gold PCR master mix (Applied Biosystems, Foster City, CA) and were cycled 45 times at 94°C for 30 s, 65°C for 30 s, and 72°C for 1 min. PCR products were cloned into pCR2.1 vectors using the TOPO cloning technique (Invitrogen, Carlsbad, CA) and sequenced on both strands as described elsewhere (20).The sequence analysis revealed high interspecific and low intraspecific variability. A TaqMan assay targeting ITS1 was then designed using Primer Express version 2.3 (Applied Biosystems, Foster City, CA). The real-time PCR assay was performed in a 20-μl total volume containing Platinum qPCR Supermix (Invitrogen, Carlsbad, CA), 0.2 μM (each) primers AcanITS1F1 (5′-TTCATGGATGGCGAACTGATAG-3′) and AcanITS1R1 (5′-GCGCCCATTGAAACATTATACTT-3′), and 0.05 μM the TaqMan probe AcanITS1P1 (5′-6-carboxyfluorescein-ATCGCATATCTACTATACGCATGTGACACCTG-BHQ-3′). The standard cycling conditions for TaqMan assays were used (i.e., 40 cycles of 95°C for 15 s and 60°C for 1 min).We evaluated the real-time PCR assay with a set of 26 Parmarion martensi slugs from Hawaii. Seventeen slugs were positive for L3 larvae as determined by pepsin digestion, and nine slugs were negative. DNA was extracted from approximately 25 mg of tissue of each slug using the DNeasy tissue and blood DNA extraction kit (Qiagen, Inc., Valencia, CA). The real-time PCR performed on this set of samples returned an identical result to the morphological analysis. The real-time PCR amplified only DNA from A. cantonensis and did not react with DNA from other nematode species (Table ​(Table1).1). The detection limit of the assay was determined by serially diluting a recombinant plasmid containing the ITS1 sequence to less than 1 copy per μl of sample. The real-time PCR reliably detected down to 10 plasmid copies in the reaction.

TABLE 1.

Comparison of conventional and real-time PCR for detection of Angiostrongylus cantonensis in mollusks and nematode samples

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

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

Growth of Arthrobacter sp. Strain JBH1 on Nitroglycerin as the Sole Source of Carbon and Nitrogen

Arthrobacter sp. strain JBH1 was isolated from nitroglycerin-contaminated soil by selective enrichment. Detection of transient intermediates and simultaneous adaptation studies with potential intermediates indicated that the degradation pathway involves the conversion of nitroglycerin to glycerol via 1,2-dinitroglycerin and 1-mononitroglycerin, with concomitant release of nitrite. Glycerol then serves as the source of carbon and energy.Nitroglycerin (NG) is manufactured widely for use as an explosive and a pharmaceutical vasodilator. It has been found as a contaminant in soil and groundwater (7, 9). Due to NG''s health effects as well as its highly explosive nature, NG contamination in soils and groundwater poses a concern that requires remedial action (3). Natural attenuation and in situ bioremediation have been used for remediation in soils contaminated with certain other explosives (16), but the mineralization of NG in soil and groundwater has not been reported.To date, no pure cultures able to grow on NG as the sole carbon, energy, and nitrogen source have been isolated. Accashian et al. (1) observed growth associated with the degradation of NG under aerobic conditions by a mixed culture originating from activated sludge. The use of NG as a source of nitrogen has been studied in mixed and pure cultures during growth on alternative sources of carbon and energy (3, 9, 11, 20). Under such conditions, NG undergoes a sequential denitration pathway in which NG is transformed to 1,2-dinitroglycerin (1,2DNG) or 1,3DNG followed by 1-mononitroglycerin (1MNG) or 2MNG and then glycerol, under both aerobic and anaerobic conditions (3, 6, 9, 11, 20), and the enzymes involved in denitration have been characterized in some detail (4, 8, 15, 21). Pure cultures capable of completely denitrating NG as a source of nitrogen when provided additional sources of carbon include Bacillus thuringiensis/cereus and Enterobacter agglomerans (11) and a Rhodococcus species (8, 9). Cultures capable of incomplete denitration to MNG in the presence of additional carbon sources were identified as Pseudomonas putida, Pseudomonas fluorescens (4), an Arthobacter species, a Klebsiella species (8, 9), and Agrobacterium radiobacter (20).Here we describe the isolation of bacteria able to degrade NG as the sole source of carbon, nitrogen, and energy. The inoculum for selective enrichment was soil historically contaminated with NG obtained at a facility that formerly manufactured explosives located in the northeastern United States. The enrichment medium consisted of minimal medium prepared as previously described (17) supplemented with NG (0.26 mM), which was synthesized as previously described (18). During enrichment, samples of the inoculum (optical density at 600 nm [OD₆₀₀] ∼ 0.03) were diluted 1/16 in fresh enrichment medium every 2 to 3 weeks. Isolates were obtained by dilution to extinction in NG-supplemented minimal medium. Cultures were grown under aerobic conditions in minimal medium at pH 7.2 and 23°C or in tryptic soy agar (TSA; 1/4 strength).Early stages of enrichment cultures required extended incubation with lag phases of over 200 h and exhibited slow degradation of NG (less than 1 μmol substrate/mg protein/h). After a number of transfers over 8 months, the degradation rates increased substantially (2.2 μmol substrate/mg protein/h). A pure culture capable of growth on NG was identified based on 16S rRNA gene analysis (504 bp) as an Arthrobacter species with 99.5% similarity to Arthrobacter pascens (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"GU246730","term_id":"289474572","term_text":"GU246730"}}GU246730). Purity of the cultures was confirmed microscopically and by formation of a single colony type on TSA plates. 16S gene sequencing and identification were done by MIDI Labs (Newark, DE) and SeqWright DNA Technology Services (Houston, TX). The Arthrobacter cells stained primarily as Gram-negative rods with a small number of Gram-positive cocci (data not shown); Gram variability is also a characteristic of the closely related Arthrobacter globiformis (2, 19). The optimum growth temperature is 30°C, and the optimum pH is 7.2. Higher pH values were not investigated because NG begins to undergo hydrolysis above pH 7.5 (data not shown). The isolated culture can grow on glycerol, acetate, succinate, citrate, and lactate, with nitrite as the nitrogen source. Previous authors described an Arthrobacter species able to use NG as a nitrogen source in the presence of additional sources of carbon. However, only dinitroesters were formed, and complete mineralization was not achieved (9).To determine the degradation pathway, cultures of the isolated strain (5 ml of inoculum grown on NG to an OD₆₀₀ of 0.3) were grown in minimal medium (100 ml) supplemented with NG at a final concentration of 0.27 mM. Inoculated bottles and abiotic controls were continuously mixed, and NG, 1,2DNG, 1,3DNG, 1MNG, 2MNG, nitrite, nitrate, CO₂, total protein, and optical density were measured at appropriate intervals. Nitroesters were analyzed with an Agilent high-performance liquid chromatograph (HPLC) equipped with an LC-18 column (250 by 4.6 mm, 5 μm; Supelco) and a UV detector at a wavelength of 214 nm (13). Methanol-water (50%, vol/vol) was used as the mobile phase at a flow rate of 1 ml/min. Nitrite and nitrate were analyzed with an ion chromatograph (IC) equipped with an IonPac AS14A anion-exchange column (Dionex, CA) at a flow rate of 1 ml/min. Carbon dioxide production was measured with a Micro Oxymax respirometer (Columbus Instruments, OH), and total protein was quantified using the Micro BCA protein assay kit (Pierce Biotechnology, IL) according to manufacturer''s instructions. During the degradation of NG the 1,2DNG concentration was relatively high at 46 and 72 h (Fig. ​(Fig.1).1). 1,3DNG, detected only at time zero, resulted from trace impurities in the NG stock solution. Trace amounts of 1MNG appeared transiently, and trace amounts of 2MNG accumulated and did not disappear. Traces of nitrite at time zero were from the inoculum. The concentration of NG in the abiotic control did not change during the experiment (data not shown).Open in a separate windowFIG. 1.Growth of strain JBH1 on NG. ×, NG; ▵, 1,2DNG; ⋄, 1MNG; □, 2MNG; ○, protein.Results from the experiment described above were used to calculate nitrogen and carbon mass balances (Tables ​(Tables11 and ​and2).2). Nitrogen content in protein was approximated using the formula C₅H₇O₂N (14). Because all nitrogen was accounted for throughout, we conclude that the only nitrogen-containing intermediate compounds are 1,2DNG and 1MNG, which is consistent with previous studies (6, 9, 20). The fact that most of the nitrogen was released as nitrite is consistent with previous reports of denitration catalyzed by reductase enzymes (4, 8, 21). The minor amounts of nitrate observed could be from abiotic hydrolysis (5, 12) or from oxidation of nitrite. Cultures supplemented with glycerol or other carbon sources assimilated all of the nitrite (data not shown).

TABLE 1.

Nitrogen mass balance

Nonpathogenic Simian Immunodeficiency Virus Infection of Sooty Mangabeys Is Not Associated with High Levels of Autologous Neutralizing Antibodies

Simian immunodeficiency virus (SIV) infection of natural-host species, such as sooty mangabeys (SMs), is characterized by a high level of viral replication and a low level of generalized immune activation, despite evidence of an adaptive immune response. Here the ability of SIV-infected SMs to mount neutralizing antibodies (Nab) against autologous virus was compared to that of human immunodeficiency virus type 1 (HIV-1) subtype C-infected subjects. While high levels of Nab were observed in HIV-1 infection, samples obtained at comparable time points from SM exhibited relatively low titers of autologous Nab. Nevertheless, SM plasma with higher Nab titers also contained elevated peripheral CD4⁺ T-cell levels, suggesting a potential immunologic benefit for SMs. These data indicate that AIDS resistance in these primates is not due to high Nab titers and raise the possibility that low levels of Nab might be an inherent feature of natural-host SIV infections.More than 40 species of African nonhuman primates (NHPs) naturally harbor CD4⁺-tropic lentiviruses that are collectively known as simian immunodeficiency viruses (SIVs) and represent the ancestors of the human pathogens human immunodeficiency virus type 1 (HIV-1) and HIV-2. Interestingly, African NHPs infected with their cognate SIV generally do not progress to AIDS, despite high levels of sustained virus replication, with the only known exception being chimpanzee SIV (SIVcpz)-infected chimpanzees (16). Among the natural hosts for SIV infection, the sooty mangabey ([SM] Cercocebus atys) is of particular interest, because cross-species transmission of SM SIV (SIVsm) from this natural host into humans initiated the HIV-2 epidemic in West Africa (17). In addition, SIVsm (herein referred to as SIV) is the ancestor of the rhesus macaque SIV (SIVmac) viruses that are used in disease pathogenesis and vaccination studies in the rhesus macaque model (17). Both naturally infected and experimentally inoculated SMs remain healthy, maintain CD4⁺ T cells, and do not progress to AIDS-like disease, despite sustained high levels of virus replication (31).Nonpathogenic infection of SMs is characterized by low levels of immune activation during the chronic phase of infection, which are reached after a transient immune activation that occurs during primary infection (reviewed in reference 31). These findings have led to the hypothesis that the absence of generalized immune activation in SIV-infected SMs during the chronic phase of infection is an important feature that favors the preservation of CD4⁺ T-cell homeostasis, thereby avoiding disease progression (31). However, most of these earlier studies focused on T cells and innate immune cells, with a significant gap existing in our understanding of whether humoral immunity might also differ between pathogenic and nonpathogenic infections. In HIV-1-infected patients, B cells produce neutralizing antibodies against the infecting (autologous) virus, which drives viral escape, continuous de novo antibody production (26-28, 32), and B-cell dysfunction (24). The striking differences in both the clinical outcomes of infection and the levels of immune activation between SIV-infected SMs and HIV-1-infected humans prompted us to compare the neutralizing antibody (Nab) response against the autologous virus in these two populations. To this end, we utilized a pseudovirus assay that has been used extensively by our group and others to evaluate Nab against HIV-1 and SIV envelope (Env) glycoproteins (15, 19, 22, 26, 28, 32, 33; also unpublished data). All SMs were housed at the Yerkes National Primate Research Center (Atlanta, GA) and maintained in accordance with National Institutes of Health guidelines. The Emory University Animal Care and Use Committee approved these studies. Details of the Zambia Emory HIV Research Project (ZEHRP) have been described elsewhere (2, 10, 21). The Emory University Institutional Review Board and the University of Zambia School of Medicine Research Ethics Committee approved informed-consent and human subject protocols. None of the subjects received antiretroviral therapy during the evaluation period.In HIV-1 infection, autologous Nabs develop to relatively high titers against the newly transmitted virus within the first few months (15, 19, 26-28, 32). Here we sought to test whether a similar increase in Nab titer occurs during nonpathogenic SIV infection of SMs. Samples were obtained from five animals that were inoculated intravenously with plasma from a naturally infected SM as part of a previous study (30). Multiple, biologically functional Envs were cloned from plasma collected at day 14 postinoculation (Table ​(Table1),1), and Nab activity was evaluated in plasma collected at 6 months postinoculation. To facilitate comparison with early HIV-1 infection, Nab activity in plasma was also evaluated between 2 and 9 months against Envs that were cloned between 31 and 88 estimated days after infection from four subtype C HIV-1-infected seroconverters in Zambia (Table ​(Table1).1). Figure ​Figure1A1A demonstrates that Nab activity in plasma diluted 1:100 was readily detectable in all HIV-1-infected subjects at levels approaching 100% neutralization. However, Nab activity in the SM plasma was significantly lower than in the human subjects (median, 10% versus 93%, respectively; P = 0.02). Binding antibody was detected in all five SMs at titers greater than 1:51,200 by enzyme-linked immunosorbent assay (ELISA), demonstrating that all monkeys had seroconverted by 6 months and maintained high titers of binding antibody throughout the evaluation period (Fig. ​(Fig.1B).1B). Thus, the low level of Nab was not due to a diminished humoral immune response.Open in a separate windowFIG. 1.Autologous Nab activity and B-cell proliferation during experimental infection of SMs. (A) Neutralization activity levels in plasma from five SMs (filled black circles), which were experimentally inoculated with plasma from a naturally SIV-infected SM, and four HIV-1-infected Zambian subjects (half-filled squares), who were recently infected through heterosexual contact, are shown. The horizontal bars represent the median for each group. To assess neutralizing activity, pseudoviruses were created by expressing each cloned Env with an HIV-1 env-deficient backbone (ΔSG3). JC53-BL (Tzm-bl) cells were infected with each pseudovirus in the presence or absence of serially diluted autologous plasma. Each point represents the average level of neutralization at a 1:100 dilution of plasma for at least two Env clones (see Table ​Table11 for number of Envs tested). Each neutralization assay was performed twice independently, using duplicate wells. Statistical significance between the groups was determined by a Mann-Whitney test, using GraphPad Prism 5. Longitudinal measurements of endpoint antibody ELISA titers in plasma (filled green circles) (23) (B), autologous neutralization activity in plasma (filled blue diamonds) (C), percentages of Ki-67⁺ CD20⁺ cells in blood (filled black triangles) (D), and percentages of CD20⁺ cells in blood (filled red squares) (E) are shown for the five experimentally inoculated SMs combined. In panel C, each point represents average neutralization at a 1:100 dilution of plasma over time for at least two day 14 Env clones from each SM. For panels D and E, PBMCs were gated by forward and side scatter, and the CD3⁻ CD20⁺ population was assessed for Ki-67 staining (D) by flow cytometry. SP34-2 was used to stain CD3, L27 was used for CD20, and B56 was used for Ki-67 (all from BD Biosciences). Error bars represent the standard errors of the means (SEMs). Plasma viral load peaked at day 14 (data not shown). Filled symbols in panels A through E indicate data generated from experimentally infected SMs.

TABLE 1.

Autologous Nab activity in experimentally SIV-infected SM and acutely HIV-1-infected humans