key: cord-0295043-o6hjt2ip
authors: Carvalho, Marilia Sa; Bastos, Leonardo Soares; Fuller, Trevon; Cruz, Owaldo Gonçalves; Damasceno, Luana; Calvet, Guilherme; Resende, Paola Cristina; Smith, Chris; Whitworth, Jimmy; Siqueira, Marilda; Brasil, Patricia
title: Incidence of SARS-CoV-2 over four epidemic waves in a low-resource community in Rio de Janeiro, Brazil: A prospective cohort study
date: 2022-08-31
journal: The Lancet Regional Health - Americas
DOI: 10.1016/j.lana.2022.100283
sha: 0c5302b66dd15a95cb8fc7924910293ca7421202
doc_id: 295043
cord_uid: o6hjt2ip

Background Incidence rates of SARS-CoV-2 infections in low-resource communities can inform vaccination strategies and non-pharmaceutical interventions (NPIs). Our objective was to estimate incidence over four epidemic waves in a slum in Rio de Janeiro, a proxy for economically deprived areas in the Global South. Methods Prospective cohort of children and household contacts screened for SARS-CoV-2 by PCR and serology (IgG). The incidence density of PCR positive infections estimated for each wave - the first wave, Zeta, Gamma and Delta - was compared to an index combining NPIs and vaccination coverage. Findings 718 families and 2501 individuals were enrolled, from May 2020 to November 2021. The incidence density of SARS-CoV-2 infection due to the first wave was 2, 3 times that of the other waves. The incidence among children was lower than that of older participants, except in later waves, when vaccination of the elderly reached 90%. Household agglomeration was significantly associated with incidence only during the first wave. Interpretation The incidence of infection greatly exceeded rates reported in similar cohorts. The observed reduction in incidence in the elderly during the Delta variant wave, in spite of the rollback of NPIs, can be attributed to increased vaccine coverage. The high incidence in young people reinforces the importance of vaccination in this age group, a policy that has yet to receive the full support of some sectors of society. Funding UK Medical Research Council, Foundation for the Advancement of Science of the State of Rio de Janeiro, National Council for Scientific and Technological Development.

Incidence of SARS-CoV-2 over four epidemic waves in a low-resource community in Rio de Janeiro, Brazil: A prospective cohort study Introduction Epidemiological studies of SARS-CoV-2 incidence are valuable from a public health standpoint. Evaluating rates of new infections in a cohort can provide data to inform vaccination strategies and non-pharmaceutical interventions (NPIs). As discussed below, while many studies have estimated the incidence of SARS-CoV-2 in cohorts in Europe and North America 1−4 relatively little is known about incident cases of SARS-CoV-2 in the Global South.

An important development in the SARS-CoV-2 pandemic has been the emergence of variants of concern (VOCs) that may result in higher rates of transmission or reinfection than the first wave lineages. 5, 6 A number of studies have estimated SARS-CoV-2 incidence over one particular epidemic wave. 1−4 However, there is no data about rates of new SARS-CoV-2 cases in lowincome, high population density areas during the period of circulation of different SARS-CoV-2 lineages. Little is known about the extent to which, if any, the incidence of SARS-CoV-2 infection changes from one SARS-CoV-2 epidemic wave to another in such areas, and if so, why.

The first SARS-CoV-2 lineage that became predominant in Rio de Janeiro in the austral summer of 2020 was B.1.1.33. 7 Subsequently, the city experienced waves associated with the Variant of Interest (VOI) Zeta (P2) 8 ; the VOCs Gamma (P.1/P.1.*), Delta (B.1.617.2/ AY.*) and Omicron (B.1.1.529) (genomic surveillance data available from: http://www.genomahcov.fiocruz.br/ dashboard-en/; epidemiologic surveillance available from: https://doi.org/10.5281/zenodo.6335536). The present study examines from the first wave through the Delta wave, prior to the introduction of the Omicron VOC.

To control COVID-19, Rio de Janeiro implemented NPIs including the closure of schools, cancellation of public events during the New Year's Eve and Carnival seasons, and mandatory use of face masks, among others. Economic support for low-income families was limited to approximately US$100 per month from April 2020 to August 2020, and half this value for the next four months. This support was made available again beginning in April 2021 for the following five months. Schools for children ages 5−15 reopened partially and gradually beginning in August 2021 and fully in November 2021. Vaccination began on 20 January 2021 for the elderly and health professionals and gradually expanded to other age groups. Vaccination of adolescents 12−17 years old began in August 2021.

Like 60 other countries in the Global South, Brazil's economy is characterized by elevated indices of income inequality. 9 Nevertheless, when formulating NPIs, governments often adopt a "one size fits all" approach without accounting for the effect of income on the capacity to comply with COVID-19 control measures. 10 There may be far fewer resources available to sustain households in low-income communities, affecting the population's willingness to comply with NPIs. Understanding the main drivers of SARS-CoV-2 incidence in low-income communities and how they may have changed over the two years of the pandemic may provide valuable information to decision-makers for formulating future NPIs and immunization priorities. Such policies are likely to be necessary for the foreseeable future, as SARS-COV-2 may persist in human populations for a prolonged period.

The objective of this study was to estimate the incidence of SARS-CoV-2 infection in a cohort of children and household contacts in a slum in Rio de Janeiro. We calculated the incidence density for each of the four pandemic waves between April 2020 and September 2021, accounting for age and household clustering.

The present investigation, which was built upon our accompaniment of previous cohorts in this community, 11−14 is part of an ongoing prospective SARS-CoV-2 study. 15 We offered study participation to all household contacts of children aged 13 and under who were seen at a primary health care centre for any reason including routine vaccinations and preventative care. The centre provides primary care to 4300 children in this age group per month. After enrolment, we collected samples in each household on study days 1 (first visit), 14 and 28, every three months during the first year, and twice in the second year. Enrolment took place from the

Evidence before the study We searched PubMed and Web of Science for publications on SARS-CoV-2 incidence density, published between 1 March 2020 and 26 February 2022. We used the keywords "incidence density" and "SARS-CoV-2" and searched for articles in English, French, Spanish and Portuguese. Cohort studies have estimated the incidence density of SARS-CoV-2 infection in families in the US and health care workers in India and Ethiopia, for the most part during the first wave of the COVID-19 pandemic. However, little is known about the incidence density of SARS-CoV-2 infection in families in lowresource, high household density communities, where compliance with non-pharmaceutical interventions (NPIs) such as lockdowns is unfeasible due to food insecurity.

The objective of this study was to estimate the incidence of SARS-CoV-2 infection in a cohort of children and household contacts for each of the four waves of the pandemic between April 2020 and September 2021, accounting for age and household clustering. Our study participants reside in a community in Rio de Janeiro with high levels of poverty and violence. Ours is one of the few studies to characterize the incidence of SARS-CoV-2 in this setting. We evaluated incidence in relation to the stringency indexa scale for evaluating the adoption of NPIs in the cityand vaccination coverage by age group.

The incidence of SARS-CoV-2 in families in a lowresource community in Rio de Janeiro was higher than in household cohorts in the US and among healthcare workers in India and Ethiopia. While the incidence of SARS-CoV-2 in children in our cohort was generally lower than that of older participants, the incidence in children increased later in the epidemic. In Brazil, immunization of children against SARS-CoV-2 continues to be opposed by some sectors of civil society and the government. Our findings highlight how important it is for the Brazilian government to stand behind the policy of vaccinating children despite political opposition.

beginning of the pandemic on May 15, 2020, to November 14, 2021.

Nasopharyngeal swabs and saliva were screened for SARS-CoV-2 by real-time RT-PCR to amplify the E gene and the RdRp region of the Orf1ab gene of SARS-CoV-2 using the SARS-CoV-2 E/RP molecular kit (Bio-Manguinhos, Rio de Janeiro, Brazil). Cycle thresholds (Ct) less than 40 were classified as positive. Serum tests were performed to screen for IgG antibodies by a Microparticle Chemiluminescent Immunoassay (CMIA), targeting the S gene (Abbott Laboratories, Abbott Park, IL, USA). SARS-CoV-2 RT-PCR and IgG serology assays were performed on all study participants at every visit according to the manufacturer's instructions.

The study was approved by institutional review boards at the Oswaldo Cruz Foundation and the London School of Tropical Medicine and Hygiene.

A new SARS-CoV-2 infection was defined as a positive RT-PCR test. Seropositive results were not considered as new infections, as the precise date of seroconversion was not available. For the first three months of the study, however, as there was no previous contact with SARS-CoV-2, seropositive cases were also considered incident infections.

Incidence density was defined as the number of new SARS-CoV-2 cases per 1000 person weeks of exposure. In each period, the disease-free observation period for each participant was defined as the duration of observation during the period in question plus two additional weeks. The two extra weeks were included to account for the maximum interval for SARS-CoV-2 transmission, and the detection of infection by the study team. The at-risk time was truncated at the beginning and the end of each period if the follow-up spanned more than one period.

Households in which no individual was examined during an entire wave were considered lost to follow-up after eight unsuccessful contact attempts by telephone and smartphone messaging apps. Individuals in each household who missed one or more study visits were excluded only if the study team was unable to contact the participant during an entire wave. During the study period, young people and adults were eligible for vaccination, following a schedule defined by the municipal health department. The first age group eligible for vaccination was older adults. Children under 12 were not eligible to be vaccinated during the study period. Therefore, the at-risk time was truncated after the first dose.

To define the start and end of each epidemic wave, we determined the number of cases per week in wave t. The week with the lowest number of cases during a wave t was defined as the beginning of the next wave, (t+1). During each wave t, multiple lineages were detected by genomic surveillance. We named each wave t based on the predominant SARS-CoV-2 lineage during the wave.

Age was included in the model divided into categories based on the age groups targeted by vaccination campaigns in Rio de Janeiro: 0−4, 5−11, 12−17, 18−59, and sixty and older. In addition, we evaluated household clustering, defined as the ratio of the number of persons living in the household and the number of rooms, as the most relevant variable to estimate the incidence in the cohort.

Stringency index. The Stringency Index was developed to track and compare policy responses around the world, combining indicators in four domains: Containment and Closure, Economic Response, Health Systems, and Vaccine Policies. 16 This indicator ranges from 0 representing no control measures to 100, which corresponds to complete lockdown. The index does not correlate specifically with incidence. Instead, it is an overall indicator of health authorities' responses to factors including local trends in incidence, mortality, hospital capacity, and economic demands.

In Brazil, the variables were available for the capitals of each state. The data set is available from https:// www.bsg.ox.ac.uk/research/research-projects/covid-19government-response-tracker. 17 Vaccination coverage. The Brazilian Health Ministry data on vaccination by age group was obtained by tabulating the open-source data from the COVID-19 vaccination information system (SI-PNI) available from https://opendatasus.saude.gov.br/dataset/covid-19-vacinacao as of February 2, 2022.

Genomic data. All genomic data from Rio de Janeiro evaluated in this study is available in the EpiCoV database at GISAID (www.gisaid.org) produced by the COVID-19 Fiocruz Genomic Surveillance Network and other genomic initiatives in Brazil. Details about genomic surveillance can be found at http://www.genomahcov.fiocruz.br/dashboard-en/.

Statistical models. As the response variable is infection vs non-infection, we used a logistic model to estimate the odds of each variable. In addition to age, which was included a priori, we also evaluated the significance of variables at the household level. Since the models were very similar according to the Watanabe−Akaike information criterion (WAIC), 18 in addition to age, we chose a variable representing household clustering (persons/ room), which we considered the most relevant indicator of housing conditions in the study site. The individuals who were included in the odds ratio calculation and who had sufficient follow-up time were retained in the incidence density calculation.

The probability of a participant being an incident infection is given by the expected value of the aforementioned logistic model. The duration of follow-up varies from participant to participant. Therefore, the estimated incidence density is the sum of the probability of being a SARS-CoV-2 positive or negative case, given the duration of follow-up of all the participants in each wave. The global incidence of infection in each wave is the expected number of infections including both SARS-CoV-2 detections by RT-PCR and seroconversion, divided by the number of people exposed in each wave. In addition, we tested the inclusion of two random effects to account for the expected dependence between individuals living in the same household and for the unstructured individual variability to account for unknown confounders (Supplemental Material).

Role of the funding source. The funders had no role in study design, data collection, data analysis, interpretation, writing of the report or decision to submit. Table 1 summarizes the sociodemographic profile of all participants according to the wave in which they were recruited. The interquartile range of age was similar across the epidemic waves. However, the median age was significantly lower in the Delta wave and higher in the Zeta wave (Supplemental Material). The age group that comprised the highest proportion of the cohort was working-age adults, which along with adolescents, represented 58% of those recruited. The majority of participants were women (60%). The total number of exams was greater during the first wave and similar for the following ones. Positive PCR ranged between 12 and 17% between the first three waves and decreased substantially in the Delta wave.

A total of 718 families were enrolled in the study; 40.5% were recruited in the first epidemic wave and approximately 20% in each of the subsequent ones. Concerning families recruited during the Zeta, Gamma and Delta waves, the proportion of families retained in subsequent waves was similar to the pattern observed among families recruited during the first wave. By the end of the study, 93.2% of the participants aged 12 and older who were eligible to be vaccinated had received at least one dose of a vaccine. The at-risk period and outcomes for each participant (RT-PCR and serology tests, and immunization) are summarized in Figure 1 and Supplementary Figure S3 . Figure 2 displays the four COVID-19 waves in Rio de Janeiro superimposed on a graph of the Stringency Index. During the initial three months, at the beginning of the epidemic, the value of the Stringency Index was indicative of lock-down. However, by the end of the first wave, these measures were relaxed, and the index oscillated between 60 and 70 points through the end of the Gamma wave. During the Delta wave, relatively few NPIs were implemented but the stringency index ticked up due to the rollout of vaccination. Four months after the start of vaccination, complete vaccination coverage among the elderly reached 60% (Figure 3 ). By the Delta wave, complete vaccination coverage reached 90% in this group. Vaccination of working-age adults coincided with the beginning of the Delta wave. By the end of the study period, complete vaccination coverage of workingage adults reached 90%. The vaccination of adolescents began in September 2021. Children 5−11 years of age were not vaccinated during the study period.

Sixteen models were evaluated (Supplementary Table  S2 ). The model with the best WAIC score included just the family random effect. Figure 4 shows the odds ratio and credible intervals (0.025 and 0.975 quantiles) for the variables age group and household clustering.

Compared to the working-age population (19 −59 years of age), the odds ratio for children 0−4 years of age had a significant protective effect during the Zeta Articles wave (odds ratio: 0.36) and there was a trend toward significance during the Delta wave (odds ratio: 0.5). A similar pattern was observed for people 60 and older. Being 60 or older had a significant protective effect during the Zeta wave (odds ratio: 0.34) and there was a trend toward a significant protective effect during the Delta wave (odds ratio: 0.25). In addition to these two age groups, being 5−11 years of age had a significant protective effect against being a SARS-CoV-2 case during the Zeta wave (odds ratio: 0.34). During the first wave and the Gamma waves, none of the age groups had a significantly increased or decreased risk of being a SARS-CoV-2 case. During the first wave, the odds of being a SARS-CoV-2 case almost tripled with the number of persons per room in the household. For the subsequent waves, the effect of household clustering decreased log-linearly with the odds of being a SARS-CoV-2 case. During the Delta wave, there was a trend toward a significant effect of household clustering. Figure 5 depicts the incidence density of SARS-CoV-2 infection adjusted by age and epidemic wave. The mean incidence density averaged across age groups during the first wave was 52.5 cases per 1000 person-weeks, which was 2.6 times greater than the incidence during the Zeta wave (19.7 per 1000 person-weeks) and 2.2 times that of the incidence during the Gamma wave (23.7 per 1000 person-weeks). The incidence density during the Delta wave was 35.6 per 1000 person-weeks. We carried out a sensitivity analysis to assess how different durations of follow-up influenced the ranking of the waves with respect to incidence. The results indicated that the aforementioned patterns -highest incidence during the first wave followed by the Delta, Gamma and Zeta waves -were robust across different durations of follow-up (Supplementary Material).

Comparing age groups, the incidence density among children was generally lower than that of older participants throughout the study. For example, during the Delta wave, the incidence density among participants 60 and older was significantly lower than in younger age groups. Furthermore, during the Delta wave incidence among children 5−11 years old was similar to that of working-age adults.

The incidence of SARS-CoV-2 infection in children and their household contacts in Rio de Janeiro greatly exceeded rates reported in similar cohorts in the Global North. The incidence of SARS-CoV-2 in our study population during the first wave was approximately eight times higher than the incidence reported for familybased cohorts in the US. 1,2 Among the factors that may have contributed to the difference in incidence rates is the geographic scope of recruitment. Since the geographic areas from which families were recruited in the US were larger than that of our study, the cohorts followed in the US may have been more diverse than ours. On the other hand, few studies have investigated a severely economically deprived region like our study site, in which participants were at risk of food insecurity, and where violence affects their daily lives. It should be noted that fieldwork under these conditions was only possible because the interviewers were residents of the region hired during our cohort studies. 11−15 Residents of low-resource communities may have higher exposure to SARS-CoV-2 than residents of highresource communities because of household crowding and lower capacity to adhere to shelter-in-place orders due to food insecurity. 19 Furthermore, due to COVID-19 denialism by prominent elected officials, there was contradictory messaging about NPIs and immunization, which further complicated compliance with COVID-19 control policies. With respect to other studies in developing countries, the incidence density in our cohort was similar to rates reported for healthcare workers in India and Ethiopia. 20,21 These frontline workers likely had higher exposure to SARS-CoV-2 than our cohort, which would have tended to increase incidence.

In our cohort, the increase in the relative risk of SARS-CoV-2 infection during the first wave and Delta periods was two to three times that of the Zeta and Gamma periods. We can conjecture that relative risk was highest during the first wave because the population was completely non-adapted to the novel virus. The association between exposure to the Delta variant and higher relative transmission compared to other VOCs observed in our cohort is supported by experimental studies. 22 The risk may have been higher during the Delta period because mutations in the SARS-CoV-2 Spike protein receptor-binding domain associated with the VOC Delta may have increased the transmissibility of the virus relative to previous lineages. The observed reduction in incidence in the elderly during the Delta variant wave, despite the rollback of NPIs, may have been attributable to increased vaccine coverage. A number NPIs were announced during this wave but not fully implemented; for instance, bars and restaurants were supposed to require proof of complete vaccination, but most did not. Masks were also not regularly required in stores and shopping centres.

In addition, we found that the risk of infection increased with household density. When we considered the number of persons per room in the household, the risk of being infected almost tripled in the first wave compared to the subsequent periods ( Figure 3 ). One possible explanation for this difference is that residents were infected earlier in the epidemic in households with particularly high density. Subsequently, in the Zeta, Gamma, and Delta waves, the credible intervals of the variables persons/room overlapped with 1 indicating that the variable was neither associated with increased risk of being a SARS-CoV-2 case nor with decreased risk.

We evaluated two indicators estimated at the ecological level: The Stringency Index and vaccination coverage. The Stringency Index graph depicts how health authorities responded to factors such as SARS-CoV-2 incidence in the community. Since there may have been Articles a delay between a change in incidence and the implementation of health authorities' response, the index may not be perfectly synchronized with patterns of incidence, as is the case for the Gamma wave. However, the index reflects overall trends in incidence. It is well-established that compliance with NPIs and vaccination is protective against SARS-CoV-2 infection. Studies in other municipalities have demonstrated that the Stringency Index captures the population's response to NPIs at the city scale. 23, 24 The community we studied is located just 15 min from downtown by public transport and most residents commute to other neighbourhoods. Thus, the pattern we observe is very similar to other low-income communities in Rio de Janeiro. According to observational data on vaccination coverage on the scale of the city of Rio de Janeiro, there was a trend of increasing coverage during the study period. It is plausible to infer that this increase in vaccination was made possible because there was a growing acceptance of vaccination among individuals in the age groups that were eligible to be immunized. We can also reasonably surmise that with increasing vaccination coverage at the city scale, there was a decline in the risk that the virus would be introduced into any particular community or household.

In our study, the incidence of SARS-CoV-2 in children was generally lower than that of older participants, in accordance with a previous study of this population. 15 However, the incidence in children increased later in the epidemic, underscoring the importance of making vaccines fully available to the young. The pattern of SARS-CoV-2 incidence by age group observed in this cohort can provide potentially valuable information to inform health policies. The World Health Organization endorsed the vaccination of children aged 5−11 with the BNT162b2 vaccine in November 2021. 25 However, in Brazil, the vaccination of children in this age group continues to be opposed by certain groups in civil society and the government. The immunization of children aged 5−11 only began in Rio de Janeiro on 17 January 2022, after our study ended. It is reasonable to conclude that the high incidence of SARS-CoV-2 in children aged 5−11 in our cohort during the Delta wave (July to November 2021) was attributable at least in part to the fact that these children were not eligible for vaccination. Our findings highlight that it is crucial that the Brazilian government stand behind the policy of vaccinating children despite political opposition.

The main strengths of this study are the longitudinal design, which followed participants prospectively, and the eligibility criteria, which included participants regardless of symptoms, providing an accurate measure of incidence. These inclusion criteria and active follow up could provide an unbiased estimate of incidence in other similar populations in Latin America. Although communities with such high levels of poverty and violence are common in lower-income countries, the incidence of SARS-CoV-2 in this setting has seldom been studied due to the challenges of conducting fieldwork in this context. Ours is one of the few studies to characterize the incidence of SARS-CoV-2 in this setting. Among the limitations of the study is that households without children were not included.

During the final epidemic wave analysed, the Delta VOC, we identified a proportional increase in incidence density in the young. This finding strongly suggests that younger age groups need to be vaccinated, which runs counter to antivaccine messaging by some government officials. As Brazil has a long tradition of mass vaccination campaigns and Brazilians largely trust the public healthcare system, 26 reducing SARS-CoV-2 incidence appears to be an attainable goal through increased access to vaccines and information campaigns publicizing the benefits of vaccinating children. 

Interested researchers should contact the corresponding author. Requests will be assessed in accordance with best practices in cohort studies to ensure that data sharing would preserve confidentiality and be compatible with the terms of the Informed Consent Form used in this study.

Incidence rates, household infection risk, and clinical characteristics of SARS-CoV-2 infection among children and adults in utah

Incidence, clinical characteristics, and risk factors of SARS-CoV-2 infection among pregnant individuals in the United States

Incidence of SARS-CoV-2 infection among health care personnel, first responders, and other essential workers during a prevaccination COVID-19 surge in Arizona

SARS-CoV-2 Seropositivity among US marine recruits attending basic training, United States, spring-fall 2020

SARS-CoV-2 variants of interest and concern naming scheme conducive for global discourse

Variant analysis of SARS-CoV-2 genomes

Evolutionary dynamics and dissemination pattern of the SARS-CoV-2 lineage B.1.1.33 during the early pandemic phase in Brazil

Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil

Poverty and Shared Prosperity 2020: Reversals of Fortune. Poverty and Shared Prosperity 2020: Reversals of Fortune. Washington

Physical activity and sociodemographic profile of Brazilian people during COVID-19 outbreak: an online and cross-sectional survey

A populational-based birth cohort study in a low-income urban area in Rio de Janeiro, Brazil: implementation and description of the characteristics of the study. Cad Sa ude P ublica

Zika virus infection in pregnant women in Rio de Janeiro

Why did ZIKV perinatal outcomes differ in distinct regions of Brazil? An exploratory study of two cohorts

A prospective, multicentre, cohort study to assess the incidence of dengue illness in households from selected communities in Brazil

SARS-CoV-2 infection dynamics in children and household contacts in a slum in Rio de Janeiro

Variation in Government Responses to COVID-19 (Version 12.0). Blavatnik School of Government Working Paper

Government responses and COVID-19 deaths: global evidence across multiple pandemic waves

Asymptotic equivalence of bayes cross validation and widely applicable information criterion in singular learning theory

High prevalence of food insecurity, the adverse impact of COVID-19 in Brazilian favela

Effectiveness of BBV152 Vaccine against SARS-CoV-2 Infections, hospitalizations and deaths among healthcare workers in the setting of high delta variant transmission

Seroepidemiology and modelbased prediction of SARS-CoV-2 in Ethiopia: longitudinal cohort study among front-line hospital workers and communities

SARS-CoV-2 B.1.617.2 Delta variant replication and immune evasion

Assessing the nationwide impact of COVID-19 mitigation policies on the transmission rate of SARS-CoV-2 in Brazil

Stringency of containment and closures on the growth of SARS-CoV-2 in Canada prior to accelerated vaccine roll-out

World Health Organization. Interim recommendations for use of the Pfizer−BioNTech COVID-19 vaccine, BNT162b2, under Emergency Use Listing

Running away from the jab: factors associated with COVID-19 vaccine hesitancy in Brazil

None declared.

Supplementary material associated with this article can be found in the online version at doi:10.1016/j. lana.2022.100283.