Combining genomics and epidemiology to track mumps virus transmission in the United States.

Unusually large outbreaks of mumps across the United States in 2016 and 2017 raised questions about the extent of mumps circulation and the relationship between these and prior outbreaks. We paired epidemiological data from public health investigations with analysis of mumps virus whole genome sequences from 201 infected individuals, focusing on Massachusetts university communities. Our analysis suggests continuous, undetected circulation of mumps locally and nationally, including multiple independent introductions into Massachusetts and into individual communities. Despite the presence of these multiple mumps virus lineages, the genomic data show that one lineage has dominated in the US since at least 2006. Widespread transmission was surprising given high vaccination rates, but we found no genetic evidence that variants arising during this outbreak contributed to vaccine escape. Viral genomic data allowed us to reconstruct mumps transmission links not evident from epidemiological data or standard single-gene surveillance efforts and also revealed connections between apparently unrelated mumps outbreaks

Direct Identification of Hundreds of Expression-Modulating Variants using a Multiplexed Reporter Assay

Although studies have identified hundreds of loci associated with human traits and diseases, pinpointing causal alleles remains difficult, particularly for non-coding variants. To address this challenge, we adapted the massively parallel reporter assay (MPRA) to identify variants that directly modulate gene expression. We applied it to 32,373 variants from 3,642 cis-expression quantitative trait loci and control regions. Detection by MPRA was strongly correlated with measures of regulatory function. We demonstrate MPRA’s capabilities for pinpointing causal alleles, using it to identify 842 variants showing differential expression between alleles, including 53 well-annotated variants associated with diseases and traits. We investigated one in detail, a risk allele for ankylosing spondylitis, and provide direct evidence of a non-coding variant that alters expression of the prostaglandin EP4 receptor. These results create a resource of concrete leads and illustrate the promise of this approach for comprehensively interrogating how non-coding polymorphism shapes human biology.National Institutes of Health (U.S.) (grant DP2OD006514)National Institutes of Health (U.S.) (grant K99HG0081)National Institutes of Health (U.S.) (grant R01HG006785

The Case for Selection at CCR5-Δ32

The C-C chemokine receptor 5, 32 base-pair deletion (CCR5-Δ32) allele confers strong resistance to infection by the AIDS virus HIV. Previous studies have suggested that CCR5-Δ32 arose within the past 1,000 y and rose to its present high frequency (5%–14%) in Europe as a result of strong positive selection, perhaps by such selective agents as the bubonic plague or smallpox during the Middle Ages. This hypothesis was based on several lines of evidence, including the absence of the allele outside of Europe and long-range linkage disequilibrium at the locus. We reevaluated this evidence with the benefit of much denser genetic maps and extensive control data. We find that the pattern of genetic variation at CCR5-Δ32 does not stand out as exceptional relative to other loci across the genome. Moreover using newer genetic maps, we estimated that the CCR5-Δ32 allele is likely to have arisen more than 5,000 y ago. While such results can not rule out the possibility that some selection may have occurred at C-C chemokine receptor 5 (CCR5), they imply that the pattern of genetic variation seen atCCR5-Δ32 is consistent with neutral evolution. More broadly, the results have general implications for the design of future studies to detect the signs of positive selection in the human genome

Reporter Assays for Ebola Virus Nucleoprotein Oligomerization, Virion-Like Particle Budding, and Minigenome Activity Reveal the Importance of Nucleoprotein Amino Acid Position 111

For highly pathogenic viruses, reporter assays that can be rapidly performed are critically needed to identify potentially functional mutations for further study under maximal containment (e.g., biosafety level 4 [BSL-4]). The Ebola virus nucleoprotein (NP) plays multiple essential roles during the viral life cycle, yet few tools exist to study the protein under BSL-2 or equivalent containment. Therefore, we adapted reporter assays to measure NP oligomerization and virion-like particle (VLP) production in live cells and further measured transcription and replication using established minigenome assays. As a proof-of-concept, we examined the NP-R111C substitution, which emerged during the 20132016 Western African Ebola virus disease epidemic and rose to high frequency. NP-R111C slightly increased NP oligomerization and VLP budding but slightly decreased transcription and replication. By contrast, a synthetic charge-reversal mutant, NP-R111E, greatly increased oligomerization but abrogated transcription and replication. These results are intriguing in light of recent structures of NP oligomers, which reveal that the neighboring residue, K110, forms a salt bridge with E349 on adjacent NP molecules. By developing and utilizing multiple reporter assays, we find that the NP-111 position mediates a complex interplay between NP\u27s roles in protein structure, virion budding, and transcription and replication

Sequence-Based Association and Selection Scans Identify Drug Resistance Loci in the Plasmodium Falciparum Malaria Parasite

Through rapid genetic adaptation and natural selection, the Plasmodium falciparum parasite—the deadliest of those that cause malaria—is able to develop resistance to antimalarial drugs, thwarting present efforts to control it. Genome-wide association studies (GWAS) provide a critical hypothesis-generating tool for understanding how this occurs. However, in P. falciparum, the limited amount of linkage disequilibrium hinders the power of traditional array-based GWAS. Here, we demonstrate the feasibility and power improvements gained by using whole-genome sequencing for association studies. We analyzed data from 45 Senegalese parasites and identified genetic changes associated with the parasites’ in vitro response to 12 different antimalarials. To further increase statistical power, we adapted a common test for natural selection, XP-EHH (cross-population extended haplotype homozygosity), and used it to identify genomic regions associated with resistance to drugs. Using this sequence-based approach and the combination of association and selection-based tests, we detected several loci associated with drug resistance. These loci included the previously known signals at pfcrt, dhfr, and pfmdr1, as well as many genes not previously implicated in drug-resistance roles, including genes in the ubiquitination pathway. Based on the success of the analysis presented in this study, and on the demonstrated shortcomings of array-based approaches, we argue for a complete transition to sequence-based GWAS for small, low linkage-disequilibrium genomes like that of P. falciparum.Molecular and Cellular BiologyOrganismic and Evolutionary Biolog

Metagenomic Sequencing of HIV-1 in the Blood and Female Genital Tract Reveals Little Quasispecies Diversity during Acute Infection.

Heterosexual transmission of human immunodeficiency virus type 1 (HIV-1) is associated with a significant bottleneck in the viral quasispecies population, yet the timing of that bottleneck is poorly understood. We characterized HIV-1 diversity in the blood and female genital tract (FGT) within 2 weeks after detection of infection in three women enrolled in a unique prospective cohort in South Africa. We assembled full-length HIV-1 genomes from matched cervicovaginal lavage (CVL) samples and plasma. Deep sequencing allowed us to identify intrahost single-nucleotide variants (iSNVs) and to characterize within-sample HIV-1 diversity. Our results demonstrated very little HIV-1 diversity in the FGT and plasma by the time viremia was detectable. Within each subject, the consensus HIV-1 sequences were identical in plasma and CVL fluid. No iSNV was present at >6% frequency. One subject had 77 low-frequency iSNVs across both CVL fluid and plasma, another subject had 14 iSNVs in only CVL fluid from the earliest time point, and the third subject had no iSNVs in CVL fluid or plasma. Overall, the small amount of diversity that we detected was greater in the FGT than in plasma and declined over the first 2 weeks after viremia was detectable, compatible with a very early HIV-1 transmission bottleneck. To our knowledge, our study represents the earliest genomic analysis of HIV-1 in the FGT after transmission. Further, the use of metagenomic sequencing allowed us to characterize other organisms in the FGT, including commensal bacteria and sexually transmitted infections, highlighting the utility of the method to sequence both HIV-1 and its metagenomic environment.IMPORTANCE Due to error-prone replication, HIV-1 generates a diverse population of viruses within a chronically infected individual. When HIV-1 is transmitted to a new individual, one or a few viruses establish the new infection, leading to a genetic bottleneck in the virus population. Understanding the timing and nature of this bottleneck may provide insight into HIV-1 vaccine design and other preventative strategies. We examined the HIV-1 population in three women enrolled in a unique prospective cohort in South Africa who were followed closely during the earliest stages of HIV-1 infection. We found very little HIV-1 diversity in the blood and female genital tract during the first 2 weeks after virus was detected in the bloodstream. These results are compatible with a very early HIV-1 population bottleneck, suggesting the need to study the HIV-1 population in the female genital tract before virus is detectable in the bloodstream

Biomarkers to Distinguish Bacterial From Viral Pediatric Clinical Pneumonia in a Malaria-Endemic Setting.

BACKGROUND: Differential etiologies of pediatric acute febrile respiratory illness pose challenges for all populations globally, but especially in malaria-endemic settings because the pathogens responsible overlap in clinical presentation and frequently occur together. Rapid identification of bacterial pneumonia with high-quality diagnostic tools would enable appropriate, point-of-care antibiotic treatment. Current diagnostics are insufficient, and the discovery and development of new tools is needed. We report a unique biomarker signature identified in blood samples to accomplish this. METHODS: Blood samples from 195 pediatric Mozambican patients with clinical pneumonia were analyzed with an aptamer-based, high-dynamic-range, quantitative assay (~1200 proteins). We identified new biomarkers using a training set of samples from patients with established bacterial, viral, or malarial pneumonia. Proteins with significantly variable abundance across etiologies (false discovery rate 90% sensitivity and >80% specificity, regardless of number of pathogen classes. Bacterial pneumonia was strongly associated with neutrophil markers-in particular, degranulation including HP, LCN2, LTF, MPO, MMP8, PGLYRP1, RETN, SERPINA1, S100A9, and SLPI. CONCLUSIONS: Blood protein signatures highly associated with neutrophil biology reliably differentiated bacterial pneumonia from other causes. With appropriate technology, these markers could provide the basis for a rapid diagnostic for field-based triage for antibiotic treatment of pediatric pneumonia

Identification and Functional Validation of the Novel Antimalarial Resistance Locus PF10_0355 in Plasmodium falciparum

The Plasmodium falciparum parasite's ability to adapt to environmental pressures, such as the human immune system and antimalarial drugs, makes malaria an enduring burden to public health. Understanding the genetic basis of these adaptations is critical to intervening successfully against malaria. To that end, we created a high-density genotyping array that assays over 17,000 single nucleotide polymorphisms (~1 SNP/kb), and applied it to 57 culture-adapted parasites from three continents. We characterized genome-wide genetic diversity within and between populations and identified numerous loci with signals of natural selection, suggesting their role in recent adaptation. In addition, we performed a genome-wide association study (GWAS), searching for loci correlated with resistance to thirteen antimalarials; we detected both known and novel resistance loci, including a new halofantrine resistance locus, PF10_0355. Through functional testing we demonstrated that PF10_0355 overexpression decreases sensitivity to halofantrine, mefloquine, and lumefantrine, but not to structurally unrelated antimalarials, and that increased gene copy number mediates resistance. Our GWAS and follow-on functional validation demonstrate the potential of genome-wide studies to elucidate functionally important loci in the malaria parasite genome.Bill & Melinda Gates FoundationEllison Medical FoundationExxon Mobil FoundationFogarty International CenterNational Institute of Allergy and Infectious Diseases (U.S.)Burroughs Wellcome FundDavid & Lucile Packard FoundationNational Science Foundation (U.S.). Graduate Research Fellowship Progra

Genome-wide detection and characterization of positive selection in human populations

With the advent of dense maps of human genetic variation, it is now possible to detect positive natural selection across the human genome. Here we report an analysis of over 3 million polymorphisms from the International HapMap Project Phase 2 (HapMap2). We used ‘long-range haplotype’ methods, which were developed to identify alleles segregating in a population that have undergone recent selection, and we also developed new methods that are based on cross-population comparisons to discover alleles that have swept to near-fixation within a population. The analysis reveals more than 300 strong candidate regions. Focusing on the strongest 22 regions, we develop a heuristic for scrutinizing these regions to identify candidate targets of selection. In a complementary analysis, we identify 26 non-synonymous, coding, single nucleotide polymorphisms showing regional evidence of positive selection. Examination of these candidates highlights three cases in which two genes in a common biological process have apparently undergone positive selection in the same population: LARGE and DMD, both related to infection by the Lassa virus3, in West Africa; SLC24A5 and SLC45A2, both involved in skin pigmentation in Europe; and EDAR and EDA2R, both involved in development of hair follicles, in Asia

Plasmepsin II–III copy number accounts for bimodal piperaquine resistance among Cambodian Plasmodium falciparum

Multidrug resistant Plasmodium falciparum in Southeast Asia endangers regional malaria elimination and threatens to spread to other malaria endemic areas. Understanding mechanisms of piperaquine (PPQ) resistance is crucial for tracking its emergence and spread, and to develop effective strategies for overcoming it. Here we analyze a mechanism of PPQ resistance in Cambodian parasites. Isolates exhibit a bimodal dose–response curve when exposed to PPQ, with the area under the curve quantifying their survival in vitro. Increased copy number for plasmepsin II and plasmepsin III appears to explain enhanced survival when exposed to PPQ in most, but not all cases. A panel of isogenic subclones reinforces the importance of plasmepsin II–III copy number to enhanced PPQ survival. We conjecture that factors producing increased parasite survival under PPQ exposure in vitro may drive clinical PPQ failures in the field