We have recently reported that Altered Schaedler Flora (ASF) can be used to durably engineer the gut microbiota to reduce ammonia production as an effective modality to reduce morbidity and mortality in the setting of liver injury. Here we investigated the effects of a low protein diet on ASF colonization and its ability to engineer the microbiota. Initially, ASF inoculation was similar between mice fed a normal protein diet or low protein diet, but the outgrowth of gut microbiota differed over the ensuing month. Notable was the inability of the dominant Parabacteroides ASF taxon to exclude other taxa belonging to the Bacteroidetes phylum in the setting of a low protein diet. Instead, a poorly classified yet highly represented Bacteroidetes family, S24-7, returned within 4 weeks of inoculation in mice fed a low protein diet, demonstrating a reduction in ASF resilience in response to dietary stress. Nevertheless, fecal ammonia levels remained significantly lower than those observed in mice on the same low protein diet that received a transplant of normal feces. No deleterious effects were observed in host physiology due to ASF inoculation into mice on a low protein diet. In total, these results demonstrate that low protein diet can have a pronounced effect on engineering the gut microbiota but modulation of ammonia is preserved.

African apes harbour at least six Plasmodium species of the subgenus Laverania, one of which gave rise to human Plasmodium falciparum. Here we use a selective amplification strategy to sequence the genome of chimpanzee parasites classified as Plasmodium reichenowi and Plasmodium gaboni based on the subgenomic fragments. Genome-wide analyses show that these parasites indeed represent distinct species, with no evidence of cross-species mating. Both P. reichenowi and P. gaboni are 10-fold more diverse than P. falciparum, indicating a very recent origin of the human parasite. We also find a remarkable Laverania-specific expansion of a multigene family involved in erythrocyte remodelling, and show that a short region on chromosome 4, which encodes two essential invasion genes, was horizontally transferred into a recent P. falciparum ancestor. Our results validate the selective amplification strategy for characterizing cryptic pathogen species, and reveal evolutionary events that likely predisposed the precursor of P. falciparum to colonize humans.

Lower respiratory tract infection (LRTI) is a major contributor to respiratory failure requiring intubation and mechanical ventilation. LRTI also occurs during mechanical ventilation, increasing the morbidity and mortality of intubated patients. We sought to understand the dynamics of respiratory tract microbiota following intubation and the relationship between microbial community structure and infection.

The vertebrate gut microbiota have been implicated in the metabolism of xenobiotic compounds, motivating studies of microbe-driven metabolism of clinically important drugs. Here, we studied interactions between the microbiota and indomethacin, a nonsteroidal anti-inflammatory drug (NSAID) that inhibits cyclooxygenases (COX) -1 and -2. Indomethacin was tested in both acute and chronic exposure models in mice at clinically relevant doses, which suppressed production of COX-1- and COX-2-derived prostaglandins and caused small intestinal (SI) damage. Deep sequencing analysis showed that indomethacin exposure was associated with alterations in the structure of the intestinal microbiota in both dosing models. Perturbation of the intestinal microbiome by antibiotic treatment altered indomethacin pharmacokinetics and pharmacodynamics, which is probably the result of reduced bacterial β-glucuronidase activity. Humans show considerable inter-individual differences in their microbiota and their responses to indomethacin - thus, the drug-microbe interactions described here provide candidate mediators of individualized drug responses.

Viruses make up a major component of the human microbiota but are poorly understood in the skin, our primary barrier to the external environment. Viral communities have the potential to modulate states of cutaneous health and disease. Bacteriophages are known to influence the structure and function of microbial communities through predation and genetic exchange. Human viruses are associated with skin cancers and a multitude of cutaneous manifestations. Despite these important roles, little is known regarding the human skin virome and its interactions with the host microbiome. Here we evaluated the human cutaneous double-stranded DNA virome by metagenomic sequencing of DNA from purified virus-like particles (VLPs). In parallel, we employed metagenomic sequencing of the total skin microbiome to assess covariation and infer interactions with the virome. Samples were collected from 16 subjects at eight body sites over 1 month. In addition to the microenviroment, which is known to partition the bacterial and fungal microbiota, natural skin occlusion was strongly associated with skin virome community composition. Viral contigs were enriched for genes indicative of a temperate phage replication style and also maintained genes encoding potential antibiotic resistance and virulence factors. CRISPR spacers identified in the bacterial DNA sequences provided a record of phage predation and suggest a mechanism to explain spatial partitioning of skin phage communities. Finally, we modeled the structure of bacterial and phage communities together to reveal a complex microbial environment with a Corynebacterium hub. These results reveal the previously underappreciated diversity, encoded functions, and viral-microbial dynamic unique to the human skin virome.

Target-site selection by retroviral integrase (IN) proteins profoundly affects viral pathogenesis. We describe the solution nuclear magnetic resonance structure of the Moloney murine leukemia virus IN (M-MLV) C-terminal domain (CTD) and a structural homology model of the catalytic core domain (CCD). In solution, the isolated MLV IN CTD adopts an SH3 domain fold flanked by a C-terminal unstructured tail. We generated a concordant MLV IN CCD structural model using SWISS-MODEL, MMM-tree and I-TASSER. Using the X-ray crystal structure of the prototype foamy virus IN target capture complex together with our MLV domain structures, residues within the CCD α2 helical region and the CTD β1-β2 loop were predicted to bind target DNA. The role of these residues was analyzed in vivo through point mutants and motif interchanges. Viable viruses with substitutions at the IN CCD α2 helical region and the CTD β1-β2 loop were tested for effects on integration target site selection. Next-generation sequencing and analysis of integration target sequences indicate that the CCD α2 helical region, in particular P187, interacts with the sequences distal to the scissile bonds whereas the CTD β1-β2 loop binds to residues proximal to it. These findings validate our structural model and disclose IN-DNA interactions relevant to target site selection.

The genomic DNAs of tailed bacteriophages are commonly modified by the attachment of chemical groups. Some forms of DNA modification are known to protect phage DNA from cleavage by restriction enzymes, but others are of unknown function. Recently, the CRISPR-Cas nuclease complexes were shown to mediate bacterial adaptive immunity by RNA-guided target recognition, raising the question of whether phage DNA modifications may also block attack by CRISPR-Cas9. We investigated phage T4 as a model system, where cytosine is replaced with glucosyl-hydroxymethylcytosine (glc-HMC). We first quantified the extent and distribution of covalent modifications in T4 DNA by single-molecule DNA sequencing and enzymatic probing. We then designed CRISPR spacer sequences targeting T4 and found that wild-type T4 containing glc-HMC was insensitive to attack by CRISPR-Cas9 but mutants with unmodified cytosine were sensitive. Phage with HMC showed only intermediate sensitivity. While this work was in progress, another group reported examples of heavily engineered CRISRP-Cas9 complexes that could, in fact, overcome the effects of T4 DNA modification, indicating that modifications can inhibit but do not always fully block attack.

Here we present APOBEC-coupled epigenetic sequencing (ACE-seq), a bisulfite-free method for localizing 5-hydroxymethylcytosine (5hmC) at single-base resolution with low DNA input. The method builds on the observation that AID/APOBEC family DNA deaminase enzymes can potently discriminate between cytosine modification states and exploits the non-destructive nature of enzymatic, rather than chemical, deamination. ACE-seq yielded high-confidence 5hmC profiles with at least 1,000-fold less DNA input than conventional methods. Applying ACE-seq to generate a base-resolution map of 5hmC in tissue-derived cortical excitatory neurons, we found that 5hmC was almost entirely confined to CG dinucleotides. The whole-genome map permitted cytosine, 5-methylcytosine (5mC) and 5hmC to be parsed and revealed genomic features that diverged from global patterns, including enhancers and imprinting control regions with high and low 5hmC/5mC ratios, respectively. Enzymatic deamination overcomes many challenges posed by bisulfite-based methods, thus expanding the scope of epigenome profiling to include scarce samples and opening new lines of inquiry regarding the role of cytosine modifications in genome biology.

The allosteric inhibitors of integrase (termed ALLINIs) interfere with HIV replication by binding to the viral-encoded integrase (IN) protein. Surprisingly, ALLINIs interfere not with DNA integration but with viral particle assembly late during HIV replication. To investigate the ALLINI inhibitory mechanism, we crystallized full-length HIV-1 IN bound to the ALLINI GSK1264 and determined the structure of the complex at 4.4 Å resolution. The structure shows GSK1264 buried between the IN C-terminal domain (CTD) and the catalytic core domain. In the crystal lattice, the interacting domains are contributed by two different dimers so that IN forms an open polymer mediated by inhibitor-bridged contacts; the N-terminal domains do not participate and are structurally disordered. Engineered amino acid substitutions at the inhibitor interface blocked ALLINI-induced multimerization. HIV escape mutants with reduced sensitivity to ALLINIs commonly altered amino acids at or near the inhibitor-bound interface, and these substitutions also diminished IN multimerization. We propose that ALLINIs inhibit particle assembly by stimulating inappropriate polymerization of IN via interactions between the catalytic core domain and the CTD and that understanding the interface involved offers new routes to inhibitor optimization.

Transposons are useful tools for creating transgenic organisms, insertional mutagenesis, and genome engineering. TcBuster, a novel hAT-family transposon system derived from the red flour beetle Tribolium castaneum, was shown to be highly active in previous studies in insect embryoes.

Genome-wide siRNA screens have identified host cell factors important for efficient HIV infection, among which are nuclear pore proteins such as RanBP2/Nup358 and the karyopherin Transportin-3/TNPO3. Analysis of the roles of these proteins in the HIV replication cycle suggested that correct trafficking through the pore may facilitate the subsequent integration step. Here we present data for coupling between these steps by demonstrating that depletion of Transportin-3 or RanBP2 altered the terminal step in early HIV replication, the selection of chromosomal sites for integration. We found that depletion of Transportin-3 and RanBP2 altered integration targeting for HIV. These knockdowns reduced HIV integration frequency in gene-dense regions and near gene-associated features, a pattern that differed from that reported for depletion of the HIV integrase binding cofactor Psip1/Ledgf/p75. MLV integration was not affected by the Transportin-3 knockdown. Using siRNA knockdowns and integration targeting analysis, we also implicated several additional nuclear proteins in proper target site selection. To map viral determinants of integration targeting, we analyzed a chimeric HIV derivative containing MLV gag, and found that the gag replacement phenocopied the Transportin-3 and RanBP2 knockdowns. Thus, our data support a model in which Gag-dependent engagement of the proper transport and nuclear pore machinery mediate trafficking of HIV complexes to sites of integration.

Antibiotic use in humans has been associated with outgrowth of fungi. Here we used a murine model to investigate the gut microbiome over 76 days of treatment with vancomycin, ampicillin, neomycin, and metronidazole and subsequent recovery. Mouse stool was studied as a surrogate for the microbiota of the lower gastrointestinal tract. The abundance of fungi and bacteria was measured using quantitative PCR, and the proportional composition of the communities quantified using 454/Roche pyrosequencing of rRNA gene tags. Prior to treatment, bacteria outnumbered fungi by >3 orders of magnitude. Upon antibiotic treatment, bacteria dropped in abundance >3 orders of magnitude, so that the predominant 16S sequences detected became transients derived from food. Upon cessation of treatment, bacterial communities mostly returned to their previous numbers and types after 8 weeks, though communities remained detectably different from untreated controls. Fungal communities varied substantially over time, even in the untreated controls. Separate cages within the same treatment group showed radical differences, but mice within a cage generally behaved similarly. Fungi increased ∼40-fold in abundance upon antibiotic treatment but declined back to their original abundance after cessation of treatment. At the last time point, Candida remained more abundant than prior to treatment. These data show that 1) gut fungal populations change radically during normal mouse husbandry, 2) fungi grow out in the gut upon suppression of bacterial communities with antibiotics, and 3) perturbations due to antibiotics persist long term in both the fungal and bacterial microbiota.

Commensal bacteria that colonize mammalian barrier surfaces are reported to influence T helper type 2 (T(H)2) cytokine-dependent inflammation and susceptibility to allergic disease, although the mechanisms that underlie these observations are poorly understood. In this report, we find that deliberate alteration of commensal bacterial populations via oral antibiotic treatment resulted in elevated serum IgE concentrations, increased steady-state circulating basophil populations and exaggerated basophil-mediated T(H)2 cell responses and allergic inflammation. Elevated serum IgE levels correlated with increased circulating basophil populations in mice and subjects with hyperimmunoglobulinemia E syndrome. Furthermore, B cell-intrinsic expression of myeloid differentiation factor 88 (MyD88) was required to limit serum IgE concentrations and circulating basophil populations in mice. Commensal-derived signals were found to influence basophil development by limiting proliferation of bone marrow-resident precursor populations. Collectively, these results identify a previously unrecognized pathway through which commensal-derived signals influence basophil hematopoiesis and susceptibility to T(H)2 cytokine-dependent inflammation and allergic disease.

Intense interest centers on the role of the human gut microbiome in health and disease, but optimal methods for analysis are still under development. Here we present a study of methods for surveying bacterial communities in human feces using 454/Roche pyrosequencing of 16S rRNA gene tags. We analyzed fecal samples from 10 individuals and compared methods for storage, DNA purification and sequence acquisition. To assess reproducibility, we compared samples one cm apart on a single stool specimen for each individual. To analyze storage methods, we compared 1) immediate freezing at -80 degrees C, 2) storage on ice for 24 or 3) 48 hours. For DNA purification methods, we tested three commercial kits and bead beating in hot phenol. Variations due to the different methodologies were compared to variation among individuals using two approaches--one based on presence-absence information for bacterial taxa (unweighted UniFrac) and the other taking into account their relative abundance (weighted UniFrac). In the unweighted analysis relatively little variation was associated with the different analytical procedures, and variation between individuals predominated. In the weighted analysis considerable variation was associated with the purification methods. Particularly notable was improved recovery of Firmicutes sequences using the hot phenol method. We also carried out surveys of the effects of different 454 sequencing methods (FLX versus Titanium) and amplification of different 16S rRNA variable gene segments. Based on our findings we present recommendations for protocols to collect, process and sequence bacterial 16S rDNA from fecal samples--some major points are 1) if feasible, bead-beating in hot phenol or use of the PSP kit improves recovery; 2) storage methods can be adjusted based on experimental convenience; 3) unweighted (presence-absence) comparisons are less affected by lysis method.

The completion of the human genome sequence has made possible genome-wide studies of retroviral DNA integration. Here we report an analysis of 3,127 integration site sequences from human cells. We compared retroviral vectors derived from human immunodeficiency virus (HIV), avian sarcoma-leukosis virus (ASLV), and murine leukemia virus (MLV). Effects of gene activity on integration targeting were assessed by transcriptional profiling of infected cells. Integration by HIV vectors, analyzed in two primary cell types and several cell lines, strongly favored active genes. An analysis of the effects of tissue-specific transcription showed that it resulted in tissue-specific integration targeting by HIV, though the effect was quantitatively modest. Chromosomal regions rich in expressed genes were favored for HIV integration, but these regions were found to be interleaved with unfavorable regions at CpG islands. MLV vectors showed a strong bias in favor of integration near transcription start sites, as reported previously. ASLV vectors showed only a weak preference for active genes and no preference for transcription start regions. Thus, each of the three retroviruses studied showed unique integration site preferences, suggesting that virus-specific binding of integration complexes to chromatin features likely guides site selection.

Genomic parasites have evolved distinctive lifestyles to optimize replication in the context of the genomes they inhabit. Here, we introduced new DNA into eukaryotic cells using bacteriophage Mu DNA transposition complexes, termed 'transpososomes'. Following electroporation of transpososomes and selection for marker gene expression, efficient integration was verified in yeast, mouse and human genomes. Although Mu has evolved in prokaryotes, strong biases were seen in the target site distributions in eukaryotic genomes, and these biases differed between yeast and mammals. In Saccharomyces cerevisiae transposons accumulated outside of genes, consistent with selection against gene disruption. In mouse and human cells, transposons accumulated within genes, which previous work suggests is a favorable location for efficient expression of selectable markers. Naturally occurring transposons and viruses in yeast and mammals show related, but more extreme, targeting biases, suggesting that they are responding to the same pressures. These data help clarify the constraints exerted by genome structure on genomic parasites, and illustrate the wide utility of the Mu transpososome technology for gene transfer in eukaryotic cells.

The vertebrate gut harbors a vast community of bacterial mutualists, the composition of which is modulated by the host immune system. Many gastrointestinal (GI) diseases are expected to be associated with disruptions of host-bacterial interactions, but relatively few comprehensive studies have been reported. We have used the rhesus macaque model to investigate forces shaping GI bacterial communities. We used DNA bar coding and pyrosequencing to characterize 141,000 sequences of 16S rRNA genes obtained from 100 uncultured GI bacterial samples, allowing quantitative analysis of community composition in health and disease. Microbial communities of macaques were distinct from those of mice and humans in both abundance and types of taxa present. The macaque communities differed among samples from intestinal mucosa, colonic contents, and stool, paralleling studies of humans. Communities also differed among animals, over time within individual animals, and between males and females. To investigate changes associated with disease, samples of colonic contents taken at necropsy were compared between healthy animals and animals with colitis and undergoing antibiotic therapy. Communities from diseased and healthy animals also differed significantly in composition. This work provides comprehensive data and improved methods for studying the role of commensal microbiota in macaque models of GI diseases and provides a model for the large-scale screening of the human gut microbiome.

Gene transfer has been used to correct inherited immunodeficiencies, but in several patients integration of therapeutic retroviral vectors activated proto-oncogenes and caused leukemia. Here, we describe improved methods for characterizing integration site populations from gene transfer studies using DNA bar coding and pyrosequencing. We characterized 160,232 integration site sequences in 28 tissue samples from eight mice, where Rag1 or Artemis deficiencies were corrected by introducing the missing gene with gamma-retroviral or lentiviral vectors. The integration sites were characterized for their genomic distributions, including proximity to proto-oncogenes. Several mice harbored abnormal lymphoproliferations following therapy--in these cases, comparison of the location and frequency of isolation of integration sites across multiple tissues helped clarify the contribution of specific proviruses to the adverse events. We also took advantage of the large number of pyrosequencing reads to show that recovery of integration sites can be highly biased by the use of restriction enzyme cleavage of genomic DNA, which is a limitation in all widely used methods, but describe improved approaches that take advantage of the power of pyrosequencing to overcome this problem. The methods described here should allow integration site populations from human gene therapy to be deeply characterized with spatial and temporal resolution.

Unequivocal demonstration of the therapeutic utility of γ-retroviral vectors for gene therapy applications targeting the hematopoietic system was accompanied by instances of insertional mutagenesis. These events stimulated the ongoing development of putatively safer integrating vector systems and analysis methods to characterize and compare integration site (IS) biosafety profiles. Continuing advances in next-generation sequencing technologies are driving the generation of ever-more complex IS datasets. Available bioinformatic tools to compare such datasets focus on the association of integration sites (ISs) with selected genomic and epigenetic features, and the choice of these features determines the ability to discriminate between datasets. We describe the scalable application of point-process coherence analysis (CA) to compare patterns produced by vector ISs across genomic intervals, uncoupled from association with genomic features. To explore the utility of CA in the context of an unresolved question, we asked whether the differing transduction conditions used in the initial Paris and London SCID-X1 gene therapy trials result in divergent genome-wide integration profiles. We tested a transduction carried out under each condition, and showed that CA could indeed resolve differences in IS distributions. Existence of these differences was confirmed by the application of established methods to compare integration datasets.

Endobronchial stents are increasingly used to treat airway complications in multiple conditions including lung transplantation but little is known about the biofilms that form on these devices. We applied deep sequencing to profile luminal biofilms of 46 endobronchial stents removed from 20 subjects primarily with lung transplantation-associated airway compromise. Microbial communities were analyzed by bacterial 16S rRNA and fungal ITS marker gene sequencing. Corynebacterium was the most common bacterial taxa across biofilm communities. Clustering analysis revealed three bacterial biofilm types: one low diversity and dominated by Corynebacterium; another was polymicrobial and characterized by Staphylococcus; and the third was polymicrobial and associated with Pseudomonas, Streptococcus, and Prevotella. Biofilm type was significantly correlated with stent material: covered metal with the Staphylococcus-type biofilm, silicone with the Corynebacterium-dominated biofilm, and uncovered metal with the polymicrobial biofilm. Subjects with sequential stents had frequent transitions between community types. Fungal analysis found Candida was most prevalent, Aspergillus was common and highly enriched in two of three stents associated with airway anastomotic dehiscence, and fungal taxa not typically considered pathogens were highly enriched in some stents. Thus, molecular analysis revealed a complex and dynamic endobronchial stent biofilm with three bacterial types that associate with stent material, a central role for Corynebacterium, and that both expected and unexpected fungi inhabit this unique niche. The current work provides a foundation for studies to investigate the relationship between stent biofilm composition and clinical outcomes, mechanisms of biofilm establishment, and strategies for improved stent technology and use in airway compromise.