Design of a population-based longitudinal cohort study of SARS-CoV-2 incidence and prevalence among adults in the San Francisco Bay Area.

PURPOSE: We describe the design of a longitudinal cohort study to determine SARS-CoV-2 incidence and prevalence among a population-based sample of adults living in six San Francisco Bay Area counties.
METHODS: Using an address-based sample, we stratified households by county and by census-tract risk. Risk strata were determined by using regression models to predict infections by geographic area using census-level sociodemographic and health characteristics. We disproportionately sampled high and medium risk strata, which had smaller population sizes, to improve precision of estimates, and calculated a desired sample size of 3400. Participants were primarily recruited by mail and were followed monthly with PCR testing of nasopharyngeal swabs, testing of venous blood samples for antibodies to SARS-CoV-2 spike and nucleocapsid antigens, and testing of the presence of neutralizing antibodies, with completion of questionnaires about socio-demographics and behavior. Estimates of incidence and prevalence will be weighted by county, risk strata and sociodemographic characteristics of non-responders, and will take into account laboratory test performance.
RESULTS: We enrolled 3842 adults from August to December 2020, and completed follow-up March 31, 2021. We reached target sample sizes within most strata.
CONCLUSIONS: Our stratified random sampling design will allow us to recruit a robust general population cohort of adults to determine the incidence of SARS-CoV-2 infection. Identifying risk strata was unique to the design and will help ensure precise estimates, and high-performance testing for presence of virus and antibodies will enable accurate ascertainment of infections.

Introduction

By the beginning of May 2021, 32.6 million people in the U.S. had been were reported as infected with SARS-CoV-2, of whom 579,634 had died [1]. These numbers under-represent the total burden of infection due to incomplete testing and the lower likelihood of asymptomatic persons coming to clinical attention. Accurate data on the extent of infection, even as vaccines are rolled out, are critical to understanding continued transmission and informing ongoing mitigation efforts.

Numerous cross-sectional studies aimed at determining population-levels of infection have been conducted in the U.S., including in Chicago, New York, Indiana, Georgia, California [2], [3], [4], [5], [6], [7], [8], [9], as well as country-wide and internationally [10], [11], [12], [13], [14], [15], [16], [17], [18]. Approaches to determining the prevalence of infection have also involved testing of remnant blood samples [19], [20], [21] including from dialysis patients [22,23]. The Centers for Disease Control and Prevention estimated an overall prevalence of infection in the U.S. of 14% based on data from community-based studies and the testing of remnant blood specimens from 10 sites nationwide, coupled with multipliers based on case reports [9]. Seroprevalence estimates have varied widely, however, due to differences in sampling approaches, the target population, and the dates during which surveys were implemented [18,24]. In addition, general population estimates may not take into account the higher rates of infection among subgroups; for example, several studies have demonstrated that Latinx communities in the U.S. are more highly affected by the pandemic, likely due to occupational hazard, higher housing density and other factors [25], [26], [27], [28], [29], [30]. Errors in prevalence estimates can also occur because of imperfect antibody test performance, which can under- or overestimate actual infections [24,31,32]. This was problematic earlier in the pandemic when rapid tests with poor test performance were used in surveys [33]. A large effort is underway to obtain nationwide estimates of prevalence and incidence by mailing home-testing kits to a household probability sample in the U.S., although results from this study are not yet published [34].

We utilized a robust epidemiologic and statistical design to enroll a representative population-based sample of adults from 6 counties in the San Francisco Bay Area into a longitudinal surveillance cohort. The original aim was to obtain regional estimates of incidence and prevalence of SARS-CoV-2 to assist local public health departments, which at the time of study conception in late March 2020 were grappling to determine the trajectory of the epidemic, to identify communities most at risk and the most effective prevention methods. Additional aims included determining the association of occupation and behaviors with infection rates, the proportion of infections that were asymptomatic, COVID-19 vaccine acceptability, and the presence of circulating viral strains. This paper describes the design and methods used for sampling, enrollment, ascertainment of infection, and analysis. As this paper is focused on methods, we do not describe study results, or the characteristics of the study sample. The project, called TrackCOVID, is a collaboration of the University of California, San Francisco (UCSF), the Stanford University Health Center, and the Zuckerberg San Francisco General Hospital, with support from the local county Departments of Public Health, and funded by the Chan-Zuckerberg Initiative.

Materials and methods

Summary

We describe the design of a longitudinal cohort study used to enroll and follow a population-based sample of adults to determine the incidence and prevalence of SARS-CoV-2 infection in the Bay Area. We used a two-stage stratification sampling scheme, based on an address-based sampling frame. We first sampled by county proportionate to the number of households (HH), and then by census tract risk strata (high, medium, low) within each county. Risk strata were determined by using regression models to predict the number of cases in each census tract. Participants were primarily recruited by mail; only one randomly selected adult from each participating household was enrolled. Participants were followed monthly with SARS-CoV-2 PCR and antibody testing, and with questionnaires. Recruitment by mail began at the end of July 2020, with enrolment between August 2020 and December 2020; follow-up was completed March 31, 2021.

Study population

The total adult population of the six counties in the Bay Area was 5,321,907 based on 2019 census data, and racially and ethnically diverse (20.4% Hispanic, 31.5% Asian, 6.0% Black, and 3.6% mixed or other) [35]. Slightly less than 1/3 (30.1%) of adults were 18–34 years of age, and 19.2% were 65 years of age or older.

The targeted sample size for this study was 3400 adults, based on estimating an incidence of 5.0 cases of infection per 100 person-years (100py) (total width of 95% confidence interval [CI] = 2.2 cases/100py). Figure 1 shows estimates of precision assuming different incidence rates and sample sizes, assuming a mean follow-up time of 6 months. The target study population included persons 18 years of age or older residing in Alameda, Contra Costa, Marin, San Francisco, San Mateo, and Santa Clara counties, who were not living in congregate settings or prisons, and who did not report a prior confirmed SARS-CoV-2 infection at screening.

Figure 1

Precision (total width of 95% confidence interval) of estimates of SARS-CoV-2 incidence as a function of sample size and different incidence rates.

Stratification. We used a stratified random sampling scheme. We first sampled by county and then by modeled risk strata within each county. Sampling was based on the number of HHs rather than number of adults. The number of HHs to be sampled within each county was determined proportionate to the number of residential HHs that were listed in the Postal Delivery Sequence file of the U.S. Postal Service [36]. (The total number of HHs within the six counties was listed as 2,442,926). We then sampled by census tract risk strata (high, medium, low) within each county; strata were based on regression model-prediction of infection (described below). Given the smaller population size within medium and high-risk census tract strata, we oversampled the number of HH adults within those strata, to ensure the precision of our incidence and prevalence estimates (Table 1 ) [37].

Table 1

Number of households listed in the US Postal Service Delivery Sequency File, by county and risk strata, and the sampling fraction.

	Risk Strata
	Low		Medium		High
County	Households N	Sampling fraction	Households N	Sampling fraction	Households N	Sampling fraction
Alameda	190,570	1.00	316,127	2.16	95,905	4.36
Contra Costa	187,079	1.00	198,771	2.05	28,299	2.62
Marin	39,767	1.00	59,580	2.29	5025	6.78
San Francisco	115,300	1.00	210,431	2.30	52,148	3.70
San Mateo	18,659	1.00	165,903	2.73	90,605	5.42
Santa Clara	347,258	1.00	281,453	1.89	40,046	2.76

Risk strata were classified based on predicting the number of SARS-CoV-2 cases that could be expected to occur within each census tract. We used predicted number of cases to identify high risk areas, rather than actual cases reported to public health departments, because at the time of study initiation, widespread access to testing was not available, particularly for communities at higher risk. In addition, residents of some communities were often hesitant to seek testing, regardless of availability [25,28]. Therefore, reported infections would not reflect the actual prevalence of infection by geographic area. Likewise, we did not rely upon reported hospitalizations and/or deaths to identify high risk areas. In the Bay Area, the majority of hospitalizations were among LatinX and Black persons, disproportionate to their representation in the general population [38,39]. Thus, reasons for increased morbidity and mortality in these populations were not only related to prevalence of infection, but to co-morbid and other conditions contributing to more severe disease. Therefore, relying on hospitalizations might overestimate the levels of infection within communities.

We classified census tracts into risk strata based on predicting the number of infections using Least Absolute Shrinkage and Selection Operator (LASSO) regression [40]. Factors potentially predictive of SARS-CoV-2 risk were selected based on existing knowledge of socioeconomic and health characteristics among persons more likely to be infected. The distribution of these factors by census-tract was abstracted from data reported in the 2018 American Community Survey (ACS) [41] and the UCSF HealthAtlas [42]. We initially included 66 census-level characteristics in the model and from these, identified 27 with the highest coefficients for predicting the cumulative numbers of cases reported by census tract and provided to us by county health departments (model R2=0.50) (Appendix A). We then applied the model using these selected factors to predict the number of SARS-CoV-2 cases that would exist within each census tract. Based on the Cochrane method [37], census tracts were grouped into strata according to the predicted cases per 100,000 adult population: high risk (>457 cases/100,000 adults), medium risk (114–457 cases/100,000 adults), and low risk (<114 cases/100,000 adults).

Recruitment and enrolment

To determine the number of households that we needed to target for recruitment, we assumed that response to recruitment letters would be 9% in low-risk, 6% in medium risk, and 4% in high-risk strata, based on previous experience with mail-based recruitment. Using these response rates and the desired sample size by strata, we determined the number of HHs to be targeted. We then purchased a stratified random sample of 60,000 HH addresses derived from the US Postal Delivery File, obtained through the Marketing Systems Group (Horsham, PA) (www.m-s-g.com). Participants were primarily recruited by mail, starting in mid-July 2020. We developed letters that described the study and encouraged enrollment, and translated them into the most prevalent languages spoken in the Bay Area (English, Spanish, Chinese, Tagalog and Vietnamese). We also developed postcards in English and Spanish. Letters and postcards were mailed in successive waves to households, with each being sent at least two letters and a postcard. We monitored response by zipcode and strata, and sent additional mailers to HHs in areas where enrolment was low.

Mailers invited the adult with the next birthday to participate, and provided a link to a study -specific website that provided more detailed information and instructions on how to enroll (trackcovidbayarea.com). Mailers listed a telephone number at a health survey research firm employed to assist with the study (The Henne Group [THG] www.thehennegroup.com), which potential participants could call to speak with someone directly. Staff fluent in five languages were available to answer questions and help with screening and enrolment.

We also attempted telephone recruitment. About one-third to one-half of addresses obtained through the Postal Delivery File are linked to telephone numbers. Starting in September 2020, THG began phoning all HHs in our target sample that had an associated phone number. Telephone recruitment was conducted in the preferred language of the prospective participant. Starting in October 2020, to increase response from residents of high-risk strata, we collaborated with local community-based organizations (CBO) working in the 6 participating counties, as well as a survey team that visited selected households to directly encourage enrollment. Prior to being deployed, teams were trained in how to guide individuals through online enrollment and scheduling; team members were bi-lingual in Spanish and English. Outreach staff also provided printed information about the study, and a gift bag with hand sanitizer and cloth masks; these items were given directly to an adult in the HH, or left outside homes at which no-one answered.

Potential participants could be screened and complete an electronic consent form directly on the study website, verbally through THG on the phone, or at their first visit. They could also schedule their first visit for testing either online or by telephone.

Laboratory testing

Enrolled participants provided samples for viral detection and for antibody testing at one of 13 testing sites that were set up throughout the 6 counties for the purposes of the study. These were co-located at existing testing sites affiliated with UCSF, Stanford Health Center, county public health departments, CBOs, and private hospitals. Sites were supported either by the UCSF or Stanford study teams. All testing platforms had received FDA emergency use authorization (EUA). We performed reverse-transcriptase polymerase chain reaction (rt-PCR) testing of nasopharyngeal swab (NP) samples to identify the presence of virus, indicative of active infection, persistent shedding, or presence of viral particles [43]. We also obtained venous blood samples for testing of antibodies to different viral antigens. Details of testing platforms and performance are provided in Appendix B.

Briefly, PCR testing of swab samples was performed using several different testing platforms, depending on whether tests were performed at the Chan-Zuckerberg BioHub, San Francisco [44], UCSF [45], [46], [47] or at Stanford Health Center laboratories [48,49]. Positive PCR samples were sent to the BioHub for whole genome sequencing [50,51].

Plasma from venous blood samples collected at UCSF-supported sites was tested for the presence of IgG antibodies to SARS CoV-2 nucleocapsid (N)-protein [52]. Blood collected at Stanford-supported sites was tested for presence of IgG antibodies to SARS-CoV-2 Spike glycoprotein (S1) and the S1 Receptor Binding Domain (RBD) [49]. Positive antibody samples were cross-tested for presence of IgG at both institutions using the above methods. All samples positive for IgG to either S1, N, or both, were assayed for the presence of neutralizing antibodies at UCSF or the Vitalant Research Institute [22]. Remnant plasma and NP eluent samples from each visit are being stored at specimen banks for confirmation testing if needed and for future research.

All participants were required to register with the electronic health record system of either UCSF or Stanford, to enable processing and reporting of laboratory tests; positive PCR test results were automatically reported to California's electronic disease reporting system [53]. Persons who had a positive PCR or a confirmed antibody test were contacted through their electronic health record system, and were also called by a study physician who counseled them on isolation guidelines, and referred them as necessary for health care and/or support services.

Questionnaire

At baseline, participants completed a detailed questionnaire (Appendix C). Socio-demographic questions included gender identity, age, race, ethnicity, education, income, occupation, household size, and numbers of hours/week working outside the home. Behaviors potentially related to the risk of infection were addressed by asking questions about the proportion of time wearing a mask outside the home in the last month, level of avoidance of people not in the home, travel outside the state, and any known exposure to someone with COVID-19. We also asked about COVID-related symptoms in the previous month and in the last 24 hours, and chronic health conditions including diabetes, obesity, immunologic compromise, among others. Questions about occupation were asked according to the Council of State and Territorial Epidemiologists Occupation Health Subcommittee recommendations [54]. Starting in December 2020, supplemental questions were added inquiring about receipt of COVID-19 vaccination, including date(s) and type of vaccine. The questionnaires were available in the five targeted languages, could be completed electronically through the study website, by phone through THG, or at a testing site with assistance from study staff.

Reimbursement

A $25 gift card was provided as reimbursement at each visit, with a one-time increase to $100 in November and December 2020 to boost enrollment and improve retention. Assistance with the cost of transportation was provided as requested.

Follow-up visits

Participants were followed monthly and were asked to complete a short questionnaire about behavior, symptoms in the last month, exposure to someone with COVID-19, and any change in health status. An NP swab and venous blood samples were also obtained. COVID-19 vaccinations were rolled out in California in a staggered fashion beginning in late December 2020. We continued to follow vaccinated individuals with PCR and antibody testing to identify vaccine breakthrough infections.

COVID-19 protections

We instituted precautions to reduce the risk of SARS-CoV-2 transmission for participants and study staff. All participants were asked to wear a face mask when arriving at the testing site, except for when an NP swab sample was being obtained. Staff collecting NP and/or venous blood specimens wore face masks, eye shields, gloves and gowns; gloves were changed between participants. Hand sanitizer was available. Most of the testing sites were outside under tents and therefore with adequate ventilation. So that participants could avoid public transportation, reimbursement was provided for travel and/or parking as requested.

Primary outcomes and statistical analysis

Our primary outcomes are prevalence and incidence of SARS-CoV-2 infection. A prevalent case is defined as someone who had either a positive PCR test and/or a confirmed antibody test at their baseline visit. A confirmed positive antibody test indicative of infection is defined as having at least 2 of 3 antibodies detected (anti-S1 anti-N, or neutralizing). An incident case is defined as someone who has a positive PCR test or a confirmed antibody test without evidence of infection at baseline or the prior visit. An infection in a vaccinated or partially vaccinated person is defined as having a positive PCR test, and/or a positive anti-nucleocapsid antibody test. Anti-spike and neutralizing antibodies can be generated in response to the vaccine and were therefore were not considered to confirm a true infection [55].

We will use a weighted binomial approach to estimate baseline prevalence with 95% CI. To estimate incidence (new SARS CoV-2 infections/100py), we will use weighted Poisson regression with person-days in the model as an offset. Persons who have evidence of a prevalent infection at their baseline visit will not be included in incidence calculations. For participants with previously negative test results, and who have a confirmed antibody test on a follow-up visit without a positive PCR test, the date of infection will be imputed as the mid-point between the last negative test and the first positive antibody test. Individuals will be censored if they meet the definition of a new infection, die, withdraw from the study, or are lost to follow-up

Weights will be estimated to account for stratification, the probability of being selected based on the number of adults in the household, and differential non-response and coverage by age, education, gender, and race/ethnicity [56]. The latter relies on raking methods [57] that will be applied after determining key differences in socio-demographic characteristics between the weighted sample and the general population based on 2019 ACS data [41,58]. The standard error used in the confidence interval estimates will be obtained via bootstrapping procedures to account for uncertainty of sample size, weight estimation, and positive percent agreement (PPA) and negative percent agreement (NPA) of testing platforms.

Sensitivity analysis

We will calculate estimates of incidence and prevalence excluding persons after a first dose of a COVID-19 vaccine, and also calculate them including vaccinated individuals; the latter will provide an estimate of the general incidence in the population during vaccine uptake. Estimates of prevalence will be adjusted for the laboratory assay performance (using bootstrapping methods according to Sempos and Tian [59]). These methods will also be applied to incidence estimates. We will use bootstrapping to estimate the variance of incidence and prevalence estimates to account for the uncertainty of weight estimation and sample size.

Ethical considerations

The study was reviewed and classified as public health surveillance by both the UCSF Office of Human Research Protections and the Stanford Medical Center institutional review board, based on the definition of surveillance in the US 2018 Revised Common Rule (45 CFR 46.102(l)(2). Official support from and engagement with the local county health departments was obtained. Participants signed separate consent forms for inclusion in the main study and for banking of remnant samples for future testing. Participants indicated at enrolment whether or not they were willing to be contacted for recruitment into future studies.

Results

We enrolled 3842 participants, continuing recruitment beyond our desired sample size of 3400 to ensure adequate numbers of enrolled adults from high-risk strata. Comparison of the desired sample sizes by county and census tract strata, and actual numbers of enrolled participants, is shown in Table 2 . Overall, we enrolled the desired number of participants except from Contra Costa County, due to delays in setting up testing sites. The proportion enrolled from high-risk strata in Santa Clara (88%) was also slightly lower than desired.

Table 2

Enrolment: desired sample size (SS), and the number and proportion of participants enrolled, by county and census tract risk strata.

	Census Tract Risk Strata
	Low		Medium		High		Total
County	SS	EnrolledN (%)	SS	EnrolleddN (%)	SS	EnrolleddN (%)	SS	EnrolleddN (%)
Alameda*	116	116 (100%)	421	521 (124%)	261	263 (101%)	798	900 (113%)
Contra Costa	152	159 (105%)	334	271 (81%)	61	37 (61%)	547	467 (85%)
Marin	57	61 (107%)	194	295 (152%)	49	76 (155%)	300	432 (144%)
San Francisco	66	74 (112%)	307	407 (133%)	130	185 (142%)	503	666 (132%)
San Mateo	7	10 (143%)	171	249 (146%)	189	251 (133%)	367	510 (139%)
Santa Clara	304	293 (96%)	475	481 (101%)	106	93 (88%)	885	867 (98%)
All	702	713 (102%)	1902	2224 (117%)	796	905 (114%)	3400	3842 (113%)

Includes the City of Berkeley, which has its own Department of Health.

The response rate, or the proportion of HHs from which a participant was enrolled, from among the number of targeted HHs, is shown in Table 3 . Our overall response rate was 6%– 9% from low-risk, 7% from medium-risk, and 4% from high-risk strata. Retention at the five-month follow-up visit, meaning completion of the questionnaire as well as providing specimens for testing (NP swab and venous blood) was 86.6%. All participants who completed a follow-up survey also agreed to be tested.

Table 3

Response rate: number of households targeted for recruitment and the response (number and proportion of participants enrolled), by county and census tract risk strata

	Risk Strata
	Low		Medium		High		Total
County	Households targetedN	EnrolledN (%)	Households targetedN	EnrolledN (%)	Households targetedN	EnrolledN (%)	Households targetedN	EnrolledN (%)
Alameda	1300	116 (9%)	7000	521 (7%)	6500	263 (4%)	14,800	900 (6%)
Contra Costa	1700	159 (9%)	5600	271 (5%)	1600	37 (2%)	8900	467 (5%)
Marin	700	61 (9%)	3300	295 (9%)	1300	76 (6%)	5300	432 (8%)
San Francisco	1019	74 (7%)	5329	407 (8%)	3247	185 (6%)	9595	666 (7%)
San Mateo	159	10 (6%)	3080	249 (8%)	4678	251 (5%)	7917	510 (6%)
Santa Clara	3122	293 (9%)	7491	481 (6%)	2875	93 (3%)	13,488	867 (6%)
All Counties	8000	713 (9%)	31,800	2224 (7%)	20,200	905 (4%)	60,000	3842 (6%)

THG attempted phone calls to 21,918 residences for which we had associated telephone numbers (36.5% of the 60,000 HH sample). Among the 9258 persons who were reached, 1390 (15.0%) were not associated with the address listed in the sample, 6196 (66.9%) refused participation, and 1014 (10.9%) enrolled on the phone or on the website. Among those who refused, 2095 (33.8%) hung up the phone before indicating why they were not interested, 1413 (18.4%) said they didn't want to participate in a study, and 1697 (27.4%) did not provide a reason. Only 135 refused because they didn't want to be tested; 200 were uninterested because they felt the study required too much time.

CBOs and a survey team approached 1590 HHs in selected high risk census tracts in 5 of the 6 counties, from which 119 (7.5%) eligible adults were enrolled at the time of canvassing. This is nearly twice the response to mailers from persons in high-risk strata, and is likely an underestimate of response, as we could not track the number of persons from these HHs who decided to enroll later.

Discussion

We designed and implemented a longitudinal cohort study to recruit a probability sample of adults in the San Francisco Bay Area to estimate the population-level incidence and prevalence of SARS-CoV-2 infection. One of the main strengths of the study was the use of stratified random sampling that relied on an address-based sampling frame. The U.S. Postal Service Delivery Sequence File provides a nearly complete list of all addresses in the country, and its use in defining our sampling frame will reduce bias compared to other non-representative, but easier to implement, sampling schemes. In our study, we did not enroll persons without housing, and excluded those living in nursing homes, homeless shelters and prisons, where rates of infection were extremely high [25,[60], [61], [62], [63], [64]]. Thus, our estimates will not represent infection among these groups.

We used regression models to predict the number of infections within census tracts, as a means of identifying strata for sampling, which was a unique feature of our study. The goal of using models was not to estimate the prevalence or incidence in particular regions, but rather to identify correlates of infection to categorize strata such that a weighted stratified sample would provide more precise estimates than a strictly random sample. Predicting cases, or the ‘risk’ within geographic areas, had the advantage of not being sensitive to short-term fluctuations in the local pandemic, such as a contained outbreak of infections. On the other hand, predicted risk strata would not reflect overall shifts in infection rates among different communities as the pandemic expanded. Use of a prediction model was based on several assumptions, however, one of which was that the risk level of all HH adults living within a census tract was the same. To evaluate this, we estimated precision based on different probabilities of misclassifying HH risk, and confirmed that even with moderate misclassification, stratification would improve precision. We also assumed that the risk, or at least the comparative risk between strata, would remain constant during the study period. Finally, we assumed that the socio-demographic characteristics we included in the model were reflective of risk. Several of the predictors included in the model have empirically been shown to be associated with higher rates of infection, including being LatinX and having low income [25,30,38,39,60].

Our overall goal was to estimate incidence and prevalence among the ‘general adult population’. The Bay Area, however, is highly heterogeneous, and includes many first- and second-generation immigrants from around the globe. Due to logistical constraints and available funding, the study was not designed to determine the incidence or prevalence by risk strata, county, or race/ethnicity with precision; therefore, outcome estimates will reflect an average across communities. We will also not be able to determine precise outcomes at specific points in time; this limits the interpretation of results, as the trajectory of the local pandemic changed during the study period, with a significant surge in reported cases in November and December 2020 [65,66]. Rates of infection have also been influenced by masking, social-distancing requirements, and the roll-out of COVID-19 vaccines.

A limitation of this study, as well as of other similar surveys, is non-response bias. Although weighting can be used to account for socio-demographic differences between the enrolled sample and the general population, the validity of results relies on the assumption that those who respond are similar to those who do not. Evaluating characteristics of non-responders requires reaching and surveying them, which is often impractical. Our overall response rate was 6%, which is what we assumed when designing the study. We also attempted recruitment by phone, but only one-third of HH addresses in our sampling frame were linked to a phone number. And although phone calls increased enrolment slightly, this approach required significant staff effort. Finally, we collaborated with CBOs to increase inclusion of participants from communities with the most barriers to participation. Other study design features, however, may have negatively affected response, such as the requirement to visit sites for sample collection. We tried to reduce this barrier by placing study sites throughout the 6 counties, by including reimbursement for transportation, and by making evening and weekend appointments available. We also increased reimbursement from $25 to $100 during the last 2 months of enrolment. The combination of these methods allowed us to reach our sample size goals. Despite these efforts, determining the sampling scheme, and recruiting and following a population-based cohort were logistically complicated, time intensive and costly. We began designing the study in April 2020, and recruited our first participants in August of that year. Enrolment of the cohort itself took 5 months, which was longer than we anticipated.

An additional strength of our study was the use of multiple antibody tests and viral PCR detection which will increase our ability to identify SARS-CoV-2 infections. We used tests able to detect antibodies to both nucleocapsid and spike RBD proteins, and positive samples were additionally evaluated for neutralizing and ACE-2 receptor-binding antibodies [48]. A variety of antibodies can be generated in response to SARS-CoV-2 that may not be detected by testing for antibodies to only one antigen [67,68]. In addition, the antibody tests we used were of relatively high test performance so that even with a low population prevalence, the likelihood that a positive test indicated a true infection will be improved, and the possibility of missing an infection reduced. Because COVID-19 vaccination began at the end of the enrolment period and during follow-up, tests that detect the presence of anti-nucleocapsid antibodies can help identify vaccine breakthrough infections or those that occur before vaccine immunity has developed, whereas antibodies to spike-protein can develop in response to immunization and therefore may not indicate a true infection [55,69]. Viral detection in combination with antibody testing and monthly specimen collection will allow us to assess the relationship of antibody production to viral shedding, the frequency of asymptomatic infections, and short-term persistence of antibodies. Neutralizing antibody tests provide additional information about humoral immunity in response to infection [22].

Although undergoing repeated NP swabs and venous blood draws can be uncomfortable, we chose these sample collection methods because at the start of the study, PCR testing of other specimens (such as anterior nasal swabs), as well as rapid tests for antigen and antibody detection had not been fully developed [70,71]. Since then, rapid antigen testing and PCR testing of self-collected nasal swabs [72] and saliva [73] have been shown to be fairly accurate, and are being used in various settings; additionally, some studies are using finger-prick capillary blood samples [34,74] to test for presence of antibodies. Self-collection of samples at our testing sites or by using mail-in home test kits would likely have increased response. However, we decided not to change our testing platforms and algorithm midway through the study, to avoid accounting for potential differences in test performance. And despite the discomfort of testing, those who enrolled in the study continued with follow-up, enhanced by the personal and ongoing interaction with site staff, including physicians and nurses, and telephone calls from staff or THG whenever a participant missed a visit.

Assembly of a longitudinal general population cohort such as TrackCOVID can be used as a platform for evaluating a variety of questions. Almost all participants agreed to be contacted for further studies. We administered a supplemental questionnaire in December 2020, just prior to vaccine roll-out, that inquired about attitudes, beliefs, and willingness to receive a COVID-19 immunization[75]. The response was high and results indicated disparities in vaccine intention by race/ethnicity, even among persons working in health care. In addition, participants are being enrolled in a follow-up study to identify breakthrough infections among those who have been vaccinated, and re-infections among whose had previously had COVID-19.

One of the aims of the study was to inform and collaborate with the public health departments in participating counties. We developed a real-time dashboard of study results that was available to counties [76]. The dashboard contained information on study recruitment, incidence and prevalence of infection, retention, sociodemographic characteristics of infected participants compared to the overall cohort, and vaccine uptake. Data were presented for the overall cohort as well as by county. Monthly meetings with county health departments were instituted to obtain their feedback and share information that could inform policy.

Conclusions

The design of this study can provide guidance for other surveys, while acknowledging the inherent difficulties in recruiting a population-based sample and the restrictions on interpretation of results. The project was enabled by collaboration with public health departments that were significantly invested in our findings and provided ongoing resources and feedback during study planning and implementation. Overall, designing and implementing a study to enroll a representative sample of the general population is challenging and requires a strong multi-disciplinary team. Employing multiple methods of recruitment, including through involvement of CBOs trusted by the local population, can also be helpful.

Acknowledgments

We would like to acknowledge the following for their invaluable assistance with the study: all of the study participants; the Departments of Public Health in the counties of San Francisco, Contra Costa, Alameda, Marin, San Mateo and Santa Clara, and in the City of Berkeley; the TrackCOVID clinical research coordinators; the staff of the The Henne Group; Mansour Fahimi, PhD, of the Marketing Systems Group for assistance with determining the sampling frame; community-based organizations for their help with recruitment (the Community Health Partnership, the Canal Alliance, and The Bay Area Community Health Advisory Council); Joe DeRisi, PhD and the Chan-Zuckerberg BioHub for PCR testing, viral sequencing and storage of samples; the clinical research lab at the Stanford Health Center; the UCSF Clinical Laboratory; Jing Jin and Graham Simmons (Vitalant Research Institute), and Rodolfo Villa (UCSF) for assistance with pseudovirus neutralization assays; David Glidden, PhD, UCSF for assistance with initial study design; Sara Covin, UCSF for assistance with references.

This study was funded by a grant from the Chan-Zuckerberg Initiative (CZI). The content of this article is solely the responsibility of the authors. CZI was not involved in the study design, the collection, analysis or reporting of results, the writing of the manuscript or the decision to submit the article for publication.

Family	Variable Name	LASSO coef.
Race/Ethnicity
	% Hispanic (overall)	0.23
	% Central American	0.08
	% Mexican	0.02
	% Black
	% Foreign-born	0.05
	% Native American	-0.01
	% Native Hawaiian/Pacific Islander	0.04
	% Southeast Asian
	% South American	0.02
	% Asian (overall)
	% East Asian	-0.04
	% South Asian	-0.14
	% White, non-Hispanic	-0.12
Age / Gender
	% 18 - 40 years old	-0.02
	% Male	0.02
	% < 5 years of age
	% < 18 years old
	% Households with a resident younger than 18
	% Households with a resident older that 65
	% > 65 years old
Education
	% Less than high school
	% College-educated
Socio-Economic
	Teen-birth rate (% of women who gave birth before 20)*	0.07
	% Households with more occupants than bedrooms
	% Households on food stamp / SNAP benefits in the last year
	% Households without internet
	% Households classified as “extremely low income” (making less than 30% of the HUD Area Median Family Income)*
	% Not fluent in English
	% Households that spend > 50% of income on rent*
	% Households that are single-family homes
	% Households earning below 1.25 the poverty line	0.06
	Incarceration rate* (% of children who grew up in this census tract who were in jail on April 1, 2010)	0.01
	% Households without vehicle access
	Overall population density (per square mile)
	Average number of occupants/household
	Unemployment rate (% of 16+ population without a job)
	Food desert (binary variable: is there grocery store access within 0.5 miles for urban areas and 10 miles for rural ones?)*	-0.03
	Eviction-filing rate* (% of renter occupied housing units that have evictions filed)	-0.02
	Gini index (measure of income inequality)
	Traffic density (vehicle-kms / hour / road length within 150 m of census tract boundary: percentile)*

Family	Variable Name	LASSO coef.
	% Families moved in the last year	-0.02
	% with limited public transit (no stops within half a mile)*
	% Households that own (vs rent)
	Median rent
	Median house price	0.01
	Median household income	0.06
Job
	% Employed population working in service
	% Employed population working in production / transportation	0.04
	% Employed population working in construction / natural resources
	% Employed population working in sales / office work
	% Employed population working in military
	% Employed population working in management
Commute
	% Commute by carpool
	% Commute by public transit
	% Commute by bike or walk
	% Commute lasts <15min	-0.02
	% Commute lasts > 1hr
	Avg commute time	-0.07
	% Commute by car (solo)
	% Work from home	-0.03
Health
	% Without health insurance
	ER visits for asthma/capita*	0.11
	% Adults with poor physical health*	0.06
	% Population with a disability*	0.04
	% Adults with poor mental health*
	% Adults who get annual checkup*

*Data for variable obtained from the UCSF Health Atlas (36). Data for all other variables were taken from the ACS 2018 (37).