Gary P. Wang

Transfusion independence and HMGA2 activation after gene therapy of human β-thalassaemia

The β-haemoglobinopathies are the most prevalent inherited disorders worldwide. Gene therapy of β-thalassaemia is particularly challenging given the requirement for massive haemoglobin production in a lineage-specific manner and the lack of selective advantage for corrected haematopoietic stem cells. Compound βE/β0-thalassaemia is the most common form of severe thalassaemia in southeast Asian countries and their diasporas. The βE-globin allele bears a point mutation that causes alternative splicing. The abnormally spliced form is non-coding, whereas the correctly spliced messenger RNA expresses a mutated βE-globin with partial instability. When this is compounded with a non-functional β0 allele, a profound decrease in β-globin synthesis results, and approximately half of βE/β0-thalassaemia patients are transfusion-dependent. The only available curative therapy is allogeneic haematopoietic stem cell transplantation, although most patients do not have a human-leukocyte-antigen-matched, geno-identical donor, and those who do still risk rejection or graft-versus-host disease. Here we show that, 33 months after lentiviral β-globin gene transfer, an adult patient with severe βE/β0-thalassaemia dependent on monthly transfusions since early childhood has become transfusion independent for the past 21 months. Blood haemoglobin is maintained between 9 and 10 g dl−1, of which one-third contains vector-encoded β-globin. Most of the therapeutic benefit results from a dominant, myeloid-biased cell clone, in which the integrated vector causes transcriptional activation of HMGA2 in erythroid cells with further increased expression of a truncated HMGA2 mRNA insensitive to degradation by let-7 microRNAs. The clonal dominance that accompanies therapeutic efficacy may be coincidental and stochastic or result from a hitherto benign cell expansion caused by dysregulation of the HMGA2 gene in stem/progenitor cells.

DNA bar coding and pyrosequencing to identify rare HIV drug resistance mutations

Treatment of HIV-infected individuals with antiretroviral agents selects for drug-resistant mutants, resulting in frequent treatment failures. Although the major antiretroviral resistance mutations are routinely characterized by DNA sequencing, treatment failures are still common, probably in part because undetected rare resistance mutations facilitate viral escape. Here we combined DNA bar coding and massively parallel pyrosequencing to quantify rare drug resistance mutations. Using DNA bar coding, we were able to analyze seven viral populations in parallel, overall characterizing 118 093 sequence reads of average length 103 bp. Analysis of a control HIV mixture showed that resistance mutations present as 5% of the population could be readily detected without false positive calls. In three samples of multidrug-resistant HIV populations from patients, all the drug-resistant mutations called by conventional analysis were identified, as well as four additional low abundance drug resistance mutations, some of which would be expected to influence the response to antiretroviral therapy. Methods for sensitive characterization of HIV resistance alleles have been reported, but only the pyrosequencing method allows all the positions at risk for drug resistance mutations to be interrogated deeply for many HIV populations in a single experiment.

Intestinal Dysbiosis and Depletion of Butyrogenic Bacteria in Clostridium difficile Infection and Nosocomial Diarrhea

ABSTRACT Clostridium difficile infection (CDI) causes nearly half a million cases of diarrhea and colitis in the United States each year. Although the importance of the gut microbiota in C. difficile pathogenesis is well recognized, components of the human gut flora critical for colonization resistance are not known. Culture-independent high-density Roche 454 pyrosequencing was used to survey the distal gut microbiota for 39 individuals with CDI, 36 subjects with C. difficile-negative nosocomial diarrhea (CDN), and 40 healthy control subjects. A total of 526,071 partial 16S rRNA sequence reads of the V1 to V3 regions were aligned with 16S databases, identifying 3,531 bacterial phylotypes from 115 fecal samples. Genomic analysis revealed significant alterations of organism lineages in both the CDI and CDN groups, which were accompanied by marked decreases in microbial diversity and species richness driven primarily by a paucity of phylotypes within the Firmicutes phylum. Normally abundant gut commensal organisms, including the Ruminococcaceae and Lachnospiraceae families and butyrate-producing C2 to C4 anaerobic fermenters, were significantly depleted in the CDI and CDN groups. These data demonstrate associations between the depletion of Ruminococcaceae, Lachnospiraceae, and butyrogenic bacteria in the gut microbiota and nosocomial diarrhea, including C. difficile infection. Mechanistic studies focusing on the functional roles of these organisms in diarrheal diseases and resistance against C. difficile colonization are warranted.

Hepatitis C Virus Transmission Bottlenecks Analyzed by Deep Sequencing

ABSTRACT Hepatitis C virus (HCV) replication in infected patients produces large and diverse viral populations, which give rise to drug-resistant and immune escape variants. Here, we analyzed HCV populations during transmission and diversification in longitudinal and cross-sectional samples using 454/Roche pyrosequencing, in total analyzing 174,185 sequence reads. To sample diversity, four locations in the HCV genome were analyzed, ranging from high diversity (the envelope hypervariable region 1 [HVR1]) to almost no diversity (the 5′ untranslated region [UTR]). For three longitudinal samples for which early time points were available, we found that only 1 to 4 viral variants were present, suggesting that productive infection was initiated by a very small number of HCV particles. Sequence diversity accumulated subsequently, with the 5′ UTR showing almost no diversification while the envelope HVR1 showed >100 variants in some subjects. Calculation of the transmission probability for only a single variant, taking into account the measured population structure within patients, confirmed initial infection by one or a few viral particles. These findings provide the most detailed sequence-based analysis of HCV transmission bottlenecks to date. The analytical methods described here are broadly applicable to studies of viral diversity using deep sequencing.

Dynamics of gene-modified progenitor cells analyzed by tracking retroviral integration sites in a human SCID-X1 gene therapy trial

X-linked severe-combined immunodeficiency (SCID-X1) has been treated by therapeutic gene transfer using gammaretroviral vectors, but insertional activation of proto-oncogenes contributed to leukemia in some patients. Here we report a longitudinal study of gene-corrected progenitor cell populations from 8 patients using 454 pyrosequencing to map vector integration sites, and extensive resampling to allow quantification of clonal abundance. The number of transduced cells infused into patients initially predicted the subsequent diversity of circulating cells. A capture-recapture analysis was used to estimate the size of the gene-corrected cell pool, revealing that less than 1/100th of the infused cells had long-term repopulating activity. Integration sites were clustered even at early time points, often near genes involved in growth control, and several patients harbored expanded cell clones with vectors integrated near the cancer-implicated genes CCND2 and HMGA2, but remain healthy. Integration site tracking also documented that chemotherapy for adverse events resulted in successful control. The longitudinal analysis emphasizes that key features of transduced cell populations--including diversity, integration site clustering, and expansion of some clones--were established early after transplantation. The approaches to sequencing and bioinformatics analysis reported here should be widely useful in assessing the outcome of gene therapy trials.

DNA bar coding and pyrosequencing to analyze adverse events in therapeutic gene transfer

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.

Massively parallel pyrosequencing in HIV research

1411–1415Keywords: DNA sequence, homogeneous PCR products, pyrosequencingThe new massively parallel sequencing methods are soastonishing that one wonders whether space aliens aresecretly behind them. One technician, running a singleinstrument,canobtainuptoapproximately1billionbasesof DNA sequence in a few days. Here we describe thenew sequencing methods, briefly present a few appli-cations in HIV research, and then speculate on futuredirections.

A method to sequence and quantify DNA integration for monitoring outcome in gene therapy

Human genetic diseases have been successfully corrected by integration of functional copies of the defective genes into human cells, but in some cases integration of therapeutic vectors has activated proto-oncogenes and contributed to leukemia. For this reason, extensive efforts have focused on analyzing integration site populations from patient samples, but the most commonly used methods for recovering newly integrated DNA suffer from severe recovery biases. Here, we show that a new method based on phage Mu transposition in vitro allows convenient and consistent recovery of integration site sequences in a form that can be analyzed directly using DNA barcoding and pyrosequencing. The method also allows simple estimation of the relative abundance of gene-modified cells from human gene therapy subjects, which has previously been lacking but is crucial for detecting expansion of cell clones that may be a prelude to adverse events.

Dysbiosis and Alterations in Predicted Functions of the Subgingival Microbiome in Chronic Periodontitis

ABSTRACT Chronic periodontitis is an inflammatory disease of the periodontium affecting nearly 65 million adults in the United States. Changes in subgingival microbiota have long been associated with chronic periodontitis. Recent culture-independent molecular studies have revealed the immense richness and complexity of oral microbial communities. However, data sets across studies have not been directly compared, and whether the observed microbial variations are consistent across different studies is not known. Here, we used 16S rRNA sequencing to survey the subgingival microbiota in 25 subjects with chronic periodontal disease and 25 healthy controls and compared our data sets with those of three previously reported microbiome studies. Consistent with data from previous studies, our results demonstrate a significantly altered microbial community structure with decreased heterogeneity in periodontal disease. Comparison with data from three previously reported studies revealed that subgingival microbiota clustered by study. However, differences between periodontal health and disease were larger than the technical variations across studies. Using a prediction score and applying five different distance metrics, we observed two predominant clusters. One cluster was driven by Fusobacterium and Porphyromonas and was associated with clinically apparent periodontitis, and the second cluster was dominated by Rothia and Streptococcus in the majority of healthy sites. The predicted functional capabilities of the periodontitis microbiome were significantly altered. Genes involved in bacterial motility, energy metabolism, and lipopolysaccharide biosynthesis were overrepresented in periodontal disease, whereas genes associated with transporters, the phosphotransferase system, transcription factors, amino acid biosynthesis, and glycolysis/gluconeogenesis were enriched in healthy controls. These results demonstrate significant alterations in microbial composition and function in periodontitis and suggest genes and metabolic pathways associated with periodontal disease.

Methods for integration site distribution analyses in animal cell genomes

The question of where retroviral DNA becomes integrated in chromosomes is important for understanding (i) the mechanisms of viral growth, (ii) devising new anti-retroviral therapy, (iii) understanding how genomes evolve, and (iv) developing safer methods for gene therapy. With the completion of genome sequences for many organisms, it has become possible to study integration targeting by cloning and sequencing large numbers of host-virus DNA junctions, then mapping the host DNA segments back onto the genomic sequence. This allows statistical analysis of the distribution of integration sites relative to the myriad types of genomic features that are also being mapped onto the sequence scaffold. Here we present methods for recovering and analyzing integration site sequences.