The role of tunneling nanotubes during early stages of HIV infection and reactivation: implications in HIV cure.

Tunneling nanotubes (TNTs), also called cytonemes or tumor microtubes, correspond to cellular processes that enable long-range communication. TNTs are plasma membrane extensions that form tubular processes that connect the cytoplasm of two or more cells. TNTs are mostly expressed during the early stages of development and poorly expressed in adulthood. However, in disease conditions such as stroke, cancer, and viral infections such as HIV, TNTs proliferate, but their role is poorly understood. TNTs function has been associated with signaling coordination, organelle sharing, and the transfer of infectious agents such as HIV. Here, we describe the critical role and function of TNTs during HIV infection and reactivation, as well as the use of TNTs for cure strategies.

Mixed selling of different poultry species facilitates emergence of public-health-threating avian influenza viruses.

ABSTRACTLive poultry markets (LPMs) are regarded as hubs for avian influenza virus (AIV) transmission in poultry and are a major risk factor in human AIV infections. We performed an AIV surveillance study at a wholesale LPM, where different poultry species were sold in separate stalls, and nine retail LPMs, which received poultry from the wholesale LPM but where different poultry species were sold in one stall, in Guangdong province from 2017 to 2019. A higher AIV isolation rate was observed at the retail LPMs than the wholesale LPM. H9N2 was the dominant AIV subtype and was mainly present in chickens and quails. The genetic diversity of H9N2 viruses was greater at the retail LPMs, where a complex system of two-way transmission between different poultry species had formed. The isolated H9N2 viruses could be classed into four genotypes: G57 and the three novel genotypes, NG164, NG165, and NG166. The H9N2 AIVs isolated from chickens and quails at the wholesale LPM only belonged to the G57 and NG164 genotypes, respectively. However, the G57, NG164, and NG165 genotypes were identified in both chickens and quails at the retail LPMs. We found that the replication and transmission of the NG165 genotype were more adaptive to both poultry and mammalian models than those of its precursor genotype, NG164. Our findings revealed that mixed poultry selling at retail LPMs has increased the genetic diversity of AIVs, which might facilitate the emergence of novel viruses that threaten public health.

HIV infection dynamics and viral rebound: Modeling results from humanized mice.

Despite years of combined antiretroviral therapy (cART), HIV persists in infected individuals. The virus also rebounds after the cessation of cART. The sources contributing to viral persistence and rebound are not fully understood. When viral rebound occurs, what affects the time to rebound and how to delay the rebound remain unclear. In this paper, we started with the data fitting of an HIV infection model to the viral load data in treated and untreated humanized myeloid-only mice (MoM) in which macrophages serve as the target of HIV infection. By fixing the parameter values for macrophages from the MoM fitting, we fit a mathematical model including the infection of two target cell populations to the viral load data from humanized bone marrow/liver/thymus (BLT) mice, in which both CD4+ T cells and macrophages are the target of HIV infection. Data fitting suggests that the viral load decay in BLT mice under treatment has three phases. The loss of infected CD4+ T cells and macrophages is a major contributor to the first two phases of viral decay, and the last phase may be due to the latent infection of CD4+ T cells. Numerical simulations using parameter estimates from the data fitting show that the pre-ART viral load and the latent reservoir size at treatment cessation can affect viral growth rate and predict the time to viral rebound. Model simulations further reveal that early and prolonged cART can delay the viral rebound after cessation of treatment, which may have implications in the search for functional control of HIV infection.

A multi-strain model with asymptomatic transmission: Application to COVID-19 in the US.

COVID-19, induced by the SARS-CoV-2 infection, has caused an unprecedented pandemic in the world. New variants of the virus have emerged and dominated the virus population. In this paper, we develop a multi-strain model with asymptomatic transmission to study how the asymptomatic or pre-symptomatic infection influences the transmission between different strains and control strategies that aim to mitigate the pandemic. Both analytical and numerical results reveal that the competitive exclusion principle still holds for the model with the asymptomatic transmission. By fitting the model to the COVID-19 case and viral variant data in the US, we show that the omicron variants are more transmissible but less fatal than the previously circulating variants. The basic reproduction number for the omicron variants is estimated to be 11.15, larger than that for the previous variants. Using mask mandate as an example of non-pharmaceutical interventions, we show that implementing it before the prevalence peak can significantly lower and postpone the peak. The time of lifting the mask mandate can affect the emergence and frequency of subsequent waves. Lifting before the peak will result in an earlier and much higher subsequent wave. Caution should also be taken to lift the restriction when a large portion of the population remains susceptible. The methods and results obtained her e may be applied to the study of the dynamics of other infectious diseases with asymptomatic transmission using other control measures.

The impact of attrition on the transmission of HIV and drug resistance.

BACKGROUND: Attrition due to loss to follow-up or termination of antiretroviral therapy (ART) among HIV-infected patients in care may increase the risk of emergence and transmission of drug resistance (TDR), diminish benefit of treatment, and increase morbidity and mortality. Understanding the impact of attrition on the epidemic is essential to provide interventions for improving retention in care. METHODS: We developed a comprehensive HIV transmission dynamics model by considering CD4 + cell count dependent diagnosis, treatment, and attrition involving TDR and acquired drug resistance. The model was calibrated by 11 groups HIV/AIDS surveillance data during 2008-2018 from Guangxi, China, and validated by the prevalence of TDR among diagnosed treatment-naive individuals. We aimed to investigate how attrition would affect the transmission of HIV and drug-resistance when expanding ART. RESULTS: In the base case with CD4 + cell count dependent per capita attrition rates 0.025∼0.15 and treatment rates 0.23∼0.42, we projected cumulative total new infections, new drug-resistant infections, and HIV-related deaths over 2022-2030 would be 145 391, 7637, and 51 965, respectively. Increasing treatment rates by 0.1∼0.2 can decrease the above total new infections (deaths) by 1.63∼2.93% (3.52∼6.16%). However, even 0.0114∼0.0220 (0.0352∼0.0695) increase in attrition rates would offset this benefit of decreasing infections (deaths). Increasing treatment rates (attrition rates) by 0.05∼0.1 would increase the above drug-resistant infections by 0.16∼0.30% (22.18∼41.15%). CONCLUSION: A minor increase in attrition can offset the benefit of treatment expansion and increase the transmission of HIV drug resistance. Reducing attrition rates for patients already in treatment may be as important as expanding treatment for untreated patients.

Modeling the progression of Type 2 diabetes with underlying obesity.

Environmentally induced or epigenetic-related beta-cell dysfunction and insulin resistance play a critical role in the progression to diabetes. We developed a mathematical modeling framework capable of studying the progression to diabetes incorporating various diabetogenic factors. Considering the heightened risk of beta-cell defects induced by obesity, we focused on the obesity-diabetes model to further investigate the influence of obesity on beta-cell function and glucose regulation. The model characterizes individualized glucose and insulin dynamics over the span of a lifetime. We then fit the model to the longitudinal data of the Pima Indian population, which captures both the fluctuations and long-term trends of glucose levels. As predicted, controlling or eradicating the obesity-related factor can alleviate, postpone, or even reverse diabetes. Furthermore, our results reveal that distinct abnormalities of beta-cell function and levels of insulin resistance among individuals contribute to different risks of diabetes. This study may encourage precise interventions to prevent diabetes and facilitate individualized patient treatment.

The impact of supplementary immunization activities on measles transmission dynamics and implications for measles elimination goals: A mathematical modelling study.

BACKGROUND: Measles has re-emerged globally due to the accumulation of susceptible individuals and immunity gap, which causes challenges in eliminating measles. Routine vaccination and supplementary immunization activities (SIAs) have greatly improved measles control, but the impact of SIAs on the measles transmission dynamics remains unclear as the vaccine-induced immunity wanes. METHODS: We developed a comprehensive measles transmission dynamics model by taking into account population demographics, age-specific contact patterns, seasonality, routine vaccination, SIAs, and the waning vaccine-induced immunity. The model was calibrated by the monthly age-specific cases data from 2005 to 2018 in Jiangsu Province, China, and validated by the dynamic sero-prevalence data. We aimed to investigate the short-term and long-term impact of three-time SIAs during 2009-2012 (9.68 million and 4.25 million children aged 8 months-14 years in March 2009 and September 2010, respectively, and 140,000 children aged 8 months-6 years in March 2012) on the measles disease burden and explored whether additional SIAs could accelerate the measles elimination. RESULTS: We estimated that the cumulative numbers of measles cases from March 2009 to December 2012 (in the short run) and to December 2018 (in the long run) after three-time SIAs (base case) were 6,699 (95% confidence interval [CI]: 2,928-10,469), and 22,411 (15,146-29,675), which averted 45.0% (42.9%-47.0%) and 34.3% (30.7%-37.9%) of 12,226 (4,916-19,537) and 34,274 (21,350-47,199) cases without SIAs, respectively. The fraction of susceptibles for children aged 8-23 months and 2-14 years decreased from 8.3% and 10.8% in March 2009 to 5.8% and 5.8% in April 2012, respectively. However, the fraction of susceptibles aged 15-49 years and above 50 years increased gradually to about 15% in 2018 irrespective of SIAs due to the waning immunity. The measles elimination goal would be reached in 2028, and administrating additional one-off SIAs in September 2022 to children aged 8-23 months, or young adolescents aged 15-19 years could accelerate the elimination one year earlier. CONCLUSIONS: SIAs have greatly reduced the measles incidence and the fraction of susceptibles, but the benefit may wane over time. Under the current interventions, Jiangsu province would reach the measles elimination goal in 2028. Additional SIAs may accelerate the measles elimination one year earlier.

A Multiscale Model of COVID-19 Dynamics.

COVID-19, caused by the infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been a global pandemic and created unprecedented public health challenges throughout the world. Despite significant progresses in understanding the disease pathogenesis and progression, the epidemiological triad of pathogen, host, and environment remains unclear. In this paper, we develop a multiscale model to study the coupled within-host and between-host dynamics of COVID-19. The model includes multiple transmission routes (both human-to-human and environment-to-human) and connects multiple scales (both the population and individual levels). A detailed analysis on the local and global dynamics of the fast system, slow system and full system shows that rich dynamics, including both forward and backward bifurcations, emerge with the coupling of viral infection and epidemiological models. Model fitting to both virological and epidemiological data facilitates the evaluation of the influence of a few infection characteristics and antiviral treatment on the spread of the disease. Our work underlines the potential role that the environment can play in the transmission of COVID-19. Antiviral treatment of infected individuals can delay but cannot prevent the emergence of disease outbreaks. These results highlight the implementation of comprehensive intervention measures such as social distancing and wearing masks that aim to stop airborne transmission, combined with surface disinfection and hand hygiene that can prevent environmental transmission. The model also provides a multiscale modeling framework to study other infectious diseases when the environment can serve as a reservoir of pathogens.

Modelling the dynamics of Trypanosoma rangeli and triatomine bug with logistic growth of vector and systemic transmission.

In this paper, an insect-parasite-host model with logistic growth of triatomine bugs is formulated to study the transmission between hosts and vectors of the Chagas disease by using dynamical system approach. We derive the basic reproduction numbers for triatomine bugs and Trypanosoma rangeli as two thresholds. The local and global stability of the vector-free equilibrium, parasite-free equilibrium and parasite-positive equilibrium is investigated through the derived two thresholds. Forward bifurcation, saddle-node bifurcation and Hopf bifurcation are proved analytically and illustrated numerically. We show that the model can lose the stability of the vector-free equilibrium and exhibit a supercritical Hopf bifurcation, indicating the occurrence of a stable limit cycle. We also find it unlikely to have backward bifurcation and Bogdanov-Takens bifurcation of the parasite-positive equilibrium. However, the sustained oscillations of infected vector population suggest that Trypanosoma rangeli will persist in all the populations, posing a significant challenge for the prevention and control of Chagas disease.

The impact of vaccination on human papillomavirus infection with disassortative geographical mixing: a two-patch modeling study.

Human papillomavirus (HPV) infection can spread between regions. What is the impact of disassortative geographical mixing on the dynamics of HPV transmission? Vaccination is effective in preventing HPV infection. How to allocate HPV vaccines between genders within each region and between regions to reduce the total infection? Here we develop a two-patch two-sex model to address these questions. The control reproduction number [Formula: see text] under vaccination is obtained and shown to provide a critical threshold for disease elimination. Both analytical and numerical results reveal that disassortative geographical mixing does not affect [Formula: see text] and only has a minor impact on the disease prevalence in the total population given the vaccine uptake proportional to the population size for each gender in the two patches. When the vaccine uptake is not proportional to the population size, sexual mixing between the two patches can reduce [Formula: see text] and mitigate the consequence of disproportionate vaccine coverage. Using parameters calibrated from the data of a case study, we find that if the two patches have the same or similar sex ratios, allocating vaccines proportionally according to the new recruits in two patches and giving priority to the gender with a smaller recruit rate within each patch will bring the maximum benefit in reducing the total prevalence. We also show that a time-variable vaccination strategy between the two patches can further reduce the disease prevalence. This study provides some quantitative information that may help to develop vaccine distribution strategies in multiple regions with disassortative mixing.

Stochastic investigation of HIV infection and the emergence of drug resistance.

Drug-resistant HIV-1 has caused a growing concern in clinic and public health. Although combination antiretroviral therapy can contribute massively to the suppression of viral loads in patients with HIV-1, it cannot lead to viral eradication. Continuing viral replication during sub-optimal therapy (due to poor adherence or other reasons) may lead to the accumulation of drug resistance mutations, resulting in an increased risk of disease progression. Many studies also suggest that events occurring during the early stage of HIV-1 infection (i.e., the first few hours to days following HIV exposure) may determine whether the infection can be successfully established. However, the numbers of infected cells and viruses during the early stage are extremely low and stochasticity may play a critical role in dictating the fate of infection. In this paper, we use stochastic models to investigate viral infection and the emergence of drug resistance of HIV-1. The stochastic model is formulated by a continuous-time Markov chain (CTMC), which is derived based on an ordinary differential equation model proposed by Kitayimbwa et al. that includes both forward and backward mutations. An analytic estimate of the probability of the clearance of HIV infection of the CTMC model near the infection-free equilibrium is obtained by a multitype branching process approximation. The analytical predictions are validated by numerical simulations. Unlike the deterministic dynamics where the basic reproduction number R0 serves as a sharp threshold parameter (i.e., the disease dies out if R0<1 and persists if R0>1), the stochastic models indicate that there is always a positive probability for HIV infection to be eradicated in patients. In the presence of antiretroviral therapy, our results show that the chance of clearance of the infection tends to increase although drug resistance is likely to emerge.

A two-sex model of human papillomavirus infection: Vaccination strategies and a case study.

Vaccination is effective in preventing human papillomavirus (HPV) infection. It still remains debatable whether males should be included in a vaccination program and unclear how to allocate the vaccine in genders to achieve the maximum benefits. In this paper, we use a two-sex model to assess HPV vaccination strategies and use the data from Guangxi Province in China as a case study. Both mathematical analysis and numerical simulations show that the basic reproduction number, an important indicator of the transmission potential of the infection, achieves its minimum when the priority of vaccination is given to the gender with a smaller recruit rate. Given a fixed amount of vaccine, splitting the vaccine evenly usually leads to a larger basic reproduction number and a higher prevalence of infection. Vaccination becomes less effective in reducing the infection once the vaccine amount exceeds the smaller recruit rate of the two genders. In the case study, we estimate the basic reproduction number is 1.0333 for HPV 16/18 in people aged 15-55. The minimal bivalent HPV vaccine needed for the disease prevalence to be below 0.05% is 24050 per year, which should be given to females. However, with this vaccination strategy it would require a very long time and a large amount of vaccine to achieve the goal. In contrast with allocating the same vaccine amount every year, we find that a variable vaccination strategy with more vaccine given in the beginning followed by less vaccine in later years can save time and total vaccine amount. The variable vaccination strategy illustrated in this study can help to better distribute the vaccine to reduce the HPV prevalence. Although this work is for HPV infection and the case study is for a province in China, the model, analysis and conclusions may be applicable to other sexually transmitted diseases in other regions or countries.

The progression of secondary diabetes: A review of modeling studies.

Mathematical modeling has provided quantitative information consistent with experimental data, greatly improving our understanding of the progression of type 1 and type 2 diabetes. However, diabetes is a complex metabolic disease and has been found to be involved in crosstalk interactions with diverse endocrine diseases. Mathematical models have also been developed to investigate the quantitative impact of various hormonal disorders on glucose imbalance, advancing the precision treatment for secondary diabetes. Here we review the models established for the study of dysglycemia induced by hormonal disorders, such as excessive glucocorticoids, epinephrine, and growth hormone. To investigate the influence of hyperthyroidism on the glucose regulatory system, we also propose a hyperthyroid-diabetes progression model. Model simulations indicate that timely thyroid treatment can halt the progression of hyperglycemia and prevent beta-cell failure. This highlights the diagnosis of hormonal disorders, together withblood sugar tests, as significant measures for the early diagnosis and treatment of diabetes. The work recapitulates updated biological research on the interactions between the glucose regulatory system and other endocrine axes. Further mathematical modeling of secondary diabetes is desired to promote the quantitative study of the disease and the development of individualized diabetic therapies.

Investigating the relationship between reopening the economy and implementing control measures during the COVID-19 pandemic.

OBJECTIVES: The COVID-19 pandemic has resulted in an enormous burden on population health and the economy around the world. Although most cities in the United States have reopened their economies from previous lockdowns, it was not clear how the magnitude of different control measures-such as face mask use and social distancing-may affect the timing of reopening the economy for a local region. This study aimed to investigate the relationship between reopening dates and control measures and identify the conditions under which a city can be reopened safely. STUDY DESIGN: This was a mathematical modeling study. METHODS: We developed a dynamic compartment model to capture the transmission dynamics of COVID-19 in New York City. We estimated model parameters from local COVID-19 data. We conducted three sets of policy simulations to investigate how different reopening dates and magnitudes of control measures would affect the COVID-19 epidemic. RESULTS: The model estimated that maintaining social contact at 80% of the prepandemic level and a 50% face mask usage would prevent a major surge of COVID-19 after reopening. If social distancing were completely relaxed after reopening, face mask usage would need to be maintained at nearly 80% to prevent a major surge. CONCLUSIONS: Adherence to social distancing and increased face mask usage are keys to prevent a major surge after a city reopens its economy. The findings from our study can help policymakers identify the conditions under which a city can be reopened safely.

The risk of future waves of COVID-19: modeling and data analysis.

After a major outbreak of the coronavirus disease (COVID-19) starting in late December 2019, there were no new cases reported in mainland China for the first time on March 18, 2020, and no new cases reported in Hong Kong Special Administrative Region on April 20, 2020. However, these places had reported new cases and experienced a second wave since June 11, 2020. Here we develop a stochastic discrete-time epidemic model to evaluate the risk of COVID-19 resurgence by analyzing the data from the beginning of the outbreak to the second wave in these three places. In the model, we use an input parameter to represent a few potential risks that may cause a second wave, including asymptomatic infection, imported cases from other places, and virus from the environment such as frozen food packages. The effect of physical distancing restrictions imposed at different stages of the outbreak is also included in the model. Model simulations show that the magnitude of the input and the time between the initial entry and subsequent case confirmation significantly affect the probability of the second wave occurrence. Although the susceptible population size does not change the probability of resurgence, it can influence the severity of the outbreak when a second wave occurs. Therefore, to prevent the occurrence of a future wave, timely screening and detection are needed to identify infected cases in the early stage of infection. When infected cases appear, various measures such as contact tracing and quarantine should be followed to reduce the size of susceptible population in order to mitigate the COVID-19 outbreak.

A unified mathematical model of thyroid hormone regulation and implication for personalized treatment of thyroid disorders.

Current clinician practice for thyroid hormone regulation of patients is based upon guesswork and experience rather than quantified analysis, which exposes patients under longer risk and discomfort. To quantitatively analyze the thyroid regulation for patients of different thyroid states, we develop a two-dimensional mathematical model that can be applied to analyze the dynamic behaviors of thyroid hormones with or without drug intervention. The unified model can be employed to study the regulation of TSH (thyroid-stimulating hormone) and FT4 (free thyroxine) for euthyroid (normal thyroid) subjects, Hashimoto's thyroiditis, and Graves' disease patients, respectively. The results suggest that the level of TPOAb (thyroid peroxidase antibody) may be a factor determining whether the patient would progress from euthyroid state to subclinical or clinical hypothyroidism, and that increased TRAb (TSH receptor antibody) may lead Graves' disease to deteriorate from the early stage to overt hyperthyroidism. Given the early blood-test data, we demonstrate the feasibility for healthcare professionals to apply our model in choosing an appropriate dosage regimen for patients to achieve the desired TSH and FT4 levels within a specified time frame. This proposed model has the potential to optimize personalized treatment and shorten the therapeutic time for patients suffering from Hashimoto's thyroiditis and Graves' disease.

Dynamics of a new HIV model with the activation status of infected cells.

The activation status can dictate the fate of an HIV-infected CD4+ T cell. Infected cells with a low level of activation remain latent and do not produce virus, while cells with a higher level of activation are more productive and thus likely to transfer more virions to uninfected cells during cell-to-cell transmission. How the activation status of infected cells affects HIV dynamics under antiretroviral therapy remains unclear. We develop a new mathematical model that structures the population of infected cells continuously according to their activation status. The effectiveness of antiretroviral drugs in blocking cell-to-cell viral transmission decreases as the level of activation of infected cells increases because the more virions are transferred from infected to uninfected cells during cell-to-cell transmission, the less effectively the treatment is able to inhibit the transmission. The basic reproduction number [Formula: see text] of the model is shown to determine the existence and stability of the equilibria. Using the principal spectral theory and comparison principle, we show that the infection-free equilibrium is locally and globally asymptotically stable when [Formula: see text] is less than one. By constructing Lyapunov functional, we prove that the infected equilibrium is globally asymptotically stable when [Formula: see text] is greater than one. Numerical investigation shows that even when treatment can completely block cell-free virus infection, virus can still persist due to cell-to-cell transmission. The random switch between infected cells with different activation levels can also contribute to the replenishment of the latent reservoir, which is considered as a major barrier to viral eradication. This study provides a new modeling framework to study the observations, such as the low viral load persistence, extremely slow decay of latently infected cells and transient viral load measurements above the detection limit, in HIV-infected patients during suppressive antiretroviral therapy.

Effects of New York's Executive Order on Face Mask Use on COVID-19 Infections and Mortality: A Modeling Study.

There is growing evidence on the effect of face mask use in controlling the spread of COVID-19. However, few studies have examined the effect of local face mask policies on the pandemic. In this study, we developed a dynamic compartmental model of COVID-19 transmission in New York City (NYC), which was the epicenter of the COVID-19 pandemic in the USA. We used data on daily and cumulative COVID-19 infections and deaths from the NYC Department of Health and Mental Hygiene to calibrate and validate our model. We then used the model to assess the effect of the executive order on face mask use on infections and deaths due to COVID-19 in NYC. Our results showed that the executive order on face mask use was estimated to avert 99,517 (95% CIs 72,723-126,312) COVID-19 infections and 7978 (5692-10,265) deaths in NYC. If the executive order was implemented 1 week earlier (on April 10), the averted infections and deaths would be 111,475 (81,593-141,356) and 9017 (6446-11,589), respectively. If the executive order was implemented 2 weeks earlier (on April 3 when the Centers for Disease Control and Prevention recommended face mask use), the averted infections and deaths would be 128,598 (94,373-162,824) and 10,515 (7540-13,489), respectively. Our study provides public health practitioners and policymakers with evidence on the importance of implementing face mask policies in local areas as early as possible to control the spread of COVID-19 and reduce mortality.

Projected COVID-19 epidemic in the United States in the context of the effectiveness of a potential vaccine and implications for social distancing and face mask use.

BACKGROUND: Multiple candidates of COVID-19 vaccines have entered Phase III clinical trials in the United States (US). There is growing optimism that social distancing restrictions and face mask requirements could be eased with widespread vaccine adoption soon. METHODS: We developed a dynamic compartmental model of COVID-19 transmission for the four most severely affected states (New York, Texas, Florida, and California). We evaluated the vaccine effectiveness and coverage required to suppress the COVID-19 epidemic in scenarios when social contact was to return to pre-pandemic levels and face mask use was reduced. Daily and cumulative COVID-19 infection and death cases from 26th January to 15th September 2020 were obtained from the Johns Hopkins University Coronavirus resource center and used for model calibration. RESULTS: Without a vaccine (scenario 1), the spread of COVID-19 could be suppressed in these states by maintaining strict social distancing measures and face mask use levels. But relaxing social distancing restrictions to the pre-pandemic level without changing the current face mask use would lead to a new COVID-19 outbreak, resulting in 0.8-4 million infections and 15,000-240,000 deaths across these four states over the next 12 months. Under this circumstance, introducing a vaccine (scenario 2) would partially offset this negative impact even if the vaccine effectiveness and coverage are relatively low. However, if face mask use is reduced by 50% (scenario 3), a vaccine that is only 50% effective (weak vaccine) would require coverage of 55-94% to suppress the epidemic in these states. A vaccine that is 80% effective (moderate vaccine) would only require 32-57% coverage to suppress the epidemic. In contrast, if face mask usage stops completely (scenario 4), a weak vaccine would not suppress the epidemic, and further major outbreaks would occur. A moderate vaccine with coverage of 48-78% or a strong vaccine (100% effective) with coverage of 33-58% would be required to suppress the epidemic. Delaying vaccination rollout for 1-2 months would not substantially alter the epidemic trend if the current non-pharmaceutical interventions are maintained. CONCLUSIONS: The degree to which the US population can relax social distancing restrictions and face mask use will depend greatly on the effectiveness and coverage of a potential COVID-19 vaccine if future epidemics are to be prevented. Only a highly effective vaccine will enable the US population to return to life as it was before the pandemic.

A Dynamic Model to Assess Human Papillomavirus Vaccination Strategies in a Heterosexual Population Combined with Men Who have Sex with Men.

Vaccination is effective in preventing human papillomavirus (HPV) infection. It is imperative to investigate who should be vaccinated and what the best vaccine distribution strategy is. In this paper, we use a dynamic model to assess HPV vaccination strategies in a heterosexual population combined with gay, bisexual, and other men who have sex with men (MSM). The basic reproduction numbers for heterosexual females, heterosexual males and MSM as well as their average for the total population are obtained. We also derive a threshold parameter, based on basic reproduction numbers, for model analysis. From the analysis and numerical investigations, we have several conclusions. (1) To eliminate HPV infection, the priority of vaccination should be given to MSM, especially in countries that have already achieved high coverage in females. The heterosexual population gets great benefit but MSM only get minor benefit from vaccinating heterosexual females or males. (2) The best vaccination strategy is to vaccinate MSM firstly as many as possible, then heterosexual females, lastly heterosexual males. (3) Given a fixed vaccination coverage of MSM, distributing the remaining vaccines to only heterosexual females or males leads to a similar prevalence in the total population. This prevalence is lower than that when vaccines are distributed to both genders. The evener the distribution, the higher the prevalence in the total population. (4) Vaccination becomes less effective in reducing the prevalence as more vaccines are given. It is more effective to allocate vaccines to a region with lower vaccination coverage. This study provides information that may help policymakers formulate guidelines for vaccine distribution to reduce HPV prevalence on the basis of vaccine availability and prior vaccination coverage. Whether these guidelines are affected when the objective is to reduce HPV-associated cancer incidence remains to be further studied.