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

New reference genomes to distinguish the sympatric malaria parasites, Plasmodium ovale curtisi and Plasmodium ovale wallikeri

Despite Plasmodium ovale curtisi (Poc) and wallikeri (Pow) being important human-infecting malaria parasites that are widespread across Africa and Asia, little is known about their genome diversity. Morphologically identical, Poc and Pow are indistinguishable and commonly misidentified. Recent rises in the incidence of Poc/Pow infections have renewed efforts to address fundamental knowledge gaps in their biology, and to develop diagnostic tools to understand their epidemiological dynamics and malaria burden. A major roadblock has been the incompleteness of available reference assemblies (PocGH01, PowCR01; ~ 33.5 Mbp). Here, we applied multiple sequencing platforms and advanced bioinformatics tools to generate new reference genomes, Poc221 (South Sudan; 36.0 Mbp) and Pow222 (Nigeria; 34.3 Mbp), with improved nuclear genome contiguity (> 4.2 Mbp), annotation and completeness (> 99% Plasmodium spp., single copy orthologs). Subsequent sequencing of 6 Poc and 15 Pow isolates from Africa revealed a total of 22,517 and 43,855 high-quality core genome SNPs, respectively. Genome-wide levels of nucleotide diversity were determined to be 2.98 × 10-4 (Poc) and 3.43 × 10-4 (Pow), comparable to estimates for other Plasmodium species. Overall, the new reference genomes provide a robust foundation for dissecting the biology of Poc/Pow, their population structure and evolution, and will contribute to uncovering the recombination barrier separating these species.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Whole-Genome Resequencing Reveals Selection Signal Related to Sheep Wool Fineness

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

Whole Genome Sequencing Contributions and Challenges in Disease Reduction Focused on Malaria

Malaria elimination remains an important goal that requires the adoption of sophisticated science and management strategies in the era of the COVID-19 pandemic. The advent of next generation sequencing (NGS) is making whole genome sequencing (WGS) a standard today in the field of life sciences, as PCR genotyping and targeted sequencing provide insufficient information compared to the whole genome. Thus, adapting WGS approaches to malaria parasites is pertinent to studying the epidemiology of the disease, as different regions are at different phases in their malaria elimination agenda. Therefore, this review highlights the applications of WGS in disease management, challenges of WGS in controlling malaria parasites, and in furtherance, provides the roles of WGS in pursuit of malaria reduction and elimination. WGS has invaluable impacts in malaria research and has helped countries to reach elimination phase rapidly by providing required information needed to thwart transmission, pathology, and drug resistance. However, to eliminate malaria in sub-Saharan Africa (SSA), with high malaria transmission, we recommend that WGS machines should be readily available and affordable in the region.

Genome-Wide Selective Signatures Reveal Candidate Genes Associated with Hair Follicle Development and Wool Shedding in Sheep

Hair follicle development and wool shedding in sheep are poorly understood. This study investigated the population structures and genetic differences between sheep with different wool types to identify candidate genes related to these traits. We used Illumina ovine SNP 50K chip genotyping data of 795 sheep populations comprising 27 breeds with two wool types, measuring the population differentiation index (Fst), nucleotide diversity (θπ ratio), and extended haplotype homozygosity among populations (XP-EHH) to detect the selective signatures of hair sheep and fine-wool sheep. The top 5% of the Fst and θπ ratio values, and values of XP-EHH < −2 were considered strongly selected SNP sites. Annotation showed that the PRX, SOX18, TGM3, and TCF3 genes related to hair follicle development and wool shedding were strongly selected. Our results indicated that these methods identified important genes related to hair follicle formation, epidermal differentiation, and hair follicle stem cell development, and provide a meaningful reference for further study on the molecular mechanisms of economically important traits in sheep.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Infection, genetics, and evolution of Trichinella: Historical insights and applications to molecular epidemiology

Genetic variation in pathogen populations provides the means to answer questions in disease ecology and transmission, illuminating interactions between genetic traits, environmental exposures, and disease. Such studies elucidate the phylogeny, evolution, transmission and pathogenesis of viruses, bacteria and parasites. Here, we review how such studies have fostered understanding of the biology and epidemiology of zoonotic nematode parasites in the genus Trichinella spp., which impose considerable economic and health burdens by infecting wildlife, livestock, and people. To use such data to define ongoing chains of local transmission and source traceback, researchers first must understand the extent and distribution of genetic variation resident in regional parasite populations. Thus, genetic variability illuminates a population's past as well as its present. Here we review how such data have helped define population dynamics of Trichinella spp. in wild and domesticated hosts, creating opportunities to harness genetic variation in the quest to prevent, track, and contain future outbreaks.

Isolation of Mutants With Reduced Susceptibility to Piperaquine From a Mutator of the Rodent Malaria Parasite Plasmodium berghei

Elucidation of the mechanisms of drug resistance in malaria parasites is crucial for combatting the emergence and spread of resistant parasites, which can be achieved by tracing resistance-associated mutations and providing useful information for drug development. Previously, we produced a novel genetic tool, a Plasmodium berghei mutator (PbMut), whose base substitution rate is 36.5 times higher than that of wild-type parasites. Here, we report the isolation of a mutant with reduced susceptibility to piperaquine (PPQ) from PbMut under PPQ pressure by sequential nine-cycle screening and named it PbMut-PPQ-R-P9. The ED50 of PbMut-PPQ-R-P9 was 1.79 times higher than that of wild-type parasites, suggesting that its PPQ resistance is weak. In the 1st screen, recrudescence occurred in the mice infected with PbMut but not in those infected with wild-type parasites, suggesting earlier emergence of PPQ-resistant parasites from PbMut. Whole-genome sequence analysis of PbMut-PPQ-R-P9 clones revealed that eight nonsynonymous mutations were conserved in all clones, including N331I in PbCRT, the gene encoding chloroquine resistance transporter (CRT). The PbCRT(N331I) mutation already existed in the parasite population after the 2nd screen and was predominant in the population after the 8th screen. An artificially inserted PbCRT(N331I) mutation gave rise to reduced PPQ susceptibility in genome-edited parasites (PbCRT-N331I). The PPQ susceptibility and growth rates of PbCRT-N331I parasites were significantly lower than those of PbMut-PPQ-R-P9, implying that additional mutations in the PbMut-PPQ-R9 parasites could compensate for the fitness cost of the PbCRT(N331I) mutation and contribute to reduced PPQ susceptibility. In summary, PbMut could serve as a novel genetic tool for predicting gene mutations responsible for drug resistance. Further study on PbMut-PPQ-R-P9 could identify genetic changes that compensate for fitness costs owing to drug resistance acquisition.

Genomic analysis reveals independent evolution of Plasmodium falciparum populations in Ethiopia

Background

Plasmodium falciparum parasite populations in Ethiopia have been experiencing local selective pressures from drugs and immunity, leading to evolutionary adaptation. However, there was a paucity of data on genomic characterization and evolutionary adaptations of P. falciparum isolates from the central area of Ethiopia.

Methods

Whole-genome analysis of 25 P. falciparum isolates from central Ethiopia, specifically from West Arsi, were studied to determine their genetic diversity, population structures, and signatures of selection in known drug resistance alleles against global isolates from Cambodia, Thailand, DR Congo, and Malawi.

Results

A total of 18,517 high-quality single-nucleotide polymorphisms (SNPs) were identified in Ethiopian P. falciparum isolates. About 84% of the Ethiopian P. falciparum isolates had a F_WS value > 0.95 showing a dominant single genotype infection in most isolates at the time of collection with little potential for out-crossing as expected in areas with low transmission intensity. Within-host diversity of Ethiopian infections was significantly different from East African (p < 0.001), but not Southeast Asian infections (P > 0.05). A significant population structure has been observed by PCA and population differentiation between Ethiopian parasites and East African (Fst ~ 10%) and Southeast Asian populations (Fst ~ 18%), suggesting limited gene flow and the independent evolution of the Ethiopian parasite population. Moreover, a total of 125 genes under balancing selection was found that include ama1, trap, eba175, and lsa3, previously identified as targets of human host immunity. Recent directional selection analysis using integrated standardized haplotype score (IHS) did not detect any selection signatures in the Pfcrt, Pfdhfr, Pfdhps, Pfmdr1, and PfK13 genes. However, known drug resistance-conferring mutations analysis showed that at least one SNP marker was fixed in these genes, but not in Pfdhps and PfK13.

Conclusion

Plasmodium falciparum populations in the central region of Ethiopia was structurally diverged from both Southeast Asian and other East African populations. Malaria infections in Ethiopia had low within-host diversity, and parasites carry fixed chloroquine resistance markers despite the withdrawal of this drug for the treatment of P. falciparum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03660-y.

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.

Genome-Wide Analysis of the Malaria Parasite Plasmodium falciparum Isolates From Togo Reveals Selective Signals in Immune Selection-Related Antigen Genes

Malaria is a public health concern worldwide, and Togo has proven to be no exception. Effective approaches to provide information on biological insights for disease elimination are therefore a research priority. Local selection on malaria pathogens is due to multiple factors including host immunity. We undertook genome-wide analysis of sequence variation on a sample of 10 Plasmodium falciparum (Pf) clinical isolates from Togo to identify local-specific signals of selection. Paired-end short-read sequences were mapped and aligned onto > 95% of the 3D7 Pf reference genome sequence in high fold coverage. Data on 266 963 single nucleotide polymorphisms were obtained, with average nucleotide diversity π = 1.79 × 10⁻³. Both principal component and neighbor-joining tree analyses showed that the Togo parasites clustered according to their geographic (Africa) origin. In addition, the average genome-wide diversity of Pf from Togo was much higher than that from other African samples. Tajima’s D value of the Togo isolates was −0.56, suggesting evidence of directional selection and/or recent population expansion. Against this background, within-population analyses identifying loci of balancing and recent positive selections evidenced that host immunity has been the major selective agent. Importantly, 87 and 296 parasite antigen genes with Tajima’s D values > 1 and in the top 1% haplotype scores, respectively, include a significant representation of membrane proteins at the merozoite stage that invaded red blood cells (RBCs) and parasitized RBCs surface proteins that play roles in immunoevasion, adhesion, or rosetting. This is consistent with expectations that elevated signals of selection due to allele-specific acquired immunity are likely to operate on antigenic targets. Collectively, our data suggest a recent expansion of Pf population in Togo and evidence strong host immune selection on membrane/surface antigens reflected in signals of balancing/positive selection of important gene loci. Findings from this study provide a fundamental basis to engage studies for effective malaria control in Togo.

Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil

BACKGROUND

Plasmodium vivax is a neglected human malaria parasite that causes significant morbidity in the Americas, the Middle East, Asia, and the Western Pacific. Population genomic approaches remain little explored to map local and regional transmission pathways of P. vivax across the main endemic sites in the Americas, where great progress has been made towards malaria elimination over the past decades.

METHODOLOGY/PRINCIPAL FINDINGS

We analyze 38 patient-derived P. vivax genome sequences from Mâncio Lima (ML)-the Amazonian malaria hotspot next to the Brazil-Peru border-and 24 sequences from two other sites in Acre State, Brazil, a country that contributes 23% of malaria cases in the Americas. We show that the P. vivax population of ML is genetically diverse (π = 4.7 × 10-4), with a high polymorphism particularly in genes encoding proteins putatively involved in red blood cell invasion. Paradoxically, however, parasites display strong genome-wide linkage disequilibrium, being fragmented into discrete lineages that are remarkably stable across time and space, with only occasional recombination between them. Using identity-by-descent approaches, we identified a large cluster of closely related sequences that comprises 16 of 38 genomes sampled in ML over 26 months. Importantly, we found significant ancestry sharing between parasites at a large geographic distance, consistent with substantial gene flow between regional P. vivax populations.

CONCLUSIONS/SIGNIFICANCE

We have characterized the sustained expansion of highly inbred P. vivax lineages in a malaria hotspot that can seed regional transmission. Potential source populations in hotspots represent a priority target for malaria elimination in the Amazon.

Experimentally Engineered Mutations in a Ubiquitin Hydrolase, UBP-1, Modulate In Vivo Susceptibility to Artemisinin and Chloroquine in Plasmodium berghei

As resistance to artemisinins (current frontline drugs in malaria treatment) emerges in Southeast Asia, there is an urgent need to identify the genetic determinants and understand the molecular mechanisms underpinning such resistance. Such insights could lead to prospective interventions to contain resistance and prevent the eventual spread to other regions where malaria is endemic. Reduced susceptibility to artemisinin in Southeast Asia has been primarily linked to mutations in the Plasmodium falciparum Kelch-13 gene, which is currently widely recognized as a molecular marker of artemisinin resistance.

As resistance to artemisinins (current frontline drugs in malaria treatment) emerges in Southeast Asia, there is an urgent need to identify the genetic determinants and understand the molecular mechanisms underpinning such resistance. Such insights could lead to prospective interventions to contain resistance and prevent the eventual spread to other regions where malaria is endemic. Reduced susceptibility to artemisinin in Southeast Asia has been primarily linked to mutations in the Plasmodium falciparum Kelch-13 gene, which is currently widely recognized as a molecular marker of artemisinin resistance. However, two mutations in a ubiquitin hydrolase, UBP-1, have been previously associated with reduced artemisinin susceptibility in a rodent model of malaria, and some cases of UBP-1 mutation variants associated with artemisinin treatment failure have been reported in Africa and SEA. In this study, we employed CRISPR-Cas9 genome editing and preemptive drug pressures to test these artemisinin susceptibility-associated mutations in UBP-1 in Plasmodium berghei sensitive lines in vivo. Using these approaches, we show that the V2721F UBP-1 mutation results in reduced artemisinin susceptibility, while the V2752F mutation results in resistance to chloroquine (CQ) and moderately impacts tolerance to artemisinins. Genetic reversal of the V2752F mutation restored chloroquine sensitivity in these mutant lines, whereas simultaneous introduction of both mutations could not be achieved and appears to be lethal. Interestingly, these mutations carry a detrimental growth defect, which would possibly explain their lack of expansion in natural infection settings. Our work provides independent experimental evidence on the role of UBP-1 in modulating parasite responses to artemisinin and chloroquine under in vivo conditions.