Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential.

BACKGROUND: Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions. METHODS: Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia. RESULTS: While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region. CONCLUSIONS: These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity.

Background: Amplified copy number in the plasmepsin II/III genes within Plasmodium falciparum has been associated with decreased sensitivity to piperaquine. To examine this association and test whether additional loci might also contribute, we performed a genome-wide association study of ex vivo P. falciparum susceptibility to piperaquine. Methods: Plasmodium falciparum DNA from 183 samples collected primarily from Cambodia was genotyped at 33716 genome-wide single nucleotide polymorphisms (SNPs). Linear mixed models and random forests were used to estimate associations between parasite genotypes and piperaquine susceptibility. Candidate polymorphisms were evaluated for their association with dihydroartemisinin-piperaquine treatment outcomes in an independent dataset. Results: Single nucleotide polymorphisms on multiple chromosomes were associated with piperaquine 90% inhibitory concentrations (IC90) in a genome-wide analysis. Fine-mapping of genomic regions implicated in genome-wide analyses identified multiple SNPs in linkage disequilibrium with each other that were significantly associated with piperaquine IC90, including a novel mutation within the gene encoding the P. falciparum chloroquine resistance transporter, PfCRT. This mutation (F145I) was associated with dihydroartemisinin-piperaquine treatment failure after adjusting for the presence of amplified plasmepsin II/III, which was also associated with decreased piperaquine sensitivity. Conclusions: Our data suggest that, in addition to plasmepsin II/III copy number, other loci, including pfcrt, may also be involved in piperaquine resistance.

A new method for sequencing the hypervariable Plasmodium falciparum gene var2csa from clinical samples.

BACKGROUND: VAR2CSA, a member of the Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) family, mediates the binding of P. falciparum-infected erythrocytes to chondroitin sulfate A, a surface-associated molecule expressed in placental cells, and plays a central role in the pathogenesis of placental malaria. VAR2CSA is a target of naturally acquired immunity and, as such, is a leading vaccine candidate against placental malaria. This protein is very polymorphic and technically challenging to sequence. Published var2csa sequences, mostly limited to specific domains, have been generated through the sequencing of cloned PCR amplicons using capillary electrophoresis, a method that is both time consuming and costly, and that performs poorly when applied to clinical samples that are commonly polyclonal. A next-generation sequencing platform, Pacific Biosciences (PacBio), offers an alternative approach to overcome these issues. METHODS: PCR primers were designed that target a 5 kb segment in the 5' end of var2csa and the resulting amplicons were sequenced using PacBio sequencing. The primers were optimized using two laboratory strains and were validated on DNA from 43 clinical samples, extracted from dried blood spots on filter paper or from cryopreserved P. falciparum-infected erythrocytes. Sequence reads were assembled using the SMRT-analysis ConsensusTools module. RESULTS: Here, a PacBio sequencing-based approach for recovering a segment encoding the majority of VAR2CSA's extracellular region is described; this segment includes the totality of the first four domains in the 5' end of var2csa (~5 kb), from clinical malaria samples. The feasibility of the method is demonstrated, showing a high success rate from cryopreserved samples and more limited success from dried blood spots stored at room temperature, and characterized the genetic variation of the var2csa locus. CONCLUSIONS: This method will facilitate a detailed analysis of var2csa genetic variation and can be adapted to sequence other hypervariable P. falciparum genes.

New var reconstruction algorithm exposes high var sequence diversity in a single geographic location in Mali.

BACKGROUND: Encoded by the var gene family, highly variable Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) proteins mediate tissue-specific cytoadherence of infected erythrocytes, resulting in immune evasion and severe malaria disease. Sequencing and assembling the 40-60 var gene complement for individual infections has been notoriously difficult, impeding molecular epidemiological studies and the assessment of particular var elements as subunit vaccine candidates. METHODS: We developed and validated a novel algorithm, Exon-Targeted Hybrid Assembly (ETHA), to perform targeted assembly of var gene sequences, based on a combination of Pacific Biosciences and Illumina data. RESULTS: Using ETHA, we characterized the repertoire of var genes in 12 samples from uncomplicated malaria infections in children from a single Malian village and showed them to be as genetically diverse as vars from isolates from around the globe. The gene var2csa, a member of the var family associated with placental malaria pathogenesis, was present in each genome, as were vars previously associated with severe malaria. CONCLUSION: ETHA, a tool to discover novel var sequences from clinical samples, will aid the understanding of malaria pathogenesis and inform the design of malaria vaccines based on PfEMP1. ETHA is available at: https://sourceforge.net/projects/etha/ .

The genomes of three stocks comprising the most widely utilized live sporozoite Theileria parva vaccine exhibit very different degrees and patterns of sequence divergence.

BACKGROUND: There are no commercially available vaccines against human protozoan parasitic diseases, despite the success of vaccination-induced long-term protection against infectious diseases. East Coast fever, caused by the protist Theileria parva, kills one million cattle each year in sub-Saharan Africa, and contributes significantly to hunger and poverty in the region. A highly effective, live, multi-isolate vaccine against T. parva exists, but its component isolates have not been characterized. Here we sequence and compare the three component T. parva stocks within this vaccine, the Muguga Cocktail, namely Muguga, Kiambu5 and Serengeti-transformed, aiming to identify genomic features that contribute to vaccine efficacy. RESULTS: We find that Serengeti-transformed, originally isolated from the wildlife carrier, the African Cape buffalo, is remarkably and unexpectedly similar to the Muguga isolate. The 420 detectable non-synonymous SNPs were distributed among only 53 genes, primarily subtelomeric antigens and antigenic families. The Kiambu5 isolate is considerably more divergent, with close to 40,000 SNPs relative to Muguga, including >8,500 non-synonymous mutations distributed among >1,700 (42.5 %) of the predicted genes. These genetic markers of the component stocks can be used to characterize the composition of new batches of the Muguga Cocktail. CONCLUSIONS: Differences among these three isolates, while extensive, represent only a small proportion of the genetic variation in the entire species. Given the efficacy of the Muguga Cocktail in inducing long-lasting protection against infections in the field, our results suggest that whole-organism vaccines against parasitic diseases can be highly efficacious despite considerable genome-wide differences relative to the isolates against which they protect.

Functional dynamics of the gut microbiome in elderly people during probiotic consumption.

UNLABELLED: A mechanistic understanding of the purported health benefits conferred by consumption of probiotic bacteria has been limited by our knowledge of the resident gut microbiota and its interaction with the host. Here, we detail the impact of a single-organism probiotic, Lactobacillus rhamnosus GG ATCC 53103 (LGG), on the structure and functional dynamics (gene expression) of the gut microbiota in a study of 12 healthy individuals, 65 to 80 years old. The analysis revealed that while the overall community composition was stable as assessed by 16S rRNA profiling, the transcriptional response of the gut microbiota was modulated by probiotic treatment. Comparison of transcriptional profiles based on taxonomic composition yielded three distinct transcriptome groups that displayed considerable differences in functional dynamics. The transcriptional profile of LGG in vivo was remarkably concordant across study subjects despite the considerable interindividual nature of the gut microbiota. However, we identified genes involved in flagellar motility, chemotaxis, and adhesion from Bifidobacterium and the dominant butyrate producers Roseburia and Eubacterium whose expression was increased during probiotic consumption, suggesting that LGG may promote interactions between key constituents of the microbiota and the host epithelium. These results provide evidence for the discrete functional effects imparted by a specific single-organism probiotic and challenge the prevailing notion that probiotics substantially modify the resident microbiota within nondiseased individuals in an appreciable fashion. IMPORTANCE: Probiotic bacteria have been used for over a century to promote digestive health. Many individuals report that probiotics alleviate a number of digestive issues, yet little evidence links how probiotic microbes influence human health. Here, we show how the resident microbes that inhabit the healthy human gut respond to a probiotic. The well-studied probiotic Lactobacillus rhamnosus GG ATCC 53103 (LGG) was administered in a clinical trial, and a suite of measurements of the resident microbes were taken to evaluate potential changes over the course of probiotic consumption. We found that LGG transiently enriches for functions to potentially promote anti-inflammatory pathways in the resident microbes.

Transcriptional attenuation controls macrolide inducible efflux and resistance in Streptococcus pneumoniae and in other Gram-positive bacteria containing mef/mel(msr(D)) elements.

Macrolide resistance, emerging in Streptococcus pneumoniae and other Gram-positive bacteria, is increasingly due to efflux pumps encoded by mef/mel(msr) operons found on discrete mobile genetic elements. The regulation of mef/mel(msr) in these elements is not well understood. We identified the mef(E)/mel transcriptional start, localized the mef(E)/mel promoter, and demonstrated attenuation of transcription as a mechanism of regulation of macrolide-inducible mef-mediated macrolide resistance in S. pneumoniae. The mef(E)/mel transcriptional start site was a guanine 327 bp upstream of mef(E). Consensus pneumococcal promoter -10 (5'-TATACT-3') and -35 (5'-TTGAAC-3') boxes separated by 17 bp were identified 7 bp upstream of the start site. Analysis of the predicted secondary structure of the 327 5' region identified four pairs of inverted repeats R1-R8 predicted to fold into stem-loops, a small leader peptide [MTASMRLR, (Mef(E)L)] required for macrolide induction and a Rho-independent transcription terminator. RNA-seq analyses provided confirmation of transcriptional attenuation. In addition, expression of mef(E)L was also influenced by mef(E)L-dependent mRNA stability. The regulatory region 5' of mef(E) was highly conserved in other mef/mel(msr)-containing elements including Tn1207.1 and the 5612IQ complex in pneumococci and Tn1207.3 in Group A streptococci, indicating a regulatory mechanism common to a wide variety of Gram-positive bacteria containing mef/mel(msr) elements.

Simultaneous transcriptional profiling of bacteria and their host cells.

We developed an RNA-Seq-based method to simultaneously capture prokaryotic and eukaryotic expression profiles of cells infected with intracellular bacteria. As proof of principle, this method was applied to Chlamydia trachomatis-infected epithelial cell monolayers in vitro, successfully obtaining transcriptomes of both C. trachomatis and the host cells at 1 and 24 hours post-infection. Chlamydiae are obligate intracellular bacterial pathogens that cause a range of mammalian diseases. In humans chlamydiae are responsible for the most common sexually transmitted bacterial infections and trachoma (infectious blindness). Disease arises by adverse host inflammatory reactions that induce tissue damage & scarring. However, little is known about the mechanisms underlying these outcomes. Chlamydia are genetically intractable as replication outside of the host cell is not yet possible and there are no practical tools for routine genetic manipulation, making genome-scale approaches critical. The early timeframe of infection is poorly understood and the host transcriptional response to chlamydial infection is not well defined. Our simultaneous RNA-Seq method was applied to a simplified in vitro model of chlamydial infection. We discovered a possible chlamydial strategy for early iron acquisition, putative immune dampening effects of chlamydial infection on the host cell, and present a hypothesis for Chlamydia-induced fibrotic scarring through runaway positive feedback loops. In general, simultaneous RNA-Seq helps to reveal the complex interplay between invading bacterial pathogens and their host mammalian cells and is immediately applicable to any bacteria/host cell interaction.

Impact of oral typhoid vaccination on the human gut microbiota and correlations with s. Typhi-specific immunological responses.

The resident microbial consortia of the human gastrointestinal tract play an integral role in modulating immune responses both locally and systemically. However, detailed information regarding the effector immune responses after vaccine administration in relation to the gastrointestinal microbiota is absent. In this study, the licensed oral live-attenuated typhoid vaccine Ty21a was administered in a clinical study to investigate whether oral immunization resulted in alterations of the microbiota and to identify whether a given microbiota composition, or subsets of the community, are associated with defined S. Typhi-specific immunological responses. The fecal microbiota composition and temporal dynamics were characterized using bacterial 16S rRNA pyrosequencing from individuals who were either immunized with the Ty21a typhoid vaccine (n = 13) or served as unvaccinated controls (n = 4). The analysis revealed considerable inter- and intra-individual variability, yet no discernible perturbations of the bacterial assemblage related to vaccine administration were observed. S. Typhi-specific cell mediated immune (CMI) responses were evaluated by measurement of intracellular cytokine production using multiparametric flow cytometry, and humoral responses were evaluated by measurement of serum anti-LPS IgA and IgG titers. Volunteers were categorized according to the kinetics and magnitude of their responses. While differences in microbial composition, diversity, or temporal stability were not observed among individuals able to mount a positive humoral response, individuals displaying multiphasic CMI responses harbored more diverse, complex communities. In line with this preliminary observation, over two hundred operational taxonomic units (OTUs) were found to differentiate multiphasic and late CMI responders, the vast majority of which classified within the order Clostridiales. These results provide an unprecedented view into the dramatic temporal heterogeneity of both the gut microbiota and host immune responses.

Analysis of the gut microbiota in the old order Amish and its relation to the metabolic syndrome.

Obesity has been linked to the human gut microbiota; however, the contribution of gut bacterial species to the obese phenotype remains controversial because of conflicting results from studies in different populations. To explore the possible dysbiosis of gut microbiota in obesity and its metabolic complications, we studied men and women over a range of body mass indices from the Old Order Amish sect, a culturally homogeneous Caucasian population of Central European ancestry. We characterized the gut microbiota in 310 subjects by deep pyrosequencing of bar-coded PCR amplicons from the V1-V3 region of the 16S rRNA gene. Three communities of interacting bacteria were identified in the gut microbiota, analogous to previously identified gut enterotypes. Neither BMI nor any metabolic syndrome trait was associated with a particular gut community. Network analysis identified twenty-two bacterial species and four OTUs that were either positively or inversely correlated with metabolic syndrome traits, suggesting that certain members of the gut microbiota may play a role in these metabolic derangements.

Genome sequences of the zoonotic pathogens Chlamydia psittaci 6BC and Cal10.

Chlamydia psittaci is a highly prevalent avian pathogen and the cause of a potentially lethal zoonosis, causing life-threatening pneumonia in humans. We report the genome sequences of C. psittaci 6BC, the prototype strain of the species, and C. psittaci Cal10, a widely used laboratory strain.