Evaluation of the population structure and genetic diversity of Plasmodium falciparum in southern China

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

Population Genetic Analysis of the Theileria annulata Parasites Identified Limited Diversity and Multiplicity of Infection in the Vaccine From India

Background: Apicomplexan parasite Theileria annulata causes significant economic loss to the livestock industry in India and other tropical countries. In India, parasite control is mainly dependent on the live attenuated schizont vaccine and the drug buparvaquone. For effective disease control, it is essential to study the population structure and genetic diversity of the Theileria annulata field isolates and vaccine currently used in India. Methodology/Results: A total of 125 T. annulata isolates were genotyped using 10 microsatellite markers from four states belonging to different geographical locations of India. Limited genetic diversity was observed in the vaccine isolates when compared to the parasites in the field; a level of geographical substructuring was evident in India. The number of genotypes observed per infection was highest in India when compared to other endemic countries, suggesting high transmission intensity and abundance of ticks in the country. A reduced panel of four markers can be used for future studies in these for surveillance of the T. annulata parasites in India. Conclusion: High genetic variation between the parasite populations in the country suggests their successful spread in the field and could hamper the disease control programs. Our findings provide the baseline data for the diversity and population structure of T. annulata parasites from India. The low diversity in the vaccine advocates improving the current vaccine, possibly by increasing its heterozygosity. The reduced panel of the markers identified in this study will be helpful in monitoring parasite and its reintroduction after Theileria eradication.

Genomic Epidemiology of Antimalarial Drug Resistance in Plasmodium falciparum in Southern China

Emerging artemisinin resistance in Southeast Asia poses a significant risk to malaria control and eradication goals, including China’s plan to eliminate malaria nationwide by 2020. Plasmodium falciparum was endemic in China, especially in Southern China. Parasites from this region have shown decreased susceptibility to artemisinin and delayed parasite clearance after artemisinin treatment. Understanding the genetic basis of artemisinin resistance and identifying specific genetic loci associated with this phenotype is crucial for surveillance and containment of resistance. In this study, parasites were collected from clinical patients from Yunnan province and Hainan island. The parasites were genotyped using a P. falciparum-specific single nucleotide polymorphism (SNP) microarray. The SNP profiles examined included a total of 27 validated and candidate molecular markers of drug resistance. The structure of the parasite population was evaluated by principal component analysis by using the EIGENSOFT program, and ADMIXTURE was used to calculate maximum likelihood estimates for the substructure analysis. Parasites showed a high prevalence of resistance haplotypes of pfdhfr and pfdhps and moderate prevalence of pfcrt. There was no mutation identified on pfmdr1. Candidate SNPs on chromosomes 10, 13, and 14 that were associated with delayed parasite clearance showed a low prevalence of mutants. Parasites from Southern China were clustered and separated from those from Southeast Asia. Parasites from Yunnan province were substructured from parasites from Hainan island. This study provides evidence for a genomic population with drug resistance in Southern China and also illustrates the utility of SNP microarrays for large-scale parasite molecular epidemiology.

Genetic evidence for imported malaria and local transmission in Richard Toll, Senegal

BACKGROUND

Malaria elimination efforts can be undermined by imported malaria infections. Imported infections are classified based on travel history.

METHODS

A genetic strategy was applied to better understand the contribution of imported infections and to test for local transmission in the very low prevalence region of Richard Toll, Senegal.

RESULTS

Genetic relatedness analysis, based upon molecular barcode genotyping data derived from diagnostic material, provided evidence for both imported infections and ongoing local transmission in Richard Toll. Evidence for imported malaria included finding that a large proportion of Richard Toll parasites were genetically related to parasites from Thiès, Senegal, a region of moderate transmission with extensive available genotyping data. Evidence for ongoing local transmission included finding parasites of identical genotype that persisted across multiple transmission seasons as well as enrichment of highly related infections within the households of non-travellers compared to travellers.

CONCLUSIONS

These data indicate that, while a large number of infections may have been imported, there remains ongoing local malaria transmission in Richard Toll. These proof-of-concept findings underscore the value of genetic data to identify parasite relatedness and patterns of transmission to inform optimal intervention selection and placement.

Genome-Wide Analysis of Genetic Diversity in Plasmodium falciparum Isolates From China-Myanmar Border

Plasmodium falciparum isolates from China-Myanmar border (CMB) have experienced regional special selective pressures and adaptive evolution. However, the genomes of P. falciparum isolates from this region to date are poorly characterized. Herein, we performed whole-genome sequencing of 34 P. falciparum isolates from CMB and a series of genome-wide sequence analyses to reveal their genetic diversity, population structures, and comparisons with the isolates from other epidemic regions (Thai-Cambodia border, Thai-Myanmar border, and West Africa). Totally 59,720 high-quality single-nucleotide polymorphisms (SNPs) were identified in the P. falciparum isolates from CMB, with average nucleotide diversity (π = 4.59 × 10^-4) and LD decay (132 bp). The Tajima's D and Fu and Li's D values of the CMB isolates were -0.8 (p < 0.05) and -0.84 (p < 0.05), respectively, suggesting a demographic history of recent population expansion or purifying selection. Moreover, 78 genes of the parasite were identified that could be under positive selection, including those genes conferring drug resistance such as pfubp1. In addition, 33 SNPs were identified for tracing the source of the parasites with a high accuracy by analysis of the most differential SNPs among the four epidemic regions. Collectively, our data demonstrated high diversity of the CMB isolates' genomes forming a distinct population, and the identification of 33-SNP barcode provides a valuable surveillance of parasite migration among the regions.

Genetic diversity and drug resistance surveillance of Plasmodium falciparum for malaria elimination: is there an ideal tool for resource-limited sub-Saharan Africa?

The intensification of malaria control interventions has resulted in its global decline, but it remains a significant public health burden especially in sub-Saharan Africa (sSA). Knowledge on the parasite diversity, its transmission dynamics, mechanisms of adaptation to environmental and interventional pressures could help refine or develop new control and elimination strategies. Critical to this is the accurate assessment of the parasite’s genetic diversity and monitoring of genetic markers of anti-malarial resistance across all susceptible populations. Such wide molecular surveillance will require selected tools and approaches from a variety of ever evolving advancements in technology and the changing epidemiology of malaria. The choice of an effective approach for specific endemic settings remains challenging, particularly for countries in sSA with limited access to advanced technologies. This article examines the current strategies and tools for Plasmodium falciparum genetic diversity typing and resistance monitoring and proposes how the different tools could be employed in resource-poor settings. Advanced approaches enabling targeted deep sequencing is valued as a sensitive method for assessing drug resistance and parasite diversity but remains out of the reach of most laboratories in sSA due to the high cost of development and maintenance. It is, however, feasible to equip a limited number of laboratories as Centres of Excellence in Africa (CEA), which will receive and process samples from a network of peripheral laboratories in the continent. Cheaper, sensitive and portable real-time PCR methods can be used in peripheral laboratories to pre-screen and select samples for targeted deep sequence or genome wide analyses at these CEAs.

Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns

Estimates of Plasmodium falciparum migration may inform strategies for malaria elimination. Here we elucidate fine-scale parasite population structure and infer recent migration across Southeast Asia using identity-by-descent (IBD) approaches based on genome-wide single nucleotide polymorphisms called in 1722 samples from 54 districts. IBD estimates are consistent with isolation-by-distance. We observe greater sharing of larger IBD segments between artemisinin-resistant parasites versus sensitive parasites, which is consistent with the recent spread of drug resistance. Our IBD analyses reveal actionable patterns, including isolated parasite populations, which may be prioritized for malaria elimination, as well as asymmetrical migration identifying potential sources and sinks of migrating parasites.

Comparison of msp genotyping and a 24 SNP molecular assay for differentiating Plasmodium falciparum recrudescence from reinfection

Background

Current World Health Organization guidelines for conducting anti-malarial drug efficacy clinical trials recommend genotyping Plasmodium falciparum genes msp1 and msp2 to distinguish recrudescence from reinfection. A more recently developed potential alternative to this method is a molecular genotyping assay based on a panel of 24 single nucleotide polymorphism (SNP) markers.

Methods

Performance parameters of these two genotyping methods were compared using data from two recently completed drug efficacy trials. Blood samples from two anti-malarial therapeutic trials were analysed by both msp genotyping and the 24 SNP assay. Additionally, to conserve time and resources, the statistical program R was used to select the most informative SNPs for a set of unrelated Malawian samples to develop a truncated SNP-based assay for the region surrounding Blantyre, Malawi. The ability of this truncated assay to distinguish reinfection from recrudescence when compared to the full 24 SNP assay was then analysed using data from the therapeutic trials.

Results

A total of 360 samples were analysed; 66 for concordance of msp and SNP barcoding methodologies, and 294 for assessing the most informative of the 24 SNP markers. SNP genotyping performed comparably to msp genotyping, with only one case of disagreement among the 50 interpretable results, where the SNP assay identified the sample as reinfection and the msp typing as recrudescence. Furthermore, SNP typing was more robust; only 6% of samples were uninterpretable by SNP typing, compared to 19.7% when msp genotyping was used. For discriminating reinfection from recrudescence, a truncated 6 SNP assay was found to perform at 95.1% the accuracy of the full 24 SNP bar code.

Conclusions

The use of SNP analysis has similar sensitivity to the standard msp genotyping in determining recrudescence from reinfection. Although more expensive, SNP typing is faster and less work intensive. Limiting the assay to those SNPs most informative in the geographical region of interest may further decrease the workload and the cost, making this technique a feasible and affordable alternative in drug efficacy trials.

Electronic supplementary material

The online version of this article (10.1186/s12936-019-2695-0) contains supplementary material, which is available to authorized users.

Population Genetic Analysis of Theileria annulata from Six Geographical Regions in China, Determined on the Basis of Micro- and Mini-satellite Markers

Theileria annulata, a tick-borne apicomplexan protozoan, causes a lymphoproliferative disease of cattle with high prevalence in tropical and sub-tropical regions. Understanding the genetic diversity and structure of local populations will provide more fundamental knowledge for the population genetics and epidemics of protozoa. In this study, 78 samples of T. annulata collected from cattle/yaks representing 6 different geographic populations in China were genotyped using eight micro- and mini-satellite markers. High genetic variation within population, moderate genetic differentiation, and high level of diversity co-occurring with significant linkage disequilibrium were observed, which indicates there is gene flow between these populations in spite of the existence of reproductive and geographical barriers among populations. Furthermore, some degree of genetic differentiation was also found between samples from China and Oman. These findings provide a first glimpse of the genetic diversity of the T. annulata populations in China, and might contribute to the knowledge of distribution, dynamics, and epidemiology of T. annulata populations and optimize the management strategies for control.

Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes

Ethiopia is one of the few African countries where Plasmodium vivax is co-endemic with P. falciparum. Malaria transmission is seasonal and transmission intensity varies mainly by landscape and climate. Although the recent emergence of drug resistant parasites presents a major issue to malaria control in Ethiopia, little is known about the transmission pathways of parasite species and prevalence of resistant markers. This study used microsatellites to determine population diversity and gene flow patterns of P. falciparum (N = 226) and P. vivax (N = 205), as well as prevalence of drug resistant markers to infer the impact of gene flow and existing malaria treatment regimes. Plasmodium falciparum indicated a higher rate of polyclonal infections than P. vivax. Both species revealed moderate genetic diversity and similar population structure. Populations in the northern highlands were closely related to the eastern Rift Valley, but slightly distinct from the southern basin area. Gene flow via human migrations between the northern and eastern populations were frequent and mostly bidirectional. Landscape genetic analyses indicated that environmental heterogeneity and geographical distance did not constrain parasite gene flow. This may partly explain similar patterns of resistant marker prevalence. In P. falciparum, a high prevalence of mutant alleles was detected in codons related to chloroquine (pfcrt and pfmdr1) and sulfadoxine-pyrimethamine (pfdhps and pfdhfr) resistance. Over 60% of the samples showed pfmdr1 duplications. Nevertheless, no mutation was detected in pfK13 that relates to artemisinin resistance. In P. vivax, while sequences of pvcrt-o were highly conserved and less than 5% of the samples showed pvmdr duplications, over 50% of the samples had pvmdr1 976F mutation. It remains to be tested if this mutation relates to chloroquine resistance. Monitoring the extent of malaria spread and markers of drug resistance is imperative to inform policy for evidence-based antimalarial choice and interventions. To effectively reduce malaria burden in Ethiopia, control efforts should focus on seasonal migrant populations.

Controlling Aedes albopictus and Culex pipiens pallens using silver nanoparticles synthesized from aqueous extract of Cassia fistula fruit pulp and its mode of action

Mosquitoes act as key vector for transmission of devastating parasites and pathogens which affect millions of people globally. In this research, the green synthesis of silver nanoparticles of Cassia fistula fruit pulp as an innovative and operative tool against vector mosquitoes is presented. Silver nanoparticles were characterized by a series of techniques including Fourier transform infrared spectroscopy, Transmission Electron Microscope and confirmed by Scanning Electron Microscope, UV-Vis spectrophotometry and X-ray diffraction. Silver nanoparticles were highly effective against the larvae (I-IV instar) and pupae of Aedes albopictus and Culex pipiens pallens after 24, 48 and 72 h of treatment. Ae. albopictus had LC50 values ranging from 8.3 mg/L (I instar) to 17.3 mg/L (pupae) and LC50 ranging from 1.1 mg/L (I instar) to 19.0 mg/L (pupae) against Cx. pipiens pallens. The systemic effect of AgNPs was further assessed in the fourth instar larvae of Ae. albopictus and Cx. pipiens pallens by measuring the levels of total proteins and activity of two important marker enzymes: Acetylcholinesterase and α- and β-carboxylesterase. Overall, the findings of the study suggest that the use of Cassia fistula-fruit pulp extract mediated synthesis of silver nanoparticles can be used for controlling vector mosquitoes. This is the first report on the mosquito larvicidal and pupicidal activity of AgNPs synthesized by Cassia fistula fruit pulp and its possible mechanism of action.

Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution

To determine whether the major human malaria parasite Plasmodium falciparum exhibits fragmented population structure or local adaptation at the northern limit of its African distribution where the dry Sahel zone meets the Sahara, samples were collected from diverse locations within Mauritania over a range of ~1000 km. Microsatellite genotypes were obtained for 203 clinical infection samples from eight locations, and Illumina paired‐end sequences were obtained to yield high coverage genomewide single nucleotide polymorphism (SNP) data for 65 clinical infection samples from four locations. Most infections contained single parasite genotypes, reflecting low rates of transmission and superinfection locally, in contrast to the situation seen in population samples from countries further south. A minority of infections shared related or identical genotypes locally, indicating some repeated transmission of parasite clones without recombination. This caused some multilocus linkage disequilibrium and local divergence, but aside from the effect of repeated genotypes there was minimal differentiation between locations. Several chromosomal regions had elevated integrated haplotype scores (|iHS|) indicating recent selection, including those containing drug resistance genes. A genomewide F_ST scan comparison with previous sequence data from an area in West Africa with higher infection endemicity indicates that regional gene flow prevents genetic isolation, but revealed allele frequency differentiation at three drug resistance loci and an erythrocyte invasion ligand gene. Contrast of extended haplotype signatures revealed none to be unique to Mauritania. Discrete foci of infection on the edge of the Sahara are genetically highly connected to the wider continental parasite population, and local elimination would be difficult to achieve without very substantial reduction in malaria throughout the region.

THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites

As many malaria-endemic countries move towards elimination of Plasmodium falciparum, the most virulent human malaria parasite, effective tools for monitoring malaria epidemiology are urgent priorities. P. falciparum population genetic approaches offer promising tools for understanding transmission and spread of the disease, but a high prevalence of multi-clone or polygenomic infections can render estimation of even the most basic parameters, such as allele frequencies, challenging. A previous method, COIL, was developed to estimate complexity of infection (COI) from single nucleotide polymorphism (SNP) data, but relies on monogenomic infections to estimate allele frequencies or requires external allele frequency data which may not available. Estimates limited to monogenomic infections may not be representative, however, and when the average COI is high, they can be difficult or impossible to obtain. Therefore, we developed THE REAL McCOIL, Turning HEterozygous SNP data into Robust Estimates of ALelle frequency, via Markov chain Monte Carlo, and Complexity Of Infection using Likelihood, to incorporate polygenomic samples and simultaneously estimate allele frequency and COI. This approach was tested via simulations then applied to SNP data from cross-sectional surveys performed in three Ugandan sites with varying malaria transmission. We show that THE REAL McCOIL consistently outperforms COIL on simulated data, particularly when most infections are polygenomic. Using field data we show that, unlike with COIL, we can distinguish epidemiologically relevant differences in COI between and within these sites. Surprisingly, for example, we estimated high average COI in a peri-urban subregion with lower transmission intensity, suggesting that many of these cases were imported from surrounding regions with higher transmission intensity. THE REAL McCOIL therefore provides a robust tool for understanding the molecular epidemiology of malaria across transmission settings.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.