Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region.

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

Genomic Analysis of Plasmodium vivax in Southern Ethiopia Reveals Selective Pressures in Multiple Parasite Mechanisms.

The Horn of Africa harbors the largest reservoir of Plasmodium vivax in the continent. Most of sub-Saharan Africa has remained relatively vivax-free due to a high prevalence of the human Duffy-negative trait, but the emergence of strains able to invade Duffy-negative reticulocytes poses a major public health threat. We undertook the first population genomic investigation of P. vivax from the region, comparing the genomes of 24 Ethiopian isolates against data from Southeast Asia to identify important local adaptions. The prevalence of the Duffy binding protein amplification in Ethiopia was 79%, potentially reflecting adaptation to Duffy negativity. There was also evidence of selection in a region upstream of the chloroquine resistance transporter, a putative chloroquine-resistance determinant. Strong signals of selection were observed in genes involved in immune evasion and regulation of gene expression, highlighting the need for a multifaceted intervention approach to combat P. vivax in the region.

Genomic analysis of dengue virus serotype 1 (DENV-1) genotypes from Surabaya, Indonesia.

Dengue has caused a significant public health impact globally. With the diverse genetic of the causative viruses, analysis of dengue virus (DENV) genomes is important to supplement epidemiological data with information that can be used to reconstruct the history of epidemics in time and space. We have reported the clinical and virological characteristics of dengue in Surabaya, Indonesia and revealed the presence of all four DENV serotypes and the predominance of DENV-1. The further classification of Surabaya DENV-1 into two different genotypes warrants in-depth genomic analysis to study the dynamics of both genotypes and their contribution to virus evolution, virus transmission, and disease. We performed full-length genome sequencing to nine isolates' representatives from DENV-1 Genotype I and Genotype IV. Phylogenetic and evolutionary analyses suggested the more recent introduction of Genotype I viruses compared to the more endemic Genotype IV. Comparative analysis of Surabaya DENV-1 genomes and other sequences available publicly revealed that the majority of the DENV-1 codons were under strong purifying selection, while seven codon sites identified to be under positive selection. We highlight a unique codon site under the positive pressure in the NS1 gene of DENV-1. Our results provide additional genomic data of DENV from Indonesia that may contribute to the better understanding of dengue disease dynamics.

Transcription factor 7-like 2 single nucleotide polymorphisms are associated with lipid profile in the Balinese.

Transcription factor 7-like 2 (TCF7L2) protein plays an important role in glucose and lipid metabolisms. Single nucleotide polymorphisms (SNPs) in the TCF7L2 gene contribute to increased fasting plasma glucose (FPG) and body mass index (BMI), and altered lipid concentrations in various population. We investigated whether the TCF7L2 SNPs were associated with obesity, high FPG and altered lipid profile in the Balinese. A total of 608 Balinese from rural and urban Bali, Indonesia, were recruited. Triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), total cholesterol (TC) and FPG were measured, and BMI was calculated. Ratios of TG/HDL-C, LDL-C/HDL-C, and TC/HDL-C were determined. Genotyping of SNPs rs7903146, rs10885406, and rs12255372 were done in all samples. Genetic association analyses under a dominant model showed that the rs7903146 (OR 5.50, 95% CI 2.34-12.91, p = 8.5 × 10-5), rs12255372 (OR 4.15, 95% CI 1.66-10.33, p = 0.003) and rs10885406 (OR 2.43, 95% CI 1.39-4.25, p = 0.003) were significantly associated with high TC/HDL-C ratio. The rs10885406 also presented a significant association with high TG (OR 2.21, 95% CI 1.29-3.81, p = 0.004) and low HDL-C (OR 3.02, 95% CI 1.58-5.80, p = 0.001) concentrations, as well as high TG/HDL-C ratio (OR 1.95, 95% CI 1.16-3.27, p = 0.013). None of the SNPs exhibited significant association with obesity or high FPG. SNPs in the TCF7L2 gene are associated with altered lipid profile in the Balinese.

Production of recombinant dengue non-structural 1 (NS1) proteins from clinical virus isolates.

Dengue is a febrile disease caused by infection of dengue virus (DENV). Early diagnosis of dengue infection is important for better management of the disease. The DENV Non-Structural Protein 1 (NS1) antigen has been routinely used for the early dengue detection. In dengue epidemic countries such as Indonesia, clinicians are increasingly relying on the NS1 detection for confirmation of dengue infection. Various NS1 diagnostic tests are commercially available, however different sensitivities and specificities were observed in various settings. This study was aimed to generate dengue NS1 recombinant protein for the development of dengue diagnostic tests. Four Indonesian DENV isolates were used as the source of the NS1 gene cloning, expression, and purification in bacterial expression system. Recombinant NS1 proteins were successfully purified and their antigenicities were assessed. Immunization of mice with recombinant proteins observed the immunogenicity of the NS1 protein. The generated recombinant proteins can be potentially used in the development of NS1 diagnostic test. With minimal modifications, this method can be used for producing NS1 recombinant proteins from isolates obtained from other geographical regions.

Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow.

BACKGROUND: Plasmodium vivax malaria remains a major public health burden in Myanmar. Resistance to chloroquine (CQ), the first-line treatment for P. vivax, has been reported in the country and has potential to undermine local control efforts. METHODS: Patients over 6 years of age with uncomplicated P. vivax mono-infection were enrolled into clinical efficacy studies in Myawaddy in 2014 and Kawthoung in 2012. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. Pvmdr1 copy number (CN) and microsatellite diversity were assessed on samples from the patients enrolled in the clinical study and additional cross-sectional surveys undertaken in Myawaddy and Shwegyin in 2012. RESULTS: A total of 85 patients were enrolled in the CQ clinical studies, 25 in Myawaddy and 60 in Kawthoung. One patient in Myawaddy (1.2%) had an early treatment failure and two patients (2.3%) in Kawthoung presented with late treatment failures on day 28. The day 28 efficacy was 92.0% (95% CI 71.6-97.9) in Myawaddy and 98.3% (95% CI 88.7-99.8) in Kawthoung. By day 2, 92.2% (23/25) in Myawaddy and 85.0% (51/60) in Kawthoung were aparasitaemic. Genotyping and pvmdr1 CN assessment was undertaken on 43, 52 and 46 clinical isolates from Myawaddy, Kawthoung and Shwegyin respectively. Pvmdr1 amplification was observed in 3.2% (1/31) of isolates in Myawaddy, 0% (0/49) in Kawthoung and 2.5% (1/40) in Shwegyin. Diversity was high in all sites (H E 0.855-0.876), with low inter-population differentiation (F ST 0.016-0.026, P < 0.05). CONCLUSIONS: Treatment failures after chloroquine were observed following chloroquine monotherapy, with pvmdr1 amplification present in both Myawaddy and Shwegyin. The results emphasize the importance of ongoing P. vivax drug resistance surveillance in Myanmar, particularly given the potential connectivity between parasite population at different sites.

Phylogenetic and evolutionary analyses of dengue viruses isolated in Jakarta, Indonesia.

Dengue has affected Indonesia for the last five decades and become a major health problem in many cities in the country. Jakarta, the capital of Indonesia, reports dengue cases annually, with several outbreaks documented. To gain information on the dynamic and evolutionary history of dengue virus (DENV) in Jakarta, we conducted phylogenetic and evolutionary analyses of DENV isolated in 2009. Three hundred thirty-three dengue-suspected patients were recruited. Our data revealed that dengue predominantly affected young adults, and the majority of cases were due to secondary infection. A total of 171 virus isolates were successfully serotyped. All four DENV serotypes were circulating in the city, and DENV-1 was the predominant serotype. The DENV genotyping of 17 isolates revealed the presence of Genotypes I and IV in DENV-1, while DENV-2 isolates were grouped into the Cosmopolitan genotype. The grouping of isolates into Genotype I and II was seen for DENV-3 and DENV-4, respectively. Evolutionary analysis revealed the relatedness of Jakarta isolates with other isolates from other cities in Indonesia and isolates from imported cases in other countries. We revealed the endemicity of DENV and the role of Jakarta as the potential source of imported dengue cases in other countries. Our study provides genetic information regarding DENV from Jakarta, which will be useful for upstream applications, such as the study of DENV epidemiology and evolution and transmission dynamics.

A zoonotic human infection with simian malaria, Plasmodium knowlesi, in Central Kalimantan, Indonesia.

BACKGROUND: The Indonesian archipelago is endemic for malaria. Although Plasmodium falciparum and P. vivax are the most common causes for malaria cases, P. malariae and P. ovale are also present in certain regions. Zoonotic case of malaria had just became the attention of public health communities after the Serawak study in 2004. However, zoonotic case in Indonesia is still under reported; only one published report of knowlesi malaria in South Kalimantan in 2010. CASE PRESENTATION: A case of Plasmodium knowlesi infection in a worker from a charcoal mining company in Central Kalimantan, Indonesia was described. The worker suffered from fever following his visit to a lowland forest being cut and converted into a new mining location. CONCLUSION: This study confirmed a zoonotic infection using polymerase chain reaction amplification and Sanger sequencing of plasmodial DNA encoding the mitochondrial cytochrome c oxidase subunit I (mtCOI).

Where chloroquine still works: the genetic make-up and susceptibility of Plasmodium vivax to chloroquine plus primaquine in Bhutan.

BACKGROUND: Bhutan has made substantial progress in reducing malaria incidence. The national guidelines recommend chloroquine (CQ) and primaquine (PQ) for radical cure of uncomplicated Plasmodium vivax, but the local efficacy has not been assessed. The impact of cases imported from India on the genetic make-up of the local vivax populations is currently unknown. METHODS: Patients over 4 years of age with uncomplicated P. vivax mono-infection were enrolled into a clinical efficacy study and molecular survey. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. On day 28 a 14-day regimen of PQ (0.25 mg/kg/day) was commenced under direct observation. After day 42, patients were followed up monthly for a year. The primary and secondary endpoints were risk of treatment failure at day 28 and at 1 year. Parasite genotyping was undertaken at nine tandem repeat markers, and standard population genetic metrics were applied to examine population diversity and structure in infections thought to be acquired inside or outside of Bhutan. RESULTS: A total of 24 patients were enrolled in the clinical study between April 2013 and October 2015. Eight patients (33.3 %) were lost to follow-up in the first 6 months and another eight patients lost between 6 and 12 months. No (0/24) treatment failures occurred by day 28 and no (0/8) parasitaemia was detected following PQ treatment. Some 95.8 % (23/24) of patients were aparasitaemic by day 2. There were no haemolytic or serious events. Genotyping was undertaken on parasites from 12 autochthonous cases and 16 suspected imported cases. Diversity was high (H E 0.87 and 0.90) in both populations. There was no notable differentiation between the autochthonous and imported populations. CONCLUSIONS: CQ and PQ remains effective for radical cure of P. vivax in Bhutan. The genetic analyses indicate that imported infections are sustaining the local vivax population, with concomitant risk of introducing drug-resistant strains.

Genetic diversity and population structure of Plasmodium vivax in Central China.

BACKGROUND: In Central China the declining incidence of Plasmodium vivax has been interrupted by epidemic expansions and imported cases. The impact of these changes on the local parasite population, and concurrent risks of future resurgence, was assessed. METHODS: Plasmodium vivax isolates collected from Anhui and Jiangsu provinces, Central China between 2007 and 2010 were genotyped using capillary electrophoresis at seven polymorphic short tandem repeat markers. Spatial and temporal analyses of within-host and population diversity, population structure, and relatedness were conducted on these isolates. RESULTS: Polyclonal infections were infrequent in the 94 isolates from Anhui (4%) and 25 from Jiangsu (12%), with a trend for increasing frequency from 2008 to 2010 (2 to 19%) when combined. Population diversity was high in both provinces and across the years tested (H(E) = 0.8 - 0.85). Differentiation between Anhui and Jiangsu was modest (F'(ST) = 0.1). Several clusters of isolates with identical multi-locus haplotypes were observed across both Anhui and Jiangsu. Linkage disequilibrium was strong in both populations and in each year tested (I(A)(S) = 0.2 - 0.4), but declined two- to four-fold when identical haplotypes were accounted for, indicative of occasional epidemic transmission dynamics. None of five imported isolates shared identical haplotypes to any of the central Chinese isolates. CONCLUSIONS: The population genetic structure of P. vivax in Central China highlights unstable transmission, with limited barriers to gene flow between the central provinces. Despite low endemicity, population diversity remained high, but the reservoirs sustaining this diversity remain unclear. The challenge of imported cases and risks of resurgence emphasize the need for continued surveillance to detect early warning signals. Although parasite genotyping has potential to inform the management of outbreaks, further studies are required to identify suitable marker panels for resolving local from imported P. vivax isolates.

Whole-Genome Sequencing and Mutation Analyses of SARS-CoV-2 Isolates from Indonesia.

The SARS-CoV-2 infection that caused the COVID-19 pandemic has become a significant public health concern. New variants with distinct mutations have emerged, potentially impacting its infectivity, immune evasion capacity, and vaccine response. A whole-genome sequencing study of 292 SARS-CoV-2 isolates collected from selected regions of Indonesia between January and October 2021 was performed to identify the distribution of SARS-CoV-2 variants and common mutations in Indonesia. During January-April 2021, Indonesian lineages B.1.466.2 and B.1.470 dominated, but from May 2021, Delta's AY.23 lineage outcompeted them. An analysis of 7515 published sequences from January 2021 to June 2022 revealed a decline in Delta in November 2021, followed by the emergence of Omicron variants in December 2021. We identified C241T (5'UTR), P314L (NSP12b), F106F (NSP3), and D614G (Spike) mutations in all sequences. The other common substitutions included P681R (76.4%) and T478K (60%) in Spike, D377Y in Nucleocapsid (61%), and I82T in Membrane (60%) proteins. Breakthrough infection and prolonged viral shedding cases were associated with Delta variants carrying the Spike T19R, G142D, L452R, T478K, D614G, P681R, D950N, and V1264L mutations. The dynamic of SARS-CoV-2 variants in Indonesia highlights the importance of continuous genomic surveillance in monitoring and identifying potential strains leading to disease outbreaks.

Performance of commercial dengue NS1 ELISA and molecular analysis of NS1 gene of dengue viruses obtained during surveillance in Indonesia.

BACKGROUND: Early diagnosis of dengue infection is crucial for better management of the disease. Diagnostic tests based on the detection of dengue virus (DENV) Non Structural Protein 1 (NS1) antigen are commercially available with different sensitivities and specificities observed in various settings. Dengue is endemic in Indonesia and clinicians are increasingly using the NS1 detection for dengue confirmation. This study described the performance of Panbio Dengue Early NS1 and IgM Capture ELISA assays for dengue detection during our surveillance in eight cities in Indonesia as well as the genetic diversity of DENV NS1 genes and its relationship with the NS1 detection. METHODS: The NS1 and IgM/IgG ELISA assays were used for screening and confirmation of dengue infection during surveillance in 2010-2012. Collected serum samples (n = 440) were subjected to RT-PCR and virus isolation, in which 188 samples were confirmed for dengue infection. The positivity of the ELISA assays were correlated with the RT-PCR results to obtain the sensitivity of the assays. The NS1 genes of 48 Indonesian virus isolates were sequenced and their genetic characteristics were studied. RESULTS: Using molecular data as gold standard, the sensitivity of NS1 ELISA assay for samples from Indonesia was 56.4% while IgM ELISA was 73.7%. When both NS1 and IgM results were combined, the sensitivity increased to 89.4%. The NS1 sensitivity varied when correlated with city/geographical origins and DENV serotype, in which the lowest sensitivity was observed for DENV-4 (19.0%). NS1 sensitivity was higher in primary (67.6%) compared to secondary infection (48.2%). The specificity of NS1 assay for non-dengue samples were 100%. The NS1 gene sequence analysis of 48 isolates revealed the presence of polymorphisms of the NS1 genes which apparently did not influence the NS1 sensitivity. CONCLUSIONS: We observed a relatively low sensitivity of NS1 ELISA for dengue detection on RT-PCR-positive dengue samples. The detection rate increased significantly when NS1 data was combined with IgM. In our study, the low sensitivity of NS1 antigen detection did not relate to NS1 genetic diversity. Rather, the performance of the NS1 antigen test was affected by the infection status of patients and geographical origin of samples.

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas.

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

Lineage-informative microhaplotypes for spatio-temporal surveillance of Plasmodium vivax malaria parasites.

Challenges in understanding the origin of recurrent Plasmodium vivax infections constrains the surveillance of antimalarial efficacy and transmission of this neglected parasite. Recurrent infections within an individual may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or new inoculations (reinfection). Molecular inference of familial relatedness (identity-by-descent or IBD) based on whole genome sequence data, together with analysis of the intervals between parasitaemic episodes ("time-to-event" analysis), can help resolve the probable origin of recurrences. Whole genome sequencing of predominantly low-density P. vivax infections is challenging, so an accurate and scalable genotyping method to determine the origins of recurrent parasitaemia would be of significant benefit. We have developed a P. vivax genome-wide informatics pipeline to select specific microhaplotype panels that can capture IBD within small, amplifiable segments of the genome. Using a global set of 615 P. vivax genomes, we derived a panel of 100 microhaplotypes, each comprising 3-10 high frequency SNPs within <200 bp sequence windows. This panel exhibits high diversity in regions of the Asia-Pacific, Latin America and the horn of Africa (median HE = 0.70-0.81) and it captured 89% (273/307) of the polyclonal infections detected with genome-wide datasets. Using data simulations, we demonstrate lower error in estimating pairwise IBD using microhaplotypes, relative to traditional biallelic SNP barcodes. Our panel exhibited high accuracy in predicting the country of origin (median Matthew's correlation coefficient >0.9 in 90% countries tested) and it also captured local infection outbreak and bottlenecking events. The informatics pipeline is available open-source and yields microhaplotypes that can be readily transferred to high-throughput amplicon sequencing assays for surveillance in malaria-endemic regions.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia.

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Uncoupling protein 2 gene polymorphisms are associated with obesity.

BACKGROUND: Uncoupling protein 2 (UCP2) gene polymorphisms have been reported as genetic risk factors for obesity and type 2 diabetes mellitus (T2DM). We examined the association of commonly observed UCP2 G(-866)A (rs659366) and Ala55Val (C > T) (rs660339) single nucleotide polymorphisms (SNPs) with obesity, high fasting plasma glucose, and serum lipids in a Balinese population. METHODS: A total of 603 participants (278 urban and 325 rural subjects) were recruited from Bali Island, Indonesia. Fasting plasma glucose (FPG), triglyceride (TG), high density lipoprotein cholesterol (HDL-C), low density lipoprotein cholesterol (LDL-C) and total cholesterol (TC) were measured. Obesity was determined based on WHO classifications for adult Asians. Participants were genotyped for G(-866)A and Ala55Val polymorphisms of the UCP2 gene. RESULTS: Obesity prevalence was higher in urban subjects (51%) as compared to rural subjects (23%). The genotype, minor allele (MAF), and heterozygosity frequencies were similar between urban and rural subjects for both SNPs. All genotype frequencies were in Hardy-Weinberg equilibrium. A combined analysis of genotypes and environment revealed that the urban subjects carrying the A/A genotype of the G(-866)A SNP have higher BMI than the rural subjects with the same genotype. Since the two SNPs showed strong linkage disequilibrium (D' = 0.946, r2 = 0.657), a haplotype analysis was performed. We found that the AT haplotype was associated with high BMI only when the urban environment was taken into account. CONCLUSIONS: We have demonstrated the importance of environmental settings in studying the influence of the common UCP2 gene polymorphisms in the development of obesity in a Balinese population.

Genomic profiles of Indonesian colorectal cancer patients.

Background: Colorectal cancer (CRC) is one of the most commonly diagnosed cancers worldwide and genetic mutation plays a vital role in CRC development. A previous study has suggested that genetic alterations among Indonesian patients with CRC might differ from those known in developed countries. This study aimed to describe the genomic profiles of Indonesian patients with CRC. Methods: A total of 13 patients were recruited for this study from May to July 2019. Tissue samples were collected, and genomic DNA was extracted from the samples. AmpliSeq for Illumina Cancer HotSpot Panel v2 Next-generation sequencing was used for DNA sequencing and a genome analysis toolkit was used for local realignment around the discovered variants. Results: A total of 45 genes comprising 391 single nucleotide variants (SNVs) with a depth >10 were observed. The genes with the most variants were STK11, SMAD4, EGFR, and ERBB4 and the genes with the most non-synonymous variants were SMAD4, TP53, FGFR3, CDKN2A, and STK11. Genes and SNVs in at least 90% of all samples consisted of 43 genes comprising 286 variants. Genes with the most non-synonymous SNVs were EGFR, SMO, FGFR3, TP53, STK11, CDKN2A. Genes related to the chromosomal instability pathway, such as TP53, SMAD4, KRAS, and APC, are also found in the analysis. Conclusions: Our findings showed that all patients with CRC in this study had genetic mutations in the chromosomal instability pathway. Analysis of genetic mutation of Indonesian patients with CRC might be crucial for advanced targeted therapy and for better clinical outcomes.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales.

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications ('use cases') have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria.

Traditionally, patient travel history has been used to distinguish imported from autochthonous malaria cases, but the dormant liver stages of Plasmodium vivax confound this approach. Molecular tools offer an alternative method to identify, and map imported cases. Using machine learning approaches incorporating hierarchical fixation index and decision tree analyses applied to 799 P. vivax genomes from 21 countries, we identified 33-SNP, 50-SNP and 55-SNP barcodes (GEO33, GEO50 and GEO55), with high capacity to predict the infection's country of origin. The Matthews correlation coefficient (MCC) for an existing, commonly applied 38-SNP barcode (BR38) exceeded 0.80 in 62% countries. The GEO panels outperformed BR38, with median MCCs > 0.80 in 90% countries at GEO33, and 95% at GEO50 and GEO55. An online, open-access, likelihood-based classifier framework was established to support data analysis (vivaxGEN-geo). The SNP selection and classifier methods can be readily amended for other use cases to support malaria control programs.

Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia.

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.