Identifying recombinant sequences in an era of large genomic databases is challenging as it requires an efficient algorithm to identify candidate recombinants and parents, as well as appropriate statistical methods to correct for the large number of comparisons performed. In 2007, a computation was introduced for an exact nonparametric mosaicism statistic that gave high-precision P values for putative recombinants. This exact computation meant that multiple-comparisons corrected P values also had high precision, which is crucial when performing millions or billions of tests in large databases. Here, we introduce an improvement to the algorithmic complexity of this computation from O(mn3) to O(mn2), where m and n are the numbers of recombination-informative sites in the candidate recombinant. This new computation allows for recombination analysis to be performed in alignments with thousands of polymorphic sites. Benchmark runs are presented on viral genome sequence alignments, new features are introduced, and applications outside recombination analysis are discussed.

SummaryBackgroundArtemisinin combination therapies (ACTs) are used worldwide as first-line treatment against confirmed or suspected Plasmodium falciparum malaria. Despite the success of ACTs at reducing the global burden of malaria, emerging resistance to artemisinin threatens these gains. Countering onset of resistance might need deliberate tactics aimed at slowing the reduction in ACT effectiveness. We assessed optimum use of ACTs at the population level, specifically focusing on a strategy of multiple first-line therapies (MFT), and comparing it with strategies of cycling or sequential use of single first-line ACTs.MethodsWith an individual-based microsimulation of regional malaria transmission, we looked at how to apply a therapy as widely as possible without accelerating reduction of efficacy by drug resistance. We compared simultaneous distribution of artemether–lumefantrine, artesunate–amodiaquine, and dihydroartemisinin–piperaquine (ie, MFT) against strategies in which these ACTs would be cycled or used sequentially, either on a fixed schedule or when population-level efficacy reaches the WHO threshold of 10% treatment failure. The main assessment criterion was total number of treatment failures per 100 people per year. Additionally, we analysed the benefits of including a single non-ACT therapy in an MFT strategy, and did sensitivity analyses in which we varied transmission setting, treatment coverage, partner-drug half-life, fitness cost of drug resistance, and the relation between drug concentration and resistance evolution.FindingsUse of MFT was predicted to reduce the long-term number of treatment failures compared with strategies in which a single first-line ACT is recommended. This result was robust to various epidemiological, pharmacological, and evolutionary features of malaria transmission. Inclusion of a single non-ACT therapy in an MFT strategy would have substantial benefits in reduction of pressure on artemisinin resistance evolution, delaying its emergence and slowing its spread.InterpretationAdjusting national antimalarial treatment guidelines to encourage simultaneous use of MFT is likely to extend the useful therapeutic life of available antimalarial drugs, resulting in long-term beneficial outcomes for patients.FundingWellcome Trust, UK Medical Research Council, Li Ka Shing Foundation.

Human respiratory syncytial virus (RSV) is the major cause of lower respiratory tract infections in children ,2 years of age. Little is known about RSV intra-host genetic diversity over the course of infection or about the immune pressures that drive RSV molecular evolution. We performed whole-genome deep-sequencing on 53 RSV-positive samples (37 RSV subgroup A and 16 RSV subgroup B) collected from the upper airways of hospitalized children in southern Vietnam over two consecutive seasons. RSV A NA1 and RSV B BA9 were the predominant genotypes found in our samples, consistent with other reports on global RSV circulation during the same period. For both RSV A and B, the M gene was the most conserved, confirming its potential as a target for novel therapeutics. The G gene was the most variable and was the only gene under detectable positive selection. Further, positively selected sites in G were found in close proximity to and in some cases overlapped with predicted glycosylation motifs, suggesting that selection on amino acid glycosylation may drive viral genetic diversity. We further identified hotspots and coldspots of intra-host genetic diversity in the RSV genome, some of which may highlight previously unknown regions of functional importance.

Avian influenza outbreaks have been occurring on smallholder poultry farms in Asia for two decades. Farmer responses to these outbreaks can slow down or accelerate virus transmission. We used a longitudinal survey of 53 small-scale chicken farms in southern Vietnam to investigate the impact of outbreaks with disease-induced mortality on harvest rate, vaccination, and disinfection behaviors. We found that in small broiler flocks (≤16 birds/flock) the estimated probability of harvest was 56% higher when an outbreak occurred, and 214% higher if an outbreak with sudden deaths occurred in the same month. Vaccination and disinfection were strongly and positively correlated with the number of birds. Small-scale farmers – the overwhelming majority of poultry producers in low-income countries – tend to rely on rapid sale of birds to mitigate losses from diseases. As depopulated birds are sent to markets or trading networks, this reactive behavior has the potential to enhance onward transmission.

Understanding global influenza migration and persistence is crucial for vaccine strain selection. Using 240 new human influenza A virus whole genomes collected in Vietnam during 2001–2008, we looked for persistence patterns and migratory connections between Vietnam and other countries. We found that viruses in Vietnam migrate to and from China, Hong Kong, Taiwan, Cambodia, Japan, South Korea, and the United States. We attempted to reduce geographic bias by generating phylogenies subsampled at the year and country levels. However, migration events in these phylogenies were still driven by the presence or absence of sequence data, indicating that an epidemiologic study design that controls for prevalence is required for robust migration analysis. With whole-genome data, most migration events are not detectable from the phylogeny of the hemagglutinin segment alone, although general migratory relationships between Vietnam and other countries are visible in the hemagglutinin phylogeny. It is possible that virus lineages in Vietnam persisted for >1 year.

Owing to the finding that Dengvaxia ® (the only licensed dengue vaccine to date) increases the risk of severe illness among seronegative recipients, the World Health Organization has recommended screening individuals for their serostatus prior to vaccination. To decide whether and how to carry out screening, it is necessary to estimate the transmission intensity of dengue and to understand the performance of the screening method. In this study, we inferred the annual force of infection (FOI; a measurement of transmission intensity) of dengue virus in three locations in Vietnam: An Giang (FOI = 0.04 for the below 10 years age group and FOI = 0.20 for the above 10 years age group), Ho Chi Minh City (FOI = 0.12) and Quang Ngai (FOI = 0.05). In addition, we show that using a quantitative approach to immunoglobulin G (IgG) levels (measured by indirect enzyme-linked immunosorbent assays) can help to distinguish individuals with primary exposures (primary seropositive) from those with secondary exposures (secondary seropositive). We found that primary-seropositive individuals—the main targets of the vaccine—tend to have a lower IgG level, and, thus, they have a higher chance of being misclassified as seronegative than secondary-seropositive cases. However, screening performance can be improved by incorporating patient age and transmission intensity into the interpretation of IgG levels.

BackgroundIn temperate and subtropical climates, respiratory diseases exhibit seasonal peaks in winter. In the tropics, with no winter, peak timings are irregular.MethodsTo obtain a detailed picture of influenza‐like illness (ILI) patterns in the tropics, we established an mHealth study in community clinics in Ho Chi Minh City (HCMC). During 2009‐2015, clinics reported daily case numbers via SMS, with a subset performing molecular diagnostics for influenza virus. This real‐time epidemiology network absorbs 6000 ILI reports annually, one or two orders of magnitude more than typical surveillance systems. A real‐time online ILI indicator was developed to inform clinicians of the daily ILI activity in HCMC.ResultsFrom August 2009 to December 2015, 63 clinics were enrolled and 36 920 SMS reports were received, covering approximately 1.7M outpatient visits. Approximately 10.6% of outpatients met the ILI case definition. ILI activity in HCMC exhibited strong nonannual dynamics with a dominant periodicity of 206 days. This was confirmed by time series decomposition, stepwise regression, and a forecasting exercise showing that median forecasting errors are 30%‐40% lower when using a 206‐day cycle. In ILI patients from whom nasopharyngeal swabs were taken, 31.2% were positive for influenza. There was no correlation between the ILI time series and the time series of influenza, influenza A, or influenza B (all P > 0.15).ConclusionThis suggests, for the first time, that a nonannual cycle may be an essential driver of respiratory disease dynamics in the tropics. An immunological interference hypothesis is discussed as a potential underlying mechanism.

Seroepidemiological studies aim to understand population-level exposure and immunity to infectious diseases. Their results are normally presented as binary outcomes describing the presence or absence of pathogen-specific antibody, despite the fact that many assays measure continuous quantities. A population's natural distribution of antibody titers to an endemic infectious disease may include information on multiple serological states -naiveté, recent infection, non-recent infection, childhood infection -depending on the disease in question and the acquisition and waning patterns of immunity. In this study, we investigate 20,152 general-population serum samples from southern Vietnam collected between 2009 and 2013 from which we report antibody titers to the influenza virus HA1 protein using a continuous titer measurement from a protein microarray assay. We describe the distributions of antibody titers to subtypes 2009 H1N1 and H3N2. Using a model selection approach to fit mixture distributions, we show that 2009 H1N1 antibody titers fall into four titer subgroups and that H3N2 titers fall into three subgroups. For H1N1, our interpretation is that the two highest-titer subgroups correspond to recent and historical infection, which is consistent with 2009 pandemic attack rates. Similar interpretations are available for H3N2, but right-censoring of titers makes these interpretations difficult to validate.The distribution of antibodies in a human population is a fossil imprint of the population's past exposure to infectious disease. If individuals' antibody concentrations can be measured accurately, they can be used to infer both the size and timing of past epidemics. The two key post-epidemic processes that need to be measured to make this inference possible are the rate of antibody acquisition and the rate of antibody waning. The rate of antibody acquisition post-infection is rapid (weeks) for most viral pathogens, but more difficult to measure for more complex pathogens that present the immune system with a diverse set of antigens. The rate of antibody waning, however, is rarely measured even for viral pathogens. To correctly translate a population's antibody titer distribution to its epidemic history, accurate measures of both these rates are necessary. To validate that this reconstruction has been done correctly, a large cohort with long-term follow-up and precise antibody measurements would be required. Studies like these are difficult to run and difficult to find in the scientific literature -both in methodological