Joshua P. Earl

Generation of Genic Diversity among Streptococcus pneumoniae Strains via Horizontal Gene Transfer during a Chronic Polyclonal Pediatric Infection

Although there is tremendous interest in understanding the evolutionary roles of horizontal gene transfer (HGT) processes that occur during chronic polyclonal infections, to date there have been few studies that directly address this topic. We have characterized multiple HGT events that most likely occurred during polyclonal infection among nasopharyngeal strains of Streptococcus pneumoniae recovered from a child suffering from chronic upper respiratory and middle-ear infections. Whole genome sequencing and comparative genomics were performed on six isolates collected during symptomatic episodes over a period of seven months. From these comparisons we determined that five of the isolates were genetically highly similar and likely represented a dominant lineage. We analyzed all genic and allelic differences among all six isolates and found that all differences tended to occur within contiguous genomic blocks, suggestive of strain evolution by homologous recombination. From these analyses we identified three strains (two of which were recovered on two different occasions) that appear to have been derived sequentially, one from the next, each by multiple recombination events. We also identified a fourth strain that contains many of the genomic segments that differentiate the three highly related strains from one another, and have hypothesized that this fourth strain may have served as a donor multiple times in the evolution of the dominant strain line. The variations among the parent, daughter, and grand-daughter recombinant strains collectively cover greater than seven percent of the genome and are grouped into 23 chromosomal clusters. While capturing in vivo HGT, these data support the distributed genome hypothesis and suggest that a single competence event in pneumococci can result in the replacement of DNA at multiple non-adjacent loci.

Comparative Genomic Analyses of 17 Clinical Isolates of Gardnerella vaginalis Provide Evidence of Multiple Genetically Isolated Clades Consistent with Subspeciation into Genovars

Gardnerella vaginalis is associated with a spectrum of clinical conditions, suggesting high degrees of genetic heterogeneity among stains. Seventeen G. vaginalis isolates were subjected to a battery of comparative genomic analyses to determine their level of relatedness. For each measure, the degree of difference among the G. vaginalis strains was the highest observed among 23 pathogenic bacterial species for which at least eight genomes are available. Genome sizes ranged from 1.491 to 1.716 Mb; GC contents ranged from 41.18% to 43.40%; and the core genome, consisting of only 746 genes, makes up only 51.6% of each strains genome on average and accounts for only 27% of the species supragenome. Neighbor-grouping analyses, using both distributed gene possession data and core gene allelic data, each identified two major sets of strains, each of which is composed of two groups. Each of the four groups has its own characteristic genome size, GC ratio, and greatly expanded core gene content, making the genomic diversity of each group within the range for other bacterial species. To test whether these 4 groups corresponded to genetically isolated clades, we inferred the phylogeny of each distributed gene that was present in at least two strains and absent in at least two strains; this analysis identified frequent homologous recombination within groups but not between groups or sets. G. vaginalis appears to include four nonrecombining groups/clades of organisms with distinct gene pools and genomic properties, which may confer distinct ecological properties. Consequently, it may be appropriate to treat these four groups as separate species.

Differences in genotype and virulence among four multidrug-resistant Streptococcus pneumoniae isolates belonging to the PMEN1 clone.

We report on the comparative genomics and characterization of the virulence phenotypes of four S. pneumoniae strains that belong to the multidrug resistant clone PMEN1 (Spain23F ST81). Strains SV35-T23 and SV36-T3 were recovered in 1996 from the nasopharynx of patients at an AIDS hospice in New York. Strain SV36-T3 expressed capsule type 3 which is unusual for this clone and represents the product of an in vivo capsular switch event. A third PMEN1 isolate – PN4595-T23 – was recovered in 1996 from the nasopharynx of a child attending day care in Portugal, and a fourth strain – ATCC700669 – was originally isolated from a patient with pneumococcal disease in Spain in 1984. We compared the genomes among four PMEN1 strains and 47 previously sequenced pneumococcal isolates for gene possession differences and allelic variations within core genes. In contrast to the 47 strains – representing a variety of clonal types – the four PMEN1 strains grouped closely together, demonstrating high genomic conservation within this lineage relative to the rest of the species. In the four PMEN1 strains allelic and gene possession differences were clustered into 18 genomic regions including the capsule, the blp bacteriocins, erythromycin resistance, the MM1-2008 prophage and multiple cell wall anchored proteins. In spite of their genomic similarity, the high resolution chinchilla model was able to detect variations in virulence properties of the PMEN1 strains highlighting how small genic or allelic variation can lead to significant changes in pathogenicity and making this set of strains ideal for the identification of novel virulence determinants.

Design and validation of a supragenome array for determination of the genomic content of Haemophilus influenzae isolates

BackgroundHaemophilus influenzae colonizes the human nasopharynx as a commensal, and is etiologically associated with numerous opportunistic infections of the airway; it is also less commonly associated with invasive disease. Clinical isolates of H. influenzae display extensive genomic diversity and plasticity. The development of strategies to successfully prevent, diagnose and treat H. influenzae infections depends on tools to ascertain the gene content of individual isolates.ResultsWe describe and validate a Haemophilus influenzae supragenome hybridization (SGH) array that can be used to characterize the full genic complement of any strain within the species, as well as strains from several highly related species. The array contains 31,307 probes that collectively cover essentially all alleles of the 2890 gene clusters identified from the whole genome sequencing of 24 clinical H. influenzae strains. The finite supragenome model predicts that these data include greater than 85% of all non-rare genes (where rare genes are defined as those present in less than 10% of sequenced strains). The veracity of the array was tested by comparing the whole genome sequences of eight strains with their hybridization data obtained using the supragenome array. The array predictions were correct and reproducible for ~ 98% of the gene content of all of the sequenced strains. This technology was then applied to an investigation of the gene content of 193 geographically and clinically diverse H. influenzae clinical strains. These strains came from multiple locations from five different continents and Papua New Guinea and include isolates from: the middle ears of persons with otitis media and otorrhea; lung aspirates and sputum samples from pneumonia and COPD patients, blood specimens from patients with sepsis; cerebrospinal fluid from patients with meningitis, as well as from pharyngeal specimens from healthy persons.ConclusionsThese analyses provided the most comprehensive and detailed genomic/phylogenetic look at this species to date, and identified a subset of highly divergent strains that form a separate lineage within the species. This array provides a cost-effective and high-throughput tool to determine the gene content of any H. influenzae isolate or lineage. Furthermore, the method for probe selection can be applied to any species, given a group of available whole genome sequences.

Virulence Potential and Genome-Wide Characterization of Drug Resistant Streptococcus pneumoniae Clones Selected In Vivo by the 7-Valent Pneumococcal Conjugate Vaccine

We used mouse models of pneumococcal colonization and disease combined with full genome sequencing to characterize three major drug resistant clones of S. pneumoniae that were recovered from the nasopharynx of PCV7-immunized children in Portugal. The three clones – serotype 6A (ST2191), serotype 15A (ST63) and serotype 19A (ST276) carried some of the same drug resistance determinants already identified in nasopharyngeal isolates from the pre-PCV7 era. The three clones were able to colonize efficiently the mouse nasopharyngeal mucosa where populations of these pneumococci were retained for as long as 21 days. During this period, the three clones were able to asymptomatically invade the olfactory bulbs, brain, lungs and the middle ear mucosa and established populations in these tissues. The virulence potential of the three clones was poor even at high inoculum (105 CFU per mouse) concentrations in the mouse septicemia model and was undetectable in the pneumonia model. Capsular type 3 transformants of clones 6A and 19A prepared in the laboratory produced lethal infection at low cell concentration (103 CFU per mouse) but the same transformants became impaired in their potential to colonize, indicating the importance of the capsular polysaccharide in both disease and colonization. The three clones were compared to the genomes of 56 S. pneumoniae strains for which sequence information was available in the public databank. Clone 15A (ST63) only differed from the serotype 19F clone G54 in a very few genes including serotype so that this clone may be considered the product of a capsular switch. While no strain with comparable degree of similarity to clone 19A (ST276) was found among the sequenced isolates, by MLST this clone is a single locust variant (SLV) of Denmark14-ST230 international clone. Clone 6A (ST2191) was most similar to the penicillin resistant Hungarian serotype 19A clone.

HIV-1 Genetic Variation Resulting in the Development of New Quasispecies Continues to Be Encountered in the Peripheral Blood of Well- Suppressed Patients.

As a result of antiretroviral therapeutic strategies, human immunodeficiency virus type 1 (HIV-1) infection has become a long-term clinically manageable chronic disease for many infected individuals. However, despite this progress in therapeutic control, including undetectable viral loads and CD4+ T-cell counts in the normal range, viral mutations continue to accumulate in the peripheral blood compartment over time, indicating either low level reactivation and/or replication. Using patients from the Drexel Medicine CNS AIDS Research and Eradication Study (CARES) Cohort, whom have been sampled longitudinally for more than 7 years, genetic change was modeled against to the dominant integrated proviral quasispecies with respect to selection pressures such as therapeutic interventions, AIDS defining illnesses, and other factors. Phylogenetic methods based on the sequences of the LTR and tat exon 1 of the HIV-1 proviral DNA quasispecies were used to obtain an estimate of an average mutation rate of 5.3 nucleotides (nt)/kilobasepair (kb)/year (yr) prior to initiation of antiretroviral therapy (ART). Following ART the baseline mutation rate was reduced to an average of 1.02 nt/kb/yr. The post-ART baseline rate of genetic change, however, appears to be unique for each patient. These studies represent our initial steps in quantifying rates of genetic change among HIV-1 quasispecies using longitudinally sampled sequences from patients at different stages of disease both before and after initiation of combination ART. Notably, while long-term ART reduced the estimated mutation rates in the vast majority of patients studied, there was still measurable HIV-1 mutation even in patients with no detectable virus by standard quantitative assays. Determining the factors that affect HIV-1 mutation rates in the peripheral blood may lead to elucidation of the mechanisms associated with changes in HIV-1 disease severity.

Development and Validation of an Haemophilus influenzae Supragenome Hybridization (SGH) Array for Transcriptomic Analyses

We previously carried out the design and testing of a custom-built Haemophilus influenzae supragenome hybridization (SGH) array that contains probe sequences to 2,890 gene clusters identified by whole genome sequencing of 24 strains of H. influenzae. The array was originally designed as a tool to interrogate the gene content of large numbers of clinical isolates without the need for sequencing, however, the data obtained is quantitative and is thus suitable for transcriptomic analyses. In the current study RNA was extracted from H. influenzae strain CZ4126/02 (which was not included in the design of the array) converted to cDNA, and labelled and hybridized to the SGH arrays to assess the quality and reproducibility of data obtained from these custom-designed chips to serve as a tool for transcriptomics. Three types of experimental replicates were analyzed with all showing very high degrees of correlation, thus validating both the array and the methods used for RNA profiling. A custom filtering pipeline for two-condition unpaired data using five metrics was developed to minimize variability within replicates and to maximize the identification of the most significant true transcriptional differences between two samples. These methods can be extended to transcriptional analysis of other bacterial species utilizing supragenome-based arrays.

Abstract 4925: Microbiome-TP53 gene interaction in human lung cancer

Lung cancer is the leading cancer diagnosis worldwide and the number one cause of cancer deaths. Exposure to cigarette smoke, the primary risk factor in lung cancer, reduces epithelial barrier function and increases susceptibility to infections. Herein, we hypothesized that somatic mutations together with cigarette smoke create a dysbiotic microbiota that is associated with lung carcinogenesis. Using lung tissue from controls (n=33) or cancer cases (n=143), we conducted 16S rRNA gene sequencing (MiSeq), with RNA-seq data from lung cancer cases in The Cancer Genome Atlas (n=1112) serving as the validation cohort. We demonstrate a lower alpha diversity in normal lung as compared to non-tumor adjacent or tumor tissue, indicating a shift in the overall microbial community in lung cancer patients as compared to those without cancer. Lung cancer cases were classified by the relative abundance of two taxa, Variovorax and Streptococcus, with an increase in Variovorax abundance in tumors as compared to non-tumor adjacent paired lung tissue (FDR corrected P=0.013). The species of Variovorax was identified histologically, and also by two additional 16S rRNA strategies (Resphera Insight analysis and PacBio sequencing). A group of taxa were associated with squamous cell carcinoma (SCC), of which Acidovorax were enriched in smokers (P =0.0013). Further, these taxa, including Acidovorax, exhibited higher abundance among the subset of SCC cases with TP53 mutations, an association not seen in adenocarcinomas (AD). Therefore, we observed a microbiome-gene and a microbiome-exposure interaction in SCC lung cancer tissue. Together, these results open the door for future biomarker research that could be used to improve screening and direct mechanistic studies of lung cancer therapy. Citation Format: Leigh Greathouse, James White, Valery Bliskovsky, Ashley Vargas, Eric Polley, Elise Bowman, Mohammed Khan, Ana Robles, Brid Ryan, Amiran Dzutsev, Giorgio Trinchieri, Marbin Pineda, Sven Bilke, Paul Meltzer, Marina Walther-Antonio, Garth Ehrlich, Joshua Mell, Joshua Earl, Sergey Balashov, Archana Bhat, Alex Valm, Clayton Deming, Sean Conlan, Julia Oh, Julie Segre, Curtis Harris. Microbiome-TP53 gene interaction in human lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4925. doi:10.1158/1538-7445.AM2017-4925