A geometric framework for evaluating rare variant tests of association.

The wave of next-generation sequencing data has arrived. However, many questions still remain about how to best analyze sequence data, particularly the contribution of rare genetic variants to human disease. Numerous statistical methods have been proposed to aggregate association signals across multiple rare variant sites in an effort to increase statistical power; however, the precise relation between the tests is often not well understood. We present a geometric representation for rare variant data in which rare allele counts in case and control samples are treated as vectors in Euclidean space. The geometric framework facilitates a rigorous classification of existing rare variant tests into two broad categories: tests for a difference in the lengths of the case and control vectors, and joint tests for a difference in either the lengths or angles of the two vectors. We demonstrate that genetic architecture of a trait, including the number and frequency of risk alleles, directly relates to the behavior of the length and joint tests. Hence, the geometric framework allows prediction of which tests will perform best under different disease models. Furthermore, the structure of the geometric framework immediately suggests additional classes and types of rare variant tests. We consider two general classes of tests which show robustness to noncausal and protective variants. The geometric framework introduces a novel and unique method to assess current rare variant methodology and provides guidelines for both applied and theoretical researchers.

Incorporating media data into a model of infectious disease transmission.

Understanding the effect of media on disease spread can help improve epidemic forecasting and uncover preventive measures to slow the spread of disease. Most previously introduced models have approximated media effect through disease incidence, making media influence dependent on the size of epidemic. We propose an alternative approach, which relies on real data about disease coverage in the news, allowing us to model low incidence/high interest diseases, such as SARS, Ebola or H1N1. We introduce a network-based model, in which disease is transmitted through local interactions between individuals and the probability of transmission is affected by media coverage. We assume that media attention increases self-protection (e.g. hand washing and compliance with social distancing), which, in turn, decreases disease model. We apply the model to the case of H1N1 transmission in Mexico City in 2009 and show how media influence-measured by the time series of the weekly count of news articles published on the outbreak-helps to explain the observed transmission dynamics. We show that incorporating the media attention based on the observed media coverage of the outbreak better estimates the disease dynamics from what would be predicted by using media function that approximate the media impact using the number of cases and rate of spread. Finally, we apply the model to a typical influenza season in Washington, DC and estimate how the transmission pattern would have changed given different levels of media coverage.

Predicting social response to infectious disease outbreaks from internet-based news streams.

Infectious disease outbreaks often have consequences beyond human health, including concern among the population, economic instability, and sometimes violence. A warning system capable of anticipating social disruptions resulting from disease outbreaks is urgently needed to help decision makers prepare appropriately. We designed a system that operates in near real-time to identify and predict social response. Over 150,000 Internet-based news articles related to outbreaks of 16 diseases in 72 countries and territories were provided by HealthMap. These articles were automatically tagged with indicators of the disease activity and population reaction. An anomaly detection algorithm was implemented on the population reaction indicators to identify periods of unusually severe social response. Then a model was developed to predict the probability of these periods of unusually severe social response occurring in the coming week, 2 and 3 weeks. This model exhibited remarkably strong performance for diseases with substantial media coverage. For country-disease pairs with a median of 20 or more articles per year, the onset of social response in the next week was correctly predicted over 60% of the time, and 87% of weeks were correctly predicted. Performance was weaker for diseases with little media coverage, and, for these diseases, the main utility of our system is in identifying social response when it occurs, rather than predicting when it will happen in the future. Overall, the developed near real-time prediction approach is a promising step toward developing predictive models to inform responders of the likely social consequences of disease spread.

Cost-Effective Control of Infectious Disease Outbreaks Accounting for Societal Reaction.

BACKGROUND: Studies of cost-effective disease prevention have typically focused on the tradeoff between the cost of disease transmission and the cost of applying control measures. We present a novel approach that also accounts for the cost of social disruptions resulting from the spread of disease. These disruptions, which we call social response, can include heightened anxiety, strain on healthcare infrastructure, economic losses, or violence. METHODOLOGY: The spread of disease and social response are simulated under several different intervention strategies. The modeled social response depends upon the perceived risk of the disease, the extent of disease spread, and the media involvement. Using Monte Carlo simulation, we estimate the total number of infections and total social response for each strategy. We then identify the strategy that minimizes the expected total cost of the disease, which includes the cost of the disease itself, the cost of control measures, and the cost of social response. CONCLUSIONS: The model-based simulations suggest that the least-cost disease control strategy depends upon the perceived risk of the disease, as well as media intervention. The most cost-effective solution for diseases with low perceived risk was to implement moderate control measures. For diseases with higher perceived severity, such as SARS or Ebola, the most cost-effective strategy shifted toward intervening earlier in the outbreak, with greater resources. When intervention elicited increased media involvement, it remained important to control high severity diseases quickly. For moderate severity diseases, however, it became most cost-effective to implement no intervention and allow the disease to run its course. Our simulation results imply that, when diseases are perceived as severe, the costs of social response have a significant influence on selecting the most cost-effective strategy.

Modelling the propagation of social response during a disease outbreak.

Epidemic trajectories and associated social responses vary widely between populations, with severe reactions sometimes observed. When confronted with fatal or novel pathogens, people exhibit a variety of behaviours from anxiety to hoarding of medical supplies, overwhelming medical infrastructure and rioting. We developed a coupled network approach to understanding and predicting social response. We couple the disease spread and panic spread processes and model them through local interactions between agents. The social contagion process depends on the prevalence of the disease, its perceived risk and a global media signal. We verify the model by analysing the spread of disease and social response during the 2009 H1N1 outbreak in Mexico City and 2003 severe acute respiratory syndrome and 2009 H1N1 outbreaks in Hong Kong, accurately predicting population-level behaviour. This kind of empirically validated model is critical to exploring strategies for public health intervention, increasing our ability to anticipate the response to infectious disease outbreaks.

The Role of Social Mobilization in Controlling Ebola Virus in Lofa County, Liberia.

The West Africa Ebola virus epidemic now appears to be coming to an end. In the proposed model, we simulate changes in population behavior that help to explain the observed transmission dynamics. We introduce an EVD transmission model accompanied by a model of social mobilization. The model was fit to Lofa County, Liberia through October 2014, using weekly counts of new cases reported by the US CDC. In simulation studies, we analyze the dynamics of the disease transmission with and without population behavior change, given the availability of beds in Ebola treatment units (ETUs) estimated from observed data. Only the model scenario that included individuals' behavioral change achieved a good fit to the observed case counts. Although the capacity of the Lofa County ETUs greatly increased in mid-August, our simulations show that the expansion was insufficient to alone control the outbreak. Modeling the entire outbreak without considering behavior change fit the data poorly, and extrapolating from early data without taking behavioral changes into account led to a prediction of exponential outbreak growth, contrary to the observed decline.  Education and awareness-induced behavior change in the population was instrumental in curtailing the Ebola outbreak in Lofa County and is likely playing an important role in stopping the West Africa epidemic altogether.

Assessing the impact of differential genotyping errors on rare variant tests of association.

Genotyping errors are well-known to impact the power and type I error rate in single marker tests of association. Genotyping errors that happen according to the same process in cases and controls are known as non-differential genotyping errors, whereas genotyping errors that occur with different processes in the cases and controls are known as differential genotype errors. For single marker tests, non-differential genotyping errors reduce power, while differential genotyping errors increase the type I error rate. However, little is known about the behavior of the new generation of rare variant tests of association in the presence of genotyping errors. In this manuscript we use a comprehensive simulation study to explore the effects of numerous factors on the type I error rate of rare variant tests of association in the presence of differential genotyping error. We find that increased sample size, decreased minor allele frequency, and an increased number of single nucleotide variants (SNVs) included in the test all increase the type I error rate in the presence of differential genotyping errors. We also find that the greater the relative difference in case-control genotyping error rates the larger the type I error rate. Lastly, as is the case for single marker tests, genotyping errors classifying the common homozygote as the heterozygote inflate the type I error rate significantly more than errors classifying the heterozygote as the common homozygote. In general, our findings are in line with results from single marker tests. To ensure that type I error inflation does not occur when analyzing next-generation sequencing data careful consideration of study design (e.g. use of randomization), caution in meta-analysis and using publicly available controls, and the use of standard quality control metrics is critical.