Browsing by Author "Childs, Lauren M."
Now showing 1 - 20 of 24
Results Per Page
Sort Options
- COVID-19 Seroprevalence in Canada Modelling Waning and Boosting COVID-19 Immunity in Canada a Canadian Immunization Research Network StudyDick, David W.; Childs, Lauren M.; Feng, Zhilan; Li, Jing; Röst, Gergely; Buckeridge, David L.; Ogden, Nick H.; Heffernan, Jane M. (MDPI, 2021-12-23)COVID-19 seroprevalence changes over time, with infection, vaccination, and waning immunity. Seroprevalence estimates are needed to determine when increased COVID-19 vaccination coverage is needed, and when booster doses should be considered, to reduce the spread and disease severity of COVID-19 infection. We use an age-structured model including infection, vaccination and waning immunity to estimate the distribution of immunity to COVID-19 in the Canadian population. This is the first mathematical model to do so. We estimate that 60–80% of the Canadian population has some immunity to COVID-19 by late Summer 2021, depending on specific characteristics of the vaccine and the waning rate of immunity. Models results indicate that increased vaccination uptake in age groups 12–29, and booster doses in age group 50+ are needed to reduce the severity COVID-19 Fall 2021 resurgence.
- Disrupting Mosquito Reproduction and Parasite Development for Malaria ControlChilds, Lauren M.; Cai, Francisco Y.; Kakani, Evdoxia G.; Mitchell, Sara N.; Paton, Douglas G.; Gabrieli, Paolo; Buckee, Caroline O.; Catteruccia, Flaminia (PLOS, 2016-12-15)The control of mosquito populations with insecticide treated bed nets and indoor residual sprays remains the cornerstone of malaria reduction and elimination programs. In light of widespread insecticide resistance in mosquitoes, however, alternative strategies for reducing transmission by the mosquito vector are urgently needed, including the identification of safe compounds that affect vectorial capacity via mechanisms that differ from fast-acting insecticides. Here, we show that compounds targeting steroid hormone signaling disrupt multiple biological processes that are key to the ability of mosquitoes to transmit malaria. When an agonist of the steroid hormone 20-hydroxyecdysone (20E) is applied to Anopheles gambiae females, which are the dominant malaria mosquito vector in Sub Saharan Africa, it substantially shortens lifespan, prevents insemination and egg production, and significantly blocks Plasmodium falciparum development, three components that are crucial to malaria transmission. Modeling the impact of these effects on Anopheles population dynamics and Plasmodium transmission predicts that disrupting steroid hormone signaling using 20E agonists would affect malaria transmission to a similar extent as insecticides. Manipulating 20E pathways therefore provides a powerful new approach to tackle malaria transmission by the mosquito vector, particularly in areas affected by the spread of insecticide resistance.
- Evolutionary consequences of feedbacks between within-host competition and and disease controlGreischar, Megan A.; Alexander, Helen K.; Bashey, Farrah; Bento, Ana, I.; Bhattacharya, Amrita; Bushman, Mary; Childs, Lauren M.; Daversa, David; Day, Troy; Faust, Christina L.; Gallagher, Molly E.; Gandon, Sylvain; Glidden, Caroline; Halliday, Fletcher; Hanley, Kathryn A.; Kamiya, Tsukushi; Read, Andrew F.; Schwabl, Philipp; Sweeny, Amy R.; Tate, Ann T.; Thompson, Robin N.; Wale, Nina; Wearing, Helen J.; Yeh, Pamela J.; Mideo, Nicole (2020-02-04)
- Extraordinary parasite multiplication rates in human malaria infectionsGreischar, Megan A.; Childs, Lauren M. (Cell Press, 2023-08)For pathogenic organisms, faster rates of multiplication promote transmission success, the potential to harm hosts, and the evolution of drug resistance. Parasite multiplication rates (PMRs) are often quantified in malaria infections, given the relative ease of sampling. Using modern and historical human infection data, we show that established methods return extraordinarily – and implausibly – large PMRs. We illustrate how inflated PMRs arise from two facets of malaria biology that are far from unique: (i) some developmental ages are easier to sample than others; (ii) the distribution of developmental ages changes over the course of infection. The difficulty of accurately quantifying PMRs demonstrates a need for robust methods and a subsequent re-evaluation of what is known even in the well-studied system of malaria.
- Immunoepidemiological Modeling of Dengue Viral InfectionNikin-Beers, Ryan Patrick (Virginia Tech, 2018-04-25)Dengue viral infection is a mosquito-borne disease with four distinct strains, where the interactions between these strains have implications on the severity of the disease outcomes. The two competing hypotheses for the increased severity during secondary infections are antibody dependent enhancement and original antigenic sin. Antibody dependent enhancement suggests that long-lived antibodies from primary infection remain during secondary infection but do not neutralize the virus. Original antigenic sin proposes that T cells specific to primary infection dominate cellular immune responses during secondary infections, but are inefficient at clearing cells infected with non-specific strains. To analyze these hypotheses, we developed within-host mathematical models. In previous work, we predicted a decreased non-neutralizing antibody effect during secondary infection. Since this effect accounts for decreased viral clearance and the virus is in quasi-equilibrium with infected cells, we could be accounting for reduced cell killing and the original antigenic sin hypothesis. To further understand these interactions, we develop a model of T cell responses to primary and secondary dengue virus infections that considers the effect of T cell cross-reactivity in disease enhancement. We fit the models to published patient data and show that the overall infected cell killing is similar in dengue heterologous infections, resulting in dengue fever and dengue hemorrhagic fever. The contribution to overall killing, however, is dominated by non-specific T cell responses during the majority of secondary dengue hemorrhagic fever cases. By contrast, more than half of secondary dengue fever cases have predominant strain-specific T cell responses. These results support the hypothesis that cross-reactive T cell responses occur mainly during severe disease cases of heterologous dengue virus infections. Finally, using the results from our within-host models, we develop a multiscale model of dengue viral infection which couples the within-host virus dynamics to the population level dynamics through a system of partial differential equations. We analytically determine the relationship between the model parameters and the characteristics of the solutions, and find thresholds under which infections persist in the population. Furthermore, we develop and implement a full numerical scheme for our model.
- The impact of within-vector parasite development on the extrinsic incubation periodChilds, Lauren M.; Prosper, Olivia F. (2020-10-07)Mosquito-borne diseases, in particular malaria, have a significant burden worldwide leading to nearly half a million deaths each year. The malaria parasite requires a vertebrate host, such as a human, and a vector host, the Anopheles mosquito, to complete its full life cycle. Here, we focus on the parasite dynamics within the vector to examine the first appearance of sporozoites in the salivary glands, which indicates a first time of infectiousness of mosquitoes. The timing of this period of pathogen development in the mosquito until transmissibility, known as the extrinsic incubation period, remains poorly understood. We develop compartmental models of within-mosquito parasite dynamics fitted with experimental data on oocyst and sporozoite counts. We find that only a fraction of oocysts burst to release sporozoites and bursting must be delayed either via a time-dependent function or a gamma-distributed set of compartments. We use Bayesian inference to estimate distributions of parameters and determine that bursting rate is a key epidemiological parameter. A better understanding of the factors impacting the extrinsic incubation period will aid in the development of interventions to slow or stop the spread of malaria.
- An introduction to compartmental modeling for the budding infectious disease modelerBlackwood, Julie C.; Childs, Lauren M. (Taylor & Francis, 2018-08-16)Mathematical models are ubiquitous in the study of the transmission dynamics of infectious diseases, In particular, the classic ‘susceptible-infectious-recovered’ (SIR) paradigm provides a modeling framework that can be adapted to describe the core transmission dynamics of a range of human and wildlife diseases. These models provide an important tool for uncovering the mechanisms generating observed disease dynamics, evaluating potential control strategies, and predicting future outbreaks. With ongoing advances in computational tools as well as access to disease incidence data, the use of such models continues to increase. Here, we provide a basic introduction to disease modeling that is primarily intended for individuals who are new to developing SIR-type models. In particular, we highlight several common issues encountered when structuring and analyzing these models.
- Linked within-host and between-host models and data for infectious diseases: a systematic reviewChilds, Lauren M.; El Moustaid, Fadoua; Gajewski, Zachary J.; Kadelka, Sarah; Nikin-Beers, Ryan; Smith, John W. Jr.; Walker, Melody; Johnson, Leah R. (PeerJ, 2019-06-19)The observed dynamics of infectious diseases are driven by processes across multiple scales. Here we focus on two: within-host, that is, how an infection progresses inside a single individual (for instance viral and immune dynamics), and between-host, that is, how the infection is transmitted between multiple individuals of a host population. The dynamics of each of these may be influenced by the other, particularly across evolutionary time. Thus understanding each of these scales, and the links between them, is necessary for a holistic understanding of the spread of infectious diseases. One approach to combining these scales is through mathematical modeling. We conducted a systematic review of the published literature on multi-scale mathematical models of disease transmission (as defined by combining within-host and between-host scales) to determine the extent to which mathematical models are being used to understand across-scale transmission, and the extent to which these models are being confronted with data. Following the PRISMA guidelines for systematic reviews, we identified 24 of 197 qualifying papers across 30 years that include both linked models at the within and between host scales and that used data to parameterize/calibrate models. We find that the approach that incorporates both modeling with data is under-utilized, if increasing. This highlights the need for better communication and collaboration between modelers and empiricists to build well-calibrated models that both improve understanding and may be used for prediction.
- Mathematical and Numerical Investigation of Immune System Development and FunctionKadelka, Mirjam Sarah (Virginia Tech, 2020-04-14)Mathematical models have long been used to describe complex biological interactions with the aim of predicting mechanistic interactions hard to distinguish from data. This dissertation uses modeling, mathematical analyses, and data fitting techniques to provide hypotheses on the mechanisms of immune response formation and function. The immune system, comprised of the innate and adaptive immune responses, is responsible for protecting the body against invading pathogens, with disease or vaccine induced immune memory leading to fast responses to subsequent infections. While there is some agreement about the underlying mechanisms of adaptive immune memory, innate immune memory is poorly understood. Stimulation with lipopolysaccharide induces differential phenotypes in innate immune cells depending on the strength of the stimulus, such that a secondary lipopolysaccharide encounter of a constant dose results in either strong or weak inflammatory cytokine expression. We model the biochemical kinetics of three molecules involved in macrophages responses to lipopolysaccharide and find that once a macrophage is programed to show a weak inflammatory response this cannot be reverted. Contrarily, a secondary lipopolysaccharide stimulus of a very high dose or applied prior to waning of the effects of the primary stimulus can induce a phenotype switch in macrophages initially programed to show strong inflammatory responses. Some pathogens, such as the hepatitis B virus, have developed strategies that hinder an efficient innate immune response. Hepatitis B virus infection is a worldwide pandemic with approximately 257 million chronically infected people. One beneficial event in disease progression is the seroclearance of hepatitis B e antigen often in combination with hepatitis B antibody formation. We propose mathematical models of within-host interactions and use them to predict that hepatitis B e antibody formation causes hepatitis B e antigen seroclearance and the subsequent reactivation of cytotoxic T cell immune responses. We use the model to quantify the time between antibody formation and antigen clearance and the average monthly hepatocyte turnover during that time. We further expand the study of hepatitis B infection, by investigating the kinetics of the virus under an experimental drug administered during a clinical trial. Available drugs usually fail to induce hepatitis B s antigen clearance, defined as the functional cure point of chronic hepatitis B infections. Drug therapy clinical trials that combined RNA interference drug ARC-520 with entecavir have shown promising results in reducing hepatitis B s antigen titers. We develop pharmacokinetic-pharmacodynamic models describing the mechanistic interactions of the drugs, hepatitis B virus DNA, and virus proteins. We fit the model to clinical trial data and predict that ARC-520 alone is responsible for the reduction of hepatitis B s and e antigens, while entecavir is the driving force behind viral reduction. This work was supported by Simons Foundation, Grant No. 427115, and National Science Foundation, Grant No. 1813011.
- Mathematical model of broadly reactive plasma cell productionErwin, Samantha; Childs, Lauren M.; Ciupe, Stanca M. (Nature Research, 2020)Strain-specific plasma cells are capable of producing neutralizing antibodies that are essential for clearance of challenging pathogens. These neutralizing antibodies also function as a main defense against disease establishment in a host. However, when a rapidly mutating pathogen infects a host, successful control of the invasion requires shifting the production of plasma cells from strain-specific to broadly reactive. In this study, we develop a mathematical model of germinal center dynamics and use it to predict the events that lead to improved breadth of the plasma cell response. We examine scenarios that lead to germinal centers that are composed of B-cells that come from a single strain-specific clone, a single broadly reactive clone or both clones. We find that the initial B-cell clonal composition, T-follicular helper cell signaling, increased rounds of productive somatic hypermutation, and B-cell selection strength are among the mechanisms differentiating between strain-specific and broadly reactive plasma cell production during infections. Understanding the contribution of these factors to emergence of breadth may assist in boosting broadly reactive plasma cells production.
- Mathematical Models of Immune Responses to Infectious DiseasesErwin, Samantha H. (Virginia Tech, 2017-04-04)In this dissertation, we investigate the mechanisms behind diseases and the immune responses required for successful disease resolution in three projects: i) A study of HIV and HPV co-infection, ii) A germinal center dynamics model, iii) A study of monoclonal antibody therapy. We predict that the condition leading to HPV persistence during HIV/HPV co-infection is the permissive immune environment created by HIV, rather than the direct HIV/HPV interaction. In the second project, we develop a germinal center model to understand the mechanisms that lead to the formation of potent long-lived plasma. We predict that the T follicular helper cells are a limiting resource and present possible mechanisms that can revert this limitation in the presence of non-mutating and mutating antigen. Finally, we develop a pharmacokinetic model of 3BNC117 antibody dynamics and HIV viral dynamics following antibody therapy. We fit the models to clinical trial data and conclude that antibody binding is delayed and that the combined effects of initial CD4 T cell count, initial HIV levels, and virus production are strong indicators of a good response to antibody immunotherapy.
- A Mathematical-Experimental Strategy to Decode the Complex Molecular Basis for Neutrophil Migratory Decision-MakingBoribong, Brittany Phatana (Virginia Tech, 2020-07-08)Neutrophils are the innate immune system's first line of defense in response to an infection. During an infection in the tissue, chemical cues called chemoattractants are released, which signal neutrophils to exit circulation and enter the tissue. Once in the tissue, neutrophils directionally migrate in response to the chemoattractant and toward the site of infection in a process called chemotaxis. At the site of infection, they initiate antimicrobial responses to clear the infection and resolve inflammation, restoring homeostasis. However, neutrophils are exposed to multiple chemoattractants and must prioritize these signals in order to correctly migrate to the appropriate site. The ability of neutrophils to properly undergo chemotaxis in the presence of infection and inflammation is crucial for resolution of inflammation and pathogen clearance. It has been recently shown that when pre-conditioned with bacterial endotoxin (LPS), innate immune function can become dysregulated. Neutrophils start to display altered antimicrobial response as well as dysfunctional migration patterns. This behavior has been seen in patients with sepsis, where a person's immune system overreacts to an infection, leading to systemic inflammation throughout the body, causing tissue damage, multiple organ failure, and in many cases, death. We explore the effects of inflammation on neutrophil migratory patterns and decision-making within chemotaxis. Additionally, to understand how inflammation within disease impacts chemotaxis, we measure the difference between neutrophils from healthy individuals and those from septic patients. We approached this using a combination of experimental and computational techniques. We developed a microfluidic assay to measure neutrophil decision-making in a competitive chemoattractant environment between an end-target (fMLP) and intermediary (LTB4) chemoattractant. Additionally, we probed for the expression level of molecules related to neutrophil chemotaxis. We also built a system of ordinary differential equations to model the dynamics of the molecular interactions underlying neutrophil chemotaxis. Our results showed that when neutrophils were induced into a highly inflammatory state, they prioritized pro-inflammatory signals over pro-resolution signals and displayed dysfunctional migration patterns. Similarly, neutrophils from patients with sepsis also displayed dysregulated migration patterns. This aberrant neutrophil chemotaxis may be implicated in the pathogenesis of sepsis, where accumulation of neutrophils in off-target organs is often seen. These results shed light onto the directional migratory decision-making of neutrophils exposed to inflammatory signals. Understanding these mechanisms may lead to the development of pro-resolution therapies that correct the neutrophil compass and reduce off-target organ damage.
- Modeling Temperature Effects on Vector-Borne Disease DynamicsEl Moustaid, Fadoua (Virginia Tech, 2019-09-09)Vector-borne diseases (VBDs) cause significant harm to humans, plants, and animals worldwide. For instance, VBDs are very difficult to manage, as they are governed by complex interactions. VBD transmission depends on the pathogen itself, vector-host movement, and environmental conditions. Mosquito-borne diseases are a perfect example of how all these factors contribute to changes in VBD dynamics. Although vectors are highly sensitive to climate, modeling studies tend to ignore climate effects. Here, I am interested in the arthropod small vectors that are sensitive to climate factors such as temperature, precipitation, and drought. In particular, I am looking at the effect of temperature on vector traits for two VBDs, namely, dengue, caused by a virus that infects humans and bluetongue disease, caused by a virus that infects ruminants. First, I collect data on mosquito traits' response to temperature changes, this includes adult traits as well as juvenile traits. Next, I use these traits to model mosquito density, and then I incorporate the density into our mathematical models to investigate the effect it has on the basic reproductive ratio R0, a measure of how contagious the disease is. I use R0 to determine disease risk. For dengue, my results show that using mosquito life stage traits response to temperature improves our vector density approximation and disease risk estimates. For bluetongue, I use midge traits response to temperature to show that the suitable temperature for bluetongue risk is between 21.5 °C and 30.7 °C. These results can inform future control and prevention strategies.
- Modeling the effects of Aedes aegypti’s larval environment on adult body mass at emergenceWalker, Melody; Chandrasegaran, Karthikeyan; Vinauger, Clément; Robert, Michael A.; Childs, Lauren M. (PLoS, 2021-11-22)Mosquitoes vector harmful pathogens that infect millions of people every year, and developing approaches to effectively control mosquitoes is a topic of great interest. However, the success of many control measures is highly dependent upon ecological, physiological, and life history traits of mosquito species. The behavior of mosquitoes and their potential to vector pathogens can also be impacted by these traits. One trait of interest is mosquito body mass, which depends upon many factors associated with the environment in which juvenile mosquitoes develop. Our experiments examined the impact of larval density on the body mass of Aedes aegypti mosquitoes, which are important vectors of dengue, Zika, yellow fever, and other pathogens. To investigate the interactions between the larval environment and mosquito body mass, we built a discrete time mathematical model that incorporates body mass, larval density, and food availability and fit the model to our experimental data. We considered three categories of model complexity informed by data, and selected the best model within each category using Akaike’s Information Criterion. We found that the larval environment is an important determinant of the body mass of mosquitoes upon emergence. Furthermore, we found that larval density has greater impact on body mass of adults at emergence than on development time, and that inclusion of density dependence in the survival of female aquatic stages in models is important. We discuss the implications of our results for the control of Aedes mosquitoes and on their potential to spread disease.
- Modeling the waning and boosting of immunity from infection or vaccinationCarlsson, Rose-Marie; Childs, Lauren M.; Feng, Zhilan; Glasser, John W.; Heffernan, Jane M.; Li, Jing; Röst, Gergely (2020-07-21)Immunity following natural infection or immunization may wane, increasing susceptibility to infection with time since infection or vaccination. Symptoms, and concomitantly infectiousness, depend on residual immunity. We quantify these phenomena in a model population composed of individuals whose susceptibility, infectiousness, and symptoms all vary with immune status. We also model age, which affects contact, vaccination and possibly waning rates. The resurgences of pertussis that have been observed wherever effective vaccination programs have reduced typical disease among young children follow from these processes. As one example, we compare simulations with the experience of Sweden following resumption of pertussis vaccination after the hiatus from 1979 to 1996, reproducing the observations leading health authorities to introduce booster doses among school-aged children and adolescents in 2007 and 2014, respectively. Because pertussis comprises a spectrum of symptoms, only the most severe of which are medically attended, accurate models are needed to design optimal vaccination programs where surveillance is less effective. (C) 2020 The Authors. Published by Elsevier Ltd.
- Modelling Allee effects in a transgenic mosquito population during range expansionWalker, Melody (Virginia Tech, 2018-06-20)Mosquitoes are vectors for many diseases that cause significant mortality and morbidity across the globe such as malaria, dengue fever and Zika. As mosquito populations expand their range into new areas, they may undergo mate-finding Allee effects such that their ability to successfully reproduce becomes difficult at low population densities. With new technology, creating target specific gene modification may now be a viable method for mosquito population control. We develop a mathematical model to investigate the effects of releasing transgenic mosquitoes into newly established low-density mosquito populations. Our model consists of two life stages (aquatic and adult), which are further divided into three genetically distinct groups: heterogeneous and homogeneous transgenic alleles that cause female infertility and a homogeneous wild type. We perform analytical and numerical analyses on the equilibria to determine the level of saturation needed to eliminate mosquitoes in a given area. This model demonstrates the potential for a gene drive system to reduce the spread of invading mosquito populations.
- Modelling Allee effects in a transgenic mosquito population during range expansionWalker, Melody; Blackwood, Julie C.; Brown, Vicki; Childs, Lauren M. (Taylor & Francis, 2018-04-27)Mosquitoes are vectors for many diseases that cause significant mortality and morbidity. As mosquito populations expand their range, they may undergo mate-finding Allee effects such that their ability to successfully reproduce becomes difficult at low population density. With new technology, creating target specific gene modification may be a viable method for mosquito population control. We develop a mathematical model to investigate the effects of releasing transgenic mosquitoes into newly established, low-density mosquito populations. Our model consists of two life stages (aquatic and adults), which are divided into three genetically distinct groups: heterogeneous and homogeneous transgenic that cause female infertility and a homogeneous wild type. We perform analytical and numerical analyses on the equilibria to determine the level of saturation needed to eliminate mosquitoes in a given area. This model demonstrates the potential for a gene drive system to reduce the spread of invading mosquito populations.
- Multiple blood feeding in mosquitoes shortens the Plasmodium falciparum incubation period and increases malaria transmission potentialShaw, W. Robert; Holmdahl, Inga E.; Itoe, Maurice A.; Werling, Kristine; Marquette, Meghan; Paton, Douglas G.; Singh, Naresh; Buckee, Caroline O.; Childs, Lauren M.; Catteruccia, Flaminia (PLOS, 2020-12-31)Many mosquito species, including the major malaria vector Anopheles gambiae, naturally undergo multiple reproductive cycles of blood feeding, egg development and egg laying in their lifespan. Such complex mosquito behavior is regularly overlooked when mosquitoes are experimentally infected with malaria parasites, limiting our ability to accurately describe potential effects on transmission. Here, we examine how Plasmodium falciparum development and transmission potential is impacted when infected mosquitoes feed an additional time. We measured P. falciparum oocyst size and performed sporozoite time course analyses to determine the parasite’s extrinsic incubation period (EIP), i.e. the time required by parasites to reach infectious sporozoite stages, in An. gambiae females blood fed either once or twice. An additional blood feed at 3 days post infection drastically accelerates oocyst growth rates, causing earlier sporozoite accumulation in the salivary glands, thereby shortening the EIP (reduction of 2.3 ± 0.4 days). Moreover, parasite growth is further accelerated in transgenic mosquitoes with reduced reproductive capacity, which mimic genetic modifications currently proposed in population suppression gene drives. We incorporate our shortened EIP values into a measure of transmission potential, the basic reproduction number R0, and find the average R0 is higher (range: 10.1%–12.1% increase) across sub-Saharan Africa than when using traditional EIP measurements. These data suggest that malaria elimination may be substantially more challenging and that younger mosquitoes or those with reduced reproductive ability may provide a larger contribution to infection than currently believed. Our findings have profound implications for current and future mosquito control interventions.
- Nonsmooth Bifurcations and the Role of Density Dependence in a Chaotic Infectious Disease ModelHughes, Ryan Patrick (Virginia Tech, 2020-01-23)Discrete dynamical systems can exhibit rich and interesting dynamics at lower dimensions (and co-dimensions) than that of ODE models. Classically, the minimal dimension to observe chaotic behavior in an ODE model is three; whereas it can be achieved in a one-dimensional discrete map. It is often the choice of mathematical biologists to use discrete systems as it fills many roles such as sparse data, incorporation of life cycle stages and noisy measurements. This work is analyzes a discrete time model of an infected salmon population. It provides an in-depth analysis of non-smooth bifurcations for alternate functional forms for density dependence in the growth function of a given model. These demonstrate interesting structures and chaotic behaviors with biologically feasible interpretations such as intrinsic growth rate and probability of death. The choice of density dependence function, as well as parameterization, leads to whether chaos occurs or not.
- Parasitism of Aedes albopictus by Ascogregarina taiwanensis lowers its competitive ability against Aedes triseriatusStump, Emma; Childs, Lauren M.; Walker, Melody (2021-01-25)Background Mosquitoes are vectors for diseases such as dengue, malaria and La Crosse virus that significantly impact the human population. When multiple mosquito species are present, the competition between species may alter population dynamics as well as disease spread. Two mosquito species, Aedes albopictus and Aedes triseriatus, both inhabit areas where La Crosse virus is found. Infection of Aedes albopictus by the parasite Ascogregarina taiwanensis and Aedes triseriatus by the parasite Ascogregarina barretti can decrease a mosquito’s fitness, respectively. In particular, the decrease in fitness of Aedes albopictus occurs through the impact of Ascogregarina taiwanensis on female fecundity, larval development rate, and larval mortality and may impact its initial competitive advantage over Aedes triseriatus during invasion. Methods We examine the effects of parasitism of gregarine parasites on Aedes albopictus and triseriatus population dynamics and competition with a focus on when Aedes albopictus is new to an area. We build a compartmental model including competition between Aedes albopictus and triseriatus while under parasitism of the gregarine parasites. Using parameters based on the literature, we simulate the dynamics and analyze the equilibrium population proportion of the two species. We consider the presence of both parasites and potential dilution effects. Results We show that increased levels of parasitism in Aedes albopictus will decrease the initial competitive advantage of the species over Aedes triseriatus and increase the survivorship of Aedes triseriatus. We find Aedes albopictus is better able to invade when there is more extreme parasitism of Aedes triseriatus. Furthermore, although the transient dynamics differ, dilution of the parasite density through uptake by both species does not alter the equilibrium population sizes of either species. Conclusions Mosquito population dynamics are affected by many factors, such as abiotic factors (e.g. temperature and humidity) and competition between mosquito species. This is especially true when multiple mosquito species are vying to live in the same area. Knowledge of how population dynamics are affected by gregarine parasites among competing species can inform future mosquito control efforts and help prevent the spread of vector-borne disease.