Exploring Climate-Disease Connections in Geopolitical Versus Ecological Regions: The Case of West Nile Virus in the United States.

Many infectious disease forecasting models in the United States (US) are built with data partitioned into geopolitical regions centered on human activity as opposed to regions defined by natural ecosystems; although useful for data collection and intervention, this has the potential to mask biological relationships between the environment and disease. We explored this concept by analyzing the correlations between climate and West Nile virus (WNV) case data aggregated to geopolitical and ecological regions. We compared correlations between minimum, maximum, and mean annual temperature; precipitation; and annual WNV neuroinvasive disease (WNND) case data from 2005 to 2019 when partitioned into (a) climate regions defined by the National Oceanic and Atmospheric Administration (NOAA) and (b) Level I ecoregions defined by the Environmental Protection Agency (EPA). We found that correlations between climate and WNND in NOAA climate regions and EPA ecoregions were often contradictory in both direction and magnitude, with EPA ecoregions more often supporting previously established biological hypotheses and environmental dynamics underlying vector-borne disease transmission. Using ecological regions to examine the relationships between climate and disease cases can enhance the predictive power of forecasts at various scales, motivating a conceptual shift in large-scale analyses from geopolitical frameworks to more ecologically meaningful regions.

Correlating transcription and protein expression profiles of immune biomarkers following lipopolysaccharide exposure in lung epithelial cells.

Universal and early recognition of pathogens occurs through recognition of evolutionarily conserved pathogen associated molecular patterns (PAMPs) by innate immune receptors and the consequent secretion of cytokines and chemokines. The intrinsic complexity of innate immune signaling and associated signal transduction challenges our ability to obtain physiologically relevant, reproducible and accurate data from experimental systems. One of the reasons for the discrepancy in observed data is the choice of measurement strategy. Immune signaling is regulated by the interplay between pathogen-derived molecules with host cells resulting in cellular expression changes. However, these cellular processes are often studied by the independent assessment of either the transcriptome or the proteome. Correlation between transcription and protein analysis is lacking in a variety of studies. In order to methodically evaluate the correlation between transcription and protein expression profiles associated with innate immune signaling, we measured cytokine and chemokine levels following exposure of human cells to the PAMP lipopolysaccharide (LPS) from the Gram-negative pathogen Pseudomonas aeruginosa. Expression of 84 messenger RNA (mRNA) transcripts and 69 proteins, including 35 overlapping targets, were measured in human lung epithelial cells. We evaluated 50 biological replicates to determine reproducibility of outcomes. Following pairwise normalization, 16 mRNA transcripts and 6 proteins were significantly upregulated following LPS exposure, while only five (CCL2, CSF3, CXCL5, CXCL8/IL8, and IL6) were upregulated in both transcriptomic and proteomic analysis. This lack of correlation between transcription and protein expression data may contribute to the discrepancy in the immune profiles reported in various studies. The use of multiomic assessments to achieve a systems-level understanding of immune signaling processes can result in the identification of host biomarker profiles for a variety of infectious diseases and facilitate countermeasure design and development.

Scoping review of Culex mosquito life history trait heterogeneity in response to temperature.

BACKGROUND: Mosquitoes in the genus Culex are primary vectors in the US for West Nile virus (WNV) and other arboviruses. Climatic drivers such as temperature have differential effects on species-specific changes in mosquito range, distribution, and abundance, posing challenges for population modeling, disease forecasting, and subsequent public health decisions. Understanding these differences in underlying biological dynamics is crucial in the face of climate change. METHODS: We collected empirical data on thermal response for immature development rate, egg viability, oviposition, survival to adulthood, and adult lifespan for Culex pipiens, Cx. quinquefasciatus, Cx. tarsalis, and Cx. restuans from existing literature according to the PRISMA scoping review guidelines. RESULTS: We observed linear relationships with temperature for development rate and lifespan, and nonlinear relationships for survival and egg viability, with underlying variation between species. Optimal ranges and critical minima and maxima also appeared varied. To illustrate how model output can change with experimental input data from individual Culex species, we applied a modified equation for temperature-dependent mosquito type reproduction number for endemic spread of WNV among mosquitoes and observed different effects. CONCLUSIONS: Current models often input theoretical parameters estimated from a single vector species; we show the need to implement the real-world heterogeneity in thermal response between species and present a useful data resource for researchers working toward that goal.

Transmission risk of Oropouche fever across the Americas.

BACKGROUND: Vector-borne diseases (VBDs) are important contributors to the global burden of infectious diseases due to their epidemic potential, which can result in significant population and economic impacts. Oropouche fever, caused by Oropouche virus (OROV), is an understudied zoonotic VBD febrile illness reported in Central and South America. The epidemic potential and areas of likely OROV spread remain unexplored, limiting capacities to improve epidemiological surveillance. METHODS: To better understand the capacity for spread of OROV, we developed spatial epidemiology models using human outbreaks as OROV transmission-locality data, coupled with high-resolution satellite-derived vegetation phenology. Data were integrated using hypervolume modeling to infer likely areas of OROV transmission and emergence across the Americas. RESULTS: Models based on one-support vector machine hypervolumes consistently predicted risk areas for OROV transmission across the tropics of Latin America despite the inclusion of different parameters such as different study areas and environmental predictors. Models estimate that up to 5 million people are at risk of exposure to OROV. Nevertheless, the limited epidemiological data available generates uncertainty in projections. For example, some outbreaks have occurred under climatic conditions outside those where most transmission events occur. The distribution models also revealed that landscape variation, expressed as vegetation loss, is linked to OROV outbreaks. CONCLUSIONS: Hotspots of OROV transmission risk were detected along the tropics of South America. Vegetation loss might be a driver of Oropouche fever emergence. Modeling based on hypervolumes in spatial epidemiology might be considered an exploratory tool for analyzing data-limited emerging infectious diseases for which little understanding exists on their sylvatic cycles. OROV transmission risk maps can be used to improve surveillance, investigate OROV ecology and epidemiology, and inform early detection.

Assessing the Influence of Climate on the Spatial Pattern of West Nile Virus Incidence in the United States.

BACKGROUND: West Nile virus (WNV) is the leading cause of mosquito-borne disease in humans in the United States. Since the introduction of the disease in 1999, incidence levels have stabilized in many regions, allowing for analysis of climate conditions that shape the spatial structure of disease incidence. OBJECTIVES: Our goal was to identify the seasonal climate variables that influence the spatial extent and magnitude of WNV incidence in humans. METHODS: We developed a predictive model of contemporary mean annual WNV incidence using U.S. county-level case reports from 2005 to 2019 and seasonally averaged climate variables. We used a random forest model that had an out-of-sample model performance of R2=0.61. RESULTS: Our model accurately captured the V-shaped area of higher WNV incidence that extends from states on the Canadian border south through the middle of the Great Plains. It also captured a region of moderate WNV incidence in the southern Mississippi Valley. The highest levels of WNV incidence were in regions with dry and cold winters and wet and mild summers. The random forest model classified counties with average winter precipitation levels <23.3mm/month as having incidence levels over 11 times greater than those of counties that are wetter. Among the climate predictors, winter precipitation, fall precipitation, and winter temperature were the three most important predictive variables. DISCUSSION: We consider which aspects of the WNV transmission cycle climate conditions may benefit the most and argued that dry and cold winters are climate conditions optimal for the mosquito species key to amplifying WNV transmission. Our statistical model may be useful in projecting shifts in WNV risk in response to climate change. https://doi.org/10.1289/EHP10986.

Assessing K-12 school reopenings under different COVID-19 Spread scenarios - United States, school year 2020/21: A retrospective modeling study.

INTRODUCTION: School-age children play a key role in the spread of airborne viruses like influenza due to the prolonged and close contacts they have in school settings. As a result, school closures and other non-pharmaceutical interventions were recommended as the first line of defense in response to the novel coronavirus pandemic (COVID-19). METHODS: We used an agent-based model that simulates communities across the United States including daycares, primary, and secondary schools to quantify the relative health outcomes of reopening schools for the period of August 15, 2020 to April 11, 2021. Our simulation was carried out in early September 2020 and was based on the latest (at the time) Centers for Disease Control and Prevention (CDC)'s Pandemic Planning Scenarios released in May 2020. We explored different reopening scenarios including virtual learning, in-person school, and several hybrid options that stratify the student population into cohorts in order to reduce exposure and pathogen spread. RESULTS: Scenarios where cohorts of students return to school in non-overlapping formats, which we refer to as hybrid scenarios, resulted in significant decreases in the percentage of symptomatic individuals with COVID-19, by as much as 75%. These hybrid scenarios have only slightly more negative health impacts of COVID-19 compared to implementing a 100% virtual learning scenario. Hybrid scenarios can significantly avert the number of COVID-19 cases at the national scale-approximately between 28 M and 60 M depending on the scenario-over the simulated eight-month period. We found the results of our simulations to be highly dependent on the number of workplaces assumed to be open for in-person business, as well as the initial level of COVID-19 incidence within the simulated community. CONCLUSION: In an evolving pandemic, while a large proportion of people remain susceptible, reducing the number of students attending school leads to better health outcomes; part-time in-classroom education substantially reduces health risks.

Distinguishing viruses responsible for influenza-like illness.

The many respiratory viruses that cause influenza-like illness (ILI) are reported and tracked as one entity, defined by the CDC as a group of symptoms that include a fever of 100 degrees Fahrenheit, a cough, and/or a sore throat. In the United States alone, ILI impacts 9-49 million people every year. While tracking ILI as a single clinical syndrome is informative in many respects, the underlying viruses differ in parameters and outbreak properties. Most existing models treat either a single respiratory virus or ILI as a whole. However, there is a need for models capable of comparing several individual viruses that cause respiratory illness, including ILI. To address this need, here we present a flexible model and simulations of epidemics for influenza, RSV, rhinovirus, seasonal coronavirus, adenovirus, and SARS/MERS, parameterized by a systematic literature review and accompanied by a global sensitivity analysis. We find that for these biological causes of ILI, their parameter values, timing, prevalence, and proportional contributions differ substantially. These results demonstrate that distinguishing the viruses that cause ILI will be an important aspect of future work on diagnostics, mitigation, modeling, and preparation for future pandemics.

Updated distribution maps of predominant Culex mosquitoes across the Americas.

BACKGROUND: Estimates of the geographical distribution of Culex mosquitoes in the Americas have been limited to state and provincial levels in the United States and Canada and based on data from the 1980s. Since these estimates were made, there have been many more documented observations of mosquitoes and new methods have been developed for species distribution modeling. Moreover, mosquito distributions are affected by environmental conditions, which have changed since the 1980s. This calls for updated estimates of these distributions to understand the risk of emerging and re-emerging mosquito-borne diseases. METHODS: We used contemporary mosquito data, environmental drivers, and a machine learning ecological niche model to create updated estimates of the geographical range of seven predominant Culex species across North America and South America: Culex erraticus, Culex nigripalpus, Culex pipiens, Culex quinquefasciatus, Culex restuans, Culex salinarius, and Culex tarsalis. RESULTS: We found that Culex mosquito species differ in their geographical range. Each Culex species is sensitive to both natural and human-influenced environmental factors, especially climate and land cover type. Some prefer urban environments instead of rural ones, and some are limited to tropical or humid areas. Many are found throughout the Central Plains of the USA. CONCLUSIONS: Our updated contemporary Culex distribution maps may be used to assess mosquito-borne disease risk. It is critical to understand the current geographical distributions of these important disease vectors and the key environmental predictors structuring their distributions not only to assess current risk, but also to understand how they will respond to climate change. Since the environmental predictors structuring the geographical distribution of mosquito species varied, we hypothesize that each species may have a different response to climate change.

Trade-offs between individual and ensemble forecasts of an emerging infectious disease.

Probabilistic forecasts play an indispensable role in answering questions about the spread of newly emerged pathogens. However, uncertainties about the epidemiology of emerging pathogens can make it difficult to choose among alternative model structures and assumptions. To assess the potential for uncertainties about emerging pathogens to affect forecasts of their spread, we evaluated the performance 16 forecasting models in the context of the 2015-2016 Zika epidemic in Colombia. Each model featured a different combination of assumptions about human mobility, spatiotemporal variation in transmission potential, and the number of virus introductions. We found that which model assumptions had the most ensemble weight changed through time. We additionally identified a trade-off whereby some individual models outperformed ensemble models early in the epidemic, but on average the ensembles outperformed all individual models. Our results suggest that multiple models spanning uncertainty across alternative assumptions are necessary to obtain robust forecasts for emerging infectious diseases.

A time-varying vulnerability index for COVID-19 in New Mexico, USA using generalized propensity scores.

The coronavirus disease (COVID-19) pandemic has highlighted systemic inequities in the United States and resulted in a larger burden of negative social outcomes for marginalized communities. New Mexico, a state in the southwestern US, has a unique population with a large racial minority population and a high rate of poverty that may make communities more vulnerable to negative social outcomes from COVID-19. To identify which communities may be at the highest relative risk, we created a county-level vulnerability index. After the first COVID-19 case was reported in New Mexico on March 11, 2020, we fit a generalized propensity score model that incorporates sociodemographic factors to predict county-level viral exposure and thus, the generic risk to negative social outcomes such as unemployment or mental health impacts. We used four static sociodemographic covariates important for the state of New Mexico-population, poverty, household size, and minority population-and weekly cumulative case counts to iteratively run our model each week and normalize the exposure score to create a time-varying vulnerability index. We found the relative vulnerability between counties varied in the first eight weeks from the initial COVID-19 case before stabilizing. This framework for creating a location-specific vulnerability index in response to an ongoing disaster may be used as a quick, deployable metric to inform health policy decisions such as allocating state resources to the county level.

Using heterogeneous data to identify signatures of dengue outbreaks at fine spatio-temporal scales across Brazil.

Dengue virus remains a significant public health challenge in Brazil, and seasonal preparation efforts are hindered by variable intra- and interseasonal dynamics. Here, we present a framework for characterizing weekly dengue activity at the Brazilian mesoregion level from 2010-2016 as time series properties that are relevant to forecasting efforts, focusing on outbreak shape, seasonal timing, and pairwise correlations in magnitude and onset. In addition, we use a combination of 18 satellite remote sensing imagery, weather, clinical, mobility, and census data streams and regression methods to identify a parsimonious set of covariates that explain each time series property. The models explained 54% of the variation in outbreak shape, 38% of seasonal onset, 34% of pairwise correlation in outbreak timing, and 11% of pairwise correlation in outbreak magnitude. Regions that have experienced longer periods of drought sensitivity, as captured by the "normalized burn ratio," experienced less intense outbreaks, while regions with regular fluctuations in relative humidity had less regular seasonal outbreaks. Both the pairwise correlations in outbreak timing and outbreak trend between mesoresgions were best predicted by distance. Our analysis also revealed the presence of distinct geographic clusters where dengue properties tend to be spatially correlated. Forecasting models aimed at predicting the dynamics of dengue activity need to identify the most salient variables capable of contributing to accurate predictions. Our findings show that successful models may need to leverage distinct variables in different locations and be catered to a specific task, such as predicting outbreak magnitude or timing characteristics, to be useful. This advocates in favor of "adaptive models" rather than "one-size-fits-all" models. The results of this study can be applied to improving spatial hierarchical or target-focused forecasting models of dengue activity across Brazil.

How New Mexico Leveraged a COVID-19 Case Forecasting Model to Preemptively Address the Health Care Needs of the State: Quantitative Analysis.

BACKGROUND: Prior to the COVID-19 pandemic, US hospitals relied on static projections of future trends for long-term planning and were only beginning to consider forecasting methods for short-term planning of staffing and other resources. With the overwhelming burden imposed by COVID-19 on the health care system, an emergent need exists to accurately forecast hospitalization needs within an actionable timeframe. OBJECTIVE: Our goal was to leverage an existing COVID-19 case and death forecasting tool to generate the expected number of concurrent hospitalizations, occupied intensive care unit (ICU) beds, and in-use ventilators 1 day to 4 weeks in the future for New Mexico and each of its five health regions. METHODS: We developed a probabilistic model that took as input the number of new COVID-19 cases for New Mexico from Los Alamos National Laboratory's COVID-19 Forecasts Using Fast Evaluations and Estimation tool, and we used the model to estimate the number of new daily hospital admissions 4 weeks into the future based on current statewide hospitalization rates. The model estimated the number of new admissions that would require an ICU bed or use of a ventilator and then projected the individual lengths of hospital stays based on the resource need. By tracking the lengths of stay through time, we captured the projected simultaneous need for inpatient beds, ICU beds, and ventilators. We used a postprocessing method to adjust the forecasts based on the differences between prior forecasts and the subsequent observed data. Thus, we ensured that our forecasts could reflect a dynamically changing situation on the ground. RESULTS: Forecasts made between September 1 and December 9, 2020, showed variable accuracy across time, health care resource needs, and forecast horizon. Forecasts made in October, when new COVID-19 cases were steadily increasing, had an average accuracy error of 20.0%, while the error in forecasts made in September, a month with low COVID-19 activity, was 39.7%. Across health care use categories, state-level forecasts were more accurate than those at the regional level. Although the accuracy declined as the forecast was projected further into the future, the stated uncertainty of the prediction improved. Forecasts were within 5% of their stated uncertainty at the 50% and 90% prediction intervals at the 3- to 4-week forecast horizon for state-level inpatient and ICU needs. However, uncertainty intervals were too narrow for forecasts of state-level ventilator need and all regional health care resource needs. CONCLUSIONS: Real-time forecasting of the burden imposed by a spreading infectious disease is a crucial component of decision support during a public health emergency. Our proposed methodology demonstrated utility in providing near-term forecasts, particularly at the state level. This tool can aid other stakeholders as they face COVID-19 population impacts now and in the future.

Google Health Trends performance reflecting dengue incidence for the Brazilian states.

BACKGROUND: Dengue fever is a mosquito-borne infection transmitted by Aedes aegypti and mainly found in tropical and subtropical regions worldwide. Since its re-introduction in 1986, Brazil has become a hotspot for dengue and has experienced yearly epidemics. As a notifiable infectious disease, Brazil uses a passive epidemiological surveillance system to collect and report cases; however, dengue burden is underestimated. Thus, Internet data streams may complement surveillance activities by providing real-time information in the face of reporting lags. METHODS: We analyzed 19 terms related to dengue using Google Health Trends (GHT), a free-Internet data-source, and compared it with weekly dengue incidence between 2011 to 2016. We correlated GHT data with dengue incidence at the national and state-level for Brazil while using the adjusted R squared statistic as primary outcome measure (0/1). We used survey data on Internet access and variables from the official census of 2010 to identify where GHT could be useful in tracking dengue dynamics. Finally, we used a standardized volatility index on dengue incidence and developed models with different variables with the same objective. RESULTS: From the 19 terms explored with GHT, only seven were able to consistently track dengue. From the 27 states, only 12 reported an adjusted R squared higher than 0.8; these states were distributed mainly in the Northeast, Southeast, and South of Brazil. The usefulness of GHT was explained by the logarithm of the number of Internet users in the last 3 months, the total population per state, and the standardized volatility index. CONCLUSIONS: The potential contribution of GHT in complementing traditional established surveillance strategies should be analyzed in the context of geographical resolutions smaller than countries. For Brazil, GHT implementation should be analyzed in a case-by-case basis. State variables including total population, Internet usage in the last 3 months, and the standardized volatility index could serve as indicators determining when GHT could complement dengue state level surveillance in other countries.

Coinfections by noninteracting pathogens are not independent and require new tests of interaction.

If pathogen species, strains, or clones do not interact, intuition suggests the proportion of coinfected hosts should be the product of the individual prevalences. Independence consequently underpins the wide range of methods for detecting pathogen interactions from cross-sectional survey data. However, the very simplest of epidemiological models challenge the underlying assumption of statistical independence. Even if pathogens do not interact, death of coinfected hosts causes net prevalences of individual pathogens to decrease simultaneously. The induced positive correlation between prevalences means the proportion of coinfected hosts is expected to be higher than multiplication would suggest. By modelling the dynamics of multiple noninteracting pathogens causing chronic infections, we develop a pair of novel tests of interaction that properly account for nonindependence between pathogens causing lifelong infection. Our tests allow us to reinterpret data from previous studies including pathogens of humans, plants, and animals. Our work demonstrates how methods to identify interactions between pathogens can be updated using simple epidemic models.

Intermittent Preventive Treatment (IPT): Its Role in Averting Disease-Induced Mortality in Children and in Promoting the Spread of Antimalarial Drug Resistance.

We develop an age-structured ODE model to investigate the role of intermittent preventive treatment (IPT) in averting malaria-induced mortality in children, and its related cost in promoting the spread of antimalarial drug resistance. IPT, a malaria control strategy in which a full curative dose of an antimalarial medication is administered to vulnerable asymptomatic individuals at specified intervals, has been shown to reduce malaria transmission and deaths in children and pregnant women. However, it can also promote drug resistance spread. Our mathematical model is used to explore IPT effects on drug resistance and deaths averted in holoendemic malaria regions. The model includes drug-sensitive and drug-resistant strains as well as human hosts and mosquitoes. The basic reproduction, and invasion reproduction numbers for both strains are derived. Numerical simulations show the individual and combined effects of IPT and treatment of symptomatic infections on the prevalence of both strains and the number of lives saved. Our results suggest that while IPT can indeed save lives, particularly in high transmission regions, certain combinations of drugs used for IPT and to treat symptomatic infection may result in more deaths when resistant parasite strains are circulating. Moreover, the half-lives of the treatment and IPT drugs used play an important role in the extent to which IPT may influence spread of the resistant strain. A sensitivity analysis indicates the model outcomes are most sensitive to the reduction factor of transmission for the resistant strain, rate of immunity loss, and the natural clearance rate of sensitive infections.

Decision Support for Mitigation of Livestock Disease: Rinderpest as a Case Study.

A versatile, interactive model to predict geographically resolved epidemic progression after pathogen introduction into a population is presented. Deterministic simulations incorporating a compartmental disease model run rapidly, facilitating the analysis of mitigations such as vaccination and transmission reduction on epidemic spread and progression. We demonstrate the simulation model using rinderpest infection of cattle, a devastating livestock disease. Rinderpest has been extinguished in the wild, but it is still a threat due to stored virus in some laboratories. Comparison of simulations to historical outbreaks provides some validation of the model. Simulations of potential outbreaks demonstrate potential consequences of rinderpest virus release for a variety of possible disease parameters and mitigations. Our results indicate that a rinderpest outbreak could result in severe social and economic consequences.

Estimating the reproductive number, total outbreak size, and reporting rates for Zika epidemics in South and Central America.

As South and Central American countries prepare for increased birth defects from Zika virus outbreaks and plan for mitigation strategies to minimize ongoing and future outbreaks, understanding important characteristics of Zika outbreaks and how they vary across regions is a challenging and important problem. We developed a mathematical model for the 2015/2016 Zika virus outbreak dynamics in Colombia, El Salvador, and Suriname. We fit the model to publicly available data provided by the Pan American Health Organization, using Approximate Bayesian Computation to estimate parameter distributions and provide uncertainty quantification. The model indicated that a country-level analysis was not appropriate for Colombia. We then estimated the basic reproduction number to range between 4 and 6 for El Salvador and Suriname with a median of 4.3 and 5.3, respectively. We estimated the reporting rate to be around 16% in El Salvador and 18% in Suriname with estimated total outbreak sizes of 73,395 and 21,647 people, respectively. The uncertainty in parameter estimates highlights a need for research and data collection that will better constrain parameter ranges.

Modeling Virus Coinfection to Inform Management of Maize Lethal Necrosis in Kenya.

Maize lethal necrosis (MLN) has emerged as a serious threat to food security in sub-Saharan Africa. MLN is caused by coinfection with two viruses, Maize chlorotic mottle virus and a potyvirus, often Sugarcane mosaic virus. To better understand the dynamics of MLN and to provide insight into disease management, we modeled the spread of the viruses causing MLN within and between growing seasons. The model allows for transmission via vectors, soil, and seed, as well as exogenous sources of infection. Following model parameterization, we predict how management affects disease prevalence and crop performance over multiple seasons. Resource-rich farmers with large holdings can achieve good control by combining clean seed and insect control. However, crop rotation is often required to effect full control. Resource-poor farmers with smaller holdings must rely on rotation and roguing, and achieve more limited control. For both types of farmer, unless management is synchronized over large areas, exogenous sources of infection can thwart control. As well as providing practical guidance, our modeling framework is potentially informative for other cropping systems in which coinfection has devastating effects. Our work also emphasizes how mathematical modeling can inform management of an emerging disease even when epidemiological information remains scanty. [Formula: see text] Copyright © 2017 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Defining the Risk of Zika and Chikungunya Virus Transmission in Human Population Centers of the Eastern United States.

The recent spread of mosquito-transmitted viruses and associated disease to the Americas motivates a new, data-driven evaluation of risk in temperate population centers. Temperate regions are generally expected to pose low risk for significant mosquito-borne disease; however, the spread of the Asian tiger mosquito (Aedes albopictus) across densely populated urban areas has established a new landscape of risk. We use a model informed by field data to assess the conditions likely to facilitate local transmission of chikungunya and Zika viruses from an infected traveler to Ae. albopictus and then to other humans in USA cities with variable human densities and seasonality. Mosquito-borne disease occurs when specific combinations of conditions maximize virus-to-mosquito and mosquito-to-human contact rates. We develop a mathematical model that captures the epidemiology and is informed by current data on vector ecology from urban sites. The model demonstrates that under specific but realistic conditions, fifty-percent of introductions by infectious travelers to a high human, high mosquito density city could initiate local transmission and 10% of the introductions could result in 100 or more people infected. Despite the propensity for Ae. albopictus to bite non-human vertebrates, we also demonstrate that local virus transmission and human outbreaks may occur when vectors feed from humans even just 40% of the time. Inclusion of human behavioral changes and mitigations were not incorporated into the models and would likely reduce predicted infections. This work demonstrates how a conditional series of non-average events can result in local arbovirus transmission and outbreaks of human disease, even in temperate cities.

Two-sex mosquito model for the persistence of Wolbachia.

We develop and analyse an ordinary differential equation model to investigate the transmission dynamics of releasing Wolbachia-infected mosquitoes to establish an endemic infection in a population of wild uninfected mosquitoes. Wolbachia is a genus of endosymbiotic bacteria that can infect mosquitoes and reduce their ability to transmit some viral mosquito-transmitted diseases, including dengue fever, chikungunya, and Zika. Although the bacterium is transmitted vertically from infected mothers to their offspring, it can be difficult to establish an endemic infection in a wild mosquito population. Our transmission model for the adult and aquatic-stage mosquitoes takes into account Wolbachia-induced fitness change and cytoplasmic incompatibility. We show that, for a wide range of realistic parameter values, the basic reproduction number, [Formula: see text], is less than one. Hence, the epidemic will die out if only a few Wolbachia-infected mosquitoes are introduced into the wild population. Even though the basic reproduction number is less than one, an endemic Wolbachia infection can be established if a sufficient number of infected mosquitoes are released. This threshold effect is created by a backward bifurcation with three coexisting equilibria: a stable zero-infection equilibrium, an intermediate-infection unstable endemic equilibrium, and a high-infection stable endemic equilibrium. We analyse the impact of reducing the wild mosquito population before introducing the infected mosquitoes and observed that the most effective approach to establish the infection in the wild is based on reducing mosquitoes in both the adult and aquatic stages.