Genomic Sequencing of Dengue Virus Strains Associated with Papua New Guinean Outbreaks in 2016 Reveals Endemic Circulation of DENV-1 and DENV-2

Over the past decade, the Pacific region has experienced many arboviral outbreaks, including dengue, chikungunya, and Zika viruses. Papua New Guinea (PNG) has a high burden of arboviral diseases, but there is a paucity of knowledge about the epidemiology and circulation of these viruses in the country. In this study, we report investigations into suspected arboviral outbreaks of febrile disease in PNG from December 2015 to June 2017. DENV-1 and DENV-2 were the mostly commonly detected viruses, and low circulation of DENV-3 and ZIKV was also detected. DENV-4 and CHIKV were not detected during this period. Full genome sequencing of selected positive samples revealed that circulation was dominated by endemic indigenous strains belonging to DENV-1 (genotype IV) and DENV-2 (genotype C) that have been present in the country for up to a decade. A DENV-2 sublineage was also identified that has been associated with outbreaks of severe dengue in both PNG and the Solomon Islands.

Diversity of Leptospira spp. in bats and rodents from Papua New Guinea

Leptospirosis is the most common bacterial zoonosis globally. The pathogen, Leptospira spp., is primarily associated with rodent reservoirs. However, a wide range of other species has been implicated as reservoirs or dead-end hosts. We conducted a survey for Leptospira spp. in bats and rodents from Papua New Guinea. Kidney samples were collected from 97 pteropodid bats (five species), 37 insectivorous bats from four different families (six species) and 188 rodents (two species). Leptospires were detected in a high proportion of pteropodid bats, including Nyctimene cf. albiventer (35%), Macroglossus minimus (34%) and Rousettus amplexicaudatus (36%). Partial sequencing of the secY gene from rodent and bat leptospires showed host species clustering, with Leptospira interrogans and L. weilii detected in rodents and L. kirschneri and a potential novel species of Leptospira detected in bats. Further research is needed in Papua New Guinea and other locales in the Pacific region to gain a better understanding of the circulation dynamics of leptospires in reservoir species and the risks to public and veterinary health.

A novel function for the CoQ10 biosynthetic complex in anti-pneumococcal macrophage function

Despite the recent introduction of pneumococcal polysaccharide and conjugate vaccines, Streptococcus pneumoniae infection remains a leading cause of illness and death worldwide. In particular, infants in Papua New Guinea are at increased risk of severe pneumococcal pneumonia compared to infants in similar countries. We sought to determine if a novel genetic variant could explain this increased susceptibility. Whole exome sequencing revealed a single nucleotide variant (D308Y) in the gene encoding COQ6 (COQ6DY), a monooxygenase required for CoQ10 biosynthesis. We have utilized both a Saccharomyces cerevisiae model and a mouse model of COQ6DY to show that despite adequate production of CoQ10, this variant directly causes increased susceptibility to S. pneumoniae. This variant represents a previously unknown function of the CoQ10 biosynthetic complex that does not exert its effects through CoQ10 deficiency but rather through alterations of mitochondrial function and metabolism. These mitochondrial deficits in COQ6DY macrophages are sufficient to abrogate macrophage killing of S. pneumoniae and alter the coordination of the downstream immune response. In conclusion, we have identified a novel susceptibility allele to S. pneumoniae infection that exerts its effects via alterations in macrophage mitochondrial function.

Supported by NIH (R21 AI142723), SLCH Children's Discovery Institute (PD-II-2018-742)

A novel variant in CoQ biosynthesis highly prevalent in Papua New Guinea children increases mortality following bacterial pneumonia

To identify immune variants predisposing to severe pneumonia, we performed whole exome sequencing in a pediatric population highly susceptible to acute lower respiratory infections, identifying a candidate novel variant in the CoQ biosynthetic pathway. To evaluate the effect of this variant on immune function during bacterial pneumonia, we generated a mouse line using CRISPR-Cas9 that expresses the homologous variant in the enzyme COQ6. Interestingly, we found that the variant does not result in CoQ deficiency, as other known variants in biosynthetic proteins do, however intra-tracheal S. pneumoniae infection leads to increased bacteremia and mortality in mutant mice. Mechanistic studies show that mutant macrophages have reduced pneumococcal killing in vitro, showing an intrinsic defect in innate immune function conferred by the COQ6 mutation. Variant macrophages have decreased mitochondrial respiratory capacity both at baseline and following stimulation with LPS, as well as an inability to induce mitochondrial reactive oxygen species (mROS) in response to stimulation despite increased mROS at baseline. Thus, the novel variant in a CoQ biosynthetic enzyme leads to changes in macrophage mitochondrial function and an intrinsic inability to kill internalized bacteria. As alveolar macrophages are the first responders in the lung to bacterial challenge, the inability of these macrophages to mount a sufficient immune response leads to the observed increase in mortality following bacterial pneumonia. This work describes a novel susceptibility locus to severe childhood pneumonia, and also represents the first known pathogenic variant in a CoQ biosynthetic protein that does not cause pathology resulting from CoQ deficiency.

Trachoma, Anti-Pgp3 Serology, and Ocular Chlamydia trachomatis Infection in Papua New Guinea

BACKGROUND

In Melanesia, the prevalence of trachomatous inflammation-follicular (TF) suggests that public health-level interventions against active trachoma are needed. However, the prevalence of trachomatous trichiasis is below the threshold for elimination as a public health problem and evidence of conjunctival infection with trachoma's causative organism (Chlamydia trachomatis [CT]) is rare. Here, we examine the prevalence of ocular infection with CT and previous exposure to CT in three evaluation units (EUs) of Papua New Guinea.

METHODS

All individuals aged 1-9 years who were examined for clinical signs of trachoma in 3 Global Trachoma Mapping Project EUs were eligible to take part in this study (N = 3181). Conjunctival swabs were collected from 349 children with TF and tested by polymerase chain reaction to assess for ocular CT infection. Dried blood spots were collected from 2572 children and tested for anti-Pgp3 antibodies using a multiplex assay.

RESULTS

The proportion of children with TF who had CT infection was low across all 3 EUs (overall 2%). Anti-Pgp3 seroprevalence was 5.2% overall and there was no association between anti-Pgp3 antibody level and presence of TF. In 2 EUs, age-specific seroprevalence did not increase significantly with increasing age in the 1- to 9-year-old population. In the third EU, there was a statistically significant change with age but the overall seroprevalence and peak age-specific seroprevalence was very low.

CONCLUSIONS

Based on these results, together with similar findings from the Solomon Islands and Vanuatu, the use of TF to guide antibiotic mass drug administration decisions in Melanesia should be reviewed.

A novel variant in ubiquinone biosynthesis highly prevalent in Papua New Guinea children increases mortality following bacterial pneumonia

To identify immune variants predisposing to severe pneumonia, we performed whole exome sequencing in a pediatric population highly susceptible to acute lower respiratory infections, identifying a candidate novel variant in the ubiquinone (CoQ10) biosynthetic pathway. To evaluate the effect of this variant on immune function during bacterial pneumonia, we generated a mouse line using CRISPR-Cas9 that expresses the homologous aspartate to tyrosine variant in the enzyme COQ6. Intra-tracheal S. pneumoniae infection leads to increased bacteremia and mortality in mice homozygous for the variant despite similar numbers of immune cells in the lung. Mechanistic studies show that macrophages expressing the variant have decreased mitochondrial activity at the ubiquinone-dependent reduction of cytochrome c by complex III, as well as decreased maximum respiratory capacity in response to acute stimulation. Variant-expressing macrophages also exhibit impaired generation of mitochondrial reactive oxygen species (mROS) causing a direct, intrinsic defect in intracellular killing of internalized bacteria. Thus, the novel variant in CoQ10 biosynthesis leads to changes in macrophage mitochondria and an intrinsic inability to kill internalized bacteria. As alveolar macrophages are the first responders in the lung to bacterial challenge, the inability of these macrophages to mount a sufficient immune response can explain the observed increase in mortality following bacterial pneumonia. Because variants in CoQ10 biosynthesis can be supplemented with CoQ10, a readily available therapy may be able to correct this defect and improve survival in children with this variant

Significant geographical differences in prevalence of mutations associated with Plasmodium falciparum and Plasmodium vivax drug resistance in two regions from Papua New Guinea

Background

Drug resistance remains a major obstacle to malaria treatment and control. It can arise and spread rapidly, and vary substantially even at sub-national level. National malaria programmes require cost-effective and timely ways of characterizing drug-resistance at multiple sites within their countries.

Methods

An improved multiplexed post-PCR ligase detection reaction—fluorescent microsphere assay (LDR-FMA) was used to simultaneously determine the presence of mutations in chloroquine resistance transporter (crt), multidrug resistance 1 (mdr1), dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) genes in Plasmodium falciparum (n = 727) and Plasmodium vivax (n = 574) isolates collected in 2006 from cross-sectional community population surveys in two geographically distinct regions (Madang and East Sepik) of Papua New Guinea (PNG) where strong regional differences in in vivo aminoquinoline and antifolate therapeutic efficacy had previously been observed. Data were compared to those of a follow-up survey conducted in 2010.

Results

Despite some very low parasite densities, the assay successfully amplified all P. falciparum and P. vivax loci in 77 and 69 % of samples, respectively. In 2006, prevalences of pfdhfr (59R-108 N) double mutation/wild type pfdhps haplotype, pfcrt SVMNT haplotype (72S-76T double mutation), and 86Y pfmdr1 mutation all exceeded 90 %. For P. vivax, 65 % carried at least two pvdhfr mutations, 97 % the 647P pvdhps mutation and 54 % the 976F pvmdr1 mutation. Prevalence of mutant haplotypes was higher in Madang than East Sepik for pfcrt SVMNT (97.4 vs 83.3 %, p = 0.001), pfdhfr (59R-108 N) (100 vs 90.6 %, p = 0.001), pvdhfr haplotypes (75.8 vs 47.6 %, p = 0.001) and pvmdr1 976F (71.2 vs 26.2 %, p < 0.001). Data from a subsequent Madang survey in 2010 showed that the prevalence of pfdhps mutations increased significantly from <5 % to >30 % (p < 0.001) as did the prevalence of pvdhfr mutant haplotypes (from 75.8 to 97.4 %, p = 0.012).

Conclusions

This LDR-FMA multiplex platform shows feasibility for low-cost, high-throughput, rapid characterization of a broad range of drug-resistance markers in low parasitaemia infections. Significant geographical differences in mutation prevalence correlate with previous genotyping surveys and in vivo trials and may reflect variable drug pressure and differences in health-care access in these two PNG populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0879-9) contains supplementary material, which is available to authorized users.

Temporal changes in Plasmodium falciparum anti-malarial drug sensitivity in vitro and resistance-associated genetic mutations in isolates from Papua New Guinea

Background

In northern Papua New Guinea (PNG), most Plasmodium falciparum isolates proved resistant to chloroquine (CQ) in vitro between 2005 and 2007, and there was near-fixation of pfcrt K76T, pfdhfr C59R/S108N and pfmdr1 N86Y. To determine whether the subsequent introduction of artemisinin combination therapy (ACT) and reduced CQ-sulphadoxine-pyrimethamine pressure had attenuated parasite drug susceptibility and resistance-associated mutations, these parameters were re-assessed between 2011 and 2013.

Methods

A validated fluorescence-based assay was used to assess growth inhibition of 52 P. falciparum isolates from children in a clinical trial in Madang Province. Responses to CQ, lumefantrine, piperaquine, naphthoquine, pyronaridine, artesunate, dihydroartemisinin, artemether were assessed. Molecular resistance markers were detected using a multiplex PCR ligase detection reaction fluorescent microsphere assay.

Results

CQ resistance (in vitro concentration required for 50% parasite growth inhibition (IC₅₀) >100 nM) was present in 19% of isolates. All piperaquine and naphthoquine IC₅₀s were <100 nM and those for lumefantrine, pyronaridine and the artemisinin derivatives were in low nM ranges. Factor analysis of IC₅₀s showed three groupings (lumefantrine; CQ, piperaquine, naphthoquine; pyronaridine, dihydroartemisinin, artemether, artesunate). Most isolates (96%) were monoclonal pfcrt K76T (SVMNT) mutants and most (86%) contained pfmdr1 N86Y (YYSND). No wild-type pfdhfr was found but most isolates contained wild-type (SAKAA) pfdhps. Compared with 2005–2007, the geometric mean (95% CI) CQ IC₅₀ was lower (87 (71–107) vs 167 (141–197) nM) and there had been no change in the prevalence of pfcrt K76T or pfmdr1 mutations. There were fewer isolates of the pfdhps (SAKAA) wild-type (60 vs 100%) and pfdhfr mutations persisted.

Conclusions

Reflecting less drug pressure, in vitro CQ sensitivity appears to be improving in Madang Province despite continued near-fixation of pfcrt K76T and pfmdr1 mutations. Temporal changes in IC₅₀s for other anti-malarial drugs were inconsistent but susceptibility was preserved. Retention or increases in pfdhfr and pfdhps mutations reflect continued use of sulphadoxine-pyrimethamine in the study area including through paediatric intermittent preventive treatment. The susceptibility of local isolates to lumefantrine may be unrelated to those of other ACT partner drugs.

Trial registration

Australian New Zealand Clinical Trials Registry ACTRN12610000913077.

How molecular epidemiology studies can support the National Malaria Control Program in Papua New Guinea

Papua New Guinea (PNG) is undertaking intensified efforts to control malaria. The National Malaria Control Program aims to reduce the burden of disease by large-scale distribution of insecticide-treated bednets, improved diagnosis and implementation of new treatments. A scientific program monitoring the effect of these interventions, including molecular epidemiology studies, closely accompanies the program. Laboratory assays have been developed in (or transferred to) PNG to measure prevalence of infection and intensity of transmission as well as potential resistance to currently used drugs. These assays help to assess the impact of the National Malaria Control Program, and they reveal a much clearer picture of malaria epidemiology in PNG. In addition, analysis of the geographical clustering of parasites aids in selecting areas where intensified control will be most successful. This paper gives an overview of current research and recently completed studies in the molecular epidemiology of malaria conducted in Papua New Guinea.

Naturally acquired immune responses to P. vivax merozoite surface protein 3α and merozoite surface protein 9 are associated with reduced risk of P. vivax malaria in young Papua New Guinean children

Background

Plasmodium vivax is the most geographically widespread human malaria parasite. Cohort studies in Papua New Guinea have identified a rapid onset of immunity against vivax-malaria in children living in highly endemic areas. Although numerous P. vivax merozoite antigens are targets of naturally acquired antibodies, the role of many of these antibodies in protective immunity is yet unknown.

Methodology/Principal Findings

In a cohort of children aged 1–3 years, antibodies to different regions of Merozoite Surface Protein 3α (PvMSP3α) and Merozoite Surface Protein 9 (PvMSP9) were measured and related to prospective risk of P. vivax malaria during 16 months of active follow-up. Overall, there was a low prevalence of antibodies to PvMSP3α and PvMSP9 proteins (9–65%). Antibodies to the PvMSP3α N-terminal, Block I and Block II regions increased significantly with age while antibodies to the PvMSP3α Block I and PvMSP9 N-terminal regions were positively associated with concurrent P. vivax infection. Independent of exposure (defined as the number of genetically distinct blood-stage infection acquired over time (_molFOB)) and age, antibodies specific to both PvMSP3α Block II (adjusted incidence ratio (aIRR) = 0.59, p = 0.011) and PvMSP9 N-terminus (aIRR = 0.68, p = 0.035) were associated with protection against clinical P. vivax malaria. This protection was most pronounced against high-density infections. For PvMSP3α Block II, the effect was stronger with higher levels of antibodies.

Conclusions

These results indicate that PvMSP3α Block II and PvMSP9 N-terminus should be further investigated for their potential as P. vivax vaccine antigens. Controlling for _molFOB assures that the observed associations are not confounded by individual differences in exposure.

Plasmodium vivax is the most geographically widespread human malaria parasite. In highly endemic areas such as Papua New Guinea, a very rapid onset of immunity against vivax-malaria is observed. Although it is known that numerous P. vivax merozoite antigens are targets of naturally acquired antibodies, the role of many of these antibodies in protective immunity is yet unknown. In a cohort of 183 children aged 1–3 years, we now show that the presence of antibodies to Merozoite Surface Protein 3α (PvMSP3α) and Merozoite Surface Protein 9 (PvMSP9) are associated with a significant reduction in the burden P. vivax malaria. Antibodies increased with age and in the presence of concurrent P. vivax infections. After adjusting for both age and individual differences in exposure, the strongest reductions in risk were seen in children with antibodies to PvMSP3α Block II (41% reduction, p = 0.001) and PvMSP9 N-terminal region. (32% reduction, p = 0.035). These results indicate that PvMSP3α Block II and PvMSP9 N-terminus should be further investigated for their potential as P. vivax vaccine antigens.

Strategies for detection of Plasmodium species gametocytes

Carriage and density of gametocytes, the transmission stages of malaria parasites, are determined for predicting the infectiousness of humans to mosquitoes. This measure is used for evaluating interventions that aim at reducing malaria transmission. Gametocytes need to be detected by amplification of stage-specific transcripts, which requires RNA-preserving blood sampling. For simultaneous, highly sensitive quantification of both, blood stages and gametocytes, we have compared and optimized different strategies for field and laboratory procedures in a cross sectional survey in 315 5-9 yr old children from Papua New Guinea. qRT-PCR was performed for gametocyte markers pfs25 and pvs25, Plasmodium species prevalence was determined by targeting both, 18S rRNA genes and transcripts. RNA-based parasite detection resulted in a P. falciparum positivity of 24.1%; of these 40.8% carried gametocytes. P. vivax positivity was 38.4%, with 38.0% of these carrying gametocytes. Sensitivity of DNA-based parasite detection was substantially lower with 14.1% for P. falciparum and 19.6% for P. vivax. Using the lower DNA-based prevalence of asexual stages as a denominator increased the percentage of gametocyte-positive infections to 59.1% for P. falciparum and 52.4% for P. vivax. For studies requiring highly sensitive and simultaneous quantification of sexual and asexual parasite stages, 18S rRNA transcript-based detection saves efforts and costs. RNA-based positivity is considerably higher than other methods. On the other hand, DNA-based parasite quantification is robust and permits comparison with other globally generated molecular prevalence data. Molecular monitoring of low density asexual and sexual parasitaemia will support the evaluation of effects of up-scaled antimalarial intervention programs and can also inform about small scale spatial variability in transmission intensity.