A GPU-accelerated compute framework for pathogen genomic variant identification to aid genomic epidemiology of infectious disease: a malaria case study.

As recently demonstrated by the COVID-19 pandemic, large-scale pathogen genomic data are crucial to characterize transmission patterns of human infectious diseases. Yet, current methods to process raw sequence data into analysis-ready variants remain slow to scale, hampering rapid surveillance efforts and epidemiological investigations for disease control. Here, we introduce an accelerated, scalable, reproducible, and cost-effective framework for pathogen genomic variant identification and present an evaluation of its performance and accuracy across benchmark datasets of Plasmodium falciparum malaria genomes. We demonstrate superior performance of the GPU framework relative to standard pipelines with mean execution time and computational costs reduced by 27× and 4.6×, respectively, while delivering 99.9% accuracy at enhanced reproducibility.

Performance evaluation of RDT, light microscopy, and PET-PCR for detecting Plasmodium falciparum malaria infections in the 2018 Zambia National Malaria Indicator Survey.

BACKGROUND: Zambia continues to advance on the path to elimination with significant reductions in malaria morbidity and mortality. Crucial components that have contributed to progress thus far and are necessary for achieving the national malaria elimination goals include properly identifying and treating all malaria cases through accurate diagnosis. This study sought to compare and assess the diagnostic performance of Rapid Diagnostic Tests (RDT) and Light Microscopy (LM) with photo-induced electron transfer polymerase chain reaction (PET-PCR) as the gold standard using 2018 Malaria Indicator Survey (MIS) data across Zambia to better understand diagnostic accuracy metrics and how these vary across a transmission gradient. METHODS: Cross-sectional samples collected in a nationally representative survey from 7 provinces in Zambia were tested for the presence of malaria parasites by light microscopy (LM), rapid diagnostic test (RDT) and the gold standard PET-PCR. Diagnostic performance was assessed including sensitivity, specificity, negative- and positive-predictive values across a wide malaria transmission spectrum. Diagnostic accuracy metrics were measured, and statistically significant differences were calculated between test methods for different outcome variables. RESULTS: From the individuals included in the MIS, the overall prevalence of Plasmodium falciparum malaria was 32.9% by RDT, 19.4% by LM, and 23.2% by PET-PCR. Herein, RDT and LM diagnostic performance was compared against gold standard PET-PCR with LM displaying a higher diagnostic accuracy than RDTs (91.3% vs. 84.6% respectively) across the transmission spectrum in Zambia. However, the performance of both diagnostics was significantly reduced in low parasitaemia samples. Consistent with previous studies, RDT diagnostic accuracy was predominantly affected by a high rate of false positives. CONCLUSIONS: RDTs and LM both perform well across a range of transmission intensities within their respective target applications, i.e., in the community, for the former, where ease of use and speed of result is critical, and at the health facility, for the latter, where accuracy is prioritized. However, the performance of both diagnostic methods is adversely affected by low parasitaemia infections. As Zambia moves towards elimination more sensitive tools may be required to identify the last cases.

Genetic Evidence of Focal Plasmodium falciparum Transmission in a Pre-elimination Setting in Southern Province, Zambia.

BACKGROUND: Southern Province, Zambia has experienced a dramatic decline in Plasmodium falciparum malaria transmission in the past decade and is targeted for elimination. Zambia's National Malaria Elimination Program recommends reactive case detection (RCD) within 140 m of index households to enhance surveillance and eliminate remaining transmission foci. METHODS: To evaluate whether RCD captures local transmission, we genotyped 26 microsatellites from 106 samples collected from index (n = 27) and secondary (n = 79) cases detected through RCD in the Macha Hospital catchment area between January 2015 and April 2016. RESULTS: Participants from the same RCD event harbored more genetically related parasites than those from different RCD events, suggesting that RCD captures, at least in part, infections related through local transmission. Related parasites clustered in space and time, up to at least 250 m from index households. Spatial analysis identified a putative focal transmission hotspot. CONCLUSIONS: The current RCD strategy detects focal transmission events, although programmatic guidelines to screen within 140 m of index households may fail to capture all secondary cases. This study highlights the utility of parasite genetic data in assessing programmatic interventions, and similar approaches may be useful to malaria elimination programs seeking to tailor intervention strategies to the underlying transmission epidemiology.

High Plasmodium falciparum genetic diversity and temporal stability despite control efforts in high transmission settings along the international border between Zambia and the Democratic Republic of the Congo.

BACKGROUND: While the utility of parasite genotyping for malaria elimination has been extensively documented in low to moderate transmission settings, it has been less well-characterized in holoendemic regions. High malaria burden settings have received renewed attention acknowledging their critical role in malaria elimination. Defining the role for parasite genomics in driving these high burden settings towards elimination will enhance future control programme planning. METHODS: Amplicon deep sequencing was used to characterize parasite population genetic diversity at polymorphic Plasmodium falciparum loci, Pfama1 and Pfcsp, at two timepoints in June-July 2016 and January-March 2017 in a high transmission region along the international border between Luapula Province, Zambia and Haut-Katanga Province, the Democratic Republic of the Congo (DRC). RESULTS: High genetic diversity was observed across both seasons and in both countries. No evidence of population structure was observed between parasite populations on either side of the border, suggesting that this region may be one contiguous transmission zone. Despite a decline in parasite prevalence at the sampling locations in Haut-Katanga Province, no genetic signatures of a population bottleneck were detected, suggesting that larger declines in transmission may be required to reduce parasite genetic diversity. Analysing rare variants may be a suitable alternative approach for detecting epidemiologically important genetic signatures in highly diverse populations; however, the challenge is distinguishing true signals from potential artifacts introduced by small sample sizes. CONCLUSIONS: Continuing to explore and document the utility of various parasite genotyping approaches for understanding malaria transmission in holoendemic settings will be valuable to future control and elimination programmes, empowering evidence-based selection of tools and methods to address pertinent questions, thus enabling more efficient resource allocation.

Vectors as Epidemiological Sentinels: Patterns of Within-Tick Borrelia burgdorferi Diversity.

Hosts including humans, other vertebrates, and arthropods, are frequently infected with heterogeneous populations of pathogens. Within-host pathogen diversity has major implications for human health, epidemiology, and pathogen evolution. However, pathogen diversity within-hosts is difficult to characterize and little is known about the levels and sources of within-host diversity maintained in natural populations of disease vectors. Here, we examine genomic variation of the Lyme disease bacteria, Borrelia burgdorferi (Bb), in 98 individual field-collected tick vectors as a model for study of within-host processes. Deep population sequencing reveals extensive and previously undocumented levels of Bb variation: the majority (~70%) of ticks harbor mixed strain infections, which we define as levels Bb diversity pre-existing in a diverse inoculum. Within-tick diversity is thus a sample of the variation present within vertebrate hosts. Within individual ticks, we detect signatures of positive selection. Genes most commonly under positive selection across ticks include those involved in dissemination in vertebrate hosts and evasion of the vertebrate immune complement. By focusing on tick-borne Bb, we show that vectors can serve as epidemiological and evolutionary sentinels: within-vector pathogen diversity can be a useful and unbiased way to survey circulating pathogen diversity and identify evolutionary processes occurring in natural transmission cycles.

Distinct parasite populations infect individuals identified through passive and active case detection in a region of declining malaria transmission in southern Zambia.

BACKGROUND: Substantial reductions in the burden of malaria have been documented in parts of sub-Saharan Africa, with elimination strategies and goals being formulated in some regions. Within this context, understanding the epidemiology of low-level malaria transmission is crucial to achieving and sustaining elimination. A 24 single-nucleotide-polymorphism Plasmodium falciparum molecular barcode was used to characterize parasite populations from infected individuals identified through passive and active case detection in an area approaching malaria elimination in southern Zambia. METHODS: The study was conducted in the catchment area of Macha Hospital in Choma District, Southern Province, Zambia, where the parasite prevalence declined over the past decade, from 9.2% in 2008 to less than 1% in 2013. Parasite haplotypes from actively detected, P. falciparum-infected participants enrolled in a serial cross-sectional, community-based cohort study from 2008 to 2013 and from passively detected, P. falciparum-infected individuals enrolled at five rural health centres from 2012 to 2015 were compared. Changes in P. falciparum genetic relatedness, diversity and complexity were analysed as malaria transmission declined. RESULTS: Actively detected cases identified in the community were most commonly rapid diagnostic test negative, asymptomatic and had submicroscopic parasitaemia. Phylogenetic reconstruction using concatenated 24 SNP barcode revealed a separation of parasite haplotypes from passively and actively detected infections, consistent with two genetically distinct parasite populations. For passively detected infections identified at health centres, the proportion of detectable polyclonal infections was consistently low in all seasons, in contrast with actively detected infections in which the proportion of polyclonal infections was high. The mean genetic divergence for passively detected infections was 34.5% for the 2012-2013 transmission season, 37.8% for the 2013-2014 season, and 30.8% for the 2014-2015 season. The mean genetic divergence for actively detected infections was 22.3% in the 2008 season and 29.0% in the 2008-2009 season and 9.9% across the 2012-2014 seasons. CONCLUSIONS: Distinct parasite populations were identified among infected individuals identified through active and passive surveillance, suggesting that infected individuals detected through active surveillance may not have contributed substantially to ongoing transmission. As parasite prevalence and diversity within these individuals declined, resource-intensive efforts to identify the chronically infected reservoir may not be necessary to eliminate malaria in this setting.

Babesia microti from humans and ticks hold a genomic signature of strong population structure in the United States.

BACKGROUND: Babesia microti is an emerging tick-borne apicomplexan parasite with increasing geographic range and incidence in the United States. The rapid expansion of B. microti into its current distribution in the northeastern USA has been due to the range expansion of the tick vector, Ixodes scapularis, upon which the causative agent is dependent for transmission to humans. RESULTS: To reconstruct the history of B. microti in the continental USA and clarify the evolutionary origin of human strains, we used multiplexed hybrid capture of 25 B. microti isolates obtained from I. scapularis and human blood. Despite low genomic variation compared with other Apicomplexa, B. microti was strongly structured into three highly differentiated genetic clusters in the northeastern USA. Bayesian analyses of the apicoplast genomes suggest that the origin of the current diversity of B. microti in northeastern USA dates back 46 thousand years with a signature of recent population expansion in the last 1000 years. Human-derived samples belonged to two rarely intermixing clusters, raising the possibility of highly divergent infectious phenotypes in humans. CONCLUSIONS: Our results validate the multiplexed hybrid capture strategy for characterizing genome-wide diversity and relatedness of B. microti from ticks and humans. We find strong population structure in B. microti samples from the Northeast indicating potential barriers to gene flow.

Mitogenomes reveal diversity of the European Lyme borreliosis vector Ixodes ricinus in Italy.

In Europe, the Ixodes ricinus tick is the most important vector of the etiological agents of Lyme borreliosis and several other emerging tick-borne diseases. Because tick-borne pathogens are dependent on their vectors for transmission, understanding the vector population structure is crucial to inform public health research of pathogen dynamics and spread. However, the population structure and dynamics of this important vector species are not well understood as most genetic studies utilize short mitochondrial and nuclear sequences with little diversity. Herein we obtained and analyzed complete mitochondrial genome (hereafter "mitogenome") sequences to better understand the genetic diversity and the population structure of I. ricinus from two long-standing tick-borne disease foci in northern Italy. Complete mitogenomes of 23 I. ricinus ticks were sequenced at high coverage. Out of 23 mitogenome sequences we identified 17 unique haplotypes composed of 244 segregating sites. Phylogenetic reconstruction using 18 complete mitogenome sequences revealed the coexistence of four highly divergent I. ricinus maternal lineages despite the narrow spatial scale over which these samples were obtained (100km). Notably, the estimated coalescence time of the 18 mitogenome haplotypes is ∼427 thousand years ago (95% HPD 330, 540). This divergence between I. ricinus lineages is consistent with the mitochondrial diversity of other arthropod vector species and indicates that long-term I. ricinus populations may have been less structured and larger than previously thought. Thus, this study suggests that a rapid and accurate retrieval of full mitochondrial genomes from this disease vector enables fine-resolution studies of tick intraspecies genetic relationships, population differentiation, and demographic history.

Whole genome capture of vector-borne pathogens from mixed DNA samples: a case study of Borrelia burgdorferi.

BACKGROUND: Rapid and accurate retrieval of whole genome sequences of human pathogens from disease vectors or animal reservoirs will enable fine-resolution studies of pathogen epidemiological and evolutionary dynamics. However, next generation sequencing technologies have not yet been fully harnessed for the study of vector-borne and zoonotic pathogens, due to the difficulty of obtaining high-quality pathogen sequence data directly from field specimens with a high ratio of host to pathogen DNA. RESULTS: We addressed this challenge by using custom probes for multiplexed hybrid capture to enrich for and sequence 30 Borrelia burgdorferi genomes from field samples of its arthropod vector. Hybrid capture enabled sequencing of nearly the complete genome (~99.5 %) of the Borrelia burgdorferi pathogen with 132-fold coverage, and identification of up to 12,291 single nucleotide polymorphisms per genome. CONCLUSIONS: The proprosed culture-independent method enables efficient whole genome capture and sequencing of pathogens directly from arthropod vectors, thus making population genomic study of vector-borne and zoonotic infectious diseases economically feasible and scalable. Furthermore, given the similarities of invertebrate field specimens to other mixed DNA templates characterized by a high ratio of host to pathogen DNA, we discuss the potential applicabilty of hybrid capture for genomic study across diverse study systems.

Investigating the rise of Omicron variant through genomic surveillance of SARS-CoV-2 infections in a highly vaccinated university population.

Novel variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continue to emerge as the coronavirus disease 2019 (COVID-19) pandemic extends into its fourth year. Understanding SARS-CoV-2 circulation in university populations is vital for effective interventions in higher education settings and will inform public health policy during pandemics. In this study, we generated 793 whole-genome sequences collected over an entire academic year in a university population in Indiana, USA. We clearly captured the rapidity with which Delta variant was wholly replaced by Omicron variant across the West Lafayette campus over the length of two academic semesters in a community with high vaccination rates. This mirrored the emergence of Omicron throughout the state of Indiana and the USA. Further, phylogenetic analyses demonstrated that there was a more diverse set of potential geographic origins for Omicron viruses introduction into campus when compared to Delta. Lastly, statistics indicated that there was a more significant role for international and out-of-state migration in the establishment of Omicron variants at Purdue. This surveillance workflow, coupled with viral genomic sequencing and phylogeographic analyses, provided critical insights into SARS-CoV-2 transmission dynamics and variant arrival.

A metagenomic analysis of the phase 2 Anopheles gambiae 1000 genomes dataset reveals a wide diversity of cobionts associated with field collected mosquitoes.

The Anopheles gambiae 1000 Genomes (Ag1000G) Consortium previously utilized deep sequencing methods to catalogue genetic diversity across African An. gambiae populations. We analyzed the complete datasets of 1142 individually sequenced mosquitoes through Microsoft Premonition's Bayesian mixture model based (BMM) metagenomics pipeline. All specimens were confirmed as either An. gambiae sensu stricto (s.s.) or An. coluzzii with a high degree of confidence ( > 98% identity to reference). Homo sapiens DNA was identified in all specimens indicating contamination may have occurred either at the time of specimen collection, preparation and/or sequencing. We found evidence of vertebrate hosts in 162 specimens. 59 specimens contained validated Plasmodium falciparum reads. Human hepatitis B and primate erythroparvovirus-1 viral sequences were identified in fifteen and three mosquito specimens, respectively. 478 of the 1,142 specimens were found to contain bacterial reads and bacteriophage-related contigs were detected in 27 specimens. This analysis demonstrates the capacity of metagenomic approaches to elucidate important vector-host-pathogen interactions of epidemiological significance.

Diversity and selection analyses identify transmission-blocking antigens as the optimal vaccine candidates in Plasmodium falciparum.

Background: A highly effective vaccine for malaria remains an elusive target, at least in part due to the under-appreciated natural parasite variation. This study aimed to investigate genetic and structural variation, and immune selection of leading malaria vaccine candidates across the Plasmodium falciparum's life cycle. Methods: We analyzed 325 P. falciparum whole genome sequences from Zambia, in addition to 791 genomes from five other African countries available in the MalariaGEN Pf3k Rdatabase. Ten vaccine antigens spanning three life-history stages were examined for genetic and structural variations, using population genetics measures, haplotype network analysis, and 3D structure selection analysis. Findings: Among the ten antigens analyzed, only three in the transmission-blocking vaccine category display P. falciparum 3D7 as the dominant haplotype. The antigens AMA1, CSP, MSP119 and CelTOS, are much more diverse than the other antigens, and their epitope regions are under moderate to strong balancing selection. In contrast, Rh5, a blood stage antigen, displays low diversity yet slightly stronger immune selection in the merozoite-blocking epitope region. Except for CelTOS, the transmission-blocking antigens Pfs25, Pfs48/45, Pfs230, Pfs47, and Pfs28 exhibit minimal diversity and no immune selection in epitopes that induce strain-transcending antibodies, suggesting potential effectiveness of 3D7-based vaccines in blocking transmission. Interpretations: These findings offer valuable insights into the selection of optimal vaccine candidates against P. falciparum. Based on our results, we recommend prioritizing conserved merozoite antigens and transmission-blocking antigens. Combining these antigens in multi-stage approaches may be particularly promising for malaria vaccine development initiatives. Funding: Purdue Department of Biological Sciences; Puskas Memorial Fellowship; National Institute of Allergy and Infectious Diseases (U19AI089680).

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries.

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia.

BACKGROUND: Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries. METHODS: We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection. RESULTS: We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure. CONCLUSIONS: Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Malaria is caused by a parasite that is spread to humans via mosquito bites. It is a leading cause of death in children under five years old in sub-Saharan Africa. Analysis of the malaria parasite's complete set of DNA (its genome) can help us to understand transmission of the disease and how this changes in response to different strategies to control the disease. We analyzed the genomes of malaria parasites from children across Zambia. Our study revealed that 77% of children harbored multiple parasite strains, which suggests that local transmission (transmission between people within the same local area) is high. Genetic evidence for long-distance transmission was rarer. Furthermore, our findings suggest parasites are evolving in response to antimalarial drugs. Our study enhances our understanding of malaria dynamics in Zambia and may help to inform strategies for improved surveillance and control.

Diversity and selection analyses identify transmission-blocking antigens as the optimal vaccine candidates in Plasmodium falciparum.

BACKGROUND: A highly effective vaccine for malaria remains an elusive target, at least in part due to the under-appreciated natural parasite variation. This study aimed to investigate genetic and structural variation, and immune selection of leading malaria vaccine candidates across the Plasmodium falciparum's life cycle. METHODS: We analysed 325 P. falciparum whole genome sequences from Zambia, in addition to 791 genomes from five other African countries available in the MalariaGEN Pf3k Database. Ten vaccine antigens spanning three life-history stages were examined for genetic and structural variations, using population genetics measures, haplotype network analysis, and 3D structure selection analysis. FINDINGS: Among the ten antigens analysed, only three in the transmission-blocking vaccine category display P. falciparum 3D7 as the dominant haplotype. The antigens AMA1, CSP, MSP119 and CelTOS, are much more diverse than the other antigens, and their epitope regions are under moderate to strong balancing selection. In contrast, Rh5, a blood stage antigen, displays low diversity yet slightly stronger immune selection in the merozoite-blocking epitope region. Except for CelTOS, the transmission-blocking antigens Pfs25, Pfs48/45, Pfs230, Pfs47, and Pfs28 exhibit minimal diversity and no immune selection in epitopes that induce strain-transcending antibodies, suggesting potential effectiveness of 3D7-based vaccines in blocking transmission. INTERPRETATION: These findings offer valuable insights into the selection of optimal vaccine candidates against P. falciparum. Based on our results, we recommend prioritising conserved merozoite antigens and transmission-blocking antigens. Combining these antigens in multi-stage approaches may be particularly promising for malaria vaccine development initiatives. FUNDING: Purdue Department of Biological Sciences; Puskas Memorial Fellowship; National Institute of Allergy and Infectious Diseases (U19AI089680).

Genomic Surveillance of SARS-CoV-2 in a University Community: Insights Into Tracking Variants, Transmission, and Spread of Gamma (P.1) Variant.

Background: Using a combination of data from routine surveillance, genomic sequencing, and phylogeographic analysis, we tracked the spread and introduction events of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants focusing on a large university community. Methods: Here, we sequenced and analyzed 677 high-quality SARS-CoV-2 genomes from positive RNA samples collected from Purdue University students, faculty, and staff who tested positive for the virus between January 2021 and May 2021, comprising an average of 32% of weekly cases across the time frame. Results: Our analysis of circulating SARS-CoV-2 variants over time revealed periods when variants of concern (VOC) Alpha (B.1.1.7) and Iota (B.1.526) reached rapid dominance and documented that VOC Gamma (P.1) was increasing in frequency as campus surveillance was ending. Phylodynamic analysis of Gamma genomes from campus alongside a subsampling of >20 000 previously published P.1 genomes revealed 10 independent introductions of this variant into the Purdue community, predominantly from elsewhere in the United States, with introductions from within the state of Indiana and from Illinois, and possibly Washington and New York, suggesting a degree of domestic spread. Conclusions: We conclude that a robust and sustained active and passive surveillance program coupled with genomic sequencing during a pandemic offers important insights into the dynamics of pathogen arrival and spread in a campus community and can help guide mitigation measures.

A genotyping assay to determine geographic origin and transmission potential of Plasmodium falciparum malaria cases.

As countries work towards malaria elimination, it is important to monitor imported cases to prevent reestablishment of local transmission. The Plasmodium falciparum Pfs47 gene has strong geographic population structure, because only those parasites with Pfs47 haplotypes compatible with the mosquito vector species in a given continent are efficiently transmitted. Analysis of 4,971 world-wide Pfs47 sequences identified two SNPs (at 707 and 725 bp) as sufficient to establish the likely continent of origin of P. falciparum isolates. Pfs47 sequences from Africa, Asia, and the New World presented more that 99% frequency of distinct combinations of the SNPs 707 and 725 genotypes. Interestingly, Papua New Guinea Pfs47 sequences have the highest diversity in SNPs 707 and 725. Accurate and reproducible High-Resolution Melting (HRM) assays were developed to genotype Pfs47 SNPs 707 and 725 in laboratory and field samples, to assess the geographic origin and risk of local transmission of imported P. falciparum malaria cases.

Phylogenetic Complexity of Morphologically Identified Anopheles squamosus in Southern Zambia.

Despite dramatic reductions in malaria cases in the catchment area of Macha Hospital, Choma District, Southern Province in Zambia, prevalence has remained near 1-2% by RDT for the past several years. To investigate residual malaria transmission in the area, this study focuses on the relative abundance, foraging behavior, and phylogenetic relationships of Anopheles squamosus specimens. In 2011, higher than expected rates of anthropophily were observed among "zoophilic" An. squamosus, a species that had sporadically been found to contain Plasmodium falciparum sporozoites. The importance of An. squamosus in the region was reaffirmed in 2016 when P. falciparum sporozoites were detected in numerous An. squamosus specimens. This study analyzed Centers for Disease Control (CDC) light trap collections of adult mosquitoes from two collection schemes: one performed as part of a reactive-test-and-treat program and the second performed along a geographical transect. Morphological identification, molecular verification of anopheline species, and blood meal source were determined on individual samples. Data from these collections supported earlier studies demonstrating An. squamosus to be primarily exophagic and zoophilic, allowing them to evade current control measures. The phylogenetic relationships generated from the specimens in this study illustrate the existence of well supported clade structure among An. squamosus specimens, which further emphasizes the importance of molecular identification of vectors. The primarily exophagic behavior of An. squamosus in these collections also highlights that indoor vector control strategies will not be sufficient for elimination of malaria in southern Zambia.

Multiple Novel Clades of Anopheline Mosquitoes Caught Outdoors in Northern Zambia.

Residual vector populations that do not come in contact with the most frequently utilized indoor-directed interventions present major challenges to global malaria eradication. Many of these residual populations are mosquito species about which little is known. As part of a study to assess the threat of outdoor exposure to malaria mosquitoes within the Southern and Central Africa International Centers of Excellence for Malaria Research, foraging female anophelines were collected outside households in Nchelenge District, northern Zambia. These anophelines proved to be more diverse than had previously been reported in the area. In order to further characterize the anopheline species, sequencing and phylogenetic approaches were utilized. Anopheline mosquitoes were collected from outdoor light traps, morphologically identified, and sent to Johns Hopkins Bloomberg School of Public Health for sequencing. Sanger sequencing from 115 field-derived samples yielded mitochondrial COI sequences, which were aligned with a homologous 488 bp gene segment from known anophelines (n = 140) retrieved from NCBI. Nuclear ITS2 sequences (n = 57) for at least one individual from each unique COI clade were generated and compared against NCBI's nucleotide BLAST database to provide additional evidence for taxonomical identity and structure. Molecular and morphological data were combined for assignment of species or higher taxonomy. Twelve phylogenetic groups were characterized from the COI and ITS2 sequence data, including the primary vector species Anopheles funestus s.s. and An. gambiae s.s. An unexpectedly large proportion of the field collections were identified as An. coustani and An. sp. 6. Six phylogenetic groups remain unidentified to species-level. Outdoor collections of anopheline mosquitoes in areas frequented by people in Nchelenge, northern Zambia, proved to be extremely diverse. Morphological misidentification and underrepresentation of some anopheline species in sequence databases confound efforts to confirm identity of potential malaria vector species. The large number of unidentified anophelines could compromise the malaria vector surveillance and malaria control efforts not only in northern Zambia but other places where surveillance and control are focused on indoor-foraging and resting anophelines. Therefore, it is critical to continue development of methodologies that allow better identification of these populations and revisiting and cleaning current genomic databases.

The evolutionary dynamics of bluetongue virus.

Bluetongue virus (BTV) is a midge-borne member of the genus Orbivirus that causes an eponymous debilitating livestock disease of great agricultural impact and which has expanded into Europe in recent decades. Reassortment among the ten segments comprising the double-stranded (ds) RNA genome of BTV has played an important role in generating the epidemic strains of this virus in Europe. In this study, we investigated the dynamics of BTV genome segment evolution utilizing time-structured data sets of complete sequences from four segments, totalling 290 sequences largely sampled from ruminant hosts. Our analysis revealed that BTV genome segments generally evolve under strong purifying selection and at substitution rates that are generally lower (mean rates of approximately 0.5-7 x 10(-4) nucleotide substitutions per site, per year) than vector-borne positive-sense viruses with single-strand (ss) RNA genomes. These also represent the most robust estimates of the nucleotide substitution rate in a dsRNA virus generated to date. Additionally, we determined that patterns of geographic structure and times to most recent common ancestor differ substantially between each segment, including a relatively recent origin for the diversity of segment 10 within the past millennium. Together, these findings demonstrate the effect of reassortment to decouple the evolutionary dynamics of BTV genome segments.