Reconstructed influenza A/H3N2 infection histories reveal variation in incidence and antibody dynamics over the life course.

Humans experience many influenza infections over their lives, resulting in complex and varied immunological histories. Although experimental and quantitative analyses have improved our understanding of the immunological processes defining an individual's antibody repertoire, how these within-host processes are linked to population-level influenza epidemiology remains unclear. Here, we used a multi-level mathematical model to jointly infer antibody dynamics and individual-level lifetime influenza A/H3N2 infection histories for 1,130 individuals in Guangzhou, China, using 67,683 haemagglutination inhibition (HI) assay measurements against 20 A/H3N2 strains from repeat serum samples collected between 2009 and 2015. These estimated infection histories allowed us to reconstruct historical seasonal influenza patterns and to investigate how influenza incidence varies over time, space and age in this population. We estimated median annual influenza infection rates to be approximately 18% from 1968 to 2015, but with substantial variation between years. 88% of individuals were estimated to have been infected at least once during the study period (2009-2015), and 20% were estimated to have three or more infections in that time. We inferred decreasing infection rates with increasing age, and found that annual attack rates were highly correlated across all locations, regardless of their distance, suggesting that age has a stronger impact than fine-scale spatial effects in determining an individual's antibody profile. Finally, we reconstructed each individual's expected antibody profile over their lifetime and inferred an age-stratified relationship between probability of infection and HI titre. Our analyses show how multi-strain serological panels provide rich information on long term, epidemiological trends, within-host processes and immunity when analyzed using appropriate inference methods, and adds to our understanding of the life course epidemiology of influenza A/H3N2.

Viral kinetics of sequential SARS-CoV-2 infections.

The impact of a prior SARS-CoV-2 infection on the progression of subsequent infections has been unclear. Using a convenience sample of 94,812 longitudinal RT-qPCR measurements from anterior nares and oropharyngeal swabs, we identified 71 individuals with two well-sampled SARS-CoV-2 infections between March 11th, 2020, and July 28th, 2022. We compared the SARS-CoV-2 viral kinetics of first vs. second infections in this group, adjusting for viral variant, vaccination status, and age. Relative to first infections, second infections usually featured a faster clearance time. Furthermore, a person's relative (rank-order) viral clearance time, compared to others infected with the same variant, was roughly conserved across first and second infections, so that individuals who had a relatively fast clearance time in their first infection also tended to have a relatively fast clearance time in their second infection (Spearman correlation coefficient: 0.30, 95% credible interval (0.12, 0.46)). These findings provide evidence that, like vaccination, immunity from a prior SARS-CoV-2 infection shortens the duration of subsequent acute SARS-CoV-2 infections principally by reducing viral clearance time. Additionally, there appears to be an inherent element of the immune response, or some other host factor, that shapes a person's relative ability to clear SARS-CoV-2 infection that persists across sequential infections.

serosim: An R package for simulating serological data arising from vaccination, epidemiological and antibody kinetics processes.

serosim is an open-source R package designed to aid inference from serological studies, by simulating data arising from user-specified vaccine and antibody kinetics processes using a random effects model. Serological data are used to assess population immunity by directly measuring individuals' antibody titers. They uncover locations and/or populations which are susceptible and provide evidence of past infection or vaccination to help inform public health measures and surveillance. Both serological data and new analytical techniques used to interpret them are increasingly widespread. This creates a need for tools to simulate serological studies and the processes underlying observed titer values, as this will enable researchers to identify best practices for serological study design, and provide a standardized framework to evaluate the performance of different inference methods. serosim allows users to specify and adjust model inputs representing underlying processes responsible for generating the observed titer values like time-varying patterns of infection and vaccination, population demography, immunity and antibody kinetics, and serological sampling design in order to best represent the population and disease system(s) of interest. This package will be useful for planning sampling design of future serological studies, understanding determinants of observed serological data, and validating the accuracy and power of new statistical methods.

Alternative epidemic indicators for COVID-19 in three settings with incomplete death registration systems.

Not all COVID-19 deaths are officially reported, and particularly in low-income and humanitarian settings, the magnitude of reporting gaps remains sparsely characterized. Alternative data sources, including burial site worker reports, satellite imagery of cemeteries, and social media-conducted surveys of infection may offer solutions. By merging these data with independently conducted, representative serological studies within a mathematical modeling framework, we aim to better understand the range of underreporting using examples from three major cities: Addis Ababa (Ethiopia), Aden (Yemen), and Khartoum (Sudan) during 2020. We estimate that 69 to 100%, 0.8 to 8.0%, and 3.0 to 6.0% of COVID-19 deaths were reported in each setting, respectively. In future epidemics, and in settings where vital registration systems are limited, using multiple alternative data sources could provide critically needed, improved estimates of epidemic impact. However, ultimately, these systems are needed to ensure that, in contrast to COVID-19, the impact of future pandemics or other drivers of mortality is reported and understood worldwide.

Long term intrinsic cycling in human life course antibody responses to influenza A(H3N2): an observational and modeling study.

Background: Over a life course, human adaptive immunity to antigenically mutable pathogens exhibits competitive and facilitative interactions. We hypothesize that such interactions may lead to cyclic dynamics in immune responses over a lifetime. Methods: To investigate the cyclic behavior, we analyzed hemagglutination inhibition titers against 21 historical influenza A(H3N2) strains spanning 47 years from a cohort in Guangzhou, China, and applied Fourier spectrum analysis. To investigate possible biological mechanisms, we simulated individual antibody profiles encompassing known feedbacks and interactions due to generally recognized immunological mechanisms. Results: We demonstrated a long-term periodicity (about 24 years) in individual antibody responses. The reported cycles were robust to analytic and sampling approaches. Simulations suggested that individual-level cross-reaction between antigenically similar strains likely explains the reported cycle. We showed that the reported cycles are predictable at both individual and birth cohort level and that cohorts show a diversity of phases of these cycles. Phase of cycle was associated with the risk of seroconversion to circulating strains, after accounting for age and pre-existing titers of the circulating strains. Conclusions: Our findings reveal the existence of long-term periodicities in individual antibody responses to A(H3N2). We hypothesize that these cycles are driven by preexisting antibody responses blunting responses to antigenically similar pathogens (by preventing infection and/or robust antibody responses upon infection), leading to reductions in antigen-specific responses over time until individual's increasing risk leads to an infection with an antigenically distant enough virus to generate a robust immune response. These findings could help disentangle cohort effects from individual-level exposure histories, improve our understanding of observed heterogeneous antibody responses to immunizations, and inform targeted vaccine strategy. Funding: This study was supported by grants from the NIH R56AG048075 (DATC, JL), NIH R01AI114703 (DATC, BY), the Wellcome Trust 200861/Z/16/Z (SR), and 200187/Z/15/Z (SR). This work was also supported by research grants from Guangdong Government HZQB-KCZYZ-2021014 and 2019B121205009 (YG and HZ). DATC, JMR and SR acknowledge support from the National Institutes of Health Fogarty Institute (R01TW0008246). JMR acknowledges support from the Medical Research Council (MR/S004793/1) and the Engineering and Physical Sciences Research Council (EP/N014499/1). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Quantifying the impact of immune history and variant on SARS-CoV-2 viral kinetics and infection rebound: A retrospective cohort study.

Background: The combined impact of immunity and SARS-CoV-2 variants on viral kinetics during infections has been unclear. Methods: We characterized 1,280 infections from the National Basketball Association occupational health cohort identified between June 2020 and January 2022 using serial RT-qPCR testing. Logistic regression and semi-mechanistic viral RNA kinetics models were used to quantify the effect of age, variant, symptom status, infection history, vaccination status and antibody titer to the founder SARS-CoV-2 strain on the duration of potential infectiousness and overall viral kinetics. The frequency of viral rebounds was quantified under multiple cycle threshold (Ct) value-based definitions. Results: Among individuals detected partway through their infection, 51.0% (95% credible interval [CrI]: 48.3-53.6%) remained potentially infectious (Ct <30) 5 days post detection, with small differences across variants and vaccination status. Only seven viral rebounds (0.7%; N=999) were observed, with rebound defined as 3+days with Ct <30 following an initial clearance of 3+days with Ct ≥30. High antibody titers against the founder SARS-CoV-2 strain predicted lower peak viral loads and shorter durations of infection. Among Omicron BA.1 infections, boosted individuals had lower pre-booster antibody titers and longer clearance times than non-boosted individuals. Conclusions: SARS-CoV-2 viral kinetics are partly determined by immunity and variant but dominated by individual-level variation. Since booster vaccination protects against infection, longer clearance times for BA.1-infected, boosted individuals may reflect a less effective immune response, more common in older individuals, that increases infection risk and reduces viral RNA clearance rate. The shifting landscape of viral kinetics underscores the need for continued monitoring to optimize isolation policies and to contextualize the health impacts of therapeutics and vaccines. Funding: Supported in part by CDC contract #200-2016-91779, a sponsored research agreement to Yale University from the National Basketball Association contract #21-003529, and the National Basketball Players Association.

Near real-time surveillance of the SARS-CoV-2 epidemic with incomplete data.

When responding to infectious disease outbreaks, rapid and accurate estimation of the epidemic trajectory is critical. However, two common data collection problems affect the reliability of the epidemiological data in real time: missing information on the time of first symptoms, and retrospective revision of historical information, including right censoring. Here, we propose an approach to construct epidemic curves in near real time that addresses these two challenges by 1) imputation of dates of symptom onset for reported cases using a dynamically-estimated "backward" reporting delay conditional distribution, and 2) adjustment for right censoring using the NobBS software package to nowcast cases by date of symptom onset. This process allows us to obtain an approximation of the time-varying reproduction number (Rt) in real time. We apply this approach to characterize the early SARS-CoV-2 outbreak in two Spanish regions between March and April 2020. We evaluate how these real-time estimates compare with more complete epidemiological data that became available later. We explore the impact of the different assumptions on the estimates, and compare our estimates with those obtained from commonly used surveillance approaches. Our framework can help improve accuracy, quantify uncertainty, and evaluate frequently unstated assumptions when recovering the epidemic curves from limited data obtained from public health systems in other locations.

Estimation of Seasonal Influenza Attack Rates and Antibody Dynamics in Children Using Cross-Sectional Serological Data.

Directly measuring evidence of influenza infections is difficult, especially in low-surveillance settings such as sub-Saharan Africa. Using a Bayesian model, we estimated unobserved infection times and underlying antibody responses to influenza A/H3N2, using cross-sectional serum antibody responses to 4 strains in children aged 24-60 months. Among the 242 individuals, we estimated a variable seasonal attack rate and found that most children had ≥1 infection before 2 years of age. Our results are consistent with previously published high attack rates in children. The modeling approach highlights how cross-sectional serological data can be used to estimate epidemiological dynamics.

Correction: Practical considerations for measuring the effective reproductive number, Rt.

[This corrects the article DOI: 10.1371/journal.pcbi.1008409.].

When intuition falters: repeated testing accuracy during an epidemic.

Widespread, repeated testing using rapid antigen tests to proactively detect asymptomatic SARS-CoV-2 infections has been a promising yet controversial topic during the COVID-19 pandemic. Concerns have been raised over whether currently authorized lateral flow tests are sufficiently sensitive and specific to detect enough infections to impact transmission whilst minimizing unnecessary isolation of false positives. These concerns have often been illustrated using simple, textbook calculations of positivity rates and positive predictive value assuming fixed values for sensitivity, specificity and prevalence. However, we argue that evaluating repeated testing strategies requires the consideration of three additional factors: new infections continue to arise depending on the incidence rate, isolating positive individuals reduces prevalence in the tested population, and each infected individual is tested multiple times during their infection course. We provide a simple mathematical model with an online interface to illustrate how these three factors impact test positivity rates and the number of isolating individuals over time. These results highlight the potential pitfalls of using inappropriate textbook-style calculations to evaluate statistics arising from repeated testing strategies during an epidemic.

Estimating epidemiologic dynamics from cross-sectional viral load distributions.

Estimating an epidemic's trajectory is crucial for developing public health responses to infectious diseases, but case data used for such estimation are confounded by variable testing practices. We show that the population distribution of viral loads observed under random or symptom-based surveillance-in the form of cycle threshold (Ct) values obtained from reverse transcription quantitative polymerase chain reaction testing-changes during an epidemic. Thus, Ct values from even limited numbers of random samples can provide improved estimates of an epidemic's trajectory. Combining data from multiple such samples improves the precision and robustness of this estimation. We apply our methods to Ct values from surveillance conducted during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in a variety of settings and offer alternative approaches for real-time estimates of epidemic trajectories for outbreak management and response.

Estimation of Transmission of COVID-19 in Simulated Nursing Homes With Frequent Testing and Immunity-Based Staffing.

Importance: Nursing homes and other long-term care facilities have been disproportionately impacted by the COVID-19 pandemic. Strategies are urgently needed to reduce transmission in these high-risk populations. Objective: To evaluate COVID-19 transmission in nursing homes associated with contact-targeted interventions and testing. Design, Setting, and Participants: This decision analytical modeling study developed an agent-based susceptible-exposed-infectious (asymptomatic/symptomatic)-recovered model between July and September 2020 to examine SARS-CoV-2 transmission in nursing homes. Residents and staff of a simulated nursing home with 100 residents and 100 staff split among 3 shifts were modeled individually; residents were split into 2 cohorts based on COVID-19 diagnosis. Data were analyzed from September to October 2020. Exposures: In the resident cohorting intervention, residents who had recovered from COVID-19 were moved back from the COVID-19 (ie, infected with SARS-CoV-2) cohort to the non-COVID-19 (ie, susceptible and uninfected with SARS-CoV-2) cohort. In the immunity-based staffing intervention, staff who had recovered from COVID-19 were assumed to have protective immunity and were assigned to work in the non-COVID-19 cohort, while susceptible staff worked in the COVID-19 cohort and were assumed to have high levels of protection from personal protective equipment. These interventions aimed to reduce the fraction of people's contacts that were presumed susceptible (and therefore potentially infected) and replaced them with recovered (immune) contacts. A secondary aim of was to evaluate cumulative incidence of SARS-CoV-2 infections associated with 2 types of screening tests (ie, rapid antigen testing and polymerase chain reaction [PCR] testing) conducted with varying frequency. Main Outcomes and Measures: Estimated cumulative incidence proportion of SARS-CoV-2 infection after 3 months. Results: Among the simulated cohort of 100 residents and 100 staff members, frequency and type of testing were associated with smaller outbreaks than the cohorting and staffing interventions. The testing strategy associated with the greatest estimated reduction in infections was daily antigen testing, which reduced the mean cumulative incidence proportion by 49% in absence of contact-targeted interventions. Under all screening testing strategies, the resident cohorting intervention and the immunity-based staffing intervention were associated with reducing the final estimated size of the outbreak among residents, with the immunity-based staffing intervention reducing it more (eg, by 19% in the absence of testing) than the resident cohorting intervention (eg, by 8% in the absence of testing). The estimated reduction in transmission associated with these interventions among staff varied by testing strategy and community prevalence. Conclusions and Relevance: These findings suggest that increasing the frequency of screening testing of all residents and staff, or even staff alone, in nursing homes may reduce outbreaks in this high-risk setting. Immunity-based staffing may further reduce spread at little or no additional cost and becomes particularly important when daily testing is not feasible.

Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil.

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness.

Genomics and epidemiology of a novel SARS-CoV-2 lineage in Manaus, Brazil.

Cases of SARS-CoV-2 infection in Manaus, Brazil, resurged in late 2020, despite high levels of previous infection there. Through genome sequencing of viruses sampled in Manaus between November 2020 and January 2021, we identified the emergence and circulation of a novel SARS-CoV-2 variant of concern, lineage P.1, that acquired 17 mutations, including a trio in the spike protein (K417T, E484K and N501Y) associated with increased binding to the human ACE2 receptor. Molecular clock analysis shows that P.1 emergence occurred around early November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.4-2.2 times more transmissible and 25-61% more likely to evade protective immunity elicited by previous infection with non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness.

Using viral load and epidemic dynamics to optimize pooled testing in resource-constrained settings.

Virological testing is central to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) containment, but many settings face severe limitations on testing. Group testing offers a way to increase throughput by testing pools of combined samples; however, most proposed designs have not yet addressed key concerns over sensitivity loss and implementation feasibility. Here, we combined a mathematical model of epidemic spread and empirically derived viral kinetics for SARS-CoV-2 infections to identify pooling designs that are robust to changes in prevalence and to ratify sensitivity losses against the time course of individual infections. We show that prevalence can be accurately estimated across a broad range, from 0.02 to 20%, using only a few dozen pooled tests and using up to 400 times fewer tests than would be needed for individual identification. We then exhaustively evaluated the ability of different pooling designs to maximize the number of detected infections under various resource constraints, finding that simple pooling designs can identify up to 20 times as many true positives as individual testing with a given budget. Crucially, we confirmed that our theoretical results can be translated into practice using pooled human nasopharyngeal specimens by accurately estimating a 1% prevalence among 2304 samples using only 48 tests and through pooled sample identification in a panel of 960 samples. Our results show that accounting for variation in sampled viral loads provides a nuanced picture of how pooling affects sensitivity to detect infections. Using simple, practical group testing designs can vastly increase surveillance capabilities in resource-limited settings.

Estimating epidemiologic dynamics from cross-sectional viral load distributions.

Estimating an epidemic's trajectory is crucial for developing public health responses to infectious diseases, but incidence data used for such estimation are confounded by variable testing practices. We show instead that the population distribution of viral loads observed under random or symptom-based surveillance, in the form of cycle threshold (Ct) values, changes during an epidemic and that Ct values from even limited numbers of random samples can provide improved estimates of an epidemic's trajectory. Combining multiple such samples and the fraction positive improves the precision and robustness of such estimation. We apply our methods to Ct values from surveillance conducted during the SARS-CoV-2 pandemic in a variety of settings and demonstrate new approaches for real-time estimates of epidemic trajectories for outbreak management and response.

Estimating internationally imported cases during the early COVID-19 pandemic.

Early in the COVID-19 pandemic, predictions of international outbreaks were largely based on imported cases from Wuhan, China, potentially missing imports from other cities. We provide a method, combining daily COVID-19 prevalence and flight passenger volume, to estimate importations from 18 Chinese cities to 43 international destinations, including 26 in Africa. Global case importations from China in early January came primarily from Wuhan, but the inferred source shifted to other cities in mid-February, especially for importations to African destinations. We estimate that 10.4 (6.2 - 27.1) COVID-19 cases were imported to these African destinations, which exhibited marked variation in their magnitude and main sources of importation. We estimate that 90% of imported cases arrived between 17 January and 7 February, prior to the first case detections. Our results highlight the dynamic role of source locations, which can help focus surveillance and response efforts.

Using viral load and epidemic dynamics to optimize pooled testing in resource constrained settings.

Extensive virological testing is central to SARS-CoV-2 containment, but many settings face severe limitations on testing. Group testing offers a way to increase throughput by testing pools of combined samples; however, most proposed designs have not yet addressed key concerns over sensitivity loss and implementation feasibility. Here, we combine a mathematical model of epidemic spread and empirically derived viral kinetics for SARS-CoV-2 infections to identify pooling designs that are robust to changes in prevalence, and to ratify losses in sensitivity against the time course of individual infections. Using this framework, we show that prevalence can be accurately estimated across four orders of magnitude using only a few dozen pooled tests without the need for individual identification. We then exhaustively evaluate the ability of different pooling designs to maximize the number of detected infections under various resource constraints, finding that simple pooling designs can identify up to 20 times as many positives compared to individual testing with a given budget. We illustrate how pooling affects sensitivity and overall detection capacity during an epidemic and on each day post infection, finding that sensitivity loss is mainly attributed to individuals sampled at the end of infection when detection for public health containment has minimal benefit. Crucially, we confirm that our theoretical results can be accurately translated into practice using pooled human nasopharyngeal specimens. Our results show that accounting for variation in sampled viral loads provides a nuanced picture of how pooling affects sensitivity to detect epidemiologically relevant infections. Using simple, practical group testing designs can vastly increase surveillance capabilities in resource-limited settings.

Test sensitivity is secondary to frequency and turnaround time for COVID-19 screening.

The COVID-19 pandemic has created a public health crisis. Because SARS-CoV-2 can spread from individuals with presymptomatic, symptomatic, and asymptomatic infections, the reopening of societies and the control of virus spread will be facilitated by robust population screening, for which virus testing will often be central. After infection, individuals undergo a period of incubation during which viral titers are too low to detect, followed by exponential viral growth, leading to peak viral load and infectiousness and ending with declining titers and clearance. Given the pattern of viral load kinetics, we model the effectiveness of repeated population screening considering test sensitivities, frequency, and sample-to-answer reporting time. These results demonstrate that effective screening depends largely on frequency of testing and speed of reporting and is only marginally improved by high test sensitivity. We therefore conclude that screening should prioritize accessibility, frequency, and sample-to-answer time; analytical limits of detection should be secondary.

Genomics and epidemiology of the P 1 SARS-CoV-2 lineage in Manaus, Brazil

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness.