Impact of the COVID-19 pandemic on inbound air travel to Canada.

Background: Commercial air travel can result in global dispersal of infectious diseases. During the coronavirus disease 2019 (COVID-19) pandemic, many countries implemented border measures, including restrictions on air travel, to reduce the importation risk of COVID-19. In the context of inbound air travel to Canada, this study aimed to: 1) characterize travel trends before and during the pandemic, and 2) statistically assess the association between travel volumes and travel restrictions during the pandemic. Methods: Monthly commercial air travel volume data from March 2017 to February 2023 were obtained from the International Air Transport Association (IATA). National and airport-level travel trends to Canada were characterized by inbound travel volumes, the number of countries contributing travellers and the ranking of the top ten countries contributing travellers across the study period, by six year-length subperiod groupings (three pre-pandemic and three pandemic). Using seasonal autoregressive integrated moving average (SARIMA) models, interrupted time series (ITS) analyses assessed the association between major travel restrictions and travel volumes by including variables to represent changes to the level and slope of the time series. Results: The pre-pandemic inbound travel volume increased by 3% to 7% between consecutive subperiods, with three seasonal peaks (July-August, December-January, March). At the onset of the pandemic, travel volume decreased by 90%, with the number of contributing countries declining from approximately 200 to 140, followed by a slow recovery in volume and seasonality. A disruption in the ranking of countries that contributed travellers was also noticeable during the pandemic. Results from the ITS analysis aligned with the timing of travel restrictions as follows: implementation in March 2020 coincided with a sharp reduction in volumes, while the easing of major restrictions, starting with the authorization of fully vaccinated travellers from the United States to enter Canada in August 2021, coincided with an increase in the slope of travel volumes. Descriptive and statistical results suggest a near-return of pre-pandemic travel patterns by the end of the study period. Conclusion: Study results suggest resilience in commercial air travel into Canada. Although the COVID-19 pandemic led to a disruption in travel trends, easing of travel restrictions appeared to enable pre-pandemic trends to re-emerge. Understanding trends in air travel volumes, as demonstrated here, can provide information that supports preparedness and response regarding importation risk of infectious pathogens.

A methodology for estimating SARS-CoV-2 importation risk by air travel into Canada between July and November 2021.

BACKGROUND: Estimating rates of disease importation by travellers is a key activity to assess both the risk to a country from an infectious disease emerging elsewhere in the world and the effectiveness of border measures. We describe a model used to estimate the number of travellers infected with SARS-CoV-2 into Canadian airports in 2021, and assess the impact of pre-departure testing requirements on importation risk. METHODS: A mathematical model estimated the number of essential and non-essential air travellers infected with SARS-CoV-2, with the latter requiring a negative pre-departure test result. The number of travellers arriving infected (i.e. imported cases) depended on air travel volumes, SARS-CoV-2 exposure risk in the departure country, prior infection or vaccine acquired immunity, and, for non-essential travellers, screening from pre-departure molecular testing. Importation risk was estimated weekly from July to November 2021 as the number of imported cases and percent positivity (PP; i.e. imported cases normalised by travel volume). The impact of pre-departure testing was assessed by comparing three scenarios: baseline (pre-departure testing of all non-essential travellers; most probable importation risk given the pre-departure testing requirements), counterfactual scenario 1 (no pre-departure testing of fully vaccinated non-essential travellers), and counterfactual scenario 2 (no pre-departure testing of non-essential travellers). RESULTS: In the baseline scenario, weekly imported cases and PP varied over time, ranging from 145 to 539 cases and 0.15 to 0.28%, respectively. Most cases arrived from the USA, Mexico, the United Kingdom, and France. While modelling suggested that essential travellers had a higher weekly PP (0.37 - 0.65%) than non-essential travellers (0.12 - 0.24%), they contributed fewer weekly cases (62 - 154) than non-essential travellers (84 - 398 per week) given their lower travel volume. Pre-departure testing was estimated to reduce imported cases by one third (counterfactual scenario 1) to one half (counterfactual scenario 2). CONCLUSIONS: The model results highlighted the weekly variation in importation by traveller group (e.g., reason for travel and country of departure) and enabled a framework for measuring the impact of pre-departure testing requirements. Quantifying the contributors of importation risk through mathematical simulation can support the design of appropriate public health policy on border measures.

Exploring the dynamics of the 2022 mpox outbreak in Canada.

The 2022 mpox outbreak predominantly impacted gay, bisexual, and other men who have sex with men (gbMSM). Two models were developed to support situational awareness and management decisions in Canada. A compartmental model characterized epidemic drivers at national/provincial levels, while an agent-based model (ABM) assessed municipal-level impacts of vaccination. The models were parameterized and calibrated using empirical case and vaccination data between 2022 and 2023. The compartmental model explored: (1) the epidemic trajectory through community transmission, (2) the potential for transmission among non-gbMSM, and (3) impacts of vaccination and the proportion of gbMSM contributing to disease transmission. The ABM incorporated sexual-contact data and modeled: (1) effects of vaccine uptake on disease dynamics, and (2) impacts of case importation on outbreak resurgence. The calibrated, compartmental model followed the trajectory of the epidemic, which peaked in July 2022, and died out in December 2022. Most cases occurred among gbMSM, and epidemic trajectories were not consistent with sustained transmission among non-gbMSM. The ABM suggested that unprioritized vaccination strategies could increase the outbreak size by 47%, and that consistent importation (≥5 cases per 10 000) is necessary for outbreak resurgence. These models can inform time-sensitive situational awareness and policy decisions for similar future outbreaks.

Population-level effectiveness of pre-exposure prophylaxis for HIV prevention among men who have sex with men in Montréal (Canada): a modelling study of surveillance and survey data.

INTRODUCTION: HIV pre-exposure prophylaxis (PrEP) has been recommended and partly subsidized in Québec, Canada, since 2013. We evaluated the population-level impact of PrEP on HIV transmission among men who have sex with men (MSM) in Montréal, Québec's largest city, over 2013-2021. METHODS: We used an agent-based mathematical model of sexual HIV transmission to estimate the fraction of HIV acquisitions averted by PrEP compared to a counterfactual scenario without PrEP. The model was calibrated to local MSM survey, surveillance, and cohort data and accounted for COVID-19 pandemic impacts on sexual activity, HIV prevention, and care. PrEP was modelled from 2013 onwards, assuming 86% individual-level effectiveness. The PrEP eligibility criteria were: any anal sex unprotected by condoms (past 6 months) and either multiple partnerships (past 6 months) or multiple uses of post-exposure prophylaxis (lifetime). To assess potential optimization strategies, we modelled hypothetical scenarios prioritizing PrEP to MSM with high sexual activity (≥11 anal sex partners annually) or aged ⩽45 years, increasing coverage to levels achieved in Vancouver, Canada (where PrEP is free-of-charge), and improving retention. RESULTS: Over 2013-2021, the estimated annual HIV incidence decreased from 0.4 (90% credible interval [CrI]: 0.3-0.6) to 0.2 (90% CrI: 0.1-0.2) per 100 person-years. PrEP coverage among HIV-negative MSM remained low until 2015 (<1%). Afterwards, coverage increased to a maximum of 10% of all HIV-negative MSM, or about 16% of the 62% PrEP-eligible HIV-negative MSM in 2020. Over 2015-2021, PrEP averted an estimated 20% (90% CrI: 11%-30%) of cumulative HIV acquisitions. The hypothetical scenarios modelled showed that, at the same coverage level, prioritizing PrEP to high sexual activity MSM could have averted 30% (90% CrI: 19%-42%) of HIV acquisitions from 2015-2021. Even larger impacts could have resulted from higher coverage. Under the provincial eligibility criteria, reaching 10% coverage among HIV-negative MSM in 2015 and 30% in 2019, like attained in Vancouver, could have averted up to 63% (90% CrI: 54%-70%) of HIV acquisitions from 2015 to 2021. CONCLUSIONS: PrEP reduced population-level HIV transmission among Montréal MSM. However, our study suggests missed prevention opportunities and adds support for public policies that reduce PrEP barriers, financial or otherwise, to MSM at risk of HIV acquisition.

Population-level impact of expanding PrEP coverage by offering long-acting injectable PrEP to MSM in three high-resource settings: a model comparison analysis.

INTRODUCTION: Long-acting injectable cabotegravir (CAB-LA) demonstrated superiority to daily tenofovir disoproxil fumarate/emtricitabine (TDF/FTC) for HIV pre-exposure prophylaxis (PrEP) in the HPTN 083/084 trials. We compared the potential impact of expanding PrEP coverage by offering CAB-LA to men who have sex with men (MSM) in Atlanta (US), Montreal (Canada) and the Netherlands, settings with different HIV epidemics. METHODS: Three risk-stratified HIV transmission models were independently parameterized and calibrated to local data. In Atlanta, Montreal and the Netherlands, the models, respectively, estimated mean TDF/FTC coverage starting at 29%, 7% and 4% in 2022, and projected HIV incidence per 100 person-years (PY), respectively, decreasing from 2.06 to 1.62, 0.08 to 0.03 and 0.07 to 0.001 by 2042. Expansion of PrEP coverage was simulated by recruiting new CAB-LA users and by switching different proportions of TDF/FTC users to CAB-LA. Population effectiveness and efficiency of PrEP expansions were evaluated over 20 years in comparison to baseline scenarios with TDF/FTC only. RESULTS: Increasing PrEP coverage by 11 percentage points (pp) from 29% to 40% by 2032 was expected to avert a median 36% of new HIV acquisitions in Atlanta. Substantially larger increases (by 33 or 26 pp) in PrEP coverage (to 40% or 30%) were needed to achieve comparable reductions in Montreal and the Netherlands, respectively. A median 17 additional PYs on PrEP were needed to prevent one acquisition in Atlanta with 40% PrEP coverage, compared to 1000+ in Montreal and 4000+ in the Netherlands. Reaching 50% PrEP coverage by 2032 by recruiting CAB-LA users among PrEP-eligible MSM could avert >45% of new HIV acquisitions in all settings. Achieving targeted coverage 5 years earlier increased the impact by 5-10 pp. In the Atlanta model, PrEP expansions achieving 40% and 50% coverage reduced differences in PrEP access between PrEP-indicated White and Black MSM from 23 to 9 pp and 4 pp, respectively. CONCLUSIONS: Achieving high PrEP coverage by offering CAB-LA can impact the HIV epidemic substantially if rolled out without delays. These PrEP expansions may be efficient in settings with high HIV incidence (like Atlanta) but not in settings with low HIV incidence (like Montreal and the Netherlands).

Measuring progress towards reaching zero new HIV acquisitions among key populations in Québec (Canada) using routine surveillance data: a mathematical modelling study.

INTRODUCTION: Men who have sex with men (MSM) and people who inject drugs (PWID) are disproportionately impacted by the HIV epidemic in Canada. Having the second-highest provincial diagnosis rate, an improved understanding of the epidemic among these populations in Québec could aid ongoing elimination efforts. We estimated HIV incidence and other epidemic indicators among MSM and PWID in Montréal and across Québec using a back-calculation model synthesizing surveillance data. METHODS: We developed a deterministic, compartmental mathematical model stratified by age, HIV status and disease progression, and clinical care stages. Using AIDS and HIV diagnoses data, including self-reported time since the last negative test and laboratory results of CD4 cell count at diagnosis, we estimated HIV incidence in each population over 1975-2020 by modelling a cubic M-spline. The prevalence, undiagnosed fraction, fraction diagnosed that started antiretroviral treatment (ART) and median time to diagnosis were also estimated. Since the COVID-19 pandemic disrupted testing, we excluded 2020 data and explored this in sensitivity analyses. RESULTS: HIV incidence in all populations peaked early in the epidemic. In 2020, an estimated 97 (95% CrI: 33-227) and 266 (95% CrI: 103-508) HIV acquisitions occurred among MSM in Montréal and Québec, respectively. Among PWID, we estimated 2 (95% CrI: 0-14) and 6 (95% CrI: 1-26) HIV acquisitions in those same regions. With 2020 data, unless testing rates were reduced by 50%, these estimates decreased, except among Québec PWID, whose increased. Among all, the median time to diagnosis shortened to <2 years before 2020 and the undiagnosed fraction decreased to <10%. This fraction was higher in younger MSM, with 22% of 15-24 year-olds living with HIV in Montréal (95% CrI: 9-39%) and 31% in Québec (95% CrI: 17-48%) undiagnosed by 2020 year-end. Finally, ART access neared 100% in all diagnosed populations. CONCLUSIONS: HIV incidence has drastically decreased in MSM and PWID across Québec, alongside significant improvements in diagnosis and treatment coverage-and the 2013 introduction of pre-exposure prophylaxis. Despite this, HIV transmission continued. Effective efforts to halt this transmission and rapidly diagnose people who acquired HIV, especially among younger MSM, are needed to achieve elimination. Further, as the impacts of the COVID-19 pandemic on HIV transmission are understood, increased efforts may be needed to overcome these.

Past dynamics of HIV transmission among men who have sex with men in Montréal, Canada: a mathematical modeling study.

BACKGROUND: Gay, bisexual, and other men who have sex with men (gbMSM) experience disproportionate risks of HIV acquisition and transmission. In 2017, Montréal became the first Canadian Fast-Track City, setting the 2030 goal of zero new HIV infections. To inform local elimination efforts, we estimate the evolving role of prevention and sexual behaviours on HIV transmission dynamics among gbMSM in Montréal between 1975 and 2019. METHODS: Data from local bio-behavioural surveys were analyzed to develop, parameterize, and calibrate an agent-based model of sexual HIV transmission. Partnership dynamics, HIV's natural history, and treatment and prevention strategies were considered. The model simulations were analyzed to estimate the fraction of HIV acquisitions and transmissions attributable to specific groups, with a focus on age, sexual partnering level, and gaps in the HIV care-continuum. RESULTS: The model-estimated HIV incidence peaked in 1985 (2.3 per 100 person years (PY); 90% CrI: 1.4-2.9 per 100 PY) and decreased to 0.1 per 100 PY (90% CrI: 0.04-0.3 per 100 PY) in 2019. Between 2000-2017, the majority of HIV acquisitions and transmissions occurred among men aged 25-44 years, and men aged 35-44 thereafter. The unmet prevention needs of men with > 10 annual anal sex partners contributed 90-93% of transmissions and 67-73% of acquisitions annually. The primary stage of HIV played an increasing role over time, contributing to 11-22% of annual transmissions over 2000-2019. In 2019, approximately 70% of transmission events occurred from men who had discontinued, or never initiated antiretroviral therapy. CONCLUSIONS: The evolving HIV landscape has contributed to the declining HIV incidence among gbMSM in Montréal. The shifting dynamics identified in this study highlight the need for continued population-level surveillance to identify gaps in the HIV care continuum and core groups on which to prioritize elimination efforts.

Accuracy of self-reported HIV-testing history and awareness of HIV-positive status in four sub-Saharan African countries.

BACKGROUND: In many countries in sub-Saharan Africa, self-reported HIV testing history and awareness of HIV-positive status from household surveys are used to estimate the percentage of people living with HIV (PLHIV) who know their HIV status. Despite widespread use, there is limited empirical information on the sensitivity of those self-reports, which can be affected by nondisclosure. METHODS: Bayesian latent class models were used to estimate the sensitivity of self-reported HIV-testing history and awareness of HIV-positive status in four Population-based HIV Impact Assessment surveys in Eswatini, Malawi, Tanzania, and Zambia. Antiretroviral (ARV) metabolite biomarkers were used to identify persons on treatment who did not accurately report their status. For those without ARV biomarkers, we used a pooled estimate of nondisclosure among untreated persons that was 1.48 higher than those on treatment. RESULTS: Among PLHIV, the model-estimated sensitivity of self-reported HIV-testing history ranged from 96% to 99% across surveys. The model-estimated sensitivity of self-reported awareness of HIV status varied from 91% to 97%. Nondisclosure was generally higher among men and those aged 15-24 years. Adjustments for imperfect sensitivity did not substantially influence estimates of PLHIV ever tested (difference <4%) but the proportion of PLHIV aware of their HIV-positive status was higher than the unadjusted proportion (difference <8%). CONCLUSION: Self-reported HIV-testing histories in four Eastern and Southern African countries are generally robust although adjustment for nondisclosure increases estimated awareness of status. These findings can contribute to further refinements in methods for monitoring progress along the HIV testing and treatment cascade.

Pre-exposure Prophylaxis Uptake Among Men Who Have Sex With Men Who Used nPEP: A Longitudinal Analysis of Attendees at a Large Sexual Health Clinic in Montréal (Canada).

BACKGROUND: Reducing HIV transmission using pre-exposure prophylaxis (PrEP) requires focussing on individuals at high acquisition risk, such as men who have sex with men with a history of nonoccupational post-exposure prophylaxis (nPEP). This study aims to characterize longitudinal trends in PrEP uptake and its determinants among nPEP users in Montréal. METHODS: Eligible attendees at Clinique médicale l'Actuel were recruited prospectively starting in October 2000 (nPEP) and January 2013 (PrEP). Linking these cohorts, we characterized the nPEP-to-PrEP cascade, examined the determinants of PrEP uptake after nPEP consultation using a Cox proportional-hazard model, and assessed whether PrEP persistence differed by nPEP history using Kaplan-Meier curves. RESULTS: As of August 2019, 31% of 2682 nPEP cohort participants had 2 or more nPEP consultations. Subsequent PrEP consultations occurred among 36% of nPEP users, of which 17% sought nPEP again afterward. Among 2718 PrEP cohort participants, 46% reported previous nPEP use. Among nPEP users, those aged 25-49 years [hazard ratio (HR) = 1.3, 95% confidence interval (CI): 1.1 to 1.7], with more nPEP episodes (HR = 1.4, 95% CI: 1.3 to 1.5), who reported chemsex (HR = 1.3, 95% CI: 1.1 to 1.7), with a sexually transmitted infection history (HR = 1.5; 95% CI: 1.3 to 1.7), and who returned for their first nPEP follow-up visit (HR = 3.4, 95% CI: 2.7 to 4.2) had higher rates of PrEP linkage. There was no difference in PrEP persistence between nPEP-to-PrEP and PrEP only participants. CONCLUSION: Over one-third of nPEP users were subsequently prescribed PrEP. However, the large proportion of men who repeatedly use nPEP calls for more efficient PrEP-linkage services and, among those who use PrEP, improved persistence should be encouraged.

Descriptive network analysis of a Standardbred horse training facility contact network: Implications for disease transmission.

Infectious respiratory disease is a common cause of morbidity among racehorses. Quantification of contact patterns in training facilities could help inform disease prevention strategies. The study objectives were to: i) describe the contact network among horses, locations, and humans at a Standardbred horse training facility in Ontario; ii) describe the characteristics of highly influential individuals; and iii) investigate how management changes alter the network metrics and discuss the potential implications for disease transmission. Proximity loggers detected contacts among horses, staff, and locations (n = 144). Network metrics and node centrality measures were described for a 2-mode and horse-only contact network. The 2-mode network density was 0.16. and the median node degree was 20 [interquartile range (IQR) = 12 to 27]. Yearlings and floating staff were most influential in the network suggesting biosecurity programs should emphasize reducing contacts in these groups. Removing highly influential staff or co-housing of age groups resulted in changes to network diameter and density.

Analyse descriptive du réseau de contacts d'un centre d'entraînement de chevaux Standardbred : Implications pour la transmission de maladies. Les maladies respiratoires infectieuses sont une cause commune de morbidité parmi les chevaux de course. Une quantification des patrons de contact dans les centres d'entraînement pourrait aider à avoir des stratégies appropriées de prévention des maladies. Les objectifs de la présente étude étaient de : i) décrire le réseau des contacts entre les chevaux, les localisations et les humains à un centre d'entraînement pour chevaux Stadardbred en Ontario; ii) décrire les caractéristiques d'individus très influents; iii) examiner comment les changements de gestion altèrent le réseau des systèmes de mesure et discuter les implications potentielles pour la transmission des maladies. Des enregistreurs de proximité détectèrent les contacts parmi les chevaux, le personnel et les localisations (n = 144). Les systèmes de mesure et les mesures de centralité des noeuds furent décrits pour un réseau à 2 modes et un réseau de contact entre chevaux uniquement. La densité du réseau à 2 modes était de 0,16 et le degré médian du noeud était 20 [écart interquartile (IQR) = 12 à 27]. Les yearlings et le personnel occasionnel étaient les plus influents dans le réseau suggérant que les programmes de biosécurité devraient mettre l'emphase sur une réduction des contacts dans ces groupes. Le retrait de personnel très influent ou cohabitation de groupes d'âge a résulté en des changements dans le diamètre et la densité du réseau.(Traduit par Dr Serge Messier).

Comparing the effects of non-homogenous mixing patterns on epidemiological outcomes in equine populations: A mathematical modelling study.

Disease transmission models often assume homogenous mixing. This assumption, however, has the potential to misrepresent the disease dynamics for populations in which contact patterns are non-random. A disease transmission model with an SEIR structure was used to compare the effect of weighted and unweighted empirical equine contact networks to weighted and unweighted theoretical networks generated using random mixing. Equine influenza was used as a case study. Incidence curves generated with the unweighted empirical networks were similar in epidemic duration (5-8 days) and peak incidence (30.8-46.4%). In contrast, the weighted empirical networks resulted in a more pronounced difference between the networks in terms of the epidemic duration (8-15 days) and the peak incidence (5-25%). The incidence curves for the empirical networks were bimodal, while the incidence curves for the theoretical networks were unimodal. The incorporation of vaccination and isolation in the model caused a decrease in the cumulative incidence for each network, however, this effect was only seen at high levels of vaccination and isolation for the complete network. This study highlights the importance of using empirical networks to describe contact patterns within populations that are unlikely to exhibit random mixing such as equine populations.

Validation of modified radio-frequency identification tag firmware, using an equine population case study.

BACKGROUND: Contact networks can be used to assess disease spread potential within a population. However, the data required to generate the networks can be challenging to collect. One method of collecting this type of data is by using radio-frequency identification (RFID) tags. The OpenBeacon RFID system generally consists of tags and readers. Communicating tags should be within 10m of the readers, which are powered by an external power source. The readers are challenging to implement in agricultural settings due to the lack of a power source and the large area needed to be covered. METHODS: OpenBeacon firmware was modified to use the tag's onboard flash memory for data storage. The tags were deployed within an equine facility for a 7-day period. Tags were attached to the horses' halters, worn by facility staff, and placed in strategic locations around the facility to monitor which participants had contact with the specified locations during the study period. When the tags came within 2m of each other, they recorded the contact event participant IDs, and start and end times. At the end of the study period, the data were downloaded to a computer and analyzed using network analysis methods. RESULTS: The resulting networks were plausible given the facility schedule as described in a survey completed by the facility manager. Furthermore, changes in the daily facility operations as described in the survey were reflected in the tag-collected data. In terms of the battery life, 88% of batteries maintained a charge for at least 6 days. Lastly, no consistent trends were evident in the horses' centrality metrics. DISCUSSION: This study demonstrates the utility of RFID tags for the collection of equine contact data. Future work should include the collection of contact data from multiple equine facilities to better characterize equine disease spread potential in Ontario.

Comparison of the dynamic networks of four equine boarding and training facilities.

Contact networks can be analyzed to assess the potential for disease spread throughout the network. The lack of Canadian facility-level equine contact data makes the characterization of the equine contact structure difficult. Therefore, the purpose of this study was to use empirical contact data to characterize and compare equine network characteristics between equine facilities in Ontario. Contact pattern data from 4 equine facilities were collected using radio-frequency identification tags. The collected data were used to form 7 static contact networks (1 for each study day) for each facility. The assumption of homogenous mixing, where each individual in a population has an equal probability of coming in contact, was assessed for each network, since homogenous mixing is often used to describe mixing patterns in disease transmission models. At the facility level, neither the day-long static networks, nor a combined, week-long network were representative of homogenous mixing. The Jaccard Similarity Index indicated that 11-62% of the contacts were repeated throughout the study period. A network generated with survey-based data enabled the prediction of 8.7-79.6% of the contacts that were recorded with the RFID tags. With respect to the node centrality, the normalized node degree ranged from 0.0 to 0.96, with a mean of 0.31. The node strength ranged from 0 to 1 with a mean of 0.38. For both the node degree and node strength, a node's centrality score relative to the other nodes' centrality scores tended to be consistent throughout the study week. A significant (p < 0.05), weak positive correlation existed between the node degree and strength (0.41 < r < 0.54). The normalized betweenness centrality ranged from 0.00 to 1.00, with a mean of 0.11. Lastly, an exponential random graph model was used to quantify the relationship between the distance between the horses' stalls and edge formation. The distance parameter was not significant for all of the facilities. To conclude, the non-homogenous nature of the contact patterns, coupled with the large range of the centrality measures indicate the importance of using empirical data to understand processes such as disease spread potential within equine populations. Although the collection of a full set of data is optimal, the study results suggest an ability to infer contact networks using observational data in situations where little-to-no data exist. This study serves as a starting point for the characterization of equine contact networks in Ontario.

HPV Screening and Vaccination Strategies in an Unscreened Population: A Mathematical Modeling Study.

Human papillomavirus (HPV), a sexually transmitted infection, is the necessary cause of cervical cancer, the third most common cancer affecting women worldwide. Prevention and control strategies include vaccination, screening, and treatment. While HPV prevention and control efforts are important worldwide, they are especially important in low-income areas with a high infection rate or high rate of cervical cancer. This study uses mathematical modeling to explore various vaccination and treatment strategies to control for HPV and cervical cancer while using Nepal as a case study. Two sets of deterministic models were created with the goal of understanding the impact of various prevention and control strategies. The first set of models examines the relative importance of screening and vaccination in an unscreened population, while the second set examines various screening scenarios. Partial rank correlation coefficients confirm the importance of screening and treatment in the reduction of HPV infections and cancer cases even when vaccination uptake is high. Results also indicate that less expensive screening technologies can achieve the same overall goal as more expensive screening technologies.