Population genetic analysis of the Plasmodium falciparum circumsporozoite protein in two distinct ecological regions in Ghana

Background

Extensive genetic diversity in the Plasmodium falciparum circumsporozoite protein (PfCSP) is a major contributing factor to the moderate efficacy of the RTS,S/AS01 vaccine. The transmission intensity and rates of recombination within and between populations influence the extent of its genetic diversity. Understanding the extent and dynamics of PfCSP genetic diversity in different transmission settings will help to interpret the results of current RTS,S efficacy and Phase IV implementation trials conducted within and between populations in malaria-endemic areas such as Ghana.

Methods

Pfcsp sequences were retrieved from the Illumina-generated paired-end short-read sequences of 101 and 131 malaria samples from children aged 6–59 months presenting with clinical malaria at health facilities in Cape Coast (in the coastal belt) and Navrongo (Guinea savannah region), respectively, in Ghana. The sequences were mapped onto the 3D7 reference strain genome to yield high-quality genome-wide coding sequence data. Following data filtering and quality checks to remove missing data, 220 sequences were retained and analysed for the allele frequency spectrum, genetic diversity both within the host and between populations and signatures of selection. Population genetics tools were used to determine the extent and dynamics of Pfcsp diversity in P. falciparum from the two geographically distinct locations in Ghana.

Results

Pfcsp showed extensive diversity at the two sites, with the higher transmission site, Navrongo, exhibiting higher within-host and population-level diversity. The vaccine strain C-terminal epitope of Pfcsp was found in only 5.9% and 45.7% of the Navrongo and Cape Coast sequences, respectively. Between 1 and 6 amino acid variations were observed in the TH2R and TH3R epitope regions of PfCSP. Tajima’s D was negatively skewed, especially for the population from Cape Coast, given the expected historical population expansion. In contrast, a positive Tajima’s D was observed for the Navrongo P. falciparum population, consistent with balancing selection acting on the immuno-dominant TH2R and TH3R vaccine epitopes.

Conclusion

The low frequencies of the Pfcsp vaccine haplotype in the analysed populations indicate a need for additional molecular and immuno-epidemiological studies with broader temporal and geographic sampling in endemic populations targeted for RTS,S application. These results have implications for the efficacy of the vaccine in Ghana and will inform the choice of alleles to be included in future multivalent or chimeric vaccines.

Within-host heterogeneity and flexibility of mcr-1 transmission in chicken gut

OBJECTIVES

To characterize the colistin-resistant bacterial population in the gut and assess diversity of mcr-1 transmission within a single individual.

METHODS

Large numbers of isolates (>100 colonies/chicken cecum sample) were collected from nine randomly selected mcr-1-positive chickens in China and used for comprehensive microbiological, molecular and comparative genomics analyses.

RESULTS

Of 1273 colonies, 968 were mcr-1 positive (962 Escherichia coli, two Escherichia fergusonii, two Klebsiella pneumoniae and two Klebsiella quasipneumoniae). One to six colistin-resistant species and three to 10 E. coli pulsed-field gel electrophoresis (PFGE) clusters could be identified from each sample. Whole-genome sequencing (WGS) analysis of the representative E. coli strains revealed three to nine sequence types observed in a single chicken host. The mcr-1 genes are located in either chromosomes or plasmids of different types, including IncI2 (n=30), IncHI2 (n=14), IncX4 (n=4), p0111(n=2) and IncHI1(n=1). Strikingly, in single cecum samples, one to five Inc type plasmids harbouring mcr-1 could be identified. Great diversity was also observed for the same IncI2 plasmid within a single chicken host. In addition, up to eight genetic contexts of the mcr-1 gene occurred within a single chicken.

CONCLUSIONS

There is extensive heterogeneity and flexibility of mcr-1 transmission in chicken gut due to bacterial species differences, distant clonal relatedness of isolates, many types and variations of mcr-positive plasmids, and the flexible genetic context of the mcr-1 gene. These compelling findings indicate that the gut is a 'melting pot' for active horizontal transfer of the mcr-1 gene.

Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing

BACKGROUND

Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/μl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/μl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated.

RESULTS

Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/μl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/μl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/μl, but only 61% at low Plasmodium concentrations (< 5 genomes/μl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases.

CONCLUSIONS

Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

Next generation amplicon sequencing improves detection of Blastocystis mixed subtype infections

Blastocystis is a highly prevalent enteric protist parasite of humans and animals. Transmission occurs via the fecal-oral route through ingestion of contaminated food or water. Genetic diversity studies have identified numerous subtypes (STs) within the genus Blastocystis based on polymorphism at the SSU rRNA gene. Although there is evidence of frequent mixed subtype infections, the extent of within-host subtype diversity remains largely unexplored. Accurate assessment of Blastocystis ST diversity is crucial to understand epidemiology and sources of Blastocystis transmission to humans. Here, we report the application of next generation sequencing (NGS) for detection and characterization of Blastocystis subtypes to investigate intra-host Blastocystis diversity. A total of 75 specimens obtained from cattle feces, previously identified as Blastocystis positive, were examined using next generation amplicon sequencing. A fragment of the SSU rRNA gene was amplified using Blastocystis-specific primers and resulting amplicons were used for NGS. Comparison of Sanger and NGS results suggest greater sensitivity using the NGS approach. Using Sanger sequencing, mixed infections were suspected in 18 specimens but only confirmed through cloning in three, while NGS identified 49 mixed infections (16 times more). In addition, NGS revealed greater diversity of subtypes with 14 detected compared to 11 by Sanger. Nine more infections with potentially zoonotic STs were detected by NGS than Sanger. Indeed, subtype 3, the most common subtype found in humans, was found in 37% (28) of specimens tested by NGS but in only four specimens using Sanger. Our findings indicate that mixed Blastocystis infections may be far more common than previously thought due to the limitations of current detection methods. This next generation amplicon sequencing strategy improves detection of mixed subtype infections and low abundance subtypes and represents a valuable resource for future Blastocystis studies to improve our understanding of its epidemiology.

Genomic exploration of sequential clinical isolates reveals a distinctive molecular signature of persistent Staphylococcus aureus bacteraemia

BACKGROUND

Large-scale genomic studies of within-host diversity in Staphylococcus aureus bacteraemia (SAB) are needed to understanding bacterial adaptation underlying persistence and thus refining the role of genomics in management of SAB. However, available comparative genomic studies of sequential SAB isolates have tended to focus on selected cases of unusually prolonged bacteraemia, where secondary antimicrobial resistance has developed.

METHODS

To understand bacterial genetic diversity during SAB more broadly, we applied whole genome sequencing to a large collection of sequential isolates obtained from patients with persistent or relapsing bacteraemia. After excluding genetically unrelated isolates, we performed an in-depth genomic analysis of point mutations and chromosome structural variants arising within individual SAB episodes.

RESULTS

We show that, while adaptation pathways are heterogenous and episode-specific, isolates from persistent bacteraemia have a distinctive molecular signature, characterised by a low mutation frequency and high proportion of non-silent mutations. Analysis of structural genomic variants revealed that these often overlooked genetic events are commonly acquired during SAB. We discovered that IS256 insertion may represent the most effective driver of within-host microevolution in selected lineages, with up to three new insertion events per isolate even in the absence of other mutations. Genetic mechanisms resulting in significant phenotypic changes, such as increases in vancomycin resistance, development of small colony phenotypes, and decreases in cytotoxicity, included mutations in key genes (rpoB, stp, agrA) and an IS256 insertion upstream of the walKR operon.

CONCLUSIONS

This study provides for the first time a large-scale analysis of within-host genomic changes during invasive S. aureus infection and describes specific patterns of adaptation that will be informative for both understanding S. aureus pathoadaptation and utilising genomics for management of complicated S. aureus infections.

Tracking carbapenemase-producing bacteria by molecular typing: Population diversity and sampling pitfall

While typing methods are increasingly refined, the sampling of bacteria to be typed in healthcare-associated infection context retains less attention. Through 2 emblematic cases of in-hospital transmission of extensively drug-resistant bacteria producing carbapenemases, we demonstrate the impact of colony sampling in typing results. Because of intra-population diversity, typing several colonies of same species and resistotype was needed to fully track the transmission among patients. Bacterial population studies could better decipher transmission routes of healthcare-associated bacteria, thereby improving outbreak control.

Multiplicity and molecular epidemiology of Plasmodium vivax and Plasmodium falciparum infections in East Africa

Background

Parasite genetic diversity and multiplicity of infection (MOI) affect clinical outcomes, response to drug treatment and naturally-acquired or vaccine-induced immunity. Traditional methods often underestimate the frequency and diversity of multiclonal infections due to technical sensitivity and specificity. Next-generation sequencing techniques provide a novel opportunity to study complexity of parasite populations and molecular epidemiology.

Methods

Symptomatic and asymptomatic Plasmodium vivax samples were collected from health centres/hospitals and schools, respectively, from 2011 to 2015 in Ethiopia. Similarly, both symptomatic and asymptomatic Plasmodium falciparum samples were collected, respectively, from hospitals and schools in 2005 and 2015 in Kenya. Finger-pricked blood samples were collected and dried on filter paper. Long amplicon (> 400 bp) deep sequencing of merozoite surface protein 1 (msp1) gene was conducted to determine multiplicity and molecular epidemiology of P. vivax and P. falciparum infections. The results were compared with those based on short amplicon (117 bp) deep sequencing.

Results

A total of 139 P. vivax and 222 P. falciparum samples were pyro-sequenced for pvmsp1 and pfmsp1, yielding a total of 21 P. vivax and 99 P. falciparum predominant haplotypes. The average MOI for P. vivax and P. falciparum were 2.16 and 2.68, respectively, which were significantly higher than that of microsatellite markers and short amplicon (117 bp) deep sequencing. Multiclonal infections were detected in 62.2% of the samples for P. vivax and 74.8% of the samples for P. falciparum. Four out of the five subjects with recurrent P. vivax malaria were found to be a relapse 44–65 days after clearance of parasites. No difference was observed in MOI among P. vivax patients of different symptoms, ages and genders. Similar patterns were also observed in P. falciparum except for one study site in Kenyan lowland areas with significantly higher MOI.

Conclusions

The study used a novel method to evaluate Plasmodium MOI and molecular epidemiological patterns by long amplicon ultra-deep sequencing. The complexity of infections were similar among age groups, symptoms, genders, transmission settings (spatial heterogeneity), as well as over years (pre- vs. post-scale-up interventions). This study demonstrated that long amplicon deep sequencing is a useful tool to investigate multiplicity and molecular epidemiology of Plasmodium parasite infections.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2337-y) contains supplementary material, which is available to authorized users.

Host specificity, pathogen exposure, and superinfections impact the distribution of Anaplasma phagocytophilum genotypes in ticks, roe deer, and livestock in a fragmented agricultural landscape

Anaplasma phagocytophilum is a bacterial pathogen mainly transmitted by Ixodes ricinus ticks in Europe. It infects wild mammals, livestock, and, occasionally, humans. Roe deer are considered to be the major reservoir, but the genotypes they carry differ from those that are found in livestock and humans. Here, we investigated whether roe deer were the main source of the A. phagocytophilum genotypes circulating in questing I. ricinus nymphs in a fragmented agricultural landscape in France. First, we assessed pathogen prevalence in 1837 I. ricinus nymphs (sampled along georeferenced transects) and 79 roe deer. Prevalence was dramatically different between ticks and roe deer: 1.9% versus 76%, respectively. Second, using high-throughput amplicon sequencing, we characterized the diversity of the A. phagocytophilum genotypes found in 22 infected ticks and 60 infected roe deer; the aim was to determine the frequency of co-infections. Only 22.7% of infected ticks carried genotypes associated with roe deer. This finding fits with others suggesting that cattle density is the major factor explaining infected tick density. To explore epidemiological scenarios capable of explaining these patterns, we constructed compartmental models that focused on how A. phagocytophilum exposure and infection dynamics affected pathogen prevalence in roe deer. At the exposure levels predicted by the results of this study and the literature, the high prevalence in roe deer was only seen in the model in which superinfections could occur during all infection phases and when the probability of infection post exposure was above 0.43. We then interpreted these results from the perspective of livestock and human health.

Deep sequencing of RSV from an adult challenge study and from naturally infected infants reveals heterogeneous diversification dynamics

As RNA virus mutation occurs during replication within host cells, we hypothesized that viral evolution during acute infections in healthy hosts reflects host immune pressure. We therefore investigated the within-host diversification of human respiratory syncytial virus (RSV), a highly prevalent cause of acute respiratory infections. We evaluated healthy adults experimentally infected with an identical inoculum and infants hospitalized with naturally acquired infections. In aggregate, viral diversification in adults peaked at day 3, with overrepresentation of diversity in the matrix protein 2 (M2) and non-structural protein 2 (NS2) genes. In one subject, delayed viral clearance was accompanied by a late peak of diversity at day 10 in known and predicted B and T cell epitopes. In contrast, infant infections showed much less viral diversity. Our findings suggest multiple overlapping mechanisms for early control of acute viral infections, which may differ between age groups and host immune responses.

Whole-genome, deep pyrosequencing analysis of a duck influenza A virus evolution in swine cells

We studied the sub-population level evolution of a duck influenza A virus isolate during passage in swine tracheal cells. The complete genomes of the A/mallard/Netherlands/10-Nmkt/1999 strain and its swine cell-passaged descendent were analysed by 454 pyrosequencing with coverage depth ranging from several hundred to several thousand reads at any point. This allowed characterization of defined minority sub-populations of gene segments 2, 3, 4, 5, 7, and 8 present in the original isolate. These minority sub-populations ranged between 9.5% (for segment 2) and 46% (for segment 4) of their respective gene segments in the parental stock. They were likely contributed by one or more viruses circulating within the same area, at the same period and in the same or a sympatric host species. The minority sub-populations of segments 3, 4, and 5 became extinct upon viral passage in swine cells, whereas the minority sub-populations of segments 2, 7 and 8 completely replaced their majority counterparts. The swine cell-passaged virus was therefore a three-segment reassortant and also harboured point mutations in segments 3 and 4. The passaged virus was more homogenous than the parental stock, with only 17 minority single nucleotide polymorphisms present above 5% frequency across the whole genome. Though limited here to one sample, this deep sequencing approach highlights the evolutionary versatility of influenza viruses whereby they exploit their genetic diversity, predilection for mixed infection and reassortment to adapt to a new host environmental niche.